%---------------------------------------------------------------------------------------
% exa010101.m,  for  example 1.1.1
% 說明三個分段相連的信號的傅里葉變換和時頻分析；
%  注：在該程序中，用到了子程序 tfstft.m, 還要用到顯示程序 tfrqview及其他程序，
%      所有這些程序請讀者在如下的網(wǎng)站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因?yàn)樵摼W(wǎng)站上的時頻分析軟件并沒有列入MATLAB的工具箱，涉及到知識產(chǎn)權(quán)問題，
%      因此不能將其列入國內(nèi)公開出版物上，但讀者可以自由下載。
%----------------------------------------------------------------------------------------
echo off
clear
n=256;
t=1:n;
   w=pi/5;
   for k=1:85,
   %    signal(k)=exp(j*k*w);
        signal(k)=sin(k*w);
   end;
   for k=86:170,
   %    signal(k)=exp(j*2*k*w);
        signal(k)=sin(2*k*w);

   end;
   for k=171:n,
   %    signal(k)=exp(j*3*k*w);
         signal(k)=sin(3*k*w);

   end;

tfrstft(signal');grid on;