目录
1.背景
关于频谱分析主要研究对应频率信号的幅度和相位,因此对信号进行FFT变换后,需要求分别画出其幅度谱和相位谱。同时我们验证一个理论:实信号的频谱是共轭对称的。同时上节提出:关于采样频率fs和采样深度如何选择?这个问题,看看会出现哪些不可思的现象。
2.程序及原理分析
clear all; clc; close all;
fs= 1000;
N=1000;
t=(0:N-1)/fs;
f1=50;f2=65.75;
x=cos(2*pi*f1*t+pi/4)+cos(2*pi*f2*t);
X=fft(x,N)/N;
Y=fftshift(X);
freq1=(-N/2:N/2-1)*fs/N;
Y1=abs(Y);
Theta=angle(Y);
%作图
subplot 211;
stem(freq1,Y1,'k');xlim([-150 150]);
xlabel('频率(Hz)');ylabel('幅值');
set(gcf,'color','w');
subplot 212;
stem(freq1,Theta,'k');xlim([-150 150]);
xlabel('频率(Hz)');ylabel('相位');
set(gcf,'color','w');
3.幅度和相位谱及栅栏效应
图 1 幅度谱和相位谱
对于信号x=cos(2*pi*50*t+pi/4)+cos(2*pi*65.75*t),包括两个单音信号分别为f1 频率为50Hz,初始相位为pi/4。 f2 频率65.75Hz,初始相位为0;对于频率f1的信号幅度都为0.5,相位为0.799(pi/4),呈现出共轭对称;对于频率f1的信号,我们发现频谱图上不存在准确的62.75Hz的信号,而是在其附近鼓起来很多信号,似乎距离越近幅度也越大,其实这也符合能量守恒定理。。最近的信号是66Hz,幅度都为0.4505,相位为-0.7864,虽然也呈现出共轭对称,但是明显小于理想的0.5,相位更不是理想的0。这就是比较典型的栅栏效应。
4.问题
问题:如何消除栅栏效应?对于信号的频率与FFT后的谱线不重合如何求解幅度和相位谱?