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

N=64;
alpha=0.05;f0=0.25a=alpha/pi;a=sqrt(a);a=sqrt(a);
t1=((1:N)-N/2-1)';y = exp(-t1.^2 * alpha/2);
y=y.*a;

y1=exp(-i*2*pi*t1*f0);
y=y.*y1;

subplot(221)
plot((-N/2:N/2-1),real(y));grid on;
xlabel(' Gauss modulated signal x(t)');

dsp=fftshift(abs(fft(y)).^2);
subplot(222)
plot((0:N/2-1)./N,dsp(1:N/2));grid on;
xlabel(' the Spectrum of x(t)');

[t0,T]=loctime(y)
t0=t0-N/2-1
[f0,B]=locfreq(y)
T*B