工程控制实验报告..doc

机械工程控制基础 实验报告 班级：09机设5班 学号：0916101179 姓名：孙 婧 玮 2011年 12 月 27 日 概述： MATLAB是Math Works公司的软件产品，是一个高级的数值分析、处理和计算的软件，其强大的矩阵运算能力和完美的图形可视化功能，使得它成为国际控制界应用最广泛的首选计算机工具。 MATLAB具有良好的的可扩展性，其函数大多数为ASCII文件，可以直接进行编辑、修改；其工具箱可以任意增加，任何人可以生成自己的MATLAB工具箱。因此，很多研究成果被直接做成MATLAB工具箱发表。 SIMULINK是基于模型化图形的动态系统仿真软件，是MATLAB的一个工具箱，它使系统分析进入一个崭新的阶段，它不需要过多地了解数值问题，而是侧重于系统的建模、分析和设计。其良好的人机界面及周到的帮助功能使得它广为科技界和工程界采用。 因此，本试验将尽可能把MATLAB和SIMULINK工具应用于控制系统的分析和计算中。 试验一：用MATLAB进行部分分式展开 1、试验目的： （1）对MATLAB进行初步的了解； （2）掌握应用MATLAB对高阶函数进行部分分式的展开。 2、试验学时：2学时 3、试验方法：MATLAB有一个命令用于求B（S）/A（S）的部分分式展开式。 设S的有理分式为 F（S）=B（S）/A(S)=num/den=(b0Sn+b1Sn-1+…+bn)/(Sn+a1Sn-1+…+an) 式中ai和bi的某些值可能是零。在MATLAB的行向量中，num和den分别表示F(S)分子和分母的系数，即 num=[b0 b1 … bn] den=[1 a1 … an] 命令 [r,p,q]=residue(num,den) MATLAB将按下式给出F（S）部分分式展开式中的留数、极点和余项： 有：r(1)、r(2)…r(n)是函数的留数；p(1)、p（2）、…p(n)是函数的极点；K（s）是函数的余项。 4、试验内容 （1）试求下列函数的部分分式的展开式，写出程序语句和部分分式的结果 程序语句：?? num=[1? 11? 39? 52? 26] ???????? den= [1? 10? 35? 50? 24] 命令?????? [r,p,k]=residue(num,den) 结果： ? [r,p,k]=residue(num,den)?? ???????? r= ?????????????? 1.0000 ?????????????? 2.5000????????? ????????????? -3.0000 ?????????????? 0.5000 ??????????? p= ????????????? -4.0000 ????????????? -3.0000 ????????????? -2.0000 ????????????? -1.0000 ??????????? k= 1 则得 如果F(s)中含重极点，则部分分式展开式将包括下列诸项 ? 式中，p(j)为一个q重极点。 程序语句：?????????? num=[0?? 1?? 4?? 6] ??????????? den =[1?? 3?? 3?? 1] 命令 ?????????????? [r,p,k]=residue(num,den) ?： [r,p,k]=residue(num,den) ?????????????? r= ???????????????? 1.0000 ???????????????? 2.0000 ???????????????? 3.0000 ????????????? ?p= ???????????????? -1.0000 ???????????????? -1.0000 ???????????????? -1.0000 ??????????????? k= ????????????????? [? ] 所以可得 注意，本例的余项k为零。 ? 求该系统的单位脉冲响应、单位阶跃响应和单位斜坡响应和单位加速度响应。 分别写出相对应的程序语句，并给出结果: A)单位脉冲响应： 程序语句：G=tf([2,20,50],[1,15,84,223,309,240,100]) impulse(G) Transfer function: 2 s^2 + 20 s + 50 --