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

    %% Synchrosqueezed and Inverse Synchrosqueezed Transform of Speech Signal
% Obtain the wavelet synchrosqueezed transform and inverse synchrosqueezed 
% transform of a speech sample using the bump wavelet. 
%%
% Load the speech signal and obtain the synchrosqueezed transform and 
% inverse synchrosqueezed transform.
load mtlb;
dt = 1/Fs;
t = 0:dt:numel(mtlb)*dt-dt;
Txmtlb = wsst(mtlb,'bump');
xrec = iwsst(Txmtlb,'bump');
%% 
% Obtain the L-infinity norm of the difference between the original 
% waveform and the reconstruction. Plot the results.    
Linf = norm(abs(mtlb-xrec),Inf);
plot(t,mtlb);
hold on;
xlabel('Seconds'); ylabel('Amplitude');
plot(t,xrec,'r');
title({'Reconstruction of Wavelet Synchrosqueezed Transform'; ...
['Largest Absolute Difference ' num2str(Linf,'%1.2f')]});