clear all;
load A4.txt;
figure(1)
%装载第二列信号
d=A4(:,4);
plot(d);title('1-3 原始信号');
xlabel('时间/s');ylabel('振幅/m');
%用db4小波对原始信号进行8层分解并提取系数
[C,L]=wavedec(d,8,'db8');
%用ddencmp()函数获得信号的默认阈值,使用wdencmp()命令函数实现去噪过程
[thr,sorh,keepapp]=ddencmp('den','wv',d);
d2=wdencmp('gbl',C,L,'db8',8,thr,sorh,keepapp);
figure(2)
plot(d2);title('2-3 默认阈值去噪信号');
xlabel('时间/s');ylabel('振幅/m');
%用db4小波对去噪信号进行8层分解
wpt=wpdec(d2,8,'db8','shannon');
figure(3)
plot(wpt);title('小波包树图');
%示意一下，重构的第8层125号节点的小波图
rcfs=wprcoef(wpt,[8 125]);
figure(4)
plot(rcfs); title('重构节点(8,125)小波包');
%第8层小波包能量
wpe=wenergy(wpt);
figure(5);
bar(wpe);
title('3-3 第8层小波包能谱');
xlabel('频带数');ylabel('能量')

高手帮我指导一下，这个程序对吗？谢谢