“控制系统数字仿真与CAD”实验指导书

合集下载

控制系统数字仿真与CAD课程设计

控制系统数字仿真与CAD课程设计

控制系统数字仿真与CAD课程设计课程背景控制系统数字仿真与CAD课程是一门涵盖了控制系统的基础理论以及控制系统CAD仿真实践的课程。

本课程主要教授学生掌握数字仿真与CAD技术在控制系统领域中的应用和实践。

控制系统数字仿真与CAD课程意在为学生开拓思路提供技术支持,同时也为学生将来走向控制系统领域提供必要的基础技能与实践经验。

课程目标•培养学生使用数字仿真软件进行控制系统仿真的能力。

•培养学生使用CAD软件进行控制系统图形设计与绘制的能力。

•培养学生将仿真与CAD技术应用于控制系统设计、分析与解决问题的能力。

•培养学生掌握控制系统相关的专业术语、技能与知识。

课程内容第一章基础知识本章主要是介绍控制系统的基础知识,包括控制系统的定义、分类、特点以及控制系统分析与设计的基础知识。

此外,还会介绍数字仿真与CAD技术在控制系统领域中的应用、主要功能与特点。

第二章仿真技术本章主要介绍数字仿真技术在控制系统中的应用,包括仿真的概念、方法、分类、技术流程以及仿真软件的选择、应用与实践。

此外,还会介绍仿真软件的主要功能与应用场景。

第三章 CAD技术本章主要介绍CAD技术在控制系统中的应用,包括CAD的概念、原理、主要技术与CAD软件的应用。

此外,还会介绍CAD技术在控制系统中的实践应用以及CAD软件的主要功能与应用场景。

第四章仿真与CAD技术在控制系统中的应用本章主要介绍仿真与CAD技术在控制系统中的应用,包括如何将仿真与CAD技术应用于控制系统设计、分析与解决问题;例如,如何绘制控制系统的传动布置图、控制系统的电气接线图等等。

此外,还会介绍仿真与CAD技术在控制系统维护、故障排查和升级改造中的应用。

第五章课程设计本章主要是针对控制系统数字仿真与CAD技术进行综合性的课程设计。

在课程设计中,学生需要综合应用数字仿真与CAD技术进行控制系统的设计与模拟仿真,从而提高问题解决能力、创新实践能力和综合应用能力。

课程评分1.课堂表现:20%2.仿真实验报告:30%3.CAD图形设计:20%4.课程设计报告:30%总结控制系统数字仿真与CAD课程设计是一门涵盖控制系统的基础理论,同时也重点介绍了CAD与数字仿真技术在控制系统领域中的应用和实践。

控制系统CAD与仿真实验指导书

控制系统CAD与仿真实验指导书

实验一MATLAB的实验环境及基本命令一实验目的:1.学习了解MA TLAB的实验环境2.在MA TLAB系统命令窗口练习有关MA TLAB命令的使用。

二实验步骤1.学习了解MA TLAB的实验环境:在Windows桌面上,用mouse双击MA TLAB图标,即可进入MA TLAB系统命令窗口:图1-1 MA TLAB系统命令窗口①在命令提示符”>>”位置键入命令:help此时显示MA T ALAB 的功能目录, 其中有“Matlab\general ”,“toolbox\control ”等;阅读目录的内容;② 键入命令:intro此时显示MA TLAB 语言的基本介绍,如矩阵输入、数值计算、曲线绘图等。

要求阅读命令平台上的注释内容,以尽快了解MA TLAB 语言的应用。

③ 键入命令:help help显示联机帮助查阅的功能,要求仔细阅读。

④ 键入命令:into显示工具箱中各种工具箱组件和开发商的联络信息。

⑤ 键入命令:demo显示MA TLAB 的各种功能演示。

2. 练习MA TLAB 系统命令的使用。

① 表达式MA TLAB 的表达式由变量、数值、函数及操作符构成。

实验前应掌握有关变量、数值、函数及操作符的有关内容及使用方法。

练习1-1: 计算下列表达式:要求计算完毕后,键入相应的变量名,查看并记录变量的值。

