时域采样与频域采样

备注：（1）、按照要求独立完成实验内容。

（2）、实验结束后，把电子版实验报告按要求格式改名（例：09号_张三_实验七.doc ）后，

实验室统一刻盘留档。

实验四 时域采样与频域采样

1、时域采样定理的要点 一、实验目的

时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。

二、实验原理

a) 对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样，形成的采样信号的频谱)(ˆΩj X

是原模拟信号频谱(

)a X j Ω以采样角频率s Ω（T s

/2π=Ω）为周期进行周期延拓。公

式为： )](ˆ[)(ˆt x

FT j X a a =Ω )(1∑∞

-∞

=Ω-Ω=n s a jn j X T b) 采样频率s Ω必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的

频谱不产生频谱混叠。

利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。

理想采样信号)(ˆt x

a 和模拟信号)(t x a 之间的关系为：∑∞

-∞

=-=n a a nT t t x t x )()()(ˆδ

对上式进行傅立叶变换，得到：

dt e nT t t x j X t j n a a Ω-∞∞

-∞

-∞

=⎰∑

-=Ω])()([)(ˆδ

dt e nT t t x t j n a Ω-∞

-∞

=∞

∞

-∑⎰

-)()( δ＝

在上式的积分号内只有当nT t =时，才有非零值，因此

∑∞

-∞

=Ω-=Ωn nT j a

a

e nT x

j X )()(ˆ

上式中，在数值上)(nT x a ＝)(n x ，再将T Ω=ω代入，得到：

∑∞

-∞

=-=Ωn n

j a

e

n x j X ω)()(ˆ

上式的右边就是序列的傅立叶变换)(ω

j e

X ，即

T j a e X j X Ω==Ωωω)()(ˆ

上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量ω用T Ω代替即可。

2、频域采样定理的要点

a) 对信号x(n)的频谱函数X(e j

ω

)在[0，2π]上等间隔采样N 点，得到

2()()

, 0,1,2,

,1j N k N

X k X e k N ωπω===-

则N 点IDFT[()N X k ]得到的序列就是原序列x(n)以N 为周期进行周期延拓后的主值区序列，公式为：

()IDFT[()][

()]()N N N N

i x n X k x n iN R

n ∞

=-∞

==+∑

b) 由上式可知，频域采样点数N 必须大于等于时域离散信号的长度M(即N≥M)，才能使

时域不产生混叠，则N 点IDFT[()N X k ]得到的序列()N x n 就是原序列x(n),即

()N x n =x(n)。如果N>M ，()N x n 比原序列尾部多N-M 个零点；如果N

短，因此。()N x n 与x(n)不相同。

在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。 对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。

三、实验内容（包括代码与产生的图形及结果分析）

1. 给定模拟信号如下：

xa(t)=Ae -αt sin(Ω0t)u(t)

式中, A=444.128，α=50 π， Ω0=50 π rad/s ，将这些参数带入上式中，对x a (t

进行傅里叶变换，它的幅频特性曲线如图1所示。

现用DFT(FFT)求该模拟信号的幅频特性，以验证时域采样理论。

按照xa(t)的幅频特性曲线，选取三种采样频率，即Fs=1 kHz ，300 Hz ，200 Hz 。观测时间选Tp=64 ms 。

要求： 编写实验程序，计算x 1(n)、 x 2(n)和x 3(n)

的幅度特性，并绘图显示。观察分析频谱混叠失真。

function tstem(xn,yn) %时域序列绘图函数

%xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名称（字符串） n=0:length(xn)-1; stem(n,xn,'.'); xlabel('n');ylabel('yn');

axis([0,n(end),min(xn),1.2*max(xn)]);

Tp=64/1000; %观察时间Tp=64微秒 %产生M 长采样序列x(n) % Fs=1000;T=1/Fs; Fs=1000;T=1/Fs; M=Tp*Fs;n=0:M-1;

A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5; xnt=A*exp(-alph*n*T).*sin(omega*n*T); Xk=T*fft(xnt,M); %M 点FFT[xnt)] yn='xa(nT)';subplot(3,2,1);

tstem(xnt,yn); %调用自编绘图函数tstem 绘制序列图

22图1 x a (t)的幅频特性曲线