Obtain the MODWPT of an ECG waveform and demonstrate perfect reconstruction using the inverse MODWPT.

loadwecg; 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'))

ans = 1.7903e-11

Obtain the MODWPT of Southern Oscillation Index data using the Daubechies extremal phase wavelet with two vanishing moments ('db2'). Reconstruct the signal using the inverse MODWPT.

loadsoi; wsoi = modwpt(soi,'db2'); xrec = imodwpt(wsoi,'db2');

Obtain the MODWPT of Southern Oscillation Index data using specified scaling and wavelets filters with the Daubechies extremal phase wavelet with two vanishing moments ('db2').

loadsoi; [lo,hi] = wfilters('db2'); wpt = modwpt(soi,lo,hi); xrec = imodwpt(wpt,lo,hi);

Plot the original SOI waveform and the reconstructed waveform.

subplot(2,1,1) plot(soi) title('Original SOI Waveform'); subplot(2,1,2) plot(xrec) title('Reconstructed SOI Waveform')