②.向量运算: )6sin(/250π=d 2/)101(+=a )sin(3.2-=e c i b 53+=n 维向量是由n 个成员组成的行或列数组。

在MA TLAB 中,由分号分隔的方括号中的元素产生一个列向量;由逗号或空号分隔的方括号中的元素产生一个列向量;同维的向量可进行加减运算,乘法须遵守特殊的原则。

练习1-2已知:X=[2 ;-4;8]求 :Y=R ';P=5*R ;E=X .*Y ;S=X '* Y练习1-3⑴产生每个元素为1的4维的行向量;⑵产生每个元素为0的4维的列向量;⑶产生一个从1到8的整数行向量,默认步长为1;⑷产生一个从π到0,间隔为π/3的行向量;③矩阵基本运算操作。

控制系统CAD及数字仿真实验指导书

控制系统CAD及数字仿真实验指导书

控制系统CAD及数字仿真实验指导书自动化系目录第一章前言............... (2)第二章控制系统CAD及数字仿真实验1.控制系统稳定性的MATLAB辅助分析 (3)2.数据处理方法的MATLAB实现 (4)3.控制系统数字仿真的实现. (5)4.基于双闭环PID控制的一阶倒立摆控制系统仿真实验 (6)5.双闭环直流调速系统设计仿真实验............ .... .... .... .7 6.水箱液位控制仿真实验........................... .. .... . .8第一章前言本实验课程是与学科基础选修课程《控制系统CAD及数字仿真》相配合的实践课程。

作为联系自动控制理论、自动控制系统设计、课程设计、毕业设计等教学环节的仿真技术类课程,其不仅可以使学生加强课程的学习效果,而且还可为学生在毕业设计中提供一个强有力的工具,有效加强教学中的实践性教学环节,提高学生的独立工作能力和创造性思维能力。

开设本课程的目的,主要是培养学生运用MATLAB语言进行编程和仿真的能力,为今后从事科研工作和与专业有关的工程技术工作打好基础。

一、上机实验要求1、要求学生熟悉MATLAB中的控制系统工具箱与SIMULINK软件包。

2、能根据有关控制算法,编写有关的MATLAB程序。

3、能对实验结果进行分析和讨论,得到相关的实验结论。

二、上机实验的基本程序:1、明确实验任务。

2、提出实验方案。

3、编制有关的MATLAB程序或利用SIMULINK工具建立系统的仿真模型。

4、进行实验操作,作好观测和记录,保存有关的实验数据。

5、整理实验数据,得出结论,撰写实验报告。

在进行上机实验时,上述程序应让学生独立完成,教师给予必要的指导,以培养学生的动手能力。

要做好各个上机实验,就应做到:实验前做准备,实验中有条理,实验后勤分析。

实验一控制系统稳定性的MATLAB辅助分析一、实验目的学会在MATLAB环境下对线性控制系统进行稳定性分析。

“控制系统数字仿真与CAD”实验指导书.

“控制系统数字仿真与CAD”实验指导书.

写,均要求用标准A4纸进行撰写,单栏排版,单面打印,并左侧装订,以便于报告最终的批阅与存档,(对于存在“逻辑混乱” 、“文字不清” 、“作图潦草” 、“排版混乱”等问题的报告,将予以退回重新撰写)。

封页:(参考最后一页的“封页”格式)正文:(小四字体)仿真实验题目 1、 2、 3、 4、 5、五、思考题 1.在系统启动过程的第 2 阶段中,理想的电流特性为:实际值小于给定/设定值,试说明为何?引言原理/建模设计/分析/论述仿真实验/结果分析结论(思考题解答) 2.动态性能中,电流/转速特性的“超调量”与理论值是否有偏差?;如有偏差,试给出分析/解释。

3.在“双闭环直流电动机调速系统”中,电流调节器与速度调节器的输出都要设置“限幅” ,试说明:你是如何选取限幅值的? 4.假设系统中的励磁电压减小/增加,试说明:系统转速将可能怎样变化?参考文献: [1] 张晓华主编《控制系统数字仿真与 CAD》第 3 版机械工业出版社 2009 [2] 张晓华主编《系统建模与仿真》清华大学出版社 2006 [3] 陈伯时主编《电力拖动自动控制系统》第 3 版机械工业出版社 2008 25
“控制系统数字仿真与CAD” 仿真实验报告姓班学名:级:号:联系电话: Email: 提交日期: 26。

《控制系统仿真与CAD》实验指导书朱品伟

《控制系统仿真与CAD》实验指导书朱品伟

《控制系统仿真与CAD》实验指导书朱品伟《控制系统仿真与CAD》实验指导书江苏理工学院电气信息工程学院朱品伟 2021.12实验一 Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法;2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容 1、帮助命令使用help命令,查找sqrt(开方)函数的使用方法; 2、矩阵运算(1)矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B。

(2)矩阵除法已知 A=[1 1 2;4 9 6;5 7 9]; B=[1 2 1;3 2 4;5 8 3]; 求A\\B,A/B。

(3)矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i]; 求A.', A'。

