信号频谱分析.pdfVIP

信号频谱分析

摘要信号频谱分析就是利用MATLAB实现快速傅里叶变换的分析的方法，通

：

过DSP数字信号处理器对实时采集到的信号进行FFT运算后实现时域与频域的

转换，频谱分析主要分析信号是由哪些频率的正弦信号叠加得到的，以及这些正

弦信号的振幅，所以FFT运算结果反映的是频域中各频率分量幅值的大小。

关键词MATLAB、傅里叶变换、频谱分析

：。

一、引言

信号的频谱分析的主要是求出时域描述相对应的频域描述，从中找出

信号频谱的变化规律，已达到特征提取的目的。不同信号的傅里叶分

析理论与方法有差别，目前常用的是FFT分析法，在特定的时间段

中队时域信号进行FFT变换，得到频域信息并获取相对于频率的幅

度、相位信息，适用于非周期信号和持续时间很短的瞬态信号的频谱

测量。

二、实验设计思想

先将实时信号的采样值送入DSP系统，DSP程序对这些采样值进行

FFT变换，经运算求出对应的信号频谱数据，并将结果送到PC机屏

幕上进行显示，使DSP硬件系统完成一台信号频谱分析仪的功能，

如下流程所示：

在A/D端加入输入信号→装在并运行FFT程序→读入AD采样值

→进行FFT运算求信号的DFT→将计算结果送入PC机→显示时域信

号及频谱图

1、实验原理说明

一个模拟信号，经过ADC采样之后，就变成了数字信号。采样定理

告诉我们，采样频率要大于信号频率的两倍。采样得到的数字信号，

然后做FFT变换。N个采样点，经过FFT之后，就可以得到N个点

的FFT结果。为了方便进行FFT运算，通常N取2的整数次方。假

设采样频率为Fs，信号频率F，采样点数为N。那么FFT之后结果

就是一个为N点的复数。每一个点就对应着一个频率点。这个点的

模值，就是该频率值下的幅度特性。若原始信号的峰值为A，那么

FFT的结果的每个点（除了第一个点直流分量之外）的模值就是A

的N/2倍。而第一个点就是直流分量，它的模值就是直流分量的N

倍。而每个点的相位就是在该频率下的信号的相位。第一个点表示直

流分量（即0Hz），而最后一个点N的再下一个点（实际上这个点是

不存在的，这里是假设的第N+1个点，也可以看做是将第一个点分

做两半分，另一半移到最后）则表示采样频率Fs，这中间被N-1个

点平均分成N等份，每个点的频率依次增加。例如某点n所表示的

频率为：Fn=(n-1)*Fs/N。由上面的公式可以看出，Fn所能分辨到频

率为Fs/N，如果采样频率Fs为1024Hz，采样点数为1024点，则可

以分辨到1Hz。1024Hz的采样率采样1024点，刚好是1秒，也就是

说，采样1秒时间的信号并做FFT，则结果可以分析到1Hz，如果采

样2秒时间的信号并做FFT，则结果可以分析到0.5Hz。如果要提高

频率分辨力，则必须增加采样点数，也即采样时间。频率分辨率和采

样时间是倒数关系。

假设FFT之后某点n用复数a+bi表示，那么这个复数的模就是

An=a*a+b*b的二次开方，相位就是Pn=atan2(b,a)。根据以上的结果，

就可以计算出n点（n≠1，且n=N/2）对应的信号的表达式为：

An/(N/2)*cos(2*pi*Fn*t+Pn)，即2*An/N*cos(2*pi*Fn*t+Pn)。对于n=1

点的信号，是直流分量，幅度即为A1/N。由于FFT结果的对称性，

通常我们只使用前半部分的结果，即小于采样频率一半的结果。

2、实验过程

（1）、先运行上位软件Signal，进入频谱分析窗口。

（2）、分析函数发生器产生的模拟信号，实验时先将函数发生器的输

出端接到数字滤波器的输入端，然后运行FFT程序，最后做实时分

析，即可在PC机屏幕上观察到输入信号的波形和快速傅里叶变换的

结果（可通过选择不同波形进行观察