- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验二-时域采样和频域采样汇编
实验二:时域采样与频域采样
一、实验目的
时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化,以及如何选择采样频率才能使采样后的信号不丢失信息;要求掌握频率域采样会引起时域周期化的概念,以及频率域采样定理及其对频域采样点数选择的指导作用。
二、实验原理及方法
1、时域采样定理的要点:
a 对模拟信号以间隔T进行时域等间隔理想采样,形成的采样信号的频谱是原模拟信号频谱以采样角频率()为周期进行周期延拓。公式为:
b 采样频率必须大于等于模拟信号最高频率的两倍以上,才能使采样信号的频谱不产生频谱混叠。利用计算机计算上式并不方便,下面我们导出另外一个公式,以便用计算机上进行实验。
理想采样信号和模拟信号之间的关系为:
对上式进行傅立叶变换,得到:
在上式的积分号内只有当时,才有非零值,因此:
上式中,在数值上=,再将代入,得到:
上式的右边就是序列的傅立叶变换,即
上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到,只要将自变量ω用代替即可。
2、频域采样定理的要点:
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,它的幅频特性曲线如图2.1 图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代表的频率为 。
要求: 编写实验程序,计算、和的幅度特性,并绘图显示
析频谱混叠失真。
Matlab源代码:
A 444.128;a 50*sqrt 2 *pi;w0 50*sqrt 2 *pi;
Tp 50/1000;F1 1000;F2 300;F3 200; %观察时间Tp 50ms
T1 1/F1;T2 1/F2;T3 1/F3; %不同的采样频率
n1 0:Tp*F1-1;n2 0:Tp*F2-1;n3 0:Tp*F3-1; %产生不同的长度区间n1,n2,n3
x1 A*exp -a*n1*T1 .*sin w0*n1*T1 ; %产生采样序列x1 n
x2 A*exp -a*n2*T2 .*sin w0*n2*T2 ; %产生采样序列x2 n
x3 A*exp -a*n3*T3 .*sin w0*n3*T3 ; %产生采样序列x3 n
f1 fft x1,length n1 ; %采样序列x1 n 的FFT变换
f2 fft x2,length n2 ; %采样序列x2 n 的FFT变换
f3 fft x3,length n3 ; %采样序列x3 n 的FFT变换
k1 0:length f1 -1;
fk1 k1/Tp; %x1 n 的频谱的横坐标的取值
k2 0:length f2 -1;
fk2 k2/Tp; %x2 n 的频谱的横坐标的取值
k3 0:length f3 -1;
fk3 k3/Tp; %x3 n 的频谱的横坐标的取值
subplot 3,2,1
stem n1,x1,.
title a Fs 1000Hz ;
xlabel n ;ylabel x1 n ;
subplot 3,2,3
stem n2,x2,.
title b Fs 300Hz ;
xlabel n ;ylabel x2 n ;
subplot 3,2,5
stem n3,x3,.
title c Fs 200Hz ;
xlabel
文档评论(0)