www.gusucode.com > wavelet 源码程序 matlab案例代码 > wavelet/PerfectReconstructionWithTheInverseMODWPTExample.m

    %% Perfect Reconstruction with the Inverse MODWPT
% Obtain the MODWPT of an ECG waveform and demonstrate perfect
% reconstruction using the inverse MODWPT.
load wecg;
wpt = modwpt(wecg);
xrec = imodwpt(wpt);
subplot(2,1,1)
plot(wecg);
title('Original ECG Waveform');
subplot(2,1,2)
plot(xrec);
title('Reconstructed ECG Waveform');
%%
% Find the largest absolute difference between the original signal and the 
% reconstruction. The difference is on the order of $10^{-11}$, which 
% demonstrates perfect reconstruction.
max(abs(wecg-xrec'))