www.pudn.com > matlabxiaobobianhuan.rar > 多尺度分解2.m, change:2009-05-11,size:1762b


s=af1522(1000:8000); 
subplot(go);plot(s);title('原始信号'); 
Xlabel('时间');Ylabel('幅值'); 
[c,l]=wavedec(s,7,'db5'); 
w=waverec(c,l,'db5'); 
ca1=appcoef(c,l,'db5',1); 
cd1=detcoef(c,l,1); 
subplot(8,2,3);plot(ca1);title('第一层低频信号'); 
Xlabel('时间');Ylabel('幅值'); 
subplot(8,2,4);plot(cd1);title('第一层高频信号'); 
Xlabel('时间');Ylabel('幅值'); 
ca2=appcoef(c,l,'db5',2); 
cd2=detcoef(c,l,2); 
subplot(8,2,5);plot(ca2);title('第二层低频信号'); 
Xlabel('时间');Ylabel('幅值'); 
subplot(8,2,6);plot(cd2);title('第二层高频信号'); 
Xlabel('时间');Ylabel('幅值'); 
ca3=appcoef(c,l,'db5',3); 
cd3=detcoef(c,l,3); 
subplot(8,2,7);plot(ca3);title('第三层低频信号'); 
Xlabel('时间');Ylabel('幅值'); 
subplot(8,2,8);plot(cd3);title('第三层高频信号'); 
Xlabel('时间');Ylabel('幅值'); 
ca4=appcoef(c,l,'db5',4); 
cd4=detcoef(c,l,4); 
subplot(8,2,9);plot(ca4);title('第四层低频信号'); 
Xlabel('时间');Ylabel('幅值'); 
subplot(8,2,10);plot(cd4);title('第四层高频信号'); 
Xlabel('时间');Ylabel('幅值'); 
ca5=appcoef(c,l,'db5',5); 
cd5=detcoef(c,l,5); 
subplot(8,2,11);plot(ca5);title('第五层低频信号'); 
Xlabel('时间');Ylabel('幅值'); 
subplot(8,2,12);plot(cd5);title('第五层高频信号'); 
Xlabel('时间');Ylabel('幅值'); 
ca6=appcoef(c,l,'db5',6); 
cd6=detcoef(c,l,6); 
subplot(8,2,13);plot(ca6);title('第六层低频信号'); 
Xlabel('时间');Ylabel('幅值'); 
subplot(8,2,14);plot(cd6);title('第六层高频信号'); 
Xlabel('时间');Ylabel('幅值'); 
ca7=appcoef(c,l,'db5',7); 
cd7=detcoef(c,l,7); 
subplot(8,2,15);plot(ca7);title('第七层低频信号'); 
Xlabel('时间');Ylabel('幅值'); 
subplot(8,2,16);plot(cd7);title('第七层高频信号'); 
Xlabel('时间');Ylabel('幅值'); 
a=waverec(c,l,'db5'); 
subplot(8,2,16);plot(a);title('重构信号'); 
Xlabel('时间');Ylabel('幅值'); 
err=norm(s-a);