%---------------------------------------------------------------------------------------
%  exa030301.m,  for  example 3.3.1 and fig 3.3.1
%  說明一矩形窗的WVD;
%
%  注：在該程序中，用到了子程序 tfrpwv.m,
%      還要用到顯示程序 tfrqview及其他程序，
%      所有這些程序請(qǐng)讀者在如下的網(wǎng)站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因?yàn)樵摼W(wǎng)站上的時(shí)頻分析軟件并沒有列入MATLAB的工具箱，涉及到知識(shí)產(chǎn)權(quán)問題，
%      因此不能將其列入國內(nèi)公開出版物上，但讀者可以自由下載。
%----------------------------------------------------------------------------------------
clear;

n=64;
sig=zeros(1,n);
   for l=n/4:3*n/4
     sig(l)=1;
   end;

tfrpwv(sig');
% 注：應(yīng)在程序tfrpwv 倒數(shù)第6行處加一句： tfr=change(tfr),則可生成如圖3.3.1的圖形，
%      否則，其Wigner分布被分成兩半。