db9小波包分解matlab,小波包分解 仿真 matlab

clc;

close all;

clear;

fs = 100000;

t = 1: 100;

x = sin(2*pi*4000* t/fs) + sin(2*pi*40000*t/fs);

%db8

[Lo_D1, Hi_D1] = wfilters('db1', 'd');%从db1到db8. 滤波器系数个数不同 2

[Lo_D2, Hi_D2] = wfilters('db2', 'd');%从db1到db8. 4

[Lo_D3, Hi_D3] = wfilters('db4', 'd');%从db1到db8. 8

[Lo_D4, Hi_D4] = wfilters('db8', 'd');%从db1到db8. 16

freqz(Lo_D1);

hold on;

freqz(Lo_D2);

hold on;

freqz(Lo_D3);

hold on;

freqz(Lo_D4);%等级越高，衰减越厉害，过渡带越窄

close all;

%cwt连续小波变换，dwt离散小波变换

%dwt仅仅对单一尺度的变化

[cA1, cD1 ]= dwt(x, 'db1');%cA1低频分量 信号的近似， cD1高频分量 信号的细节

%A1=upcoef('a',cA1,'db1');%单尺度低频系数向上一步的重构信号

%D1=upcoef('a',cA2,'db1');%高频系数 ...

X = idwt(cA1, cD1, 'db1');%

figure(10);

subplot(3, 1,1)

plot(x);

title('原始');

subplot(3, 1,2)

plot(X);

title('重构');

subplot(3, 1,3)

plot(x-X);

title('误差');

%wavedec多尺度分解 仅仅对低频分量分解

%[C,L] = WAVEDEC(X,N,Lo_D,Hi_D) 已有的滤波器系数

%wavedec 由多次调用 dwt ，将输出的低频分量，再次传入，如此

[C, L]= wavedec(x, 2, 'db1');% L对应于每个节点