www.gusucode.com > signal 案例源码程序 matlab代码 > signal/InverseFourierSynchrosqueezedTransformOfSpeechSignalExample.m

    %% Inverse Fourier Synchrosqueezed Transform of Speech Signal
% Load a speech signal sampled at $F_s=7418\;\rm Hz$. The file contains a
% recording of a female voice saying the word "MATLAB®." Compute the
% Fourier synchrosqueezed transform of the signal.

%%

load mtlb         % To hear, type sound(mtlb,Fs)

[sst,f] = fsst(mtlb,Fs);

%%
% Invert the transform to reconstruct the signal. Plot the original and
% reconstructed signals, as well as the difference between them.

xrec = ifsst(sst);
 
t = (0:length(mtlb)-1)/Fs;
plot(t,mtlb,t,xrec,t,mtlb-xrec)

xlabel('Time (s)')
legend('Original','Reconstructed','Difference')

%%
% Check the accuracy of the reconstruction by computing the $\ell_\infty$
% norm of the difference between the original signal and the inverse
% transform.

Linf = norm(abs(mtlb-xrec),Inf)

% To hear, type sound(mtlb-xrec,Fs)