Matlab程序设计实训

1、实习报告实习名称 matlab程序设计实训专业班级 *成绩评定考核内容实习表现实习报告实习考核综合评定成绩成绩削毛与信息工程学院和谐勤奋求是创新实习考核和成绩评定办法1. 实习成绩的考核由指导教师根据实习表现、实习报告、实习成果、现场操作、设计、 口试或笔试等几个方面，给出各项权重，综合评定。该实习考核教研室主任审核，主 管院长审批备案。2. 成绩评定釆用五级分制，即优、良、中、及格、不及格。3. 参加本次实习吋间不足三分之二或旷课四天以上者，不得参加本次考核，按不及格处 理。4. 实习结束一周内，指导教师提交实习成绩和实习总结。5. 实习过程考核和实习成绩在教师手册中有记载。实习报告内容实习

2、报告内容、格式各专业根据实习类别（技能实习、认识实习、生产实习、毕 业实习等）统一规范，经教研室主任审核、主管院长审批备案。注：1.实习任务书和实习指导书在实习前发给学生，实习任务书放置在实习报告封面后和正文目录前。2.为了节省纸张，保护环境，便于保管实习报告，统一采用a4纸，实习报告建议双面打印（正文采 用宋体五号字）或手写，右侧装订。2016/2017学年第二学期matlab程序设计实训计划指导教师：*班级：自动化1541、2班实习地点：1教8楼机房一、实习目的实习的目的：了解并掌握matlab软件的各种应用开发环境；了解并掌握matlab软件强 大的科学运算功能；了解并掌握利用matla

3、b语言程序设计流稈进行程序设计的方法；了解并 掌握matlab的图形处理功能；了解matlab软件的交互式仿真功能；了解并掌握matlab软 件在控制系统建模、分析以及设计屮的应用。实习的任务：通过实习可以使学生将matlab软件与控制系统理论及仿真相结合，加深对 控制系统理论知识的理解，验证理论知识的正确性，提高软件的实际操作能力，可以有效地 培养学生分析及解决实际问题的能力，不断提高自学能力，养成良好的科学态度以及实事求 是、严谨踏实的学习风气，成为真正的应用型本科人才。通过实训达到边学边练、理论与实 际操作相结合的学习目的。二、实习内容掌握matlab编程环境、基础知识：掌握矩阵生成及运

4、算、数组运算；矩阵函数；掌握常 用绘图命令及方法；掌握基本图形处理、二维图形处理方法；掌握可视化m文件程序设计及 调试；掌握simulink仿真工具的建立和基本仿真方法。了解并掌握基于matlab的控制系统数学模型及分析方法；了解并掌握基于matlab的控 制系统pid控制器的基本设计方法。综合应用matlab程序设计语言，根据设计任务要求，完 成软件程序的编制并调试运行。三、时间安排(2017.07.0102017.07.014)序号实习 内容时间（天）1matlab编程环境、基础知识、矩阵生成及运算、数组运算、 矩阵函数、常用绘图命令及方法12基本图形处理、二维图形处理方法、可视化m文件程

5、序设计 及调试、simulink仿真工具的建立和基本仿真方法13基于matlab的控制系统数学模型及分析方法14基于matlab的控制系统pid控制器的基本设计方法、综合性 应用15上机考试，撰写实习报告1合计5四、 实习要求实习期间，严格遵守课堂纪律，不许迟到、早退；实习期间，学习态度端正，认真上机 练习；熟练掌握matlab软件的应用方法；完成实习报告一份，内容充实，写出实习的体会与 收获；实习报告打印成文。目录1实习目的和意义错误!未定义书签。2 matalab/simulink简介错误!未定义书签。2matalab的发展历程错误!未定义书签。2.2 matalab的影响错误!未定义书签

6、。2.3 matalab与控制系统仿真错误!未定义书签。2.4 simulink 简介23 matalab语言基础23matalab基础知识34 simulink 仿真工具84.1 运彳亍simulik 演示程序94.2 simulink模型的建立105控制系统数学模型105.1动态过程微分方程描述105.2 高阶系统的时域分析116控制系统分析126.1 时域分析126.2 频域分析157控制器设计157比例控制器157.2积分控制器167.3比例积分控制器177.4比例微分控制器187.5比例积分微分控制器188上机考试189实习总结201实习目的和意义随着电子工业的飞速发展和电子计算机技

