时域采样理论的验证.docVIP

时域采样理论的验证

一、实验目的

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

二、实验原理及方法

时域采样定理的要点是：

(a)对模拟信号以间隔T进行时域等间隔理想采样，形成的采样信号的频谱是原模拟信号频谱以采样角频率〔〕为周期进行周期延拓。公式为：

〔b〕采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的

频谱不产生频谱混叠。

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

理想采样信号和模拟信号之间的关系为：

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

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

上式中，在数值上＝，再将代入，得到：

上式的右边就是序列的傅立叶变换，即

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

三、实验内容及步骤

给定模拟信号，

式中A=444.128，=50π，=50πrad/s，它的幅频特性曲线如图

图1的幅频特性曲线

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

安照的幅频特性曲线，选取三种采样频率，即=1kHz，300Hz，200Hz。观测时间选。为使用DFT，首先用下面公式产生时域离散信号，对三种采样频率，采样序列按顺序用，，表示。

因为采样频率不同，得到的，，的长度不同，长度〔点数〕用公式计算。选FFT的变换点数为M=64，序列长度不够64的尾部加零。

X(k)=FFT[x(n)]，k=0,1,2,3,-----,M-1式中k代表的频率为。

四、实验要求：

编写实验程序，计算、和的幅度特性，并绘图显示。观察分析频谱混叠失真。

五、程序清单和信号波形

程序清单：

Tp=64/1000;

Fs=1000;T=1/Fs;

M=Tp*Fs;n=0:M-1;

f=n*Fs/M;

A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;

xn=A*exp(-alph*n*T).*sin(omega*n*T);

Xk=T*fft(xn,M);

subplot(3,1,1);

plot(f,abs(Xk));

xlabel(f/Hz);

ylabel(|x1(jf)|);

title(x1(n)的幅度特性);

%================================================================

Tp=64/1000;

Fs=300;T=1/Fs;

M=Tp*Fs;n=0:M-1;

f=n*Fs/M;

A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;

xn=A*exp(-alph*n*T).*sin(omega*n*T);

Xk=T*fft(xn,M);

subplot(3,1,2);

plot(f,abs(Xk));

xlabel(f/Hz);

ylabel(|x2(jf)|);

title(x2(n)的幅度特性);

%================================================================

Tp=64/1000;

Fs=200;T=1/Fs;

M=Tp*Fs;n=0:M-1;

f=n*Fs/M;

A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;

xn=A*exp(-alph*n*T).*sin(omega*n*T);

Xk=T*fft(xn,M);

subplot(3,1,3);

plot(f,abs(Xk));

xlabel(f/Hz);

ylabel(|x3(jf)|);

title(x3(n)的幅度特性);

信号波形：

六、实验总结

由图可见，采样序列的频谱确实是以采样频率为周期对模拟信号频谱的周期延拓。当采样频率为1000Hz时频谱混叠很小；当采样频率为300Hz时，在折叠频率150Hz附近频谱混叠很严重；当采样频率为200Hz时，在折叠频率110Hz附近