%??程序作者：沙威，香港大學(xué)電氣電子工程學(xué)系，wsha@eee.hku.hk
%??參考文獻：小波分析理論與MATLAB?R2007實現(xiàn)，葛哲學(xué)，沙威，第20章??小波變換在矩陣方程求解中的應(yīng)用（沙威、陳明生編寫）.

%??構(gòu)造正交小波變換矩陣，圖像大小N*N，N=2^P，P是整數(shù)。

function?ww=DWT（N）

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

%?N=256;???????????????????????????%??矩陣維數(shù)（大小為2的整數(shù)冪次）
L=length（h）;???????????????????????%??濾波器長度
rank_max=log2（N）;??????????????????%??最大層數(shù)
rank_min=double（int8（log2（L）））+1;??%??最小層數(shù)
ww=1;???%??預(yù)處理矩陣

%??矩陣構(gòu)造
for?jj=rank_min:rank_max
????
????nn=2^jj;
????
????%??構(gòu)造向量
????p1_0=sparse（[hzeros（1nn-L）]）;
????p2_0=sparse（[gzeros（1nn-L）]）;
????
????%??向量圓周移位
????for?ii=1:nn/2
????????p1（ii:）=circshift（p1_0‘2*（ii-1））‘;
????????p2（ii:）=circshift（p2_0‘2*（ii-1））‘;
????end
????
????%??構(gòu)造正交矩陣
????w1=[p1;p2];
????mm=2^rank_max-length（w1）;
????w=spa

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????66614??2009-06-08?18:00??Wavelet_OMP\lena256.bmp

?????文件???????2499??2009-10-14?18:13??Wavelet_OMP\Wavelet_OMP.m

?????文件???????1065??2009-10-14?18:27??Wavelet_OMP\DWT.m

?????文件????????912??2010-04-01?16:44??Wavelet_OMP\wavelet_OMP_colour.m

?????目錄??????????0??2009-10-14?18:56??Wavelet_OMP

-----------?---------??----------?-----??----

????????????????71090????????????????????5