一、实验目的
时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。
二、实验原理与方法
时域采样定理的要点是：
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 Ω代替即可。 频域采样定理的要点是：
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 个零点；如果NM 短，因此。()N x n 与x(n)不相同。
在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。
对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。
三、 实验内容及步骤
（1）时域采样理论的验证
给定模拟信号，)()sin()(0t u t Ae t x t
a Ω=-α
式中A =444.128，α=502π，0Ω=502πrad/s ，它的幅频特性曲线如图10.2.1
图10.2.1 )(t x a 的幅频特性曲线
现用DFT(FFT)求该模拟信号的幅频特性，以验证时域采样理论。
安照)(t x a 的幅频特性曲线，选取三种采样频率，即s F =1k Hz ，300Hz ，200Hz 。观测时间选ms T p 50=。
为使用DFT ，首先用下面公式产生时域离散信号，对三种采样频率，采样序列按顺序用)(1n x ，)(2n x ，)(3n x 表示。
)()sin()()(0nT u nT Ae
nT x n x nT
a Ω==-α 因为采样频率不同，得到的)(1n x ，)(2n x ，)(3n x 的长度不同， 长度（点数）用公式s p F T N ⨯=计算。选FFT 的变换点数为M=64，序列长度不够64的尾部加零。
X (k )=FFT[x (n )] ， k =0,1,2,3,-----,M -1 式中k 代表的频率为 k M
k π
ω2=
。 要求： 编写实验程序，计算)(1n x 、)(2n x 和)(3n x 的幅度特性，并绘图显示。观察分析频谱混叠失真。 （2）频域采样理论的验证 给定信号如下：
⎪⎩
⎪
⎨⎧≤≤-≤≤+=其它02614271301)(n n n n n x
编写程序分别对频谱函数()FT[()]j X e x n ω
=在区间]2,0[π上等间隔采样32 和16点，得到)()(1632k X k X 和： 32232
()()
, 0,1,2,31j k X k X e k ωπ
ω=
==
16216
()()
, 0,1,2,15j k X k X e k ω
πω=
==
再分别对)()(1632k X k X 和进行32点和16点IFFT ，得到)()(1632n x n x 和：
323232()IFFT[()] , 0,1,2,,31x n X k n == 161616()IFFT[()] , 0,1,2,
,15x n X k n ==