(4)使用冒号选出指定元素已知: A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素;A中所有列第2,3行的元素;(5)方括号[]用magic函数生成一个4阶魔术矩阵,删除该矩阵的第四列。

3、多项式(1)求多项式的根(2)已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] 求矩阵A的特征多项式;求特征多项式中未知数为20时的值;把矩阵A作为未知数代入到多项式中; 4、基本绘图命令(1)绘制余弦曲线y=cos(t),t∈[0,2π](2)在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5),t∈[0,2π] 5、基本绘图控制绘制[0,4π]区间上的x1=10sint曲线,并要求:(1)线形为点划线、颜色为红色、数据点标记为加号;(2)坐标轴控制:显示范围、刻度线;(3)标注控制:坐标轴名称、标题、相应文本; 6、基本程序设计(1)编写命令文件:计算1+2+…+n<2000时的最大n值;(2)编写函数文件:分别用for和while循环结构编写程序,求2的0到n次幂的和。

控制系统仿真与CAD-实验报告

控制系统仿真与CAD-实验报告

《控制系统仿真与CAD》实验课程报告一、实验教学目标与基本要求上机实验是本课程重要的实践教学环节。

实验的目的不仅仅是验证理论知识,更重要的是通过上机加强学生的实验手段与实践技能,掌握应用MATLAB/Simulink 求解控制问题的方法,培养学生分析问题、解决问题、应用知识的能力和创新精神,全面提高学生的综合素质。

通过对MATLAB/Simulink进行求解,基本掌握常见控制问题的求解方法与命令调用,更深入地认识和了解MATLAB语言的强大的计算功能与其在控制领域的应用优势。

上机实验最终以书面报告的形式提交,作为期末成绩的考核内容。

