本科实验报告课程名称：信号与系统（乙）实验姓名：王超仁系：光电.doc

本科实验报告 课程名称： 信号与系统（乙）实验 姓 名： 王超仁 系： 光电信息工程学系 专 业： 信息工程（光电系） 学 号： 3110100733 指导教师： 吕 俊 张建国 2012年 5月23 日 实验报告 课程名称： 信号与系统（乙）实验 指导老师：吕俊 张建国 成绩：______ 实验名称： 线性时不变系统的时域分析 ================================================================================= 实验目的 1. 学会使用MATLAB 对线性时不变系统的时域特性进行仿真分析； 2. 熟悉LTI 系统在典型激励下的响应及特征； 3. 掌握用MATLAB 函数数值求解系统零输入响应和零状态响应的方法； 4. 学习使用Simulink 进行系统建模和仿真。 实验原理和内容 线性时不变系统（LTI）对输入信号的响应等于系统的单位冲激响应与输入信号的卷积。通过卷积可求得LTI 系统对任意输入信号的响应。如果我们从某一初始时刻（如t = 0）开始分析某LTI 系统，则通过卷积可以求得系统的零状态响应。因此卷积是LTI 系统时域分析的基本方法之一。在上一个实验中我们已学习了在MATLAB 中实现卷积的方法。很多线性LTI 系统的数学模型可以归纳为n 阶线性常系数微分方程（对于连续时间系统）和差分方程（对于离散时间系统），或者表示为一阶常微分（差分）方程组（常称为状态空间模型）。如果已知系统的输入信号和初始状态，我们或许可以通过解析的方法求解系统的响应。但对于一般的问题，往往不一定存在解析解，而且即使有解析解，对于高阶系统，求解过程也可能相当繁琐。因此本实验我们将学习用MATLAB 对LTI 系统的时域响应进行数值求解和仿真的方法。 LTI 系统的时域分析 很多信号和系统问题可归结为求解微分方程和差分方程的问题。如连续时间LTI 系统的输入信号x(t) 和响应y(t) 用线性常系数微分方程来描述：如果已知输入信号x(t) 及系统的初始条件y (0?), y′ (0?), y′′ (0?), · · · , y(n?1) (0?)，可求出系统的响应。对于离散时间LTI 系统，其输入信号x[t] 和系统响应y[t] 用线性常系数差分方程表示： 如果已知x[n] 和y[k ? 1], y[k ? 2], · · · , y[k ? N]，也可求出n ≥ k 时y[n] 的值。通常我们从某一时刻(t = 0)时刻开始分析某LTI 系统，系统的响应可表示为零输入响应零状态响应两部分之和。零输入响应是指系统初始时刻之后的输入为零、仅由系统的初始状态引起的系统的输出，零状态响应是指系统的初始状态为零、仅由系统初始时刻之后的输入引起的系统的输出。系统响应的时域解析解法的过程是先求出微分方程（或差分方程）的齐次解，再根据输入信号的形式确定方程的特解，然后根据初始条件确定解的系数，最后得到系统的响应。零输入响应和零状态响应需要根据输入信号和初始条件的不同，分别求解得到。这种计算可能相当繁琐，而且不一定存在解析解，需要通过数值方法来求解。 数值求解的基本原理 3.2式所描述的离散时间因果LTI 系统（通常也称为数字滤波器）可分为两个部分：一部分是自递归(AR, autoregressive)的IIR 滤波器：AR 滤波器的当前的输出与之前的输出有关可递归地求得输出；另一部分是滑动平均(MA, moving average)的FIR 滤波器：MA 滤波器的输出是非递归的，只和输入有关，可通过卷积计算。因此3.2式给出的IIR 滤波器也称为ARMA 滤波器。一般来说，总是可以将3.2式写为递推的形式：根据系统的初始条件和输入信号即可递推地算出y[n]。这种方法简单直观，且适合计算机求解，但通常只能求出系统输出序列的值，不能像解析解那样得到系统输出的数学表达式。 对于3.1式描述的连续时间系统，数值求解的原理是先将微分方程离散化，近似为相应的差分方程(将dt 近似为Δt)，再求解差分方程，因此其数值求解方法方法与上述的离散系统的求解方法是类似的。MATLAB 的工具箱已有求解上述问题的函数，因此只需简单的函数调用即可数值求解系统的各种响应，如冲激响应、阶跃响应、零输入响应和零状态响应等，可大大简化复杂系统的设计和分析过程。另外MathWorks 还提供了用于系统建模和仿真的Simulink 软件包，配合MATLAB 使用可省去许多代码的重复编写工作，提高求解的效率。 实验数据记录、处理和分析 编写函数如下所示 首先编写函数‘recur.m’ function y=recur(c,d,n,x,X