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

计算机与信息工程学院验证性实验报告 专业：通信工程 年级/班级：2011级节 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)。在实际中这个方程在时间上是从到朝前计算的，因此该方程的另一种形式是： 方程的解能以下面形式求得：分别为方程的齐次解跟特解部分。已知输入和差分方程的稀疏，可用filter对差分方程进行数值求解。最简单形式为：y=filter(b,a,x) 其中b=[b0,b1,…,bM];a=[a0,a1,…,aN]; 2．上面差分方程解的形式为齐次解和特解，另外还可以求零输入解和零状态解理论计算中要用到z变换，请好好掌握z变换的内容。用MATLAB实现时，若已知初始条件，则应用y=filter(b,a, x, xic)来求完全响应。这里xic是初始状态输入数组。MATLAB还提供一种filtic函数来得到xic。 xic=filtic(b,a,Y,X )其中b和