二、题目及解答第一部分:MATLAB 必备基础知识、控制系统模型与转换、线性控制系统的计算机辅助分析1.>>f=inline('[-x(2)-x(3);x(1)+a*x(2);b+(x(1)-c)*x(3)]','t','x','flag','a','b','c');[t,x]=ode45( f,[0,100],[0;0;0],[],0.2,0.2,5.7);plot3(x(:,1),x(:,2),x(:,3)),grid,figure,plot(x(:,1),x(:,2)), grid2.>>y=@(x)x(1)^2-2*x(1)+x(2);ff=optimset;rgeScale='off';ff.TolFun=1e-30;ff.Tol X=1e-15;ff.TolCon=1e-20;x0=[1;1;1];xm=[0;0;0];xM=[];A=[];B=[];Aeq=[];Beq=[];[ x,f,c,d]=fmincon(y,x0,A,B,Aeq,Beq,xm,xM,@wzhfc1,ff)Warning: Options LargeScale = 'off' and Algorithm ='trust-region-reflective' conflict.Ignoring Algorithm and running active-set algorithm. To runtrust-region-reflective, setLargeScale = 'on'. To run active-set without this warning, useAlgorithm = 'active-set'.> In fmincon at 456Local minimum possible. Constraints satisfied.fmincon stopped because the size of the current search direction is less thantwice the selected value of the step size tolerance and constraints aresatisfied to within the selected value of the constraint tolerance.<stopping criteria details>Active inequalities (to within options.TolCon = 1e-20):lower upper ineqlin ineqnonlin2x =1.00001.0000f =-1.0000c =4d =iterations: 5funcCount: 20lssteplength: 1stepsize: 3.9638e-26algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'firstorderopt: 7.4506e-09constrviolation: 0message: [1x766 char]3.(a) >> s=tf('s');G=(s^3+4*s+2)/(s^3*(s^2+2)*((s^2+1)^3+2*s+5))G =s^3 + 4 s + 2------------------------------------------------------s^11 + 5 s^9 + 9 s^7 + 2 s^6 + 12 s^5 + 4 s^4 + 12 s^3Continuous-time transfer function.(b)>> z=tf('z',0.1);H=(z^2+0.568)/((z-1)*(z^2-0.2*z+0.99))H =z^2 + 0.568-----------------------------z^3 - 1.2 z^2 + 1.19 z - 0.99Sample time: 0.1 secondsDiscrete-time transfer function.4.>> A=[0 1 0;0 0 1;-15 -4 -13];B=[0 0 2]';C=[1 00];D=0;G=ss(A,B,C,D),Gs=tf(G),Gz=zpk(G)G =a =x1 x2 x3x1 0 1 0x2 0 0 1x3 -15 -4 -13b =u1x1 0x2 0x3 2c =x1 x2 x3y1 1 0 0d =u1y1 0Continuous-time state-space model.Gs =2-----------------------s^3 + 13 s^2 + 4 s + 15 Continuous-time transfer function.Gz =2---------------------------------(s+12.78) (s^2 + 0.2212s + 1.174) Continuous-time zero/pole/gain model.5.设采样周期为0.01s>> z=tf('z',0.01);H=(z+2)/(z^2+z+0.16) H =z + 2--------------z^2 + z + 0.16Sample time: 0.01 secondsDiscrete-time transfer function.6.>> syms J Kp Ki s;G=(s+1)/(J*s^2+2*s+5);Gc=(Kp*s+Ki)/s;GG=feedback(G*Gc,1) GG =((Ki + Kp*s)*(s + 1))/(J*s^3 + (Kp + 2)*s^2 + (Ki + Kp + 5)*s + Ki)7.(a)>>s=tf('s');G=(211.87*s+317.64)/((s+20)*(s+94.34)*(s+0.1684));Gc=(169.6*s+400)/(s*(s+4));H=1/(0.01*s+1);GG=feedback(G*Gc,H),Gd=ss(GG),Gz=zpk(GG)GG =359.3 s^3 + 3.732e04 s^2 + 1.399e05 s + 127056----------------------------------------------------------------0.01 s^6 + 2.185 s^5 + 142.1 s^4 + 2444 s^3 + 4.389e04 s^2 + 1.399e05 s + 127056Continuous-time transfer function.Gd =a =x1 x2 x3 x4 x5 x6x1 -218.5 -111.1 -29.83 -16.74 -6.671 -3.029x2 128 0 0 0 0 0x3 0 64 0 0 0 0x4 0 0 32 0 0 0x5 0 0 0 8 0 0x6 0 0 0 0 2 0b =u1x1 4x2 0x3 0x4 0x5 0x6 0c =x1 x2 x3 x4 x5 x6y1 0 0 1.097 3.559 1.668 0.7573d =u1y1 0Continuous-time state-space model.Gz =35933.152 (s+100) (s+2.358) (s+1.499)----------------------------------------------------------------------(s^2 + 3.667s + 3.501) (s^2 + 11.73s + 339.1) (s^2 + 203.1s + 1.07e04) Continuous-time zero/pole/gain model.(b)设采样周期为0.1s>>z=tf('z',0.1);G=(35786.7*z^2+108444*z^3)/((1+4*z)*(1+20*z)*(1+74.04*z));Gc= z/(1-z);H=z/(0.5-z);GG=feedback(G*Gc,H),Gd=ss(GG),Gz=zpk(GG)GG =-108444 z^5 + 1.844e04 z^4 + 1.789e04 z^3----------------------------------------------------------------1.144e05 z^5 +2.876e04 z^4 + 274.2 z^3 + 782.4 z^2 + 47.52 z + 0.5Sample time: 0.1 secondsDiscrete-time transfer function.Gd =a =x1 x2 x3 x4 x5 x1 -0.2515 -0.00959 -0.1095 -0.05318 -0.01791x2 0.25 0 0 0 0x3 0 0.25 0 0 0x4 0 0 0.125 0 0x5 0 0 0 0.03125 0b =u1x1 1x2 0x3 0x4 0x5 0c =x1 x2 x3 x4 x5y1 0.3996 0.6349 0.1038 0.05043 0.01698d =u1y1 -0.9482Sample time: 0.1 secondsDiscrete-time state-space model.Gz =-0.94821 z^3 (z-0.5) (z+0.33)----------------------------------------------------------(z+0.3035) (z+0.04438) (z+0.01355) (z^2 - 0.11z + 0.02396)Sample time: 0.1 secondsDiscrete-time zero/pole/gain model.8.>>s=tf('s');g1=1/(s+1);g2=s/(s^2+2);g3=1/s^2;g4=(4*s+2)/(s+1)^2;g5=50;g6=(s^2+2) /(s^3+14);G1=feedback(g1*g2,g4);G2=feedback(g3,g5);GG=3*feedback(G1*G2,g6) GG =3 s^6 + 6 s^5 + 3 s^4 + 42 s^3 + 84 s^2 + 42 s---------------------------------------------------------------------------s^10 + 3 s^9 + 55 s^8 + 175 s^7 + 300 s^6 + 1323 s^5 + 2656 s^4 + 3715 s^3 + 7732 s^2 + 5602 s + 1400Continuous-time transfer function.9.>>s=tf('s');T0=0.01;T1=0.1;T2=1;G=(s+1)^2*(s^2+2*s+400)/((s+5)^2*(s^2+3*s+100 )*(s^2+3*s+2500));Gd1=c2d(G,T0),Gd2=c2d(G,T1),Gd3=c2d(G,T2),step(G),figure,st ep(Gd1),figure,step(Gd2),figure,step(Gd3)Gd1 =4.716e-05 z^5 - 0.0001396 z^4 + 9.596e-05 z^3 + 8.18e-05 z^2 - 0.0001289 z + 4.355e-05----------------------------------------------------------------z^6 - 5.592 z^5 + 13.26 z^4 - 17.06 z^3 + 12.58 z^2 - 5.032 z + 0.8521Sample time: 0.01 secondsDiscrete-time transfer function.Gd2 =0.0003982 z^5 - 0.0003919 z^4 - 0.000336 z^3 + 0.0007842 z^2 - 0.000766 z + 0.0003214----------------------------------------------------------------z^6 - 2.644 z^5 + 4.044 z^4 - 3.94 z^3 + 2.549 z^2 - 1.056 z + 0.2019Sample time: 0.1 secondsDiscrete-time transfer function.Gd3 =8.625e-05 z^5 - 4.48e-05 z^4 + 6.545e-06 z^3 + 1.211e -05 z^2 - 3.299e-06 z + 1.011e-07---------------------------------------------------------------z^6 - 0.0419 z^5 - 0.07092 z^4 - 0.0004549 z^3 + 0.002495 z^2 - 3.347e-05 z + 1.125e-07Sample time: 1 secondsDiscrete-time transfer function.10.(a)>> G=tf(1,[1 2 1 2]);eig(G),pzmap(G) ans =-2.0000-0.0000 + 1.0000i-0.0000 - 1.0000i系统为临界稳定。