7、术的广泛应用，促使电子设计自动化技术日新 月异，社会对电气信息类人才的需求也日益迫切、要求也更加严格。本次matlab实习旨 在让学生们熟悉matlab语言的使用、了解matlab在电子信息课程中的应用、学握 matlab矩阵输入、运算以及matlab数值的运算功能、掌握各种数据的创建、访问、扩 建及缩减、了解利用matlab计算系统响应的方法；验证信号与系统的基本概念、基本理 论，掌握信号与系统的分析方法、了解matlab软件的基本功能及实际操作方法，熟悉利 用matlab编写简单的程序，解决一般性的问题，提高学生实际仿真控制和动手操作的综 合能力，使学生在完成设计能力基础上，提升计算机使用

8、能力，培养学生的综合设计以及仿 真控制能力，为学生今后的学习和工作打下良好的基础。2 matalab/simulink 简介2.1 matalab的发展历程matalab是mathworks公司推出的一个功能强大的计算仿真软件，是目前世界上应 用最广泛的计算机仿真软件。它最早出现于1980年，美国新墨西哥大学计算机科学系主任 cleve moler教授采用fortran语言编写了集命令翻译、科学计算于一身的一套交互式软件 系统，设计初衷是为了方便学生解决“线性代数”课程的矩阵运算问题。这个软件系统被命名 为matalab,是matrix laboratory的缩写，译为“矩阵实验”，表明其基本

9、操作单元是矩阵。 这就是最初的matalab02.2 matalab 的影响在欧美各高等院校，matalab已经被正式列入研究生和本科生教学计划，成为线性代 数、数值分析、数理统计、口动控制理论、数字信号处理、动态系统仿真、图像处理等课程 的基本教学工具，是大学生必须掌握的基本技能之一。2.3 matalab与控制系统仿真（1）强大的运算功能matlab提供了向量、数组、矩阵、复数运算，以及求解高次微分方程、常微分方程的 数值积分等强大的运算功能，这些运算功能使控制理论及控制系统中经常遇到的计算问题得 以顺利解决。简单易学的编程语言matlab的编程语言是脚本语言，这种解释性的语言简单易学。m

10、atlab命令也与数 学中的符号、公式非常接近，可读性强，容易掌握。(3) 大量配套工具箱matlab具有大量与控制系统设计相关的配套工具箱,如控制系统工具箱、系统辨识工 具箱、鲁棒控制工具箱、模糊控制工具箱、神经网络工具箱、最优化工具箱、模型预测控制 工具箱和多变量频域设计工具箱等。这些工具箱使得控制系统的仿真与计算变得便捷与高效。(4) 强大的图形功能除了一般的数据显示，matlab还支持多种形式的二维/三维图形显示，丰富的绘图命令 可以随时将计算结果可视化，使数据内容清晰可见、一目了然，便于对控制系统的数据处理结 果进行分析。(5) 高效的编程效率matlab内具有丰富的库函数，从加减乘

11、除、正弦、余弦、积分、微分、方程求解和矩 阵求逆，到快速傅里叶变换等一应俱全，而且可以直接调用，不必将其子程序的命令或语句逐 一列出，大大提高了编程效率。2.4 simulink 简介1990年math works公司在matlab中加入了新的控制系统模型化图形输入与仿真工 具，并命名为simulabo该工具很快在控制工程领域获得了广泛的认可，并在1992年被正式 更名为simulinkosimulink是matlab屮用于动态系统建模和仿真的一个软件包，它的出现使得控制系 统仿真进入模型化图形组态阶段，控制系统的分析与设计变得更加便捷和直观。simulink与 matlab语言相比,区别是其

12、与用户的交互接口是基于windows的模型化图形输入,其结果 是使得用户可以把更多的时间和精力投入到系统模型的构建,而非语言的编程上。3 matalab语言基础(1) 语言简洁紧凑，运算符和库函数极其丰富，使用方便灵活,编程效率高，matlab除了 提供和c语言一样的运算符号外,还提供了大量的矩阵和向量运算符，灵活使用matlab的 运算符可使程序变得极为简短。matlab程序书写形式自由，利用丰富的库函数避开了繁杂 的子程序编程任务，压缩了一切不必要的编程工作。(2) 具有结构化的控制语句，如for循环、while循环、break语句、if语句和switch语 句等，同时又有面向对象的编程特

