matlab 小波变换

[Lo_D,Hi_D,Lo_R,Hi_R]=wfilters('haar');%获得Haar小波的尺度和小波系数
Lo_D     %打印低通分解滤波器系数
Hi_D     %打印高通分解滤波器系数
Lo_R     %打印低通重构滤波器系数
Hi_R     %打印高通重构滤波器系数

waveinfo('haar');   %打印Haar小波相关信息

[phi,psi,xval]=wavefun('haar',10);%wavefun函数产生尺度向量和小波向量
xaxis=zeros(size(xval));  %y轴
axis square;
title('Haar Scaling Function');
  subplot ( 122 ) ; plot ( xval , psi , ' k ' , xval , xaxis , ' - k ' ) ; % 绘制小波向量
 axis ( [ 0 1 -1.5 1.5 ] ) ; % y轴
axis square ;
title ( 'Haar Wavelet Function ' ) ;

f = magic (4 ) ;  % f 是4 X 4的魔术矩阵
f
% 对f 进行一次小波分解，c1是系数矩阵，s1是c1各分量的维数。
[C1,S1]=wavedec2(f,1,'haar');
C1
S1
[ C2 , S2 ] = wavedec2 ( f , 2 , ' haar ' ) ;
C2
S2

f = imread ( '你的图 ' ) ;  
imshow ( f ) ;
[C,S] = wavefast( f , 1 ,'sym4' ) ;
% 以4格灰度图像方式显示小波分解结果
figure , wave2gray ( C , S , -6 ) ;
[ nc , y ] = wavecut ( 'a' , C , S ) ;  % 去掉分解结果中的近似系数
figure , wave2gray ( nc , S , -6 ) ;  % 显示处理后的小波系数
% 小波反变换，得到边缘图像
edges = abs ( waveback ( nc , S , 'sym4' ) );
figure , imshow ( mat2gray ( edges ) ) ;   % 显示边缘图像

f = imread ( '你的图' ) ;
[ C , S ] = wavefast ( f , 4 , 'sym4' ) ;  % 4层小波分解
wave2gray ( C , S , 20 ); % 以图像形式显示小波系数
% 第一层分解的小波系数置为0 ，并显示反变换后的图像
[ C , g8 ] = wavezero ( C , S , 1 , 'sym4' ) ;  % wavezero 函数请见7.6节
% 第一，二层分解的小波系数置为0 ，并显示反变换后的图像
[ C , g8 ] = wavezero ( C , S , 2 , 'sym4' ) ;  
% 第一，二，三层分解的小波系数置为0 ，并显示反变换后的图像
[ C , g8 ] = wavezero ( C , S , 3 , 'sym4' ) ;
% 4层分解的小波系数都置为0 ，并显示反变换后的图像
[ C , g8 ] = wavezero ( C , S , 4 , 'sym4' ) ; 

f = imread ( '你的图' ) ;  
[ C , S ] = wavefast( f,4,'jpeg9.7') ;  % 4层小波分解
wave2gray ( C , S , 8 ); % 显示处理后的小波系数
f = wavecopy ('a', C , S ) ;  % f 是第4层分解后的近似系数
figure , imshow ( mat2gray ( f ) ) ;  % 显示第4层分解后的近似图像
[ C , S ] = waveback(C,S,'jpeg9.7',1) ;   % 一次重构
%  f 是第3层分解后的近似图像
f = wavecopy ('a',C,S) ;
% 显示第3层分解后的近似图像
figure , imshow ( mat2gray ( f ) ) ;
% 二次重构
[ C , S ] = waveback(C,S,'jpeg9.7',1) ;
%  f 是第2层分解后的近似系数
f = wavecopy ('a',C,S) ;
% 显示第2层分解后的近似图像
figure , imshow ( mat2gray ( f ) ) ;
[C,S]= waveback(C, S,'jpeg9.7',1);%三次重构
f = wavecopy('a', C, S);%f是第1层分解后的近似函数
figure, imshow(mat2gray(f));%显示第1层分解后的近似图像
[C, S]= waveback(C, S, 'jpeg9.7', 1); %四次重构
f = wavecopy('a', C, S);%f是最终重构后的图像
figure, imshow(mat2gray(f)); %显示重构后的最终图像

f = imread ( '你的图 ' ) ; 
[ C , S ] = wavefast( f,2,'db4') ;  % 4层小波分解
wave2gray(C,S); % 显示处理后的小波系数
figure,wave2gray(C,S,8);
figure,wave2gray(C,S,-8);