%??mallet_wavelet.m
%??此函數用于研究Mallet算法及濾波器設計
%??此函數僅用于消噪
a=pi/8;??????????????????%角度賦初值
b=pi/8;
%低通重構FIR濾波器h0（n）沖激響應賦值
h0=cos（a）*cos（b）;
h1=sin（a）*cos（b）;
h2=-sin（a）*sin（b）;
h3=cos（a）*sin（b）;
low_construct=[h0h1h2h3];
L_fre=4;?????????????????????????????%濾波器長度
low_decompose=low_construct（end:-1:1）;???%確定h0（-n），低通分解濾波器
for?i_high=1:L_fre;??????????????????????%確定h1（n）=（-1）^n高通重建濾波器
????if（mod（i_high2）==0）;
????coefficient=-1;
????else
????coefficient=1;
????end
????high_construct（1i_high）=low_decompose（1i_high）*coefficient;
end
high_decompose=high_construct（end:-1:1）;??%高通分解濾波器h1（-n）
L_signal=1578;???%信號長度
n=1:L_signal;???%信號賦值
f=10;
t=0.001;
a=load（[Name‘.txt‘]）;
y=a（st:n2:2）;?
h1=figure（1）;
plot（y）;
title（‘原信號‘）;
fileName=[Name‘.jpg‘];
saveas（h1fileName‘jpg‘）;
check1=sum（high_decompose）;?????????%h0（n）性質校驗
check2=sum（low_decompose）;
check3=norm（high_decompose）;
check4=norm（low_decompose）;
l_fre=conv（ylow_decompose）;????????%卷積
l_fre_down=dyaddown（l_fre）;?????????%抽取得低頻細節
h_fre=conv（yhigh_decompose）;
h_fre_down=dyaddown（h_fre）;?????????%信號高頻細節
h2=figure（2）;
subplot（211）
plot（l_fre_down）;
title（‘小波分解的低頻系數‘）;
subplot（212）;
plot（h_fre_down）;
title（‘小波分解的高頻系數‘）;
fileName=[Name‘_小波分解高頻和低頻系數‘‘.jpg‘];??%截取后平滑處理后的圖片
saveas（h2fileName‘jpg‘）;
l_fre_pull=dyadup（l_fre_down）;?????%0差值
h_fre_pull=dyadup（h_fre_down）;
l_fre_denoise=conv（low_constructl_fre_pull）;
h_fre_denoise=conv（high_constructh_fre_pull）;
l_fre_keep=wkeep（l_fre_denoiseL_signal）;???%取結果的中心部分消除卷積影響
h_fre_keep=wkeep（h_fre_denoiseL_signal）;
sig_denoise=l_fre_keep+h_fre_keep;??????????%信號重構
compare=sig_denoise-y;??????????????????????%與原信號比較
h3=figure（3）;
subplot（311）
plot（y）;???????????????????????????????????
ylabel（‘y‘）;?????????????%原信號
subplot（312）;
plot（sig_denoise）;??????????????????????????
ylabel（‘sig\_denoise‘）;??%重構信號
subplot（313）;
plot（compare）;
ylabel（‘compare‘）;???????%原信號與消噪后信號的比較
fileName=[Name‘__對比‘‘.jpg‘];??%截取后平滑處理后的圖片
saveas（h3fileName‘jpg‘）;

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

?????文件???????2256??2009-03-18?20:17??Untitled2.m

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

?????????????????2256????????????????????1