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

echo off
n=32;f0=0.25;
   sig1=1-abs(1-2*(0:n-1)./n);
   y=exp(j*2*pi*(0:n-1).*f0);
   sig1=y.*sig1;

subplot(221)
plot((0:n-1),real(sig1));grid
xlabel('Time');
ylabel('x(t)'); 

subplot(222)
dsp1=fftshift(abs(fft(sig1)).^2); 
plot((-n/2:n/2-1)/n,dsp1); grid;       
xlabel('Normalized frequency');
ylabel('Spectrum of x(t)');

[t0,T]=loctime(sig1')
[ff0,B]=locfreq(sig1')
qx=B/ff0
TBx=T*B

n=2*n;
   sig1=1-abs(1-2*(0:n-1)./n);
   y=exp(j*2*pi*(0:n-1).*f0/2);
   sig1=y.*sig1;

subplot(223)
plot((0:n-1),real(sig1));
xlabel('Time');
ylabel('x(t/2)'); grid

subplot(224)
dsp1=fftshift(abs(fft(sig1)).^2); 
plot((-n/2:n/2-1)/n,dsp1); grid;       
xlabel('Normalized frequency'); 
ylabel('Spectrum of x(t/2)');
[t0,T]=loctime(sig1')
[ff0,B]=locfreq(sig1')
qs=B/ff0
TBs=T*B