%產(chǎn)生信號(hào)
N=128;
sig=fmlin(N);
%計(jì)算窄帶模糊函數(shù)
amb=ambifunb(sig);
amb=amb([N/2+1:N 1:N/2],:);
%模糊函數(shù)反變換
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));
%計(jì)算Wigner-Ville分布
[wvd2,t,f]=tfrwv(sig);
%比較
diff=max(max(abs(wvd1-wvd2)))