MATLAB实验报告——运用MATLAB求解和分析线性时不变系统1.docVIP

MATLAB实 验 报 告 课程名称 MATLAB程序设计 实验日期 2015 年 05 月 11 日 学生姓名学号班级实验名称运用MATLAB求解和分析线性时不变系统 实验仪器 MATLAB7.1 Windows XP实验目的 1.熟悉线性时不变LTI系统在典型激励信号下的响应及其特征。 2.掌握线性LTI系统单位冲激响应的求解方法。 3.熟悉MATLAB的相关函数的调用格式和作用。 4.会用MATLAB对系统进行时域分析。 实验基本原理 Impulse函数 功能：计算并画出系统的冲激响应。 调用格式： Impulse（sys）：其中sys 可以是利用命令tf，zpk或ss建立的系统函数。 Impulse（sys，t）计算并画出系统在向量t定义的时间内的冲激响应。 Y=impulse（sys，t）：保存系统的输出值。 Step函数 功能：计算并画出系统的阶跃响应曲线。 调用格式： Step（sys）：其中sys可以是利用命令tf，zpk，或ss建立的系统。 Step（sys，t）：计算并画出系统在向量t定义的时间内的阶跃响应。 Lsim函数 功能：计算并画出系统在任意输入下的零状态响应。 调用格式： Lsim（sys,x,t）:其中sys可以是利用命令tf,zpk或ss建立的系统函数，x是系统的输入，t定义的是时间范围。 Lsim（sys，x,t,zi）：计算出系统在任意输入和零状态下的全响应，sys必须是状态空间形式的系统函数，zi是系统的原始状态。 roots函数 功能：计算齐次多项式的根。 调用格式： r=roots（b）：计算多项式b的根，r为多项式的根。 5．impz函数 功能: 求离散系统单位脉冲响应,并绘制其时域波形。 调用格式：impz(b ,a) :以默认方式绘出向量a , b 定义的离散系统的单位脉冲响应的离散时域波形. impz(b ,a ,n) :绘出由向量a , b定义的离散系统在0—n (n必须为整数)离散时间范围内的单位序列响应的时域波形. impz(b ,a ,n1：n2) : 绘出由向量a , b定义的离散系统在n1—n2 (n1 , n2必须为整数,且n1n2)离散时间范围内的单位序列响应的时域波形。 y=impz(b , a , n1 ：n2): 并不绘出系统单位序列响应的时域波形,而是求出向量a , b定义的离散系统在n1—n2(n1 , n2必须为在整数,且n1n2)离散时间范围内的单位序列响应的数值。 filter函数 功能:对输入数据数字滤波. 调用格式: y=filter(b , a ,x) :返回向量a , b定义的离散系统在输入为x时的零状态响应.如果x是一个矩阵,那么函数filter对矩阵x的列进行操作;如果x是一个N维数组,函数filter对数组中的一个非零量进行操作[y , zf]=filter(b , a ,x): 返回了一个状态向量的最终值zf. [y , zf]= filter(b , a , x , zi) :指定了滤波器的初始状态向量zi. [y ,zf]= filter(b ,a , x ,zi , dim ) 则是给定x中要进行滤波的维数dim.如果要使用零初始状态,则将zi设为空向量。 实验内容和步骤  实践内容： 例12.4-1 用MATLAB在时域中求解 的齐次解。 起始条件为，求系统的零输入响应。 例12.4-2 求连续时间系统，当时的特解。 3. 例12.4-3 用户MATLAB在时域中求解 ，的零响入相应。 4.例12.4-5 用MATLAB求解方程零状态响应分量，已知系统差分方程为系统的激励序列。 数据记录及分析  例12.4-1 （1）用MATLAB在时域中求解 的齐次解。 程序截图： （2） 程序截图： 运行图像截图： 例12.4-2 求连续时间系统，当时的特解。 程序截图： 运行图像截图： 3. 例12.4-3 用户MATLAB在时域中求解 ，的零响入相应。 解法一： 程序截图： 运行图像截图： 解法二： 用MATLAB的dsolve函数直接获得零状态响应的表达式 程序截图： 运行图像截图： 解法三： 程序截图： 运行图像截图： 4.例12.4-5 用MATLAB求解方程零状态响应分量，已知系统差分方程为系统的激励序列。 程序截图： 运行图像截图： 实验注意事项 每一次M-file的修改后，都要存盘。 实验心得实验例题12.4-1中最初没搞清楚roots的用法，耗费了较长时间，经过咨询老师后才知道原来本例题并不复杂，加之由于书本印刷错误