%??This?program?is?free?software;?you?can?redistribute?it?and/or?modify
%??it?under?the?terms?of?the?GNU?General?Public?License?as?published?by
%??the?Free?Software?Foundation;?either?version?2?of?the?License?or
%??（at?your?option）?any?later?version.
%
%??This?program?is?distributed?in?the?hope?that?it?will?be?useful
%??but?WITHOUT?ANY?WARRANTY;?without?even?the?implied?warranty?of
%??MERCHANTABILITY?or?FITNESS?FOR?A?PARTICULAR?PURPOSE.??See?the
%??GNU?General?Public?License?for?more?details.
%
%??You?should?have?received?a?copy?of?the?GNU?General?Public?License
%??along?with?this?program;?if?not?write?to?the?Free?Software
%??Foundation?Inc.?51?Franklin?St?Fifth?Floor?Boston?MA??02110-1301??USA
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（N2））/2+1:N?1:（N+rem（N2））/2]:）8）;?
xlabel（‘Delay‘）;?ylabel（‘Doppler‘）;
title（‘Wigner-Ville?weighting?function‘）
set（gca‘yTickLabel‘[]）
set（gca‘xTickLabel‘[]）
hold?on
[ah]=contour（WF1[1/2]‘g‘）;
set（h‘linewidth‘2）;
hold?off

subplot（322）;
Max=max（max（tfr1））;
levels=linspace（Max/10MaxNcont）;
contour（tfr1levels）;
xlabel（‘Time‘）;?ylabel（‘Frequency‘）;
title（‘Wigner-Ville?distribution‘）
set（gca‘yTickLabel‘[]）
set（gca‘xTickLabel‘[]）

subplot（323）;
contour（dlr（[（N+rem（N2））/2+1:N?1:（N+rem（N2））/2]:）8）;?
xlabel（‘Delay‘）;?ylabel（‘Doppler‘）;
title（‘Spectrogram?weighting?function‘）;
set（gca‘yTickLabel‘[]）
set（gca‘xTickLabel‘[]）
hold?on
[ah]=contour（WF2[1/2]‘g‘）;
set（h‘linewidth‘2）;
hold?off

subplot（324）;
Max=max（max（tfr2））;
levels=linspace（Max/10MaxNcont）;
contour（tfr2（1:N/2:）levels）;
xlabel（‘Time‘）;?ylabel（‘Frequency‘）;
title（‘Spectrogram‘）
set（gca‘yTickLabel‘[]）
set（gca‘xTickLabel‘[]）

subplot（325）;
contour（dlr（[（N+rem（N2））/2+1:N?1:（N+rem（N2））/2]:）8）;?
xlabel（‘Delay‘）;?ylabel（‘Doppler‘）;
title（‘SP-WV?weighting?function‘）;
set（gca‘yTickLabel‘[]）
set（gca‘xTickLabel‘[]）
hold?on
[ah]=contour（WF3[1/2]‘g‘）;
set（h‘linewidth‘2）;
hold?off

subplot（326）;
Max=max（max（tfr3））;
levels=linspace（Max/10MaxNcont）;
contour（tfr3levels）;
xlabel（‘Time‘）;?ylabel（‘Frequency‘）;
title（‘Smoothed-pseudo-WVD‘）;
set（gca‘yTickLabel‘[]）
set（gca‘xTickLabel‘[]）


figure（2）;

subplot（221）;
contour（dlr（[（N+rem（N2））/2+1:N?1:（N+rem（N2））/2]:）8）;?
xlabel（‘Delay‘）;?ylabel（‘Doppler‘）;
title（‘Born-Jordan?weighting?function‘）;
set（gca‘yTickLabel‘[]）
set（gca‘xTickLabel‘[]）
hold?on
[ah]=contour（WF4[1/2]‘g‘）;
set（h‘linewidth‘2）;
hold?off

subplot（222）;
Max=max（max（tfr4））;
lev