《控制系统仿真与CAD》课程设计指导书

《控制系统仿真与CAD》课程设计指导书

《控制系统仿真与CAD》课程设计指导书一、目的和任务配合《控制系统仿真与CAD》课程的理论教学,通过课程设计教学环节,使学生掌握当前流行的演算式MATLAB语言的基本知识,学会运用MATLAB语言进行控制系统仿真和辅助设计的基本技能,有效地提高学生实验动手能力。

基本要求:1、利用MATLAB提供的基本工具,灵活地编制和开发程序,开创新的应用。

2、熟练地掌握各种模型之间的转换,系统的时域、频域分析及根轨迹绘制。

3、熟练运用SIMULINK对系统进行仿真。

4、掌握PID控制器参数的设计。

二、设计要求1、编制相应的程序,并绘制相应的曲线。

2、对设计结果进行分析。

3、撰写和打印设计报告(包括程序、结果分析、仿真结构框图、结果曲线)。

三、设计内容1、本次设计有八个可以选择的题目,至少选择两个题目进行设计。

2、“设计报告”要按规定的格式撰写(对于存在“逻辑混乱”、“文字不清”、“作图潦草”等问题的报告,将予以退回重新撰写)。

3、无论计算机录入/打印还是手工书写,均要求用标准A4纸进行撰写,以便于报告最终的批阅与存档。

