%---------------------------------------------------------------------------------------
% exa010305.m,  for  example 1.3.5
%      說(shuō)明 幅度為三角波的幅度調(diào)制信號(hào)在長(zhǎng)度變化前后，信號(hào)的時(shí)間中心、頻率最小、時(shí)寬
%           和帶寬的變化情況；
%  注：在該程序中，用到了子程序 loctime,locfreq.m, 
%      這些程序請(qǐng)讀者在如下的網(wǎng)站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因?yàn)樵摼W(wǎng)站上的時(shí)頻分析軟件并沒(méi)有列入MATLAB的工具箱，涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，
%      因此不能將其列入國(guó)內(nèi)公開(kāi)出版物上，但讀者可以自由下載。
%----------------------------------------------------------------------------------------
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