13、性0(3) 图形功能强大。具有对二维和三维数据可视化、图像处理、动画制作等绘图命令， 也包括可以修改图形及编制完整图形界面的绘图命令。(4) 功能强大的工具箱。工具箱可分为两类:功能性工具箱和学科性工具箱。功能性工 具箱主要用来扩充其符号计算、图示建模仿真、文字处理以及与硬件实时交互等功能，对用 于多种学科;学科性工具箱专业性比较强，包括优化工具箱、统计工具箱、控制工具箱、小波 工具箱、图像处理工具箱、通信工具箱等。(5) 源程序的开放性。除内部函数外，所有matlab的核心文件和工具箱文件都是可读 可改的源文件，用户可修改源文件和加入自己的文件，它们可以与库函数一样被调用。3.1 matal

14、ab基础知识matalab数据类型主要包括：数字、字符串、矩阵、单元型数据及结构数据等。3.1.1变量与常量与常规的程序设计语言不同，matlab语言中的变量既不需要事先定义,也不需要预先 指定变量类型，matlab会自动依据所赋予变量的值或对变量所进行的操作来识别变量的 类型。在赋值过程中，如果赋值变量己存在时，matlab语言将使用新值代替旧值，并以新值 类型代替旧值类型。常量表示数值ans计算结果pi圆周率eqs浮点运算的相对精度inf正无穷大nan表示不定值rcalmax/rcalmin最大的正实数/最小的正实数i，j虚数单位3.1.2字符串字符是matlab中符号运算的基本构成单元，

15、也是文字等表达方式的基本元素。字符 串用单撤号进行输入或赋值，也可以用函数char ()来生成。字符串的每个字符(包括空格) 都是字符数组的一个元素。3.1.3矩阵生成在matlab中，矩阵的输入必须以方括号“”作为其开始与结束标志，矩阵的行与行之 间要用分号“;''或按冋车键分开，矩阵的元素之间要用逗号丁或用空格分隔。矩阵的大小可以 不必预先定义，且矩阵元素的值可以用表达式表示。建立矩阵的方法有直接输入矩阵元素、 现有矩阵基础上添删元素、读取数据文件、直接建立特殊矩阵等。3.1.4矩阵的运算矩阵的加、减运算符分别为只有维数相同的矩阵才可以进行加、减运算。两个矩 阵的加减运算是

16、对应元素的加减，而矩阵与标量的加减运算则是矩阵中的每一个元素都与该 标量进行加减运算。矩阵的乘法运算符为“杯，只有当两个矩阵屮前一个矩阵的列数和后一个矩阵的行数相 同时，才可以进行乘法运算。标量与矩阵的乘法运算是标量与矩阵中的每一个元素进行相乘 的运算。矩阵的除法有两种形式:左除和右除对于矩阵a和b,如果a矩阵是非奇异方 阵,ab是a的逆矩阵乘b,即inv(a)*b.ab运算等效于求a*x=b的解;而b/a是b乘a的 逆矩阵，即b*inv(a),b/a等效于求x*a=b的解，右除b/a也可由b/a=(a，b)左除来实现。 如果a是奇异矩阵系统将给出警告信息。矩阵的幕运算符为ua，aap意思是a

17、的p次方。如果a是一个方阵,p是一个大于1的 整数，则ap表示a的p次幕，即a自乘p次。如果p不是整数，则矩阵的乘方是计算矩阵a 的各特征值和特征向量的乘方。如果b是方阵,a是标量,ab就是一个按特征值与特征向量 的升幕排列的b次方程阵。如果a和b都是矩阵，则/b是错误的。3.1.5用ploy创建多项式»a=1 2 3;3 4 5;5 6 7a=123345567»p=ploy(a)»ezplot (ppa, -50,50)运行结杲如下>fiture i图313.1.6多项式的曲线拟合»x=0 3.3 4 5 6 64 7 7.4 8 8.6 9

18、10»y=6.016 5.608 5.359 5.360 5.470 5.420 5.350 5.338 5.260 5.205 5.110 5.000;»n=3»p=ployfit(x,y,n);»xi=linspacc(0, 10 ,1000)»z=ployval(p,xi);>>plot(x,y,'b',x,y,t,xi,z,':g')运行结果如下图3-2(1)单曲线语句格式plot(x,y)淇中x和y为长度相同的向量，分别为储存x坐标和y坐标数据。>>y=0,0.4&0.8

