小波变换matlab代码121047

小波变换代码包。 调用形式：ww=DWT(N) N为数据大小，返回变换系数矩阵。 使用举例 X=imread('lena256.bmp'); X=double(X); % 小波变换矩阵生成 ww=DWT(a); % 小波变换让图像稀疏化（注意该步骤会耗费时间，但是会增大稀疏度） X1=ww*sparse(X)*ww';

%  程序作者：沙威，香港大学电气电子工程学系，wsha@eee.hku.hk
%  参考文献：小波分析理论与MATLAB R2007实现，葛哲学，沙威，第20章  小波变换在矩阵方程求解中的应用（沙威、陈明生编写）.

%  构造正交小波变换矩阵，图像大小N*N，N=2^P，P是整数。

function ww=DWT（N）

[hg]= wfilters（‘sym8‘‘d‘）;       %  分解低通和高通滤波器

% N=256;                           %  矩阵维数（大小为2的整数幂次）
L=length（h）;                       %  滤波器长度
rank_max=log2（N）;                  %  最大层数
rank_min=double（int8（log2（L）））+1;  %  最小层数
ww=1;   %  预处理矩阵

%  矩阵构造
for jj=rank_min