paramfun.m

时间序列的时频分析工具箱

echo off
%PARAMFUN Figures representing different TFD of the Cohen's class.
%	On the left, the ambiguity plane and the weighting functions ;
%	On the right, the time-frequency distributions.

%	O. Lemoine - February, July 1996. 


if ~exist('paramfun.mat'),
 echo on
%
% PARAMFUN does not exist.
%
% You haven't created "paramfun.mat". Please run ??? to proceed.
% The next part of the demo is skipped.
%
else
 load paramfun
 Ncont=5;
 
 subplot(321);
 contour(dlr([(N+rem(N,2))/2+1:N 1:(N+rem(N,2))/2],:),8); 
 xlabel('Delay'); ylabel('Doppler');
 title('Wigner-Ville weighting function')
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 hold on
 [a,h]=contour(WF1,[1/2],'g');
 set(h,'linewidth',2);
 hold off
 
 subplot(322);
 Max=max(max(tfr1));
 levels=linspace(Max/10,Max,Ncont);
 contour(tfr1,levels);
 xlabel('Time'); ylabel('Frequency');
 title('Wigner-Ville distribution')
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 
 subplot(323);
 contour(dlr([(N+rem(N,2))/2+1:N 1:(N+rem(N,2))/2],:),8); 
 xlabel('Delay'); ylabel('Doppler');
 title('Spectrogram weighting function');
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 hold on
 [a,h]=contour(WF2,[1/2],'g');
 set(h,'linewidth',2);
 hold off
 
 subplot(324);
 Max=max(max(tfr2));
 levels=linspace(Max/10,Max,Ncont);
 contour(tfr2(1:N/2,:),levels);
 xlabel('Time'); ylabel('Frequency');
 title('Spectrogram')
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 
 subplot(325);
 contour(dlr([(N+rem(N,2))/2+1:N 1:(N+rem(N,2))/2],:),8); 
 xlabel('Delay'); ylabel('Doppler');
 title('SP-WV weighting function');
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 hold on
 [a,h]=contour(WF3,[1/2],'g');
 set(h,'linewidth',2);
 hold off
 
 subplot(326);
 Max=max(max(tfr3));
 levels=linspace(Max/10,Max,Ncont);
 contour(tfr3,levels);
 xlabel('Time'); ylabel('Frequency');
 title('Smoothed-pseudo-WVD');
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 
 
 figure(2);
 
 subplot(221);
 contour(dlr([(N+rem(N,2))/2+1:N 1:(N+rem(N,2))/2],:),8); 
 xlabel('Delay'); ylabel('Doppler');
 title('Born-Jordan weighting function');
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 hold on
 [a,h]=contour(WF4,[1/2],'g');
 set(h,'linewidth',2);
 hold off
 
 subplot(222);
 Max=max(max(tfr4));
 levels=linspace(Max/10,Max,Ncont);
 contour(tfr4,levels);
 xlabel('Time'); ylabel('Frequency');
 title('Born-Jordan distribution');
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 
 
 subplot(223);
 contour(dlr([(N+rem(N,2))/2+1:N 1:(N+rem(N,2))/2],:),8); 
 xlabel('Delay'); ylabel('Doppler');
 title('CW weighting function');
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 hold on
 [a,h]=contour(WF5,[1/2],'g');
 set(h,'linewidth',2);
 hold off
 
 subplot(224);
 Max=max(max(tfr5));
 levels=linspace(Max/10,Max,Ncont);
 contour(tfr5,levels);
 xlabel('Time'); ylabel('Frequency');
 title('Choi-Williams distribution');
 set(gca,'yticklabels',[])
 set(gca,'xticklabels',[])
 
 echo on
end;