信号频谱分析01

一、实验要求

利用FFT分析连续正弦信号的频谱：改变时间截断长度研究其对频谱泄漏的影响，并研究采用不同窗函数对信号进行时域截断时，对频谱泄漏的影响。

二、技术方案

首先绘制出连续正弦信号x=sin(2*pi*25*t)，由理论学习可知：时间截断长度为周期整数倍时不会产生频谱泄露，而时间截断长度非周期整数倍时则会产生频谱泄露，因此设定两个时间截断长度，由于该信号周期为0.04，一个为周期整数倍0.12，另一个为非周期整数倍0.14，观察比较频谱差别。此外，分别引入三角窗bartlett，矩形窗rectwin，观察频谱在加窗之后的变化。

三、实验程序

clc;
clear;
close all;

fs=1000;
Ts=1/fs;
%t=0:Ts:(0.12-Ts);
t=0:Ts:(0.14-Ts);
x=sin(2*pi*25*t);
N=length(t);
y=fft(x)/N;
z1=abs(y);

if mod(N,2)==0
    k=-N/2:(N/2-1);
else
    k=-(N-1)/2:(N-1)/2;
end;
f=k*fs/N;

subplot(4,1,1),plot(t,x,'r'),title('正弦信号');
subplot(4,1,2),stem(f,fftshift(z1),'g'),title('幅度谱');

w=bartlett(N);
x1=x.*w';
y1=fft(x1)/N;
z3=abs(y1);
subplot(4,1,3),stem(f,fftshift(z3)),title('加三角窗后幅度谱');

w1=rectwin(N);
x2=x.*w1';
y2=fft(x2)/N;
z4=abs(y2);
subplot(4,1,4),stem(f,fftshift(z4),'c'),title('加矩形窗后幅度谱');

四、实验结果

图中分别为非整周期的正弦信号、非整周期正弦信号的幅度谱、非整周期正弦信号加三角窗后的幅度谱、非整周期正弦信号加矩形窗后的幅度谱。时间截断长度为周期整数倍时不会产生频谱泄露，而时间截断长度非周期整数倍时则会产生频谱泄露，由第二张图所示。通过第二张图与第四张图对比可知，加矩形窗相当于不加窗。由第三张图可知，三角窗旁瓣小，主瓣宽度大，防泄漏效果较差。