时域采样与频域采样

实验二：时域采样与频域采样

一、实验目的：

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

二、实验原理与方法：

1、时域采样定理的要点：

1）对模拟信号)(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 2）采样频率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

在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个

采样理论的要点。

对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两

个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。因此放在一起进行实验。

三、实验内容及步骤：

1、时域采样理论的验证：

给定模拟信号，)()sin()(0t u t Ae t x t a Ω=-α

式中A =444.128，α=502π，0Ω=502πrad/s ，它的幅频特性曲线如图2.1

图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 ==

分别画出()j X e ω、)()(1632k X k X 和的幅度谱，并绘图显示x (n)、)()(1632n x n x 和的波形，进行对比和分析，验证总结频域采样理论。

提示：频域采样用以下方法容易变程序实现。

① 直接调用MATLAB 函数fft 计算3232()FFT[()]X k x n =就得到()j X e ω在

]2,0[π的32点频率域采样

② 抽取32()X k 的偶数点即可得到()j X e ω在]2,0[π的16点频率域采样

16()X k ，即1632()(2) , 0,1,2,

,15X k X k k ==。

○

3 当然也可以按照频域采样理论，先将信号x(n)以16为周期进行周期延拓，取其主值区（16点），再对其进行16点DFT(FFT),得到的就是

()j X e ω

在]2,0[π的16点频率域采样

16()

X k 。