网站大量收购闲置独家精品文档,联系QQ:2885784924

【精品】数字信号处理实验,Matlab实验数字与信号处理.doc

【精品】数字信号处理实验,Matlab实验数字与信号处理.doc

  1. 1、本文档共10页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验一 信号、系统及系统响应 实验目的 熟悉理想采样的性质,了解信号采样前后的频谱变化,加深对采样定理的理解。 熟悉离散信号和喜用的时域特性。 熟悉线性卷积的计算变成方法,利用卷积的方法,观察、分析系统响应的时域特性。 掌握序列傅氏变换的计算机实现方法,利用序列的傅氏变换对离散信号、系统及系统响应进行频域分析。 实验原理 (1)连续时间信号的采样 采样是从连续时间信号到离散时间信号的过度桥梁,对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生的变化以及信号内容部丢失的条件,而且有助于加深对拉式变化。傅氏变化、Z变化和序列傅氏变换之间的关系。 对一个连续时间信号进行理想采样的过程可以表示为该信号的一个周期冲击脉冲的乘积,即 (1--1) 其中是连续信号的理想采样,是周期冲击脉冲 (1--2) (1--3) 其中T为采样周期,是采样角频率。设是连续时间信号的双边拉氏变换,即有: (1--4) 此时理想采样信号(1--5) 作为拉氏变换的一种特例,信号理想采样的傅里叶变换 (1---6) 由式(1--5)和(1--6)可知,信号理想采样后的频谱是原信号频谱的周期延拓,其延拓周期等于采样。根据Shannon取样定理,如果原信号是带限信号,且采样频率高于原信号最高频率分量2倍,则采样以后不会发生频谱混淆现象。 在计算机处理时,不采用式(1--6)计算信号的频谱,而是利用序列的傅里叶变换计算信号的频谱,定义序列,根据Z变换的定义,可以得到序列的Z变换为: (1--7) 以代替上式中的Z,就可以得到序列的傅里叶变换 (1--8) 式(1--6)和(1--8)具有如下关系: (1--9) 由式(1--9)可知,在分析一个连续时间信号的频谱时,可以通过取样将有关的计算转化为序列傅里叶变化的计算。 有限长序列分析 一般来说,在计算机上不可能,也不必要处理连续的曲线,通常,我们只要观察、分析在某些频率点上的值。对于长度为N的有限长序列 (1--10) 一般只需要在0—2之间均匀地取M个频率点,计算这些点上的序列傅里叶变换 (1--11) 其中是一个复函数,它的模就是幅频特性曲线。 信号卷积 一个线性时不变离散系统的响应可以用它的单位冲激响应和输入信号的卷积来表示: (1--12) 根据傅里叶变换和Z变换,与式(1--12)对应应该由 (1--13) (1--14) 式(1--12)告诉我们可以通过对两个序列的移位、相乘、累加计算信号响应;而式(1--14)告诉我们卷积运算也可以在频域上用乘积实现。 实验内容及步骤 (一)编制实验用主程序及相应子程序 信号产生子程序,包括: 理想采样信号序列:对信号进行理想采样,可以得到一个理想的采样信号序列:,其中,A为幅度因子,是衰减因子,是频率。T为采样周期。 单位脉冲序列 矩形序列 ,其中N=10 2、系统单位脉冲响应序列产生子程序,本实验中用到两种FIR系统: (1) (2) 3、有限长序列线性卷积子程序,用于计算两个给定长度(分别是M和N)的序列的卷积,输出序列长度为. (二)上机实验内容 在编制以上各部分程序以后,编制主程序调用各个功能模块实现对信号、系统和系统响应的时域和频域分析,完成以下实验内容。 1.产生理想采样信号序列xa(n) ,A=444.128; a=50*sqrt(2.0)*pi;w0=50*sqrt(2.0)*pi; 编写M文件zhic1 n=0:50; A=444.128; a=50*sqrt(2)*pi; T=input(T=); w0=50*sqrt(2)*pi; x=A*exp(-a*n*T).*sin(w0*n*T); m=input(m=);subplot(3,2,m);stem(n,x); title(理想采样信号序列); k=-25:25; f=(1/25)*k*(1/T); X=x*(exp(-j*pi/12.5)).^(n*k); magX=abs(X); o=input(o=);subplot(3,2,o);stem(f,magX); 保存zhic1 zhic1 T=1/1000 m=1 o=2 zhic1 T=1/300 m

文档评论(0)

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

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

1亿VIP精品文档

相关文档