说明:
(1)做图3.2的收获,附加图1;以高斯信号为例,来验证得到频谱幅相的两种求解过程
matlab程序代码如下:
% initial matlab workspace
% figure3.a1 ; 附加图1; 以高斯信号为例,来验证得到频谱幅相的两种求解过程
clc
clear
close all
% 时间序列
T = 150;
Nt = 150;
dt = T/Nt;
t = -T/2:dt:T/2-dt;
% 高斯信号
sigma = 2;
ft = exp( - t.^2./2./sigma.^2 );
% 得到抽样频率
fs = 1/dt;
fre = (0:Nt-1)*fs/Nt - fs/2;
% 得到理论频谱
fw = sqrt(2