3、代码
%含噪声的三角波与正弦波的组合
%利用db5小波对信号进行7层分解
%生产正弦信号
clc;close all;clear all;
N=1000;
t=1:N;
sig1=sin(0.3*t);
%生成三角形波形
sig2(1:500)=((1:500)-1)/500;
sig2(501:N)=(1000-(501:1000))/500;
figure(1);
subplot(211);
plot(t,sig1,‘linewidth’,2);
xlabel(‘样本序号 N’);
ylabel(‘幅值A’);
subplot(212);
plot(t,sig2,‘linewidth’,2);
xlabel(‘样本序号 N’);
ylabel(‘幅值A’);
%叠加信号
x=sig1+sig2+randn(1,N);
figure(2);
plot(t,x,‘linewidth’,2);
xlabel(‘样本序号 N’);
ylabel(‘幅值A’);%一维小波分解
[c,l]=wavedec(x,7,‘db5’);%重构第1-7层逼近系数
a7=wrcoef(‘a’,c,l,‘db5’,7);
a6=wrcoef(‘a’,c