电气与自动化工程学院综合实验报告系统仿真综合实验报告

1、合肥工业大学电气与自动化工程学院综合实验报告实 验 名 称:系统仿真综合实验姓 名:学 号:专 业 班 级:实 验 地 点:指 导 教 师:成 绩:日 期:2012年7月实验一 matlab基本操作实验目的1熟悉matlab实验环境，练习matlab命令、m文件、simulink的基本操作。2利用matlab编写程序进行矩阵运算、图形绘制、数据处理等。3利用simulink建立系统的数学模型并仿真求解。实验原理matlab环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。matlab有3种窗口，即：命令窗口（the command window）、m-文件编辑窗口（the edit

2、window）和图形窗口（the figure window），而simulink另外又有simulink模型编辑窗口。1命令窗口（the command window）当matlab启动后，出现的最大的窗口就是命令窗口。用户可以在提示符“”后面输入交互的命令，这些命令就立即被执行。在matlab中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。在命令窗口中输入该文件名，这一连串命令就被执行了。因为这样的文件都是以“.m”为后缀，所以称为m-文件。2m-文件编辑窗口（the edit window）我们可以用m-文件编辑窗口来产生新的m-文件，或者编辑已经存在的m-文件。在m

3、atlab主界面上选择菜单“file/new/m-file”就打开了一个新的m-文件编辑窗口；选择菜单“file/open”就可以打开一个已经存在的m-文件，并且可以在这个窗口中编辑这个m-文件。3图形窗口（the figure window）图形窗口用来显示matlab程序产生的图形。图形可以是2维的、3维的数据图形，也可以是照片等。 matlab中矩阵运算、绘图、数据处理等内容参见教材自动控制系统计算机仿真的相关章节。simulink是matlab的一个部件，它为matlab用户提供了一种有效的对反馈控制系统进行建模、仿真和分析的方式。有两种方式启动simulink:1在command w

4、indow中，键入simulink，回车。2单击工具栏上simulink图标。启动simulink后，即打开了simulink库浏览器（simulink library browser）。在该浏览器的窗口中单击“create a new model（创建新模型）”图标，这样就打开一个尚未命名的模型窗口。把simulink库浏览器中的单元拖拽进入这个模型窗口，构造自己需要的模型。对各个单元部件的参数进行设定，可以双击该单元部件的图标，在弹出的对话框中设置参数。 实验内容1 用matlab可以识别的格式输入下面两个矩阵 再求出它们的乘积矩阵c，并将c矩阵的右下角23子矩阵赋给d矩阵。赋值完成后，调

5、用相应的命令查看matlab工作空间的占用情况。答案：1） 程序：a=1 2 3 3;2 3 5 7;1 3 5 7;3 2 3 9;1 8 9 4;b=1+4i 4 3 6 7 8;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;c=a*bd=c(4:5,4:6)结果：c = 1.0e+002 * columns 1 through 4 0.1400 + 0.0400i 0.5200 + 0.2100i 0.5100 0.4000 0.2500 + 0.0800i 1.0300 + 0.3500i 1.0300 0.7700 0.2400 + 0.0400

6、i 0.9900 + 0.3500i 1.0000 0.7100 0.2200 + 0.1200i 1.0800 + 0.2100i 1.1100 0.8200 0.3900 + 0.0400i 1.1400 + 0.6300i 1.0800 0.9300 columns 5 through 6 0.4100 0.3100 + 0.0400i 0.7700 0.5900 + 0.0600i 0.7000 0.5100 + 0.0600i 0.7900 0.6500 + 0.0400i 0.9900 0.7000 + 0.1600id = 82.0000 79.0000 65.0000 + 4.

7、0000i 93.0000 99.0000 70.0000 +16.0000i工作空间:name value a b c d 82.000000000000000 + 0.000000000000000i,79.000000000000000 + 0.000000000000000i,65.000000000000000 + 4.000000000000000i;93.000000000000000 + 0.000000000000000i,99.000000000000000 + 0.000000000000000i,70.000000000000000 + 16.0000000000000