19、4,1,0.91,0.6,0.14;»plot(y)生成图33所示的折线图»x=linspace(0,2*pi,30);»y=sin(x);»plot(x，y)生成图形如图34所示图33图34(2) 多重线在同一个画面上可以画出许多曲线，语句格式为plot (xl, yl, x2, y2,，xm yn)o»t=0:pi/100:6;»y=sin(t);»yl=sin(t+0.25);»y2=sin(t+0.5);»plot(t,y,t,yl,t,y2)得到图形如图35所示另一种方法是利用hold命令，在已

20、经画好的图形上，若设置hold on,matlab将把新 的plot命令产生的图形画在原來的图形上。而命令hold off将结朿这个过程。»t=0:pi/100:6;y=sin(t);plot(x,y)»hold on»y 1=sin(t+0.25);plot(t,y 1)»y2=sin(t+0.5);plot(t,y2)»hold off得到的图形如图36所示图3-5图3-6(3) 三维曲而图与三维网线图的区别是网线图的线条有颜色，空挡是黑色的，而曲面图的线条是黑色的, 空挡是有颜色的。具体图形如图37所示。fimtv q*l* <j

21、d k > s * x a u 口图3-74 simulink仿真工具在工程实际中，控制系统的结构往往很复杂，如果不借助专用的系统建模软件，则很难 准确地把一个控制系统的复杂模型输入计算机，对其进行进一步的分析与仿真。4.1运行simulik演示程序图4-1如图所示为直流无刷电机转速的pi控制系统模型，其中标注为vdc、vab、te等的模块 实际上实现的就是示波器的功能，双击该模块，即可以打开示波器观测相应的数据。在进行仿 真之前，首先设置仿真参数，图41所示的仿真采用原模型默认的参数。选择菜单栏simulation 下的start命令，或者单击simulink工具栏上的开始按钮p,系统

22、开始按照模型中设置的参数 进行仿真，仿真结果:电压、电流和转矩等曲线将显示在示波器屮。若要停止仿真，可选择 simulation菜单栏下的sop命令，或者单击simulink工具栏上的停止按钮照仿真结束后，可选 择file菜单栏下的close命令关闭模型咸者单击窗口右上角的区來关闭当前窗口。4.1 simulink模型的建立4.1.1模型窗口由于simulink是基于matlab环境之上的高性能的系统级仿真设计平台，因此启动simulink之前必须首先运行matlab,然后才能启动simulink并建立系统系统模型。 在matlab的命令窗口屮选择file->newt new model

23、菜单项。 单击simulink工具栏中的“新建模型''图标。 选中simulink菜单系统屮的file>ncw> model菜单项 还可以使用new - system命令来建立新模型。4.1.2模型举例图42图43图444.2 simulink的仿真方法建立好了 simulink模型后就可以启动仿真过程了。最简单的方法当然是按simulink工 具栏下的【启动仿真】按钮了。启动仿真过程后将以默认参数为基础进行仿真，而用户还可 以白己设置出需要的控制参数，打开simulation> simulation parameters菜单项，用八可以从 屮填写相应的数据，

24、控制仿真过程。在对话框屮有5个标签，默认的标签为微分方程求解程序 solver的设置，在该标签下的对话框主耍接受微分方程求解的算法的选择及仿真控制参数的 设置。设置仿真参数和选择解法器，选择simulation菜单下的parameters命令，就会弹出一个仿 真参数对话框，它主要用三个页面来管理仿真的参数。(1) solver页:它允许用户设置仿真的开始和结束时间，选择解法器，说明解法器参数及选择 一些输出选项。(2) workspacel/o页:其作用是管理模型从matlab工作空间的输入和对它的输出。(3) diagnostics页:允许用户选择simulink在仿真中显示的警告信息的等级

25、。仿真结果可以用数据的形式保存在文件屮，也可以用图形的方式直观地显示出来,查看和 分析结果曲线对于了解模型的内部结构，以及判断结果的准确性具有重要意义。采用以下方 法可绘制模型的输出轨迹。(1) 利用示波器模块(scope)得到输出结果。(2) 利用输出接口模块(out)得到输出结果。(3) 通过将数据传送到工作空间模块(to workspace)得到输出结果。5控制系统数学模型通常,实际控制系统的动态性能和稳态性能是无法进行准确的分析。只有对实际系统的 各个组成部分进行数学建模，通过对实际控制系统的数学抽象，将其内部变量间的关系表示出 来，才能对系统进行仿真与计算,并将仿真分析的结杲用以设计

