武汉工程大学实验报告-数字仿真

实验名称：Matlab 的基本操作与编程一、实验目的：1）熟悉MATLAB 软件的运行环境和基本操作2）掌握MATLAB 矩阵的输入方式、元素的提取与组合 3）掌握数值运算。

4）掌握MATLAB 软件的绘图功能 5）掌握M 函数的编写。

二、实验内容：1）启动MATLAB 软件，观察其界面组成及操作方法，了解各部分的功能 2）使用基本的MATLAB 命令，并观察记录执行结果帮助、查询信息类命令：Demo 、help 、who 、whos 显示、记录格式等命令：clc 、clear 、format 尝试一下其他的命令（dos 命令）3）生成一个5阶魔方矩阵，并提取其第（3、4、5）行，第（2、3、4）列构成的新的矩阵5）用命令行方式求解下式的值42cos lim22x x e x -→（提示使用syms x 定义一个符号，使用limit 函数）6）MATLAB 的绘图（1） 二维绘图命令plot ：画出,sin x y =在]2,0[π∈x 上的图形（2） 三维绘图命令plot3： 画出三维螺旋线⎪⎩⎪⎨⎧===t z t y t x cos sin ，]4,0[π∈t 的图形．mesh 命令：绘制)2(22y x e z +-=，在]5,5[-∈x ，]5,5[-∈y 区间的曲面7）编写M 函数利用程序流程控制语句编写一个函数myfactorial （n ），实现n ！（阶乘）。

要求使用help 命令可以列出相关的帮助信息。

三、实验结果及分析2、Demohelp3、（1）实验程序：a=magic(5)运行结果：a =17 24 1 8 15 23 5 7 14 164 6 13 20 2210 12 19 21 311 18 25 2 9（2）实验程序：a(3:5,2:4)运行结果：ans =6 13 2012 19 2118 25 25、实验程序：syms x;limit((cos(x)-exp(x*x/2)/2)/4)运行结果：ans =1/86、（1）实验程序：x = 0 : pi/20 : 2*pi;y = sin(x);plot(x, y);实验结果：（2）实验程序：x = 0 : pi/20 : 2*pi;y = sin(x);plot(x, y);x=sin(t);x=sin(t); y=cos(t);t=0:pi/20:2*pi;plot3(x,y,t);实验结果：实验程序：x=-5:0.5:5;y=x;[X,Y]=meshgrid(x,y);R=(X.^2+Y.^2)/((-1)*4);Z=exp(R);surf(X,Y,Z)实验结果：7、实验程序：function[x]=myfactorial(n)x=1;for a=1:n;x=x*a;end运行举例：>> myfactorial(2)ans =2实验名称：典型闭环系统的数字仿真及计算机解题一、实验目的：1）熟悉典型闭环的仿真过程2）掌握MATLAB编程实现典型闭环环节仿真3）利用典型闭环环节仿真程序解题。

4）掌握MATLAB下对控制系统进行时域、频域和根轨迹的分析二、实验内容：1）编写典型环节阶跃响应函数典型环节冲击响应函数function [yout,t] = my_step(num,den,v,t0,tf,h,R,n)输入参数：num：传递函数的分子系数向量den：传递函数的分母系数向量v：反馈比例系数t0：仿真起始时间tf：仿真终止时间h：仿真步长R:阶跃幅值n:系统阶次输出参数：yout：响应输出t：时间向量2）用上述函数分析以下系统，同时用simulink分析该系统，并比较其结果。

3）被控对象的传递函数为)20030(400)(2++=S S S s G ，用simulin 建模并分析其单位阶跃响应。

用MATLAB 命令绘出其伯德图和根轨迹图。

三、实验结果及分析 2)结果：3）仿真结果：Bodenum=[0,0,0,400];>> den=[1,30,200,0];>> bode(num,den)>>根轨迹>> num=[0,0,0,400]; >> den=[1,30,200,0]; >> rlocus(num,den)实验名称：双闭环控制直流电动机调速系统建模与仿真分析一、实验目的：1.熟悉Matlab/Simulink仿真环境；2.掌握Simulink图形化建模方法；3.验证“直流电动机转速/电流双闭环PID控制方案”的有效性。