8、00i小结：matlab通过确认下标，可以对矩阵进行插入子块、提取子块和重排子块的操作。如果提取子块时，n或m是常数，则返回指定的行列；如果n或m是向量，则返回的是指定矩阵的子块。2 分别用for和while循环结构编写程序，求出答案：程序1：s=0;for k=0:63; s=s+2k;enddisp(the sum is),s结果：the sum iss = 1.8447e+019程序2： s=0;k=0;while k nyquist(g)奈氏曲线不包围（-1，j0），故系统稳定再输入：margin(g)从图上可以看出：幅值裕度gm=12db,相角裕度pm=41.5deg小结：可以根据乃

9、奎斯特曲线和伯德图来判断闭环系统是否稳定；利用margin不仅可以得到伯德图，还可以得到幅值裕度和相角裕度4. 已知开环传递函数令k1作伯特图，应用频域稳定判据确定闭环系统的稳定性，并确定使系统获得最大相位裕度的增益k值。答案：clear g=tf(1 1,0.1 1 0 0); margin(g),grid由图：幅值裕度gmpm1; pm1=pm; endenddisp(获得最大相角裕度的k值)kdisp(最大相角裕度)pm1结果：获得最大相角裕度的k值k = 3.1600最大相角裕度pm1 = 549032小结：利用for语句以及if语句，可以判断出获得最大相角裕度的k值 5单位负反馈系统

10、的开环传递函数为绘制根轨迹曲线，并求出使闭环系统稳定的k值范围。使用rltool（sys）函数观察不同k值的阶跃响应情况，分析k值变化对系统响应有何影响。答案：num=-0.5 1;den=conv(conv(0.5 1,0.2 1),0.1 1);g=tf(num,den);rlocus(g),grid,k,p=rlocfind(g)select a point in the graphics windowselected_point = -9.9052 - 0.0155ik = 0.0063p = -10.0918 -4.8540 -2.0542再输入：k,p=rlocfind(g)sel

11、ected_point = -0.0474 + 3.7112ik = 1.3008p = -16.9195 -0.0402 + 3.6874i -0.0402 - 3.6874i在输入rltool(g)小结：通过比较，k值越大，系统地单位阶跃响应波动就越大，震荡越激烈，当k值大于临界值时，系统变得不稳定6. 非最小相位系统的开环传递函数如下 绘制频率特性曲线，并解释为什么这样的系统被称为“非最小相位”系统，试从其频率特性加以解释。在simulink环境中建立对应闭环系统的仿真模型，观察它们的阶跃响应。答案：在右半s平面没有零极点的传递函数，是最小相位系统，反之便称为“非最小相位”系统。num1

12、=-6 24;den1=0.05 0.6 1 0 0;g=tf(num1,den1); margin(g),grid在输入：num2=10 -60 110 60;den2=1 17 82 130 100;g1=tf(num2,den2); margin(g1),grid仿真：小结：通过本题可以看出，非最小相位系统可能稳定，也可能不稳定，所以不能用频域判断法来判断。7单位负反馈系统的开环传递函数为：，其中 利用simulink工具给出simulink仿真框图及阶跃响应曲线。答案：transport delay 参数设置: time delay设置成0.5transfer fcn参数设置:nume

13、rator coefficients 设置成1denominator coefficients 设置成 1 2 2仿真模型：仿真曲线：小结：通过延时环节，可以使系统延时产生响应实验三 pid控制器的设计实验目的研究pid控制器对系统的影响。实验原理1模拟pid控制器控制系统结构如图所示模拟pid控制器的表达式为：式中，、和分别为比例系数、积分系数和微分系数。 对上式进行拉普拉斯变换，整理后得到连续pid控制器的传递函数 显然，、和三个参数一旦确定（注意：，），pid控制器的性能也就确定下来。为了避免纯微分运算，通常采用近似的pid控制器，其传递函数为 2. 离散pid控制器 如果采样周期为t，