四、时间安排1、课程设计时间为一周。

2、第1天布置设计题目,讲授设计的要求。

3、第2~4天学生进行设计。

4、第5天教师验收,然后学生撰写和打印设计报告。

选题一:二阶弹簧—阻尼系统的PID 控制器设计及参数整定考虑弹簧-阻尼系统如图1所示,其被控对象为二阶环节,传递函数()G s 如下,参数为M=1kg ,b=2N.s/m ,k=25N/m ,()1F s =。

设计要求:(1)控制器为P 控制器时,改变比例系数大小,分析其对系统性能的影响并绘制相应曲线。

(2)控制器为PI 控制器时,改变积分系数大小,分析其对系统性能的影响并绘制相应曲线。

(例如当Kp=50时,改变积分系数大小)(3)设计PID 控制器,选定合适的控制器参数,使闭环系统阶跃响应曲线的超调量σ%<20%,过渡过程时间Ts<2s, 并绘制相应曲线。

控制系统CAD实验指导书2011

控制系统CAD实验指导书2011

《控制系统CAD 》上机实验指导书一、课程实验目的通过运用MATLAB 软件进行控制系统的分析和设计,对控制类课程进行理论验证实验和设计实验,使学生能够进一步熟悉控制系统的原理和设计方法,更直观的了解自动控制系统的运动过程,从而提高学生对控制系统的分析和设计能力。

二、上机实验要求(1) 在MATLAB 环境下独立完成实验内容。

(2) 根据实验结果写出分析报告,报告中附上实验中所编写的程序,同时要提交该程序的电子版。

另外,实验报告中,对每一个实验除了说明实验过程及实验结果以外,要对实验中出现的问题进行分析,并进行小结。

封面要注明姓名、学号、专业等信息。

提交程序电子版时,请以学号建立一个文件夹,将实验中的程序放入该文件夹里,并建立一个readme.txt 文件,用以说明程序和实验的对应关系,以及姓名等个人信息。

三、上机实验内容(8学时)实验一(1学时):(1)熟悉MATLAB 环境;(2)控制系统模型及其转换。

实验二(3学时):二阶控制系统在时域和频域内的分析。

实验三(3学时):编制程序,实现串联校正(超前、滞后、超前--滞后校正),要得到校正装置的元件数值的大小,并标称化。

