%MATLAB2維小波變換經(jīng)典程序?
??
%??FWT_DB.M;
%??此示意程序用DWT實(shí)現(xiàn)二維小波變換
%??編程時(shí)間2004-4-10，編程人沙威
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;
clc;
T=256;???????%??圖像維數(shù)
SUB_T=T/2;???%??子圖維數(shù)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%??1.調(diào)原始圖像矩陣
load?wbarb;??%??下載圖像
f=X;?????????%??原始圖像
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%??2.進(jìn)行二維小波分解
l=wfilters（‘db10‘‘l‘）;????%??db10（消失矩為10）低通分解濾波器沖擊響應(yīng)（長(zhǎng)度為20）
L=T-length（l）;
l_zeros=[lzeros（1L）];????%??矩陣行數(shù)與輸入圖像一致，為2的整數(shù)冪
h=wfilters（‘db10‘‘h‘）;????%??db10（消失矩為10）高通分解濾波器沖擊響應(yīng)（長(zhǎng)度為20）
h_zeros=[hzeros（1L）];????%??矩陣行數(shù)與輸入圖像一致，為2的整數(shù)冪
for?i=1:T;???%??列變換
????row（1:SUB_Ti）=dyaddown（?ifft（?fft（l_zeros）.*fft（f（:i）‘）?）?）.‘;????%??圓周卷積<->FFT
????row（SUB_T+1:Ti）=dyaddown（?ifft（?fft（h_zeros）.*fft（f（:i）‘）?）?）.‘;??%??圓周卷積<->FFT
end;
for?j=1:T;???%??行變換
????line（j1:SUB_T）=dyaddown（?ifft（?fft（l_zeros）.*fft（row（j:））?）?）;????%??圓周卷積<->FFT
????line（jSUB_T+1:T）=dyaddown（?ifft（?f

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????4962??2007-09-05?15:34??xiaobo.m

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

?????????????????4962????????????????????1