差分方程及其求解数字信号处理实验研究报告.docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 计算机与信息工程学院验证性实验报告 专业：通信工程 年级/班级：2011级 2013—2014学年第一学期 课程名称 数字信号处理 指导教师 段新涛 本组成员 学号姓名 实验地点 计科楼111 实验时间 周五7-8节 项目名称 差分方程及其求解 实验类型 验证性 一、实验目地 1．学习并掌握系统地差分方程表示方法以及差分方程地相关概念. 2．熟练使用filter函数对差分方程进行数值求解. 3．掌握差分方程地求解及MATLAB实现方法. 二、实验原理及方法 1．一LTI系统可以用一个线性常系数差分方程表示： ，任意n 如果，那么这个差分方程就是N阶地，已知系统地输入序列，用这个方程可以根据当前输入x(n)和以前M点地输入x(n-m),…,x(n-1),以及以前N点地输出y(n-N),…,y(n-1)来计算当前输出y(n).在实际中这个方程在时间上是从到朝前计算地，因此该方程地另一种形式是：b5E2RGbCAP 方程地解能以下面形式求得：分别为方程地齐次解跟特解部分.已知输入和差分方程地稀疏，可用filter对差分方程进行数值求解.最简单形式为：y=filter(b,a,x) 其中b=[b0,b1,…,bM];a=[a0,a1,…,aN];p1EanqFDPw 2．上面差分方程解地形式为齐次解和特解，另外还可以求零输入解和零状态解理论计算中要用到z变换，请好好掌握z变换地内容.用MATLAB实现时，若已知初始条件，则应用y=filter(b,a, x, xic)来求完全响应.这里xic是初始状态输入数组.MATLAB还提供一种filtic函数来得到xic.DXDiTa9E3d xic=filtic(b,a,Y,X )其中b和a是滤波器系数数组，Y和X是分别从y(n)和x(n)地初始条件得来地初始状态数组.RTCrpUDGiT 三.实验内容 1、已知下面差分方程： ;任意n 要求： a.计算并画出在n=-20,…,100地脉冲响应h(n); b.计算并画出在n=-20,…,100地单位阶跃相应s(n). 2.解以下差分方程： 其中，初始条件为y(-1)=4和y(-2)=10. 要求先用理论计算，再用MATLAB编程实现，并对比两个结果. 参考流程图： 实验内容1． 开始输入差分方程地系数矩阵 开始 输入差分方程地系数矩阵 输入一冲激或阶跃序列 结束 利用 filter 函数计算冲激响应或阶跃响应 画出结果图形 实验内容2. 开始 开始 输入输入序列 x 输入参数及初始条件 结束 计算输出序列 画出图形 四、实验报告要求 1．总结差分方程地性质及应用． 2．写出实验程序．记录实验数据并与理论计算作比较，总结结果. 五、实验程序及结果 针对问题一： 问题a地程序如下： 先建立一个生成单位抽样序列地M文件 function [x,n]=impseq(n0,n1,n2) if((n0n1)|(n0n2)|(n1n2)) error(参数必须满足n1=n0=n2) end n=n1:n2; x=[(n-n0)==0]; 主程序如下： clc;clear all b=1;a=[1 -1 0.9]; [x,n]=impseq(0,-20,100); y=filter(b,a,x); plot(n,y) title(输出脉冲响应h(n)) xlabel(n) 问题b地程序如下： 先建立一个生成单位阶跃序列地M文件 function [x,n]=stepseq(n0,n1,n2) if((n0n1)|(n0n2)|(n1n2)) error(参数必须满足n1=n0=n2) end n=[n1:n2]; x=[(n-n0)=0]; 主程序如下： clc;clear all b=1;a=[1 -1 0.9]; [x,n]= stepseq(0,-20,100); y=filter(b,a,x); plot(n,y) title(输出单位阶跃响应s(n)) xlabel(n) 针对问题二： 1、理论计算 系统函数为 输入信号地Z变换为 对给出地差分方程求单边Z变换，得： 故 把n=1，2，3…11代入上式得：y(1)=2, y(2)=1.25, y(3)= 0.9375, y(4)= 0.7969,jLBHrnAILg y(5)= 0.7305, y(6)= 0.6982,y(7)= 0.6824, y(8)= 0.6745, y(9)= 0.6706, y(10)= 0.6686, y(11)= 0.6676xHAQX74J0X 2、MATLAB编程实现 clc;clear all b=[1