① 设一单位负反馈控制系统,如果控制对象的传递函数为:)80)(4()(++=s s s K s G p 试设计一个串联超前校正装置。

要求:相角裕度≥45。

;当系统的输入信号是单位斜坡信号时,稳态误差e ss ≤0.04;取C=1μF 时,确定该串联超前校正装置的元件数据;绘制出校正后系统和未校正系统的Bode 图及其闭环系统的单位阶跃响应曲线,并进行对比。

(提示:稳态误差e ss ≤0.04 —> 取kv=1/e ss =25,k0=8000 )② 设一单位负反馈控制系统,其控制对象的传递函数为:)80)(4()(++=s s s K s G p 试设计一个串联滞后校正装置。

要求:相角裕度≥45。

;当系统的输入信号是单位斜坡信号时,稳态误差e ss ≤0.04;取C=100μF 时,确定该串联滞后校正装置的元件数据;绘制出校正后系统和未校正系统的Bode 图及其闭环系统的单位阶跃响应曲线,并进行对比。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

(17)
一般以选择 h=5 为好所以:
n 5 T n , K N
6 50 T n 2
(18)
经过如上设计,得到的双闭环控制系统从理论上讲有如下动态性能: 电动机起动 过程中电流的超调量为 4.3%,转速的超调量为 8.3%。
6
(3)ACR 和 ASR 的理论设计及结果 ①电流环的设计 电流环的设计具体设计步骤如下:
U d 0 ( s)
E ( s)
1/ R Tl s 1
I d (s)
I d (s)
R Tm s
E ( s)
a)
式(1)的结构图
1/ R I d ( s ) Tl s 1
b) 式(2)的结构图
I dL ( s )
R E (s) 1 Tm s Ce
U d 0 (s)
n( s )
c) 整个直流电动机的动态结构图 图 2 额定励磁下直流电动机的动态结构图
U d 0 (s) K Ks K s e Ts s Tsss 1 1 U ct ( s ) e 1 Ts s Ts 2 s 2 Ts 3 s 3 2! 3! 考虑到 Ts 很小,忽略其高次项,则晶闸管触发和整流装置的传递函数可近似成
一阶惯性环节
U d 0 (s) Ks U ct ( s ) Ts s 1
(20)
(d) 校验近似条件
由电流环截止频率,晶闸管装置传递函数近似条件,忽略反电势对电流环影响 的条件,小时间常数近似处理条件等考虑得 电流调节器传递函数为
WACR ( s ) 0.292 0.018s 1 0.018s 1 0.018s 0.062 s
(21)
②转速环的设计 具体设计步骤如下:
对于电力拖动控制系统,转速环通常希望具有良好的抗扰性能,因此我们要把 转速环校正成典型Ⅱ型系统。 要把转速环校正成典型Ⅱ型系统,ASR 也应该采用 PI 调节器,其传递函数为
WASR ( s ) K n
式中 Kn—电流调节器的比例系数;
ns 1 ns
(13)
n —电流调节器的超前时间常数。
其结构图如图 3 所示。
(4)
U ct ( s )
K s eTs s
U d 0 ( s)
U ct ( s)
Ks Ts s 1
U d 0 (s)
a) 准确的结构图 b)近似的结构图 图 3 晶闸管触发和整流装置的动态结构图
3
D.比例放大器、测速发电机和电流互感器的动态数学模型 比例放大器、测速发电机和电流互感器的响应都可以认为是瞬时的,因此它们 的放大系数也就是它们的传递函数,即
WASR ( s ) 19.33
③ASR 输出限幅值的确定
0.0917 s 1 0.0917 s 1 0.0917 s 0.005s
(23)
转速外环呈开环状态, 转速的变化对系统不再产 当 ASR 输出达到限幅值 U*im, 生影响。双闭环系统变成一个电流无静差的单闭环系统。稳态时
Id U im*

I dm
(24)
式中,最大电流 Idm 是由设计者选定的,取决于电机的过载能力和拖动系统允 许的最大加速度。在这里,我们选取 Idm=20A,那么 ASR 输出限幅值为
U im* I dm 0.4 20 8V
(25)
4、SIMULINK 建模
借助 SIMULINK,根据上节理论计算得到的参数,可得双闭环调速系统的动态 结构图如下所示:
U ct ( s ) Kp U n ( s )
U n (s) n( s ) U i (s) I d ( s)
(5) (6) (7)
E.双闭环控制直流电动机调速系统的动态数学模型 根据以上分析,可得双闭环控制系统的动态结构图如下
I dL
Un

WASR ( s)
Un
Ui

WACR ( s )
控制系统数字仿真与 CAD
实验指导书
王高林 编
哈尔滨工业大学电气工程系 2010 年 11 月
“双闭环控制直流电动机调速系统”数字仿真实验
一、 实验目的
1. 熟悉 Matlab/Simulink 仿真环境; 2. 掌握 Simulink 图形化建模方法; 3. 验证 “直流电动机转速/电流双闭环 PID 控制方案”的有效性。
WACR ( s ) K i
is 1 is
(8)
式中 Ki—电流调节器的比例系数;
i —电流调节器的超前时间常数。
为了让调节器零点对消掉控制对象的大时间常数(极点) ,选择
i Tl
(9)
一般情况下,希望超调量 σ%≤5%时,取阻尼比 ξ=0.707, K I T i 0.5 ,得:
5
KI
1 , ( T i Ts Toi ) 2T i
KI Ki K s iR
(10) (11) (12)
又因为 得到 (2)转速调节器的设计ห้องสมุดไป่ตู้
Ki K I
iR Tl R R Tl 0.5 K s 2 K s Ti K s Ti
二、 实验内容
1. “双闭环直流电动机调速系统”的建模 2. 电流环/调节器设计 3. 电流环动态跟随性能仿真实验 4. 转速环/调节器设计 5. 转速环动态抗扰性能仿真实验 6. 系统动态性能分析 (给出仿真实验结果与理论分析结果的对比/分析/结论)
三、 实验步骤 1、系统建模
A.控制对象的建模 建立线性系统动态数学模型的基本步骤如下: (1)根据系统中各环节的物理定律,列写描述据该环节动态过程的微分方程; (2)求出各环节的传递函数; (3)组成系统的动态结构图并求出系统的传递函数。 下面分别建立双闭环调速系统各环节的微分方程和传递函数。 B.额定励磁下的直流电动机的动态数学模型 图 1 给出了额定励磁下他励直流电机的等效电路,其中电枢回路电阻 R 和电感 L 包含整流装置内阻和平波电抗器电阻与电感在内,规定的正方向如图所示。
Ui
U ct
Ks U d 0 Ts s 1
1/ R I d Tl s 1
R Tm s
1 Ce
n


