2014080330000DSP实验报告.docx

  1. 1、本文档共21页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
2014080330000DSP实验报告

信息学院数字信号处理实验报告姓  名:班 级 学 号:专 业:指 导 教 师:一. 实验目的2二. 实验原理和方法2三. 实验内容4信号、系统及系统响应4实验结论13用FFT做谱分析13实验结论20四. 思考题20一. 实验目的熟悉连续信号经过理想抽样前后的频谱变化关系,加深对时域抽样定理的理解。熟悉时域离散系统的时域特性。利用卷积方法观察分析系统的时域特性。掌握序列傅里叶变换的计算机实验方法,利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。进一步加深DFT算法原理和基本性质的理解(因为 FFT只是DFT的一种快速算法,所以FFT的运算结果必然满足DFT的基本性质)。熟悉FFT算法原理和FFT子程序的应用。学习用FFT对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便在实际中正确应用FFT。二. 实验原理和方法抽样是连续信号数字处理的第一个关键环节。对抽样过程的研究不仅可以了解抽样前后信号时域和频域特性发生的变化以及信号信息不丢失的条件,而且可以加深对傅里叶变换、Z变换和序列傅里叶变换之间关系式的理解。我们知道,对一个连续信号xa(t)进行理想抽样的过程可用(1.1)式表示。 = δT(t) ——(1.1)其中为xa(t)的理想抽样,δT(t)为周期冲激脉冲,即——(1.2)的傅里叶变换为= ——(1.3)(1.3)式表明为的周期延拓,其延拓周期为抽样角频率(Ωs=2π/T)。抽样前后信号的频谱示意图见“参考教材图1-29”。只有满足抽样定理时,才不会发生频率混叠失真。在计算机上用高级语言编程直接按(1.3)式计算理想抽样的频谱很不方便。下面导出用序列的傅里叶变换来计算的公式。将(1.2)式代入(1.1)式并进行傅里叶变换, = = = = ——(1.4)式中的xa(nT)就是采样后得到的序列x(n),即x(n) = xa(nT)x(n)的序列傅里叶变换为X(ejω) = —— (1.5)比较(1.5)和(1.4)可知=X(ejω) |ω = ΩT——(1.6)这说明两者之间只在频率度量上差一个常数因子T。实验过程中应注意这一差别。离散信号和系统在时域均可用序列来表示。序列图形给人以形象直观的印象,它可加深我们对信号和系统的时域特征的理解。本实验还将观察分析几种信号及系统的时域特性。为了在数字计算机上观察分析各种序列的频域特性,通常对X(ejω)在[0,2π]上进行M点采样来观察分析。对长度为N的有限长序列x(n),有——(1.7)其中 , k = 0,1, …, M-1通常M应取得大一些,以便观察谱的细节变化。取模||可绘出幅频待性曲线。一个时域离散线性非移变系统的输入/输出关系为y(n) = x(n) * h(n) = ——(1.8)这里,y(n)为系统的输出序列,x(n)为输入序列。h(n)、x(n)可以是无限长,也可以是有限长。为了计算机绘图观察方便,主要讨论有限长情况。如果h(n)和x(n)的长度分别为N和M,则y(n)的长度为L = N + M - 1。这样,(1.8)式所描述的卷积运算就是序列移位、相乘和累加的过程,所以编程十分简单。上述卷积运算也可以在频域实现(即卷积定理:时域卷积,频域相乘。)Y(ejω) = X(ejω)H(ejω) —— (1.9)(1.9)式右边的相乘是在各频点{ωk}上的频谱值相乘。三. 实验内容信号、系统及系统响应1.单边指数脉冲;xa(t)=10*exp(-5t)连续信号源程序Dt=0.0001;t=0:Dt:0.5;x=exp(-10*t);Wmax=2*pi*200;K=500;k=-K:1:K;W=k*Wmax/K;X=x*exp(-j*t'*W)*Dt;X=abs(X);figure(1);subplot(2,1,1);plot(t,x,'b');grid;xlabel('t/s');ylabel('x(t)')title('连续信号')subplot(2,1,2);plot(W/(2*pi),X,'b');xlabel('f/HZ');ylabel('X(jf)');tltle('连续时间信号的傅里叶变换');离散信号:fs=50T=1/fst=0:T:0.5x=exp(-10*t);Wmax=2*pi*200;K=500;k=-K:1:K;W=k*Wmax/K;X=x*exp(-j*t'*W)*Dt;X=abs(X);figure(1);subplot(2,1,1);stem(t,x,'b');grid;xlabel('t/s');ylabel('x(t)')title('离散信号')subplot(2,1,2);plot(W/(2*pi),X,'b');xlabel(

文档评论(0)

haihang2017 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档