m=load（‘G:\myjob\1-32\動車抗蛇行全拆\sf01a_27_H.txt‘）;
t=m（1:145801）;
s=m（1:145802）;


subplot（221）;
plot（s）;
title（‘原信號‘）;

wpt=wpdec（s3‘sym2‘‘shannon‘）;
%?plot（wpt）;
%第一個節點的stein無偏似然估計閾值去噪
s0=wpcoef（wpt[30]）;
s30=wprcoef（wpt[30]）;
[s0xs0ls0]=wden（s0‘rigrsure‘‘s‘‘one‘1‘sym2‘）;
si=zeros（18251）
%中間6個節點的滑動極小極大準則閾值去噪
for?i=1:6;
????a=10
????s3_i=wpcoef（wpt[3i]）;%提取小波系數
????thr=thselect（s3_i‘minimaxi‘）;
????THR=a*thr
????sy=wthresh（s3_i‘s‘THR）;
????si=cat（2sisy）;
end
si=si（:2:7）;
si_06=cat（2s0si）;
%最后一個節點的固定形式閾值去