%---------------------------------------------------------------------------------------
% exa010303.m,  for  example 1.3.3
%      說明 高斯調(diào)制Chirp信號的時間中心、頻率最小、時寬和帶寬；
%  注：在該程序中，用到了子程序 fmlin,amgauss,loctime,locfreq.m, 
%      這些程序請讀者在如下的網(wǎng)站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因為該網(wǎng)站上的時頻分析軟件并沒有列入MATLAB的工具箱，涉及到知識產(chǎn)權(quán)問題，
%      因此不能將其列入國內(nèi)公開出版物上，但讀者可以自由下載。
%----------------------------------------------------------------------------------------
clear;

N=256;
sig1=fmlin(N).*amgauss(N);
subplot(221)
plot(real(sig1));
xlabel('Time'); 
ylabel('Real part of x(t)');
title('Linear frequency modulation'); grid

dsp1=fftshift(abs(fft(sig1)).^2); 
subplot(222)
plot((-N/2:N/2-1)/N,dsp1); grid;       

xlabel('Normalized frequency'); 
ylabel('the Spectrum of x(t)');

[t0,T]=loctime(sig1)
[f0,B]=locfreq(sig1)
T*B