%---------------------------------------------------------------------------------------
% exa010501a.m,  for  example 1.5.1 and fig 1.5.1
%      給出兩個chirp信號相加后，由瞬時頻率定義求出的瞬時頻率曲線，說明瞬時頻率的不確切性；
%  注：在該程序中，用到了子程序 loctime,locfreq.m, instfreq.m，tfrideal.m等，
%      這些程序請讀者在如下的網(wǎng)站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因為該網(wǎng)站上的時頻分析軟件并沒有列入MATLAB的工具箱，涉及到知識產(chǎn)權(quán)問題，
%      因此不能將其列入國內(nèi)公開出版物上，但讀者可以自由下載。
%----------------------------------------------------------------------------------------
clear;

N=128;
[x1,if1]=fmlin(N,0,0.3);
[x2,if2]=fmlin(N,0.2,0.5);
y=x1+x2;
t=1:N;
subplot(311)plot(t,real(y));grid on;
ylabel('real(x(t))')

subplot(312)
ifr=instfreq(y); 
plot(ifr'); grid;
axis([1 128 0 0.5]);
ylabel(' frequency'); 

subplot(313)
tfrideal([if1,if2]);