%??程序作者：沙威，香港大學電氣電子工程學系，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:rank_max
????
????nn=2^jj;
????
????%??構造向量
????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
????
????%??構造正交矩陣
????w1=[p1;p2];
????mm=2^rank_max-length（w1）;
????w=sparse（[w1zeros（length（w1）mm）;zeros（mmlength（w1））eye（mmmm）]）;
????ww=ww*w;
????
????clear?p1;clear?p2;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2009-10-14?18:56??Wavelet_OMP\
?????文件????????1065??2009-10-14?18:27??Wavelet_OMP\DWT.m
?????文件???????66614??2009-06-08?18:00??Wavelet_OMP\lena256.bmp
?????文件????????2499??2009-10-14?18:13??Wavelet_OMP\Wavelet_OMP.m