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


subplot（221）;
plot（s）;
title（‘原信號(hào)‘）;

wpt=wpdec（s3‘sym2‘‘shannon‘）;
%?plot（wpt）;
%第一個(gè)節(jié)點(diǎn)的stein無偏似然估計(jì)閾值去噪
s0=wpcoef（wpt[30]）;
s30=wprcoef（wpt[30]）;
[s0xs0ls0]=wden（s0‘rigrsure‘‘s‘‘one‘1‘sym2‘）;
si=zeros（18251）
%中間6個(gè)節(jié)點(diǎn)的滑動(dòng)極小極大準(zhǔn)則閾值去噪
for?i=1:6;
????a=10
????s3_i=wpcoef（wpt[3i]）;%提取小波系數(shù)
????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）;
%最后一個(gè)節(jié)點(diǎn)的固定形式閾值去