matlab实现s小波变换代码

suffer

输入图像矩阵，输出一幅图像，分为四块；

1. function wvlet=waveup1(x)
   
2. %function [ll,hl,lh,hh]=waveup1(x)
   
3. x=double(x);
   
4. wvlet=x;
   
5. [hp,lp]=size(x);
   
6. hc=hp/2;
   
7. lc=lp/2;
   
8. %奇偶列重排
   
9. for n=1:lc
   
10.         j=n*2-1;
    
11.         wvlet(:,n)=x(:,j);
    
12.         wvlet(:,lc+n)=x(:,j+1);
    
13. end
    
14. %行变换开始
    
15. for n=1:lc-1       
    
16.         k=floor((wvlet(:,n)+wvlet(:,n+1))/2);
    
17.         wvlet(:,lc+n)=wvlet(:,lc+n)-k;
    
18. end
    
19.         k=floor((wvlet(:,lc)+wvlet(:,lc-1))/2);   
    
20.         wvlet(:,lp)=wvlet(:,lp)-k;
    
21. 
    
22.         k=floor((wvlet(:,lc+2)+wvlet(:,lc+1))/4+0.5);
    
23.         wvlet(:,1)=wvlet(:,1)+k;
    
24. for n=2:lc
    
25.         k=floor((wvlet(:,lc+n)+wvlet(:,lc+n-1))/4+0.5);
    
26.         wvlet(:,n)=wvlet(:,n)+k;
    
27. end
    
28. %行变换结束
    
29. x=wvlet;
    
30. %奇偶行重排
    
31. for n=1:hc
    
32.         j=n*2-1;
    
33.         wvlet(n,:)=x(j,:);
    
34.         wvlet(hc+n,:)=x(j+1,:);
    
35. end
    
36. %列重排开始       
    
37. for n=1:hc-1       
    
38.         k=floor((wvlet(n,:)+wvlet(n+1,:))/2);
    
39.         wvlet(hc+n,:)=wvlet(hc+n,:)-k;
    
40. end
    
41.         k=floor((wvlet(hc,:)+wvlet(hc-1,:))/2);
    
42.         wvlet(hp,:)=wvlet(hp,:)-k;
    
43. 
    
44.         k=floor((wvlet(hc+2,:)+wvlet(hc+1,:))/4+0.5);
    
45.         wvlet(1,:)=wvlet(1,:)+k;
    
46. for n=2:hc
    
47.         k=floor((wvlet(hc+n,:)+wvlet(hc+n-1,:))/4+0.5);
    
48.         wvlet(n,:)=wvlet(n,:)+k;
    
49. end
    
50. %列重排结束
    
51. wvlet=uint8(wvlet);
    
52. ll=wvlet(1:hc,1:lc);
    
53. hl=wvlet(1:hc,lc+1:lp);
    
54. lh=wvlet(hc+1:hp,1:lc);
    
55. hh=wvlet(hc+1:hp,lc+1:lp);

复制代码

suffer

1. %function wvlet=waveup2(x)
   
2. function [ll,hl,lh,hh]=waveup2(x)
   
3. x=double(x);
   
4. wvlet=x;
   
5. [hp,lp]=size(x);
   
6. hc=hp/2;
   
7. lc=lp/2;
   
8. %奇偶列重排
   
9. for n=1:lc
   
10.         j=n*2-1;
    
11.         wvlet(:,n)=x(:,j);
    
12.         wvlet(:,lc+n)=x(:,j+1);
    
13. end
    
14. %行变换开始
    
15. for n=1:lc-1       
    
16.         k=(wvlet(:,n)+wvlet(:,n+1))/2;
    
17.         wvlet(:,lc+n)=wvlet(:,lc+n)-k;
    
18. end
    
19.         k=(wvlet(:,lc)+wvlet(:,lc-1))/2;   
    
20.         wvlet(:,lp)=wvlet(:,lp)-k;
    
21. 
    
