原创力文档- 1、本文档共5页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验2 离散系统的差分方程、冲激响应和卷积分析
实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。
实验原理：离散系统
其输入、输出关系可用以下差分方程描述：
输入信号分解为冲激信号，。记系统单位冲激响应
，则系统响应为如下的卷积计算式：
当时，h[n]是有限长度的（n：[0，M]），称系统为FIR系统；反之，称系统为IIR系统。
在MATLAB中，可以用函数y=Filter(p,d,x) 求解差分方程，也可以用函数 y=Conv(x,h)计算卷积。
实验内容：编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。
实验要求：给出理论计算结果和程序计算结果并讨论。
实验过程：
（1）单位冲激响应：
a=[1,0.6,0.08];
b=[1,-1];
N=20;
x=[1,zeros(1,N)];
y=filter(b,a,x);
stem(y);
xlabel(时间序列n);
ylabel(信号幅度);
title(单位冲激响应h(n));
（2）单位阶跃响应：
a=[1,0.6,0.08];
b=[1,-1];
N=20;
x=[ ones(1,N)];
y=filter(b,a,x);
stem(y);
xlabel(时间序号);
ylabel(信号幅度);
title(单位阶跃响应h（n）);
理论分析：
由差分方程得系统函数为：
利用分部分式法可得：，z反变换得：
h(n)即为单位冲击响应。
对于阶跃响应， ,所以阶跃响应
=
=
所以，
可见，实验结果与理论分析是一致的。
单位冲激响应
a=[1];
b=[0.2,0.2,0.2,0.2,0.2,0.2];
N=20;
x=[1,zeros(1,N)];
y=filter(b,a,x);
stem(y);
xlabel(时间序号n);
ylabel(信号幅度);
title(单位冲激响应h（n）);
单位阶跃响应
a=1;
b=[0.2,0.2,0.2,0.2,0.2,0.2];
N=20;
x=[ones(1,N)];
y=filter(b,a,x);
stem(y);
xlabel(时间序号);
ylabel(信号幅度);
title(单位阶跃响应h（n）);
理论分析：
由差分方程得系统函数为：
z反变换得：
h(n)即为单位冲击响应。
对于阶跃响应， ,所以阶跃响应:
= 所以，
可见，实验结果与理论分析是一致的。
文档评论（0）