图 4 双闭环控制系统的动态结构图
2、实验系统参数
系统中采用三相桥式晶闸管整流装置,基本参数如下: 直流电动机:220V,13.6A,1480r/min,Ce =0.131V/(r/min) , 允许过载倍数 λ=1.5。 晶闸管装置: K s 76 。 电枢回路总电阻:R=6.58Ω。 时间常数: Tl =0.018s, Tm =0.25s。 反馈系数:α=0.00337V/(r/min) ,β=0.4V/A。 反馈滤波时间常数: Toi =0.005s, Ton =0.005s。
(a) 确定时间常数
按小时间常数近似处理,取 Tn 2Ti Ton 0.01334 0.005 0.01834s 。
(b) 选择转速调节器结构
由于设计要求无静差,转速调节器必须含有积分环节;又根据动态要求,应按
7
典型Ⅱ型系统设计转速环。故 ASR 选用 PI 调节器,其传递函数为
4
3. PID 调节器参数设计
设计多闭环控制系统的一般原则是:从内环开始,一环一环地逐步向外扩展。 在这里是:先从电流环入手,首先设计好电流调节器,然后把整个电流环看作是转 速调节系统中的一个环节,再设计转速调节器。 双闭环控制系统的动态结构图绘于图 5,它增加了滤波环节,包括电流滤波、 转速滤波和两个给定滤波环节。 其中 Toi 为电流反馈滤波时间常数,Ton 为转速反馈滤波时间常数
WASR ( s ) K n
ns 1 ns
(22)
(c) 选择转速调节器参数
按典型Ⅱ型系统最佳参数的原则,取 h=5,则 ASR 的超前时间常数为
n hT n 5 0.01834s 0.0917s
转速开环增益
KN h 1 6 1/s 2 356.77s -2 2 2 2 2h Tn 2 25 0.01834
用单位阶跃函数来表示滞后,则晶闸管触发和整流装置的输入输出关系为
U d 0 K sU ct 1(t Ts )
按拉氏变换的位移定理,则传递函数为
U d 0 (s) K s eTs s U ct ( s )
(3)
由于式(3)中含有指数函数 e Ts s ,它使系统成为非最小相位系统,分析和设计 都比较麻烦。为了简化,先将 e Ts s 按台劳级数展开,则式(3)变成
(a) 确定时间常数
整流装置滞后时间常数 Ts 按表 1,三相桥式电路的平均失控时间 Ts=0.00167s。 电流滤波时间常数 Toi=0.005s。 电流环小时间常数 Ti 取 Ti Ts Toi 0.00167 0.005 0.00667s 。
(b) 选择电流调节器结构
于是,ASR 的比例系数为
Kn (h 1) CeTm 6 0.4 0.131 0.25 19.33 2h RTn 2 5 0.00337 6.58 0.01834
(d) 校验近似条件
从转速环截止频率,电流环传递函数简化条件,小时间常数近似处理条件等考 虑得: 转速调节器传递函数为
转速开环增益 按照典型Ⅱ型系统的参数选择方法, ( T n 2T i Ton ) n hT n ,
KN h 1 2h 2Tn 2
KN
K n R n CeTm
(14)
(15) (16)
考虑到式(14)和(15) ,得到 ASR 的比例系数
Kn ( h 1) CeTm 2h RTn
电流调节器选择 PI 型,其传递函数为
WACR ( s ) K i
相关文档
最新文档