%產生信號
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)))