- 1、本文档共9页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
本科学生验证性实验报告
学号 104090459 姓名 张静
学院 物电学院 专业、班级10电子
实验课程名称 数字信号处理实验
教师及职称 杨卫平 教授
开课学期 2013 至 2013 学年 下 学期
填报时间 2013 年 5 月 23 日
云南师范大学教务处编印
实验序号 2 实验名称 离散系统分析 实验时间 2013.5.21 实验室 同析3-313 一.实验预习 1.实验目的
深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义,熟练掌握利用MATLAB分析离散系统的时域响应、频响特性和零极点的方法。掌握利用DTFT和DFT确定系统特性的原理和方法。 2.实验原理、实验流程或装置示意图
MATLAB提供了许多可用于分析线性时不变连续系统的函数,主要包含有系统函数、系统时域响应、系统频域响应等分析函数。
1. 离散系统的时域响应
离散时间LTI系统可用如下的线性常系数差分方程来描述:
已知输入信号x[k]以及系统初始状态y[-1],y[-2] ...,就可以求出系统的响应。MATLAB提供了计算差分方程的数值解的函数,可以计算上述差分方程描述的离散系统的单位脉冲响应、零状态响应和完全响应等。
在调用MATLAB函数时,需要利用描述该离散系统的系数函数。对差分方程进行Z变换即可得系统函数:
在MATLAB中可使用向量a和向量b分别保存分母多项式和分子多项式的系数:
这些系数均从z0按z的降幂排列。
离散系统的单位冲激响应h[k]的计算
[h,k] = impz(b, a):计算系统的单位脉冲响应h[k]和相应的时间向量k;也可简写为:h = impz(b, a)。其中:
[h,k] = impz(b, a, n):计算n点单位脉冲响应h[k];
也可简写为:h = impz(b, a, n)。
impz(b, a):绘制单位脉冲响应h[k]的图形。
离散系统响应y[k]的计算
y = filter(b, a, x):
计算系统在输入x作用下的零状态响应y[k];
y = filter(b, a, x, zi):
计算系统在输入x和初始状态作用下的完全响应y[k]。
zi是由系统的初始状态经过filtic函数转换而得到的初始条件:zi= filtic (b, a, Y0) ,Y0为系统的初始状态,
Y0= [ y[-1],y[-2],y[-3],...]。
2.离散系统的系统函数零极点分析
离散LTI系统的系统函数H(z)可以表示为零极点形式:
使用MATLAB提供的roots函数计算离散系统的零极点;
使用zplane函数绘制离散系统的零极点分布图。
注意:在利用这些函数时,要求H(z)的分子多项式和分母多项式的系数的个数相等,若不等则需要补零。
3.离散系统的频率响应
当离散因果LTI系统的系统函数H(z)的极点全部位于z平面单位圆内时,系统的频率响应可由H(z)求出,即
[H, w]=freqz(b, a, n):
计算系统的n点频率响应H,w为频率点向量。
H=freqz(b, a, w) :计算系统在指定频率点向量w上的频响;
freqz(b,a): 绘制频率响应曲线。
其中:b和a分别为系统函数H(z)的分子分母系数矩阵;
4.利用DTFT和DFT确定离散系统的特性
在很多情况下,需要根据LTI系统的输入和输出对系统进行辨识,即通过测量系统在已知输入x[k]激励下的响应y[k]来确定系统的特性。
若系统的脉冲响应为h[k],由于存在y[k]=x[k]* h[k],因而可在时域通过解卷积方法求解h[k]。
在实际应用中,进行信号解卷积比较困难。因此,通常从频域来分析系统,这样就可以将时域的卷积转变为频域的乘积,从而通过分析系统输入序列和输出序列的频率特性来确定系统的频率特性H() ,再由 H() 得到系统的脉冲响应h[k]。
若该LTI系统输入x[k]的DTFT为 X() ,系统输出y[k]的DTFT为 Y() ,则系统的频率特性可表示为:
有限长序列的DTFT可以利用FFT计算出其在区间内的N个等间隔频率点上的样点值:
X = fft(x,N); Y = fft(y,N);
再利用 H = Y./X 和 h = ifft(H,N) 可以得到系统的单位脉冲响应h[k]
3.实验设备及材料
计算机一台及MATLA
您可能关注的文档
- 安徽省合肥市2014届高三第二次教学质量检测试题.doc
- 安徽省合肥市2015年一模试卷.doc
- 安徽省合肥市肥西县肥西中学2014—2015学年高三上学期第3轮周考物理试题及答案.doc
- 安徽省名校联盟2013届高三第一次联考文科综合试题.docx
- 安徽省太湖中学2014届高三语文复习班上学期期中试题.doc
- 安徽省安庆市十校2011-2012学年八年级语文上学期期中联考试题新人教版.doc
- 安徽省宿州市十三所重点中学2011-2012学年高一下学期期中质量检测地理试题.doc
- 安徽省怀宁中学2011届高三第二次质量检测语文.doc
- 安徽省淮南一中等“四校”2013届高三第一次联考语文试题.doc
- 安徽省滁州市高中信息技术《第三章加工信息》教案.doc
- 某某单位2024年党建工作总结及2025年工作计划.doc
- 某某市发改委关于2024年度落实党风廉政建设工作责任制情况的报告.doc
- 某某局2024年全面从严治党和党风廉政建设工作总结.doc
- 某某区财政局2024年法治政府建设总结及2025年工作谋划.doc
- 2024年党管武装工作述职报告2篇.doc
- 2024年度国企党委书记抓基层党建工作述职报告3篇.doc
- 公司党委书记2024年述职述廉报告.docx
- 2024年度乡镇党委领导班子民主生活会(四个带头)对照检查材料.doc
- 市医疗保障局关于2024年法治政府建设工作情况的报告.docx
- 市民政局党组2024年巡察整改工作情况报告.docx
文档评论(0)