[重点]对正弦信号的采样频谱分析.pdf

一、题目要求：

给定采样频率fs，两个正弦信号相加，两信号幅度不同、频率不同。要求

给定正弦信号频率的选择与采样频率成整数关系和非整数关系两种情况，信号持

续时间选择多种情况分别进行频谱分析。

二、题目原理与分析：

本题目要对正弦信号进行抽样，并使用fft对采样信号进行频谱分析。因此

首先对连续正弦信号进行离散处理。实际操作中通过对连续信号间隔相同的抽样

周期取值来达到离散化的目的。根据抽样定理，如果信号带宽小于奈奎斯特频率

（即采样频率的二分之一），那么此时这些离散的采样点能够完全表示原信号。

高于或处于奈奎斯特频率的频率分量会导致混叠现象。设抽样周期为TS(抽样角

频率为ωS)，则

可见抽样后的频谱是原信号频谱的周期性重复，当信号带宽小于奈奎斯特频率的

二分之一时不会产生频谱混叠现象。

因此，我们对采样频率的选择采取fs2fo,fs=2fo,fs2fo三种情况进行分析。

对信号采样后，使用fft函数对其进行频谱分析。为了使频谱图像更加清楚，更

能准确反映实际情况并接近理想情况，我们采用512点fft。取512点fft不仅可

以加快计算速度，而且可以使频谱图更加精确。若取的点数较少，则会造成频谱

较大的失真。

三、实验程序：

本实验采用matlab编写程序，实验中取原信号为

ft=sin(2πfXt)+2sin(10πfXt)，取频率f=1kHz，实验程序如下：

f=1000;fs=20000;Um=1;

N=512;T=1/fs;

t=0:1/fs:0.01;

ft=Um*sin(2*pi*f*t)+2*Um*sin(10*pi*f*t);

subplot(3,1,1);

plot(t,ft);gridon;

axis([00.011.1*min(ft)1.1*max(ft)]);

xlabel(t),ylabel(ft);

title(抽样信号的连续形式);

subplot(3,1,2);

stem(t,ft);gridon;

axis([00.011.1*min(ft)1.1*max(ft)]);

xlabel(t),ylabel(ft);

title(实际抽样信号);

k=0:N-1;

Fw=fft(ft,N);

subplot(3,1,3);

plot(k,abs(Fw));gridon;

axis([0550-0.265*pi]);

title(抽样信号幅度谱)

在实际操作过程中，对于信号频率与采样频率所成整数倍与非整数倍关系时，信号持续

时间不同时，只需改变程序中的相关语句即可。既t=0:1/fs:to；语句控制信号持续时

间，改变to即可。改变抽样频率只需对fs取不同的值即可。

四、实验过程及图示：

1.信号持续时间为0.01s，信号频率与采样频率成整数关系：

（1）fs2fo，取fs=20kHz，得到频谱图：

（2）fs=2fo，取fs=10kHz，得到频谱图：

（3）fs2fo，取fs=5kHz，得到频谱图：

通过比较三个图形发现当抽样信号频率大于原信号频率的二倍时抽样信号

能较好的反应原信号，并且抽样信号频谱呈现两个峰值，与正弦信号的理想频谱

既冲击函数较为接近。但是由于实际信号的持续时间是有限的，因此频谱不可能

完全表现为冲击函数的情况，会有尾部延伸。当抽样频率等于原信号频率的二倍

时，抽样信号只能表现为单个正弦信号的形式，因此频谱只能表现为单峰情况，

且幅度也较前者有较大的下降。当抽样信号频率小于原信号频率的两倍时，抽样

信号波形有较大的失真，且幅度有更大的下降，频谱的尾部所占比例更大，失真

较为严重。

2.持续时间为0.01s，信号频率与采样频率成非整数关系：

（1）fs2fo,取fs为16.5kHz，得到频谱为：

（2）fs=2fo的情况同1，省略。

（3）fs2fo，取fs为2.5kHz，得到频谱为：

通过观察频谱图发现，对抽样频率取三种情况时频谱的规律与成整数关系时

的规律基本相同，但是纵向比较时，抽样信号的波形与原信号波形有较大的失真，

这是由于抽样信号的频率不为原信号的整数倍造成的，反应到频率谱上，导致出

现的峰值下降，较为弱的趋向理想冲击函数。