DSP_03离散时间信号-差分方程例析.pptVIP

数字信号处理 陈鹏 3299340 线性常系数差分方程Linear Constant Coefficients Difference Equation Lesson 3 复习提问 如何判断周期序列？ 什么响应反映了系统的特性？ 从哪两个方面来理解h(n)？ 离散卷积的四个步骤是什么？ 线性常系数差分方程Linear Constant Coefficients Difference Equation 线性常系数差分方程Linear Constant Coefficients Difference Equation 1.零输入响应的求解 Zero Input Response 1.零输入响应的求解 Zero Input Response 1.零输入响应的求解 Zero Input Response 2.零状态响应的求解 Zero State Response 2.零状态响应的求解 Zero State Response 2.零状态响应的求解 Zero State Response 2.零状态响应的求解 Zero State Response 递推法 递推法 根据因果系统规定的初始条件递推如下 差分方程表示 差分方程表示 用MATLAB解线性常系数差分方程 MATLAB信号处理工具箱中提供的filter函数，可以实现线性常系数差分方程的递推解法，调用格式如下： y=filter(b,a,x,xi); 其中x是输入信号向量，b和a是系统差分方程的系数矩阵，即b=[b0,b1,…,bM]和a=[a0,a1,…,aN]，a0=1。xi是和初始条件有关的向量，用函数xi=filtic(b,a,ys,xs)得到，其中ys和xs是初始条件向量，即ys=[y(-1),y(-2),...]，xs=[x(-1),x(-2),...]。如果是因果序列，则xs=0。 用函数filter(b,a,x,xi)计算输出y，如果和输入信号和系统的初始状态有关，称为系统的零输入响应。如果系统的输入条件为零，就默认xi=0，调用格式为y=filter(b,a,xi)。 用MATLAB解线性常系数差分方程 例2.12 用MATLAB计算差分方程 当输入序列为x(n)=?(n) 时的输出结果y(n),0?n?30 。 用MATLAB解线性常系数差分方程 解 MATLAB程序如下： N=31; b=[0.8 -0.44 0.36 0.22]; a=[1 0.7 -0.45 -0.6]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(b,a,x); stem(k,y) xlabel(n);ylabel(输出y(n)) 用MATLAB解线性常系数差分方程 用MATLAB解线性常系数差分方程 例2.13 用MATLAB计算差分方程 当输入为x(n)=?(n)，初始条件y(-1)=1 ，求输出结果y(n),0?n?30。 用MATLAB解线性常系数差分方程 解 MATLAB程序如下： N=31; a=[1 -0.8];ys=1; b=1; x=[1 zeros(1,N-1)]; xi=filtic(b,a,ys);%finds the initial conditions k=0:1:N-1; y=filter(b,a,x,xi); stem(k,y) xlabel(n); ylabel(输出y(n)) 用MATLAB解线性常系数差分方程 线性非移变系统可以用下列常系数差分方程来描述 其中 称为差分方程或系统的阶次。 引入单位延迟算子 ，即 ，于是上述差分方程可以写为： 求系统的零输入响应，响应形式与系统的齐次解相同，由系统初始条件确定齐次解中的待定系数。 求系统的零状态响应。 系统的完全响应为零输入响应与零状态响应之和。 求解下列线性非移变系统的差分方程可以按照以下3步进行。 首先写出齐次方程 Homogeneous Equation 得到特征方程 Characteristic Equation N个特征根 Characteristic Roots 零输入响应就是下面齐次解的形式 Homogene