examp16_13.m

来自「小波分析理论及MATLAB7实现的例程,第16章小波工具箱的应用基础」· M 代码 · 共 28 行

clc;
clear;
%调入含突变点的信号
load cnoislop;
x=cnoislop;
N=length(x);
t=1:N;
figure(1);
plot(t,x,'LineWidth',2);xlabel('时间 t/s');ylabel('幅值 A');
%一维小波分解
[c,l] = wavedec(x,6,'db3');
%重构第1～6层逼近系数
a6 = wrcoef('a',c,l,'db3',6);
a5 = wrcoef('a',c,l,'db3',5);
a4 = wrcoef('a',c,l,'db3',4);
a3 = wrcoef('a',c,l,'db3',3);
a2 = wrcoef('a',c,l,'db3',2);
a1 = wrcoef('a',c,l,'db3',1);
%显示逼近系数
figure(2)
subplot(6,1,1);plot(a6,'LineWidth',2);ylabel('a6');
subplot(6,1,2);plot(a5,'LineWidth',2);ylabel('a5');
subplot(6,1,3);plot(a4,'LineWidth',2);ylabel('a4');
subplot(6,1,4);plot(a3,'LineWidth',2);ylabel('a3');
subplot(6,1,5);plot(a2,'LineWidth',2);ylabel('a2');
subplot(6,1,6);plot(a1,'LineWidth',2);ylabel('a1');
xlabel('时间 t/s');