26、系统并优化控制器，使系统的 性能指标符合实际控制系统的预期设计需要。控制系统的数学模熨分为静态数学模熨和动态 数学模型。静态数学模型是指描述各阶导数为零的变量间关系的代数方程。动态数学模型是 指描述各阶导数不为零的变量间关系的代数方程。线性系统常用的数学模型有传递函数模型、 状态方程模型、零极点模型和部分分式模型等。不同的应用需要不同的数学模型，它们之间 存在着内在联系并且可以互相转换。5.1动态过程微分方程描述分析和设计控制系统首先要对实际控制系统进行数学建模,线性定常系统的基础模型即使微分方程，它是描述事物最基本的工具。5微分方程的形式线性定常系统或元件的微分方程形式为dc(t)dt+ %

27、二° dtm 1 d严心 dt m )等式左侧表示输出变量的各阶导数，右侧表示输入变量的各阶导数。通常n>m,m>o,n>l 且an0, tuogj及bj均为实数)。5.1.2传递函数定义线性定常系统在零初始条件下，系统输出量的拉普拉斯变换与输入量拉普拉斯变换之比, 称为该系统的传递函数。定义屮的零初始条件有两方而的含义：(1) 输入作用是在t=0时刻后作用于系统，因此输入量及其各阶导数在t=0时的值为 零。(2) 输入信号作用于系统之前系统是静止的，即a0时系统的输出量及各阶导数为零。根据传递函数的定义，则系统的传递函数为：g( ；)= cg) _ 畧 + 处心

28、+ + 怙 + bm = m g)r(s) aqsm + 空心 + + d心 s + 如n(s)m(5)= bosm + 勺严 + + bm_xs + bmn(s) = aqsn + hf an_xs + anc(s) = g(s) * r(s)输入量r (s)经过传递函数g (s)的传递后，得到了输入量c (s),可以用框图形象c(s)r(s)直观地表示11!这一传递关系，箭头表示信号传递的方向。> g (s)5.2高阶系统的时域分析通常把三阶以上的系统称为高阶系统。在实际控制系统中,几乎都是由高阶微分方程來 描述，即均为高阶系统。对于高阶系统的建模及分析比较困难，通常采用突出主要因素

29、,忽略次 要因素的方法，将高阶系统近似为一个二阶系统。5.2.1高阶系统的特点(1) 高阶系统的时间响应表达式rh简单函数组成。(2) 如果闭环极点都具有负实部，高阶系统则是稳定的。(3) 高阶系统吋间响应的类型取决于闭环极点的性质和大小，响应曲线形状与闭环零点 有关。5.2.2高阶系统的闭环主导极点系统输出各动态响应分量衰减快慢取决于对应的闭环极点距离s平面虚轴的远近，其中 最靠近虚轴的闭环极点所对应的动态分量衰减得最慢，在所有各分量中起主要作用。如果高 阶系统中，所有英他极点的实部比距离虚轴最近的闭环极点的实部大5倍以上，并且在该极点 附近不存在闭环零点，则这种离虚轴最近的闭环极点将对系统

30、的动态响应起主导作用,并称其 为闭环主导极点主导极点常以共觇复数形式出现，此时可用二阶系统的动态响应指标来估计 高阶系统的性能。所以，主导极点具有重耍的实用意义。6控制系统分析分析控制系统首先要建立数学模型，英次再应用不同的方法分析控制系统各项性能。在 经典控制理论中分析方法有很多种，主要有时域分析法、根轨迹分析法、频域分析法等。不 同的分析方法具有各自的优点，它们的适用范围和对象也有所不同。时域分析法主要以拉普 拉斯变换为数学工具直接在时间域研究控制系统动态性能和稳态性能的方法。时域分析法的 优点是可以形象、直观、准确地得到系统时间响应的全部信息。根轨迹法是一种简便的图解 分析方法，尤英对于

31、多回路系统分析更加有效。频域分析法是在正眩信号作用下分析系统性 能的图解方法，不仅适用线性定常系统，还可以推广应用于非线性系统。系统的频域特性可以 通过分析法和实验法定性与定量地对各项性能进行分析。6.1时域分析时域分析法就是通过求解控制系统的时间响应來分析系统的稳定性、快速性和准确性。 这是-种在吋间域对系统进行分析的方法，具有形象、直观、准确等特点。尤其适用于对二 阶控制系统的各项性能进行分析。但是由于时域分析法计算繁琐，因此对于高阶系统此方法 并不适用。6.1.1系统传递函数如下：g=10/s4+4s3+18s2+20s+10.试着绘制其阶跃响应曲线。解：matlab程序如下：num=1