二、实验内容：（详见后面的实验指导）1.“双闭环直流电动机调速系统”的建模2.电流环／调节器设计3.电流环动态跟随性能仿真实验4.转速环／调节器设计5.转速环动态抗扰性能仿真实验6.系统动态性能分析（给出仿真实验结果与理论分析结果的对比／分析／结论）三、实验结果及分1、“双闭环直流电动机调速系统”的建模仿真波形2、电流环模型及其环动态跟随性能仿真实验伯德图[num,den]=linmod('current_loop')sys=tf(num,den)margin(sys)[mag,phase,w]=bode(sys);奈奎斯特图[num,den]=linmod('current_loop')sys=tf(num,den)margin(sys)[gm,pm,wcg,wcp]=margin(mag,phase,w) Nyquist(sys)阶跃响应曲线[num,den]=linmod('current_loop')sys=tf(num,den)margin(sys)Step(sys)运行程序得到相关参数：Transfer function:-1.164e-010 s^3 + 2.038e007 s^2 + 4.977e009 s + 2.136e011 -----------------------------------------------------------------------s^5 + 1043 s^4 + 3.253e005 s^3 + 4.531e007 s^2 + 3.22e009 s + 8.544e010gm =4.3078pm =48.4499wcg =345.6682wcp =163.79233、转速环建模及其抗干扰性能仿真模型仿真波形参考文献：[1] 张晓华主编《控制系统数字仿真与CAD》第2版机械工业出版社2005[2] 陈伯时主编《电力拖动自动控制系统》第2版机械工业出版社2001附件：“双闭环控制直流电动机调速系统”数字仿真实验实验指导四、实验目的4.熟悉Matlab/Simulink仿真环境；5.掌握Simulink图形化建模方法；6.验证“直流电动机转速/电流双闭环PID控制方案”的有效性。

五、实验内容7.“双闭环直流电动机调速系统”的建模8.电流环／调节器设计9.电流环动态跟随性能仿真实验10.转速环／调节器设计11.转速环动态抗扰性能仿真实验12.系统动态性能分析（给出仿真实验结果与理论分析结果的对比／分析／结论）六、实验步骤1、系统建模A．控制对象的建模建立线性系统动态数学模型的基本步骤如下：（1）根据系统中各环节的物理定律，列写描述据该环节动态过程的微分方程；（2）求出各环节的传递函数；（3）组成系统的动态结构图并求出系统的传递函数。

下面分别建立双闭环调速系统各环节的微分方程和传递函数。

B．额定励磁下的直流电动机的动态数学模型图1给出了额定励磁下他励直流电机的等效电路，其中电枢回路电阻R和电感L包含整流装置内阻和平波电抗器电阻与电感在内，规定的正方向如图所示。

1图1 直流电动机等效电路由图1可列出微分方程如下：0dd d dI U RI LE dt=++ （主电路，假定电流连续） e E C n = （额定励磁下的感应电动势）2375e L GD dnT T dt-=⋅ （牛顿动力学定律，忽略粘性摩擦）e m d T C I = （额定励磁下的电磁转矩）定义下列时间常数：l LT R =——电枢回路电磁时间常数，单位为s ；2375m e mGD R T C C =——电力拖动系统机电时间常数，单位为s ； 代入微分方程，并整理后得：0()dd d ldI U E R I T dt -=+ m d d L T dE I I R dt-=⋅ 式中，/dL L m I T C =——负载电流。

在零初始条件下，取等式两侧得拉氏变换，得电压与电流间的传递函数0()1/()()1d d l I s R U s E s T s =-+（1）电流与电动势间的传递函数为()()()d dL m E s RI s I s T s=- （2）d Ua) b)Uc)图2 额定励磁下直流电动机的动态结构图2a) 式（1）的结构图 b)式（2）的结构图c)整个直流电动机的动态结构图C ．晶闸管触发和整流装置的动态数学模型在分析系统时我们往往把它们当作一个环节来看待。

这一环节的输入量是触发电路的控制电压U ct ，输出量是理想空载整流电压U d0。

把它们之间的放大系数K s 看成常数，晶闸管触发与整流装置可以看成是一个具有纯滞后的放大环节，其滞后作用是由晶闸管装置的失控时间引起的。

下面列出不同整流电路的平均失控时间：表1 各种整流电路的平均失控时间（f=50Hz ）用单位阶跃函数来表示滞后，则晶闸管触发和整流装置的输入输出关系为01()d s ct s U K U t T =⋅-按拉氏变换的位移定理，则传递函数为0()()s T s d s ct U s K e U s -= （3） 由于式（3）中含有指数函数s T s e -，它使系统成为非最小相位系统，分析和设计都比较麻烦。

为了简化，先将s T s e -按台劳级数展开，则式（3）变成02233()()12!3!s s T s d s ss T s ct s s s U s K K K e U s e T s T s T s -===++++ 考虑到T s 很小，忽略其高次项，则晶闸管触发和整流装置的传递函数可近似成一阶惯性环节0()()1d sct s U s K U s T s ≈+ （4）其结构图如图3所示。

3a)b)图3 晶闸管触发和整流装置的动态结构图 a) 准确的结构图 b)近似的结构图D ．比例放大器、测速发电机和电流互感器的动态数学模型比例放大器、测速发电机和电流互感器的响应都可以认为是瞬时的，因此它们的放大系数也就是它们的传递函数，即()()ct p n U s K U s =∆ （5） ()()n U s n s α= （6） ()()i d U s I s β= （7）E ．双闭环控制直流电动机调速系统的动态数学模型 根据以上分析，可得双闭环控制系统的动态结构图如下图4 双闭环控制系统的动态结构图2、实验系统参数系统中采用三相桥式晶闸管整流装置，基本参数如下：直流电动机：220V ，13.6A ，1480r/min ，e C =0.131V/（r/min ）， 允许过载倍数λ=1.5。