说明:
(1)做图3.2的收获,附加图2;以高斯信号为例,来说明如何进行正确的频域插值-频谱的幅度和相位都与理论值一致
matlab程序代码如下:
% initial matlab workspace
% figure3.a2 ; 附加图2; 以高斯信号为例,来说明如何进行正确的频域插值
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 );
% 对频域插值,插值之后频域的抽样点数为Nt+Nz
Nz = 4*Nt;
N = Nt + Nz;
% 得到抽样频率
fs = 1/dt;
fre = (0:N-1)*fs/N - fs/2;
% 得到理论频谱
fw = sqrt(2*pi