%??程序作者：沙威，香港大學(xué)電氣電子工程學(xué)系，wsha@eee.hku.hk
%??參考文獻(xiàn)：小波分析理論與MATLAB?R2007實(shí)現(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

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

?????文件???????1065??2010-03-18?15:53??Wavelet_IRLS\DWT.m

?????文件????????747??2012-12-27?22:25??Wavelet_IRLS\irls.m

?????文件??????66614??2004-06-30?00:47??Wavelet_IRLS\lena.bmp

?????文件???????7409??2013-01-02?22:49??Wavelet_IRLS\Wavelet_IRLS.m

?????目錄??????????0??2013-01-19?12:51??Wavelet_IRLS

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

????????????????75835????????????????????5