《武汉工程大学实验报告-数字仿真_17286解析》由会员分享，可在线阅读，更多相关《武汉工程大学实验报告-数字仿真_17286解析（38页珍藏版）》请在金锄头文库上搜索。
1、实验报告专业班级：姓名学号：课程名称：控制系统数字仿真与CAD实验成绩： 指导老师： 实验名称：Matlab的基本操作与编程一、实验目的：1）熟悉MATLAB软件的运行环境和基本操作2）掌握MATLAB矩阵的输入方式、元素的提取与组合3）掌握数值运算。4）掌握MATLAB软件的绘图功能5）掌握M函数的编写。二、实验内容：1）启动MATLAB软件，观察其界面组成及操作方法，了解各部分的功能2）使用基本的MATLAB命令，并观察记录执行结果帮助、查询信息类命令：Demo、help、who、whos显示、记录格式等命令：clc、clear、format尝试一下其他的命令（dos命令）3）生成一个5阶
2、魔方矩阵，并提取其第（3、4、5）行，第（2、3、4）列构成的新的矩阵 a=magic(5)a = 17 24 1 8 15 23 5 7 14 16 4 6 13 20 22 10 12 19 21 3 11 18 25 2 9 a(3:5,2:4)ans = 6 13 20 12 19 21 18 25 25）用命令行方式求解下式的值（提示使用syms x 定义一个符号，使用limit函数） syms x limit(cos(x)-exp(x2/2)/2)/4,x,0) ans = 1/86）MATLAB的绘图（1） 二维绘图命令plot：画出在上的图形 x=0:0.01:pi*2; y=
3、sin(x); plot(x,y)（2） 三维绘图命令plot3： 画出三维螺旋线，的图形mesh命令：绘制，在，区间的曲面 t=0:0.01:pi*4; x=sin(t);y=cos(t); z=t; plot3(x,y,z) x=-5:0.1:5; y=x;X,Y=meshgrid(x,y);Z=exp(-(X.2+Y.2)/22);mesh(X,Y,Z)7）编写M函数利用程序流程控制语句编写一个函数myfactorial（n），实现n！（阶乘）。要求使用help命令可以列出相关的帮助信息。functionx=myfactorial(n)%the function is to slove
4、n! x=1; for a=1:n; x=x*a; end help myfactorial the function is to slove n! myfactorial(8)ans = 40320三、实验结果及分析武汉工程大学实验报告专业班级：楼宇自动化01班姓名：刘斐学号：1104170109课程名称：控制系统数字仿真与CAD实验成绩： 指导老师：王利恒 实验名称：典型闭环系统的数字仿真及计算机解题一、实验目的：1）熟悉典型闭环的仿真过程2）掌握MATLAB编程实现典型闭环环节仿真3）利用典型闭环环节仿真程序解题。4）掌握MATLAB下对控制系统进行时域、频域和根轨迹的分析二、实验内容：
5、1）编写典型环节阶跃响应函数典型环节冲击响应函数function yout,t = my_step(num,den,v,t0,tf,h,R,n)输入参数：num：传递函数的分子系数向量den：传递函数的分母系数向量v：反馈比例系数t0：仿真起始时间tf：仿真终止时间h：仿真步长R:阶跃幅值n:系统阶次输出参数：yout：响应输出t：时间向量function yout,t = my_step(b,a,X0,V,T0,Tf,h,R,n)b=b/a(1);a=a/a(1);A=a(2:n+1);A=rot90(rot90(eye(n-1,n);-fliplr(A);B=zeros(1,n-1),1;
6、m1=length(b);C=fliplr(b),zeros(1,n-m1);Ab=A-B*C*V;X=X0;y=0;t=T0; N=round(Tf-T0)/h);for i=1:N; K1=Ab*X+B*R; K2=Ab*(X+h*K1/2)+B*R; K3=Ab*(X+h*K2/2)+B*R; K4=Ab*(X+h*K3)+B*R; X=X+h*(K1+2*K2+2*K3+K4)/6 y=y,C*X; t=t,t(i)+h;end;t,yplot(t,y)2）用上述函数分析以下系统，同时用simulink分析该系统，并比较其结果。程序为：function yout,t = my_step
7、(b,a,X0,V,T0,Tf,h,R,n)b=b/a(1);a=a/a(1);A=a(2:n+1);A=rot90(rot90(eye(n-1,n);-fliplr(A);B=zeros(1,n-1),1;m1=length(b);C=fliplr(b),zeros(1,n-m1);Ab=A-B*C*V;X=X0;y=0;t=T0; N=round(Tf-T0)/h);for i=1:N; K1=Ab*X+B*R; K2=Ab*(X+h*K1/2)+B*R; K3=Ab*(X+h*K2/2)+B*R; K4=Ab*(X+h*K3)+B*R; X=X+h*(K1+2*K2+2*K3+K4)/6
8、 y=y,C*X; t=t,t(i)+h;end;t,yplot(t,y)my_step(1 2 3,3 4 5 1,0 0 0,2,0,10,0.01,1,3)grid仿真图形为：Simulink对比可得，用m函数编写得到的典型环节冲击响应函数和用simulink绘制得到的函数结果一样。3）被控对象的传递函数为，用simulink建模并分析其单位阶跃响应。用MATLAB命令绘出其伯德图和根轨迹图。程序为： num=400;den=1 30 200 0;sys=tf(num,den); margin(sys) rlocus(sys)仿真结果为：三、实验结果及分析武汉工程大学实验报告专业班级：姓
9、名：学号：课程名称：控制系统数字仿真与CAD实验成绩： 指导老师： 实验名称：双闭环控制直流电动机调速系统建模与仿真分析一、实验目的：1. 熟悉Matlab/Simulink仿真环境；2. 掌握Simulink图形化建模方法；3. 验证 “直流电动机转速/电流双闭环PID控制方案”的有效性。二、实验内容：（详见后面的实验指导）1. “双闭环直流电动机调速系统”的建模Simulink画图：仿真图形为：2. 电流环调节器设计3. 电流环动态跟随性能仿真实验伯德图num,den=linmod(current_loop); sys=tf(num,den); margin(sys)奈魁斯特图阶跃曲线Tr
10、ansferfunction:-1.164e-010s3+2.038e007s2+4.977e009s+2.136e011-s5+1043s4+3.253e005s3+4.531e007s2+3.22e009s+8.544e010gm=4.3078pm=48.4499wcg=345.6682wcp=163.79234. 转速环调节器设计5. 转速环动态抗扰性能仿真实验6. 系统动态性能分析（给出仿真实验结果与理论分析结果的对比分析结论）三、实验结果及分析附件：“双闭环控制直流电动机调速系统”数字仿真实验实验指导一、 实验目的4. 熟悉Matlab/Simulink仿真环境；5. 掌握Simul
11、ink图形化建模方法；6. 验证 “直流电动机转速/电流双闭环PID控制方案”的有效性。二、 实验内容7. “双闭环直流电动机调速系统”的建模8. 电流环调节器设计9. 电流环动态跟随性能仿真实验10. 转速环调节器设计11. 转速环动态抗扰性能仿真实验12. 系统动态性能分析（给出仿真实验结果与理论分析结果的对比分析结论）三、 实验步骤1、系统建模A控制对象的建模建立线性系统动态数学模型的基本步骤如下：（1）根据系统中各环节的物理定律，列写描述据该环节动态过程的微分方程；（2）求出各环节的传递函数；（3）组成系统的动态结构图并求出系统的传递函数。下面分别建立双闭环调速系统各环节的微分方程和传
12、递函数。B额定励磁下的直流电动机的动态数学模型图1给出了额定励磁下他励直流电机的等效电路，其中电枢回路电阻R和电感L包含整流装置内阻和平波电抗器电阻与电感在内，规定的正方向如图所示。图1 直流电动机等效电路由图1可列出微分方程如下： （主电路，假定电流连续） （额定励磁下的感应电动势） （牛顿动力学定律，忽略粘性摩擦） （额定励磁下的电磁转矩）定义下列时间常数：电枢回路电磁时间常数，单位为s；电力拖动系统机电时间常数，单位为s；代入微分方程，并整理后得： 式中，负载电流。在零初始条件下，取等式两侧得拉氏变换，得电压与电流间的传递函数 （1）电流与电动势间的传递函数为 （2） a) b)c)图2 额定励磁下直流电动机的动态结构图a) 式（1）的结构图 b)式（2）的结构图c)整个直流电动机的动态结构图C晶闸管触发和整流装置的动态数学模型在分析系统时我们往往把它们当作一个环节来看待。这一环节的输入量是触发电路的控制电压Uct，输出量是理想空载整流电压Ud0。把它们之间的放大系数Ks看成常数，晶闸管触发与整流装置可以看成是一个具有纯滞后的放大环节，其滞后作用是由晶闸管装置的失控时间引起的。下面列出不同整流电路的平均失控时间：表1 各种整流电路的平均失控时间（f=50Hz）整流电路形式平均失控
