数字信号处理实验一 系统响应及系统稳定性.docx

数字信号处理实验一 系统响应及系统稳定性.docx

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

实验一系统响应及系统稳定性实验目的掌握求系统响应的方法。2. 掌握时域离散系统的时域特性。3. 分析、观察及检验系统的稳定性。二、实验原理与方法在时域中,描写系统特性的方法是差分方程和单位脉冲响应,在频域中可以用系统函数描述系统特性。已知差分方程、单位脉冲响应或者系统函数求出系统对于该输入信号的响应。本实验采用matlab语言工具箱中的filter函数和conv函数计算输入信号和系统的单位脉冲响应的线性卷积,即求出系统的响应。三、实验内容1.编程,包括产生输入信号,单位脉冲响应序列的子序列,用filter函数和conv函数求解系统输出响应的子程序。程序中要有绘制信号波形的功能。2. 给定一低通滤波器的差分方程为y(n)=0.05x(n)+0.05x(n-1)+0.9y(n-1) 输入信号x1(n)=R8(n), x2(n)=u(n).(1)输入为x1(n)时系统响应。(程序及波形)设初始状态 y(-1)=1ys=1;xn=[1,zeros(1,50)];B=[0.05,0.05];A=[1,-0.9];xi=filtic(B,A,ys);hn=filter(B,A,xn,xi);n=0:length(hn)-1;stem(n,hn,.);x1=ones(1,8);yn=conv(x1,hn);n=0:length(yn)-1;stem(n,yn,.)输入为u(n)时的系统响应:ys=1;xn=[1,zeros(1,50)];B=[0.05,0.05];A=[1,-0.9];xi=filtic(B,A,ys);hn=filter(B,A,xn,xi);n=0:length(hn)-1;stem(n,hn,.);x2=ones(1,50);yn=conv(x2,hn);n=0:length(yn)-1;stem(n,yn,.)(2)求出系统的单位脉冲响应:ys=1;xn=[1,zeros(1,50)];B=[0.05,0.05];A=[1,-0.9];xi=filtic(B,A,ys);hn=filter(B,A,xn,xi);n=0:length(hn)-1;stem(n,hn,.);3.给定系统的单位脉冲响应为h1(n)=R10(n),h2(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3).用线性卷积求出x1(n)=R8(n)分别对于两系统的输出响应,并画出波形。h1=[ones(1,10),zeros(1,30)];x1=[ones(1,8),zeros(1,30)];yn1=conv(x1,h1);n=0:length(yn1)-1;stem(n,yn1,.);对h2(n)的系统响应:h2=[1,2.5,2.5,1,zeros(1,30)];x1=[ones(1,8),zeros(1,30)];yn2=conv(x1,h2);n=0:length(yn2)-1;stem(n,yn2,.);4.给定一谐振器的差分方程为y(n)=1.8237y(n-1)-0.9801y(n-2)+b0x(n)-b0x(n-2),b0=1/100.49用实验方法检查系统是否稳定。输入信号为单位阶跃序列u(n),稳态输出为一常数(包括零),则系统一定稳定,否则系统不稳定.由实验可知,该系统是稳定的。下面分别取单位阶跃序列长度为100和300的情况。设初始状态y(-1)=1,y(-2)=1ys=[1,1];xn=[1,zeros(1,100)];B=[1/100.49,0,-1/100.49];A=[1,-1.8237,0.9801];xi=filtic(B,A,ys);hn=filter(B,A,xn,xi);x1=ones(1,50);yn=conv(x1,hn);n=0:length(yn)-1;stem(n,yn,.);ys=[1,1];xn=[1,zeros(1,300)];B=[1/100.49,0,-1/100.49];A=[1,-1.8237,0.9801];xi=filtic(B,A,ys);hn=filter(B,A,xn,xi);x1=ones(1,50);yn=conv(x1,hn);n=0:length(yn)-1;stem(n,yn,.);(2).给定输入信号为x(n)=sin(0.014n)+sin(0.4n),求系统输出响应。ys=[1,1];xn=[1,zeros(1,50)];B=[1/100.49,0,-1/100.49];A=[1,-1.8237,0.9801];xi=filtic(B,A,ys);hn=filter(B,A,xn,xi);n=0:50;x2=sin(0.014*n)+sin(0.4*n);yn=conv(x2,hn);n=0:length(yn)-1;stem

文档评论(0)

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

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

1亿VIP精品文档

相关文档