小波包阈值去噪

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）;
%最后一个节点的固定形式阈值去