22.         k=(wvlet(:,lc+2)+wvlet(:,lc+1))/4;
    
23.         wvlet(:,1)=wvlet(:,1)+k;
    
24. for n=2:lc
    
25.         k=(wvlet(:,lc+n)+wvlet(:,lc+n-1))/4;
    
26.         wvlet(:,n)=wvlet(:,n)+k;
    
27. end
    
28. %行变换结束
    
29. x=wvlet;
    
30. %奇偶行重排
    
31. for n=1:hc
    
32.         j=n*2-1;
    
33.         wvlet(n,:)=x(j,:);
    
34.         wvlet(hc+n,:)=x(j+1,:);
    
35. end
    
36. %列重排开始       
    
37. for n=1:hc-1       
    
38.         k=(wvlet(n,:)+wvlet(n+1,:))/2;
    
39.         wvlet(hc+n,:)=wvlet(hc+n,:)-k;
    
40. end
    
41.         k=(wvlet(hc,:)+wvlet(hc-1,:))/2;
    
42.         wvlet(hp,:)=wvlet(hp,:)-k;
    
43. 
    
44.         k=(wvlet(hc+2,:)+wvlet(hc+1,:))/4;
    
45.         wvlet(1,:)=wvlet(1,:)+k;
    
46. for n=2:hc
    
47.         k=(wvlet(hc+n,:)+wvlet(hc+n-1,:))/4;
    
48.         wvlet(n,:)=wvlet(n,:)+k;
    
49. end
    
50. %列重排结束
    
51. %wvlet=uint8(wvlet);
    
52. ll=wvlet(1:hc,1:lc);
    
53. hl=wvlet(1:hc,lc+1:lp);
    
54. lh=wvlet(hc+1:hp,1:lc);
    
55. hh=wvlet(hc+1:hp,lc+1:lp);

复制代码

suffer

1. function wvlet=waveup3(x)
   
2. %function [ll,hl,lh,hh]=waveup3(x)
   
3. x=double(x);
   
4. wvlet=x;
   
5. [hp,lp]=size(x);
   
6. hc=hp/2;
   
7. lc=lp/2;
   
8. %奇偶列重排
   
9. for n=1:lc
   
10.         j=n*2-1;
    
11.         wvlet(:,n)=x(:,j);
    
12.         wvlet(:,lc+n)=x(:,j+1);
    
13. end
    
14. %行变换开始
    
15. for n=1:lc       
    
16.         wvlet(:,lc+n)=wvlet(:,lc+n)-wvlet(:,n);
    
17. end
    
18. for n=1:lc
    
19.     k=wvlet(:,lc+n)/2-mod(wvlet(:,lc+n)/2,1);
    
20.         wvlet(:,n)=wvlet(:,n)+k;
    
21. end
    
22. %行变换结束
    
23. x=wvlet;
    
24. %奇偶行重排
    
25. for n=1:hc
    
26.         j=n*2-1;
    
27.         wvlet(n,:)=x(j,:);
    
28.         wvlet(hc+n,:)=x(j+1,:);
    
29. end
    
30. %列重排开始       
    
31. for n=1:hc       
    
32.         wvlet(hc+n,:)=wvlet(hc+n,:)-wvlet(n,:);
    
33. end
    
34. for n=1:hc
    
35.         k=wvlet(hc+n,:)/2-mod(wvlet(hc+n,:)/2,1);
    
36.         wvlet(n,:)=wvlet(n,:)+k;
    
37. end
    
38. %列重排结束
    
39. %wvlet=uint8(wvlet);
    
40. ll=wvlet(1:hc,1:lc);
    
41. hl=wvlet(1:hc,lc+1:lp);
    
42. lh=wvlet(hc+1:hp,1:lc);
    
43. hh=wvlet(hc+1:hp,lc+1:lp);

复制代码

PSZ10

s小波是哪个小波啊，能不能详细描述一下这个程序的功能啊

suffer

该函数的功能：输入图像矩阵，输出一幅图像，分为四块

上面已经有相关说明了