clear all
close all
clc
% 加載信號(hào)
load leleccum;
% 采樣點(diǎn)
x = leleccum(300:2348);

% x = importdata(‘V1.txt’);
n = 3;
wpt=wpdec(x,n,‘dmey’);
% Plot wavelet packet tree (binary tree)
plot(wpt)
%% wavelet packet coefficients.default：use the front 4.
cfs0=wpcoef(wpt,7);
cfs1=wpcoef(wpt,8);
cfs2=wpcoef(wpt,9);
cfs3=wpcoef(wpt,10);
figure;
subplot(5,1,1);
plot(x);
title(‘原始信號(hào)’);
subplot(5,1,2);
plot(cfs0);
title([‘結(jié)點(diǎn) 7’,’ 系數(shù)’])
subplot(5,1,3);
plot(cfs1);
title([‘結(jié)點(diǎn) 8’,’ 系數(shù)’])
subplot(5,1,4);
plot(cfs2);
title([‘結(jié)點(diǎn) 9’,’ 系數(shù)’])
subplot(5,1,5);
plot(cfs3);
title([‘結(jié)點(diǎn) 10’,’ 系數(shù)’])
%% reconstruct wavelet packet coefficients.
rex0=wprcoef(wpt,7);
rex1=wprcoef(wpt,8);
rex2=wprcoef(wpt,9);
rex3=wprcoef(wpt,10);
figure;
subplot(5,1,1);
plot(x);
title(‘原始信號(hào)’);
subplot(5,1,2);
plot(rex0);
title([‘重構(gòu)結(jié)點(diǎn)7’,’ 系數(shù)’])
subplot(5,1,3);
plot(rex1);
title([‘重構(gòu)結(jié)點(diǎn) 8’,’ 系數(shù)’])
subplot(5,1,4);
plot(rex2);
title([‘重構(gòu)結(jié)點(diǎn) 9’,’ 系數(shù)’])
subplot(5,1,5);
plot(rex3);
title([‘重構(gòu)結(jié)點(diǎn)10’,’ 系數(shù)’])

總結(jié)

以上是生活随笔為你收集整理的小波包分解、重构 matlab代码的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。