N1 ,对连续时间信号进行采样,然后利用 DFT 计算出采样长度 N1 情况下的频谱,并给出时域和频域图像。一部分采用加窗的方法对时域图像进行采样,然后利用 DFT 计算出采样长度 N1 情况下的频谱,并给出时域和频域图像。这里要注意的是我们采用加窗的方法时, 设定信号的长度为 400 ,然后对长度为 400 的信号进行截短加窗,这时只有加窗处我们是对信号采样的,即得到的信号是有效的,其余位置信号都为 0,所以这时我们仍然认为采样长度为 128. (2) 实验代码: N1=128; n=0:N1-1; k=0:N1-1; xn=sin(5*2*pi*n/20)+cos(3*2*pi*n/20); WNnk=dftmtx(N1); Xk=xn*WNnk; subplot(2,2,1); stem(n,xn); m=(0:N1/2-1)/N1; subplot(2,2,2); stem(m,abs(Xk(1:N1/2))) n1=1:400; %设定信号长度为 400 xn1=(sin(n1*5*2*pi/20)+cos(n1*3*2*pi/20)).*(heaviside (n1)-heaviside(n1-N1+1)); %对数据进行截短加窗 WNnk=dftmtx(400); Xk1=xn1*WNnk; subplot(2,2,3); stem(n1,xn1); m=(0:400/2-1)/400; subplot(2,2,4); stem(m,abs(Xk1(1:400/2))); (3) 实验结果: 图3.1 为实验得到的总体图。图 3.2 为采样长度为 32 的时域图像,图 3.3 为对应的频域图像。图 2.4 为加窗后的时域图像,图 2.5 为加窗后的对应的频域图像。得到的ms =0.0078 ,ms1 =0.0078 。图3.1 图 3.2 图3.3 图3.4
全文预览
