test0805gausscmod.m

来自「信号处理实例」· M 代码 · 共 28 行

%---------------------------------------------------------------------------------------
% test0805gausscmod.m,  for  example 08.05
% gauss signal and its spectrum
% 2006.12.15 LYB
%----------------------------------------------------------------------------------------
clear;

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

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

xlabel('Normalized frequency'); 
ylabel('the Spectrum of x(t)');
text(0.32,520,'test0805gausscmod');

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