实验4时域采样理论与频域采样定理验证文档格式.doc
《实验4时域采样理论与频域采样定理验证文档格式.doc》由会员分享，可在线阅读，更多相关《实验4时域采样理论与频域采样定理验证文档格式.doc（9页珍藏版）》请在冰点文库上搜索。
（b）采样频率必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的
频谱不产生频谱混叠。
利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。
理想采样信号和模拟信号之间的关系为：
对上式进行傅立叶变换，得到：
在上式的积分号内只有当时，才有非零值，因此：
上式中，在数值上＝，再将代入，得到：
上式的右边就是序列的傅立叶变换，即
上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量ω用代替即可。
频域采样定理的要点是：
a）对信号x（n）的频谱函数X（ejω）在[0，2π]上等间隔采样N点，得到
则N点IDFT[]得到的序列就是原序列x（n）以N为周期进行周期延拓后的主值区序列，公式为：
（b）由上式可知，频域采样点数N必须大于等于时域离散信号的长度M（即N≥M），才能使时域不产生混叠，则N点IDFT[]得到的序列就是原序列x（n）,即=x（n）。
如果N>
M，比原序列尾部多N-M个零点；
如果N<
M，z则=IDFT[]发生了时域混叠失真，而且的长度N也比x（n）的长度M短，因此。
与x（n）不相同。
在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。
对比上面叙述的时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：
“时域采样频谱周期延拓，频域采样时域信号周期延拓”。
因此放在一起进行实验。
三、实验内容及步骤
（1）时域采样理论的验证。
给定模拟信号，
式中A=444.128，=50π，=50πrad/s，它的幅频特性曲线如图10.2.1
图10.2.1的幅频特性曲线
现用DFT（FFT）求该模拟信号的幅频特性，以验证时域采样理论。
安照的幅频特性曲线，选取三种采样频率，即=1kHz，300Hz，200Hz。
观测时间选。
为使用DFT，首先用下面公式产生时域离散信号，对三种采样频率，采样序列按顺序用，，表示。
因为采样频率不同，得到的，，的长度不同，长度（点数）用公式计算。
选FFT的变换点数为M=64，序列长度不够64的尾部加零。
X（k）=FFT[x（n）]，k=0,1,2,3,-----,M-1式中k代表的频率为。
要求：
编写实验程序，计算、和的幅度特性，并绘图显示。
观察分析频谱混叠失真。
（2）频域采样理论的验证。
给定信号如下：
编写程序分别对频谱函数在区间上等间隔采样32
和16点，得到：
再分别对进行32点和16点IFFT，得到：
分别画出、的幅度谱，并绘图显示x（n）、的波形，进行对比和分析，验证总结频域采样理论。
提示：
频域采样用以下方法容易变程序实现。
①直接调用MATLAB函数fft计算就得到在的32点频率域采样
②抽取的偶数点即可得到在的16点频率域采样，即。
当然也可以按照频域采样理论，先将信号x（n）以16为周期进行周期延拓，取其主值区（16点），再对其进行16点DFT（FFT）,得到的就是在的16点频率域采样。
四、思考题
如果序列x（n）的长度为M，希望得到其频谱在上的N点等间隔采样，当N<
M时，如何用一次最少点数的DFT得到该频谱采样？
五、实验报告及要求
（1）运行程序,打印要求显示的图形。
（2）分析比较实验结果，简述由实验得到的主要结论。
（3）简要回答思考题。
（4）附上程序清单和有关曲线。
六、程序清单和信号波形
1、时域采样理论的验证
程序清单：
%时域采样理论验证程序
Tp=64/1000;
%观察时间Tp=64微秒
%产生M长采样序列x（n）
%Fs=1000;
T=1/Fs;
Fs=1000;
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）;
%M点FFT[xnt）]
subplot（3,1,1）;
plot（f,abs（Xk））;
xlabel（'
f/Hz'
）;
ylabel（'
|x1（jf）|'
title（'
x1（n）的幅度特性'
%====================================================================
%Fs=300Hz
Fs=300;
%Fs=200Hz
Fs=200;
xnt=A*exp（-alph*n*T）.*sin（omega*n*T）;
Xk=T*fft（xnt,M）;
yn='
xa（nT）'
;
subplot（3,2,5）;
tstem（xnt,yn）;
boxon;
title（'
（a）Fs=1000Hz'
k=0:
fk=k/Tp;
subplot（3,2,6）;
plot（fk,abs（Xk））;
（a）T*FT[xa（nT）],Fs=1000Hz'
f（Hz）'
幅度'
axis（[0,Fs,0,1.2*max（abs（Xk））]）;
信号波形：
2、频域采样理论的验证
M=27;
N=32;
M;
%产生M长三角波序列x（n）
xa=0:
floor（M/2）;
xb=ceil（M/2）-1:
-1:
0;
xn=[xa,xb];
Xk=fft（xn,1024）;
%1024点FFT[x（n）],用于近似序列x（n）的TF
X32k=fft（xn,32） ;
%32点FFT[x（n）]
x32n=ifft（X32k）;
%32点IFFT[X32（k）]得到x32（n）
X16k=X32k（1:
2:
N）;
%隔点抽取X32k得到X16（K）
x16n=ifft（X16k,N/2）;
%16点IFFT[X16（k）]得到x16（n）
subplot（3,2,2）;
stem（n,xn,'
.'
boxon
（b）三角波序列x（n）'
n'
x（n）'
axis（[0,32,0,20]）
1023;
wk=2*k/1024;
%
subplot（3,2,1）;
plot（wk,abs（Xk））;
（a）FT[x（n）]'
\omega/\pi'
|X（e^j^\omega）|'
axis（[0,1,0,200]）
N/2-1;
subplot（3,2,3）;
stem（k,abs（X16k）,'
（c）16点频域采样'
k'
|X_1_6（k）|'
axis（[0,8,0,200]）
n1=0:
subplot（3,2,4）;
stem（n1,x16n,'
（d）16点IDFT[X_1_6（k）]'
x_1_6（n）'
N-1;
stem（k,abs（X32k）,'
（e）32点频域采样'
|X_3_2（k）|'
axis（[0,16,0,200]）
stem（n1,x32n,'
（f）32点IDFT[X_3_2（k）]'
x_3_2（n）'
思考题简答
先对原序列x（n）以N为周期进行周期延拓后取主值区序列，
再计算N点DFT则得到N点频域采样：
七、实验总结
1由图可见，采样序列的频谱的确是以采样频率为周期对模拟信号频谱的周期延拓。
当采样频率为1000Hz时频谱混叠很小；
当采样频率为300Hz时，在折叠频率150Hz附近频谱混叠很严重；
当采样频率为200Hz时，在折叠频率110Hz附近频谱混叠更很严重。
2频域采样定理的图验证了频域采样理论和频域采样定理。
对信号x（n）的频谱函数X（ejω）在[0，2π]上等间隔采样N=16时，N点IDFT[]得到的序列正是原序列x（n）以16为周期进行周期延拓后的主值区序列：
升级会员 