%% 学习目标:多层小波的重构
%% 程序1
clear all;
close all;
load noissin.mat;
x=noissin;
[C,L]=wavedec(x,3,'db1'); %小波多层分解
y=waverec(C,L,'db1'); %重构,必须小波类型一致
e=max(abs(x-y)) %重构的误差
%% 程序2
clear all;
close all;
load noissin.mat;
x=noissin;
[C,L]=wavedec(x,3,'db1'); %3层小波分解
A=wrcoef('a',C,L,'db1'); %对近似系数进行单支重构,去掉了高频部分
D1=wrcoef('d',C,L,'db1',1); %对细节系数进行单支重构
D2=wrcoef('d',C,L,'db1',2);
D3=wrcoef('d',C,L,'db1',3);
figure;
subplot(511);
plot(x);
subplot(512);
plot(A);
subplot(513);
plot(D1);
subplot(514);
plot(D2);
subplot(515);
plot(D3);
set(gcf,'position',[50,50,800,600]);