%---------------------------------------------------------------------------------------
%  exa030402.m,  for fig 3.4.2
%  說明在計算WVD時，加窗和不加窗的影響;
%
%  注：在該程序中，用到了子程序 tfrwv.m,
%      還要用到顯示程序 tfrqview及其他程序，
%      所有這些程序請讀者在如下的網站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因為該網站上的時頻分析軟件并沒有列入MATLAB的工具箱，涉及到知識產權問題，
%      因此不能將其列入國內公開出版物上，但讀者可以自由下載。
%----------------------------------------------------------------------------------------
clear;

n=256;
t=1:n;
   w=pi/8;
   for k=1:n/4,
      signal(k)=exp(j*1*k*w);
   end;
   for k=n/4:n/2,
     signal(k)=exp(j*2.5*k*w);
   end;
   for k=n/2:n,
      signal(k)=exp(j*3.5*k*w);
   end;

tfrwv(signal');