f5_9.m

matlab时频分析配套源码第3章的程序

%产生信号
N=128;
sig=fmlin(N);
%计算窄带模糊函数
amb=ambifunb(sig);
amb=amb([N/2+1:N 1:N/2],:);
%模糊函数反变换
ambi=ifft(amb).';
tdr=zeros(N);
tdr(1:N/2,:)=ambi(N/2:N-1,:);
tdr(N:-1:N/2+2,:)=ambi(N/2-1:-1:1,:);
wvd1=real(fft(tdr));
%计算Wigner-Ville分布
[wvd2,t,f]=tfrwv(sig);
%比较
diff=max(max(abs(wvd1-wvd2)))