14、在第k个采样周期的导数可近似表示为； 在k个采样周期内对的积分可近似表示为； 因此，离散pid控制器的表达式为； 离散pid控制器的表达式可简化为； 离散pid控制器的脉冲传递函数为；式中，、和分别为比例系数、积分系数和微分系数。实验内容1已知三阶对象模型，利用matlab编写程序，研究闭环系统在不同控制情况下的阶跃响应，并分析结果。(1) 时，在不同kp值下，闭环系统的阶跃响应；(2) 时，在不同值下，闭环系统的阶跃响应；(3) 时，在不同值下，闭环系统的阶跃响应；答案：（1）num=1;den=conv(conv(1 1,1 1),1 1);gk=tf(num,den); for kp=0

15、.1:0.4:0.9 sys=feedback(kp*gk,1,-1); step(sys); hold on pauseendtitle(不同kp下的阶跃响应)xlabel(时间（秒）)ylabel(幅值)gtext(kp=0.1),gtext(kp=0.5),gtext(kp=0.9)小结：当kp取值越小时系统的响应较慢，而当kp取值较大时，响应加快，但超调量增加（2）num=1;den=conv(conv(1 1,1 1),1 1);gk=tf(num,den); for ti=0.1:1:2.1 gc=tf(ti 1,ti); sys=feedback(gc*gk,1,-1); ste

16、p(sys); hold on pauseendtitle(不同ti下的阶跃响应)xlabel(时间（秒）)ylabel(幅值)gtext(ti=0.1),gtext(ti=1.1),gtext(td=2.1) 小结：当ti较小时，进入稳态的速度较慢，当ki取值较大时，进入稳态的速度较快，但相应的，超调量也增加（3）num=1;den=conv(conv(1 1,1 1),1 1);gk=tf(num,den); for td=0.1:1.5:3.1 gc=tf(td 1 1,1); sys=feedback(gc*gk,1,-1); step(sys); hold on pauseendti

17、tle(不同td下的阶跃响应)xlabel(时间（秒）)ylabel(幅值)gtext(td=0.1),gtext(td=1.6),gtext(td=3.1)小结：td取值较大时，系统对变化去世的调节变慢，超调量较大；当td取值较大时，对变化趋势的调节加强，阶跃响应的初期会出现尖脉冲*2. 已知被控对象为一电机模型，传递函数为，输入信号为，采用pid控制方法设计控制器，利用matlab编程进行仿真，绘制具有pid控制器的控制系统正弦跟踪曲线。答案：参数设置：小结：输出曲线比输入曲线的幅值更小，表明系统变的更稳定3. 完成自动控制系统计算机仿真教材第7章中的例7-4题。【例7-4】 某直流电机速

18、度控制系统如图所示，采用pid控制方案，使用期望特性法来确定 、 和 这三个参数。建立该系统的simulink模型，观察其单位阶跃响应曲线，并且分析这三个参数分别对控制性能的影响。解： 使用期望特性法来设计pid控制器。假设pid控制器的传递函数为系统闭环的传递函数为不妨假设希望闭环极点为：-300，-300，-30+30i,-30-30i则期望特征多项式为: 对应系数相等，可求得:kp=4.4156,ki=119.34,kd=0.067仿真模型：参数值如图：（2）kp=0.5: kp=5kp=20（3）ki=20ki=120ki=300(4)kd=0.01kd=0.07kd=0.24. 完成自动控制系统计算机仿真教材第7章中的例7-5题。【例7-5】 如图7-19所示的系统，被控对象为一个带有延迟的惯性环节，试用ziegler-nichols经验整定公式，计算pid控制器的参数，并且绘制其仿真系统单位阶跃响应曲线。解： 由该系统传递函数可知， 由ziegler- nichols经验整定公式， 可得：pid控制器的传递函数为： 参数设置：仿真模型：仿真曲线小结：根据ziegler-nichols经验整定公式是针对被控对象模型为带有延迟的一阶惯性传递函数提出，