32、0;den=l 4 18 20 10;t=0:0.1:10;y=step(num,den,t)plot(t,y)xlabel(4 时间 t,);ylabel(<y,);title。单位阶跃响应曲线j;运行程序，系统单位阶跃响应曲线如图所示。图616.1.2系统传递函数如下：g(s)=36/s2+13s+36.试用matlab绘制出单位阶跃响应、单位斜 坡响应、单位脉冲响应曲线图。解：matlab程序如下所75：num=36;den=l 13 36;t=0:0.1:10;u=t;y=step(num,den,t);y 1 =lsim(num,den,u,t);y2=impulsc(num,

33、dcn,t)plot(t,y,7,t,yl,",t,y2,y);xlabe(，时间 t,);ylabe(ty,)titled单位阶跃、斜坡、脉冲输入响应曲线j;legend (伸位阶跃响应曲线,，单位斜坡响应曲线,，单位脉冲响应曲线，)运行程序，系统的单位阶跃响应曲线、单位斜坡响应曲线、单位脉冲响应曲线如图所示.图6-26.2频域分析6.2.1已知一典型二阶传递函数为：g(s) = q；心2 + 2®, +研),试分别绘出其bode so 解:matlab程序如下所示：w=0,logspace(-2,2,200);zet=0.7;wn=0:0.25:1,1.8,2.5;ho

34、ld onfor i=l:length(wn)num=wn(i)a2;den= 1,2*zet* wn(i),wn(i)a2;sys=tfnum,dcn);bode(sys,w);endbode图如图所示。ficurc 1口回冈edi tinsert 1。“ qetktop imdow k«lpd j w j k 、5 ® 0 y <3 图6-36.2.2绘制出奈特图，具体程序如下:g=tf(100,conv(l 4 1,1 6);nyquist(g);绘制的图形如图所示。> fi<ar. iellpl网108八 jfitw xn««r

35、tfindov h«lp:d | a ® 农 x "j a | o 気 口nyrqukt digram68101214 1s 18real am图647控制器设计7.1比例控制器比例控制器是最简单的控制方法z。比例控制器的输入输出误差信号成止比关系，其 传递函数为：g(s)=kp, kp为比例系数，比例系数只能改变系统增益，不影响相位，仅采用 比例控制时，系统的输出存在稳态误差，增大kp可以提高系统的开坏增益，减小稳态误差， 但是会减低系统稳定性，甚至是闭环系统的不稳定。某控制系统中g(s)=l/(s+3)(s+2)(s+l),在控制单元使加比例控制，并且采用不同

36、的比例 系数kp=0.1,0.5,l,2,5,10,观察各比例系数下的阶跃响应及控制效果。解：matlab程序如下：kp珂0.1,0.5,1,2,5,10;go=tf( 1 ,conv(conv( 1,1 ,2,1 ),3,1 );for i=l:6g=fbcdback(go.*kp ,1);step(g);hold on;endgtext('kp=0.1 ,);gtest(,kp=0.5,);gtest(,kp=l,)gtest('kp=2,);gtest(,kp=5');gtest(,kp=10,);运行程序得到不同比例系数下的系统单位阶跃响应曲线，如图所示。ib1

37、.4step respxise图7-1口回冈'1監4 八 yitv insertdesktop findovk o ® d q 口从图中可以看出，随着比例系数的增大，系统响应速度加快，稳态误差减小，超调量却 在增加，调节时间变长，比例系数增大到一定程度后，系统最终会变的不稳定。积分控制器的传递函数为：g(s) = 某控制系统中g=l/(s+l)(2s+l),在控制单元使加积分控制器，观察系统稳态误差的 变化。解：matlab程序如下所示：go=tf(l,conv(l,l,2,l);gc=tf(l,l,o);subplot(2,l, 1 );step(feedback(go, 1);titled加积分控制前j;subplot(2j ,2);step(feedback(go*gc, 1); titled加积分控制后j;运行程序得到加积分控制前后系统的单位阶跃响应曲线，如图所示。t&