




已阅读5页,还剩12页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
. 目录一. 摘要.1二. 关键词.1三. Abstract.1四. Keywords.1五. 课程设计的目的.4六. 课程设计用的仪器与器件.4七. 课程设计的具体内容.4八. 心得体会.16九. 参考文献.17十. 附录(程序).17 控制系统建模、分析、设计和仿真 自动化学生 安国杰 指导老师 侯建华摘要:本次课程实践为设计两个控制器分别为最小拍无波纹和最小拍有波纹控制器。通过这次实践可以进一步对所学的计算机控制技术有进一步的了解,并对Matlab软件的操作有一定程度的熟悉,为以后的工作或研究作基础。 Matlab最初主要用于矩阵数值的计算,随着它的版本功能来越来强大,应用范围内也越来越广阔。如今,Matlab已经发展成为国际上非常好的生产潜力这与工程计算机语言之一,它使用方便、输入快捷、运算高效、内容丰富,是高等院校理工科学和科研中必不可少的工具之一,掌握Matlab已经成为相关专业课大学生,研究生和教师的必备技能。 Matlab在我国的应用已有十多年的历史上,而自动控制器则是其最重要的应用领域之一,自动控制系统的建模、分析、设计及应用等都离不开Matlab 4.0开始,就一直将应用教学和科研工作中。对于自动化专业的本科生来说,Matlab是一种必须掌握的现代计算工工具。关键词: MALAB ; 仿真 。 Computer Simulation Course Design Student majoring in automation anguojie Tutor hooujianhuaSummary:The course of practice for the design of two controllers are the minimum ripple and minimum film shot there is no ripple controller. Through this practice can be further learned the computer-controlled technology have a better understanding, and the operation of the software matlab have a certain degree of familiarity, for future work or study as a basis. Matlab was originally mainly used for numerical calculation of the matrix, with its version of the feature to increasingly powerful, increasingly broad range of applications. Today, Matlab has become a very good production potential in the international engineering computer language which one it easy to use, enter the fast, efficient computing, content-rich, scientific and research institutions of higher learning Polytechnic indispensable tool , Matlab has become a master professional courses related to college students, graduate students and teachers the necessary skills. Matlab application in China for more than ten years of history, while the automatic controller is the most important areas of application, and automatic control systems modeling, analysis, design and application are all inseparable from Matlab 4.0 starting has been the application of teaching and research work. For the automation of the undergraduates, Matlab is a need to master work of modern computing tools.Keywords: MATLAB ; simulation。1课程设计的目的(1)学习并掌握MATLAB软件的程序编写与调试,以及计算机仿真技术的应用。(2)运用MATLAB 软件处理和仿真,分析所建立的控制系统模型的可行性。(3)学会运用Simulink对动态系统进行建模、仿真和分析。2课程设计用的仪器与器件硬件: PC计算机 软件: MATLAB7.03课程设计的具体内容 2号题 控制系统建模、分析、设计和仿真设连续被控对象的实测传递函数为:用零阶保持器离散化,采样周期取0.2秒,分别设计一单位加速度信号输入时的最少拍有波纹控制器Dy(z)和一单位速度信号输入时的最少拍无波纹控制器Dw(z)。具体要求见(二)。1、求被控对象传递函数G(s)的MATLAB描述。 z=-2 -6;p=0 0 -1 -5 -8;k=668; Gs=zpk(z,p,k)Zero/pole/gain: 668 (s+2) (s+6)-s2 (s+1) (s+5) (s+8)2、求被控对象脉冲传递函数G(z)。 Gz=c2d(Gs,0.2,zoh) Zero/pole/gain:0.68436 (z+2.904) (z-0.6703) (z-0.3011) (z+0.1908)- (z-1)2 (z-0.8187) (z-0.3679) (z-0.2019) Sampling time: 0.23、转换G(z)为零极点增益模型并按z-1形式排列。 a,b,c=zpkdata(Gz)a = 4x1 double b = 5x1 double c = 0.6844 Gz=zpk(a,b,c,0.2,variable,z-1) Zero/pole/gain:0.68436 z-1 (1+2.904z-1) (1-0.6703z-1) (1-0.3011z-1) (1+0.1908z-1)- (1-z-1)2 (1-0.8187z-1) (1-0.3679z-1) (1-0.2019z-1) Sampling time: 0.24、确定误差脉冲传递函数Ge(z)形式,满足单位加速度信号输入时闭环稳态误差为零和实际闭环系统稳定的要求。 syms z b0 b1 a0 a1 a2 a3 Gez=(1-z-1)3*(1+b0*z-1) Gez = (1-1/z)3*(1+b0/z)5、确定闭环脉冲传递函数Gc(z)形式,满足控制器Dy(z)可实现、最少拍和实际闭环系统稳定的要求。 Gcz=z-1*(1+2.904*z-1)*(a0+a1*z-1+a2*z-2) Gcz = 1/z*(1+363/125/z)*(a0+a1/z+a2/z2)6、根据4、5、列写方程组,求解Gc(z)和Ge(z)中的待定系数并最终求解Gc(z)和Ge(z) 。 f1=subs(Gcz,z,1)-1 f1 = 488/125*a0+488/125*a1+488/125*a2-1 f2=subs(diff(Gcz,1),z,1) f2 = -851/125*a0-1339/125*a1-1827/125*a2 f3=subs(diff(Gcz,2),z,1) f3 = 2428/125*a0+5106/125*a1+1752/25*a2 a0j a1j a2j=solve(f1,f2,f3) a0j = 150061125/116214272 a1j = -94337625/58107136 a2j = 68382125/116214272 A=double(a0j a1j a2j)A =1.2912 -1.6235 0.5884 Gcz=subs(Gcz,a0 a1 a2,A) Gcz =1/z*(1+363/125/z)*(5815251560779019/4503599627370496-7311647450616351/4503599627370496/z+5299963719925943/9007199254740992/z2) Gez=(1-z-1)3*(1+b0*z-1) Gez = (1-1/z)3*(1+b0/z) f4=subs(Gez,z,-2.904)-1 f4 = 68382125/4783214717363069361*b0 boj=solve(f4) boj = 198581691/116214272 B=double(boj)B =1.7088 Gez=subs(Gez,1 b0,1 boj) Gez = (1-1/z)3*(1+198581691/116214272/z)7、求针对单位加速度信号输入的最少拍有波纹控制器Dy(z)并说明Dy(z)的可实现性。 Gz=0.68436*z-1*(1+2.904*z-1)*(1-0.6703*z-1)*(1-0.3011*z-1)*(1+0.1908*z-1)/(1-z-1)2/(1-0.8187*z-1)/(1-0.3679*z-1)/(1-0.2019*z-1) Gz =17109/25000/z*(1+363/125/z)*(1-6703/10000/z)*(1-3011/10000/z)*(1+477/2500/z)/(1-1/z)2/(1-8187/10000/z)/(1-3679/10000/z)/(1-2019/10000/z) Dyz=Gcz/Gez/Gz Dyz =25000/17109*(5815251560779019/4503599627370496-7311647450616351/4503599627370496/z+5299963719925943/9007199254740992/z2)/(1-1/z)/(1+198581691/116214272/z)/(1-6703/10000/z)/(1-3011/10000/z)/(1+477/2500/z)*(1-8187/10000/z)*(1-3679/10000/z)*(1-2019/10000/z) 8、用程序仿真方法分析加速度信号输入时闭环系统动态性能和稳态性能。 N,D=numden(simplify(Gcz); numc=sym2poly(N)numc =1.0e+018 * 1.4538 2.3940 -4.6458 1.9239 denc=sym2poly(D)denc = 1.0e+018 * 1.1259 0 0 0 0 t=0:0.2:10t = Columns 1 through 14 0 0.2000 0.4000 0.6000 0.8000 1.0000 1.2000 1.4000 1.6000 1.8000 2.0000 2.2000 2.4000 2.6000 Columns 15 through 28 2.8000 3.0000 3.2000 3.4000 3.6000 3.8000 4.0000 4.2000 4.4000 4.6000 4.8000 5.0000 5.2000 5.4000 Columns 29 through 42 5.6000 5.8000 6.0000 6.2000 6.4000 6.6000 6.8000 7.0000 7.2000 7.4000 7.6000 7.8000 8.0000 8.2000 Columns 43 through 51 8.4000 8.6000 8.8000 9.0000 9.2000 9.4000 9.6000 9.8000 10.0000 u=0.5*(t.2)u = Columns 1 through 14 0 0.0200 0.0800 0.1800 0.3200 0.5000 0.7200 0.9800 1.2800 1.6200 2.0000 2.4200 2.8800 3.3800 Columns 15 through 28 3.9200 4.5000 5.1200 5.7800 6.4800 7.2200 8.0000 8.8200 9.6800 10.5800 11.5200 12.5000 13.5200 14.5800 Columns 29 through 42 15.6800 16.8200 18.0000 19.2200 20.4800 21.7800 23.1200 24.5000 25.9200 27.3800 28.8800 30.4200 32.0000 33.6200 Columns 43 through 51 35.2800 36.9800 38.7200 40.5000 42.3200 44.1800 46.0800 48.0200 50.0000 dlsim(numc,denc,u)9、用图形仿真方法(Simulink)分析单位加速度信号输入时闭环系统动态性能和稳态性能。 numdy=sym2poly(N)numdy =1.0e+037 * 0.8250 -2.1827 2.2623 -1.1331 0.2664 -0.0229 dendy=sym2poly(D)dendy =1.0e+036 *4.3723 -0.3141 -9.8182 6.0515 -0.0038 -0.287710、确定误差脉冲传递函数Ge(z)形式,满足单位速度信号输入时闭环稳态误差为零和实际闭环系统稳定的要求。 syms z b1 b2 b3 b4 b5 a3 a4 a5 Gezw=(1-z-1)2*(b1+b2*z-1+b3*z-2+b4*z-3+b5*z-4) Gezw = (1-1/z)2*(b1+b2/z+b3/z2+b4/z3+b5/z4)11、确定闭环脉冲传递函数Gc(z)形式,满足控制器Dw(z)可实现、无波纹、最少拍和实际闭环系统稳定的要求。 Gczw=z-1*(1+2.904*z-1)*(1-0.6703*z-1)*(1-0.3011*z-1)*(1+0.1908*z-1)*(a3+a4*z-1) Gczw = 1/z*(1+363/125/z)*(1-6703/10000/z)*(1-3011/10000/z)*(1+477/2500/z)*(a3+a4/z)12、根据10、11、列写方程组,求解Gc(z)和Ge(z)中的待定系数并最终求解Gc(z)和Ge(z) 。 f5=subs(Gczw,z,1)-1 f5 = 4184491184601/3906250000000*a3+4184491184601/3906250000000*a4-1 f6=subs(diff(Gczw,1),z,1) f6 = 3747960650717/6250000000000*a3-14736126223223/31250000000000*a4 a3j a4j=solve(f5,f6) a3j = 7195374132433105468750000/17509966474003480259529201a4j = 9150294557414550781250000/17509966474003480259529201 A=double(a3j a4j)A = 0.4109 0.5226 Gczw=subs(Gczw,a3 a4,A) Gczw =1/z*(1+363/125/z)*(1-6703/10000/z)*(1-3011/10000/z)*(1+477/2500/z)*(925332561593647/2251799813685248+2353474703694299/4503599627370496/z) f7=limit(Gezw,z,inf)-1 f7 = b1-1 f8=subs(Gezw,z,-2.904)-1 f8 =238144/131769*b1-29768000/47832147*b2+3721000000b3-465125000000/6302794178043*b4+58140625000000/2287914286629609*b5-1 f9=subs(Gezw,z,0.6703)-1 f9 =10870209/44930209*b1+108702090000/301167190927*b2+1087020900000000/2018723680783681*b3+10870209000000000000/13531504832293013743*b4+108702090000000000000000/90701676890860071119329*b5-1 f10=subs(Gezw,z,0.3011)-1 f10 =48846121/9066121*b1+488461210000/27298090331*b2+4884612100000000/82194549986641*b3+48846121000000000000/247487790009776051*b4+488461210000000000000000/745185735719435689561*b5-1 f11=subs(Gezw,z,-0.1908)-1 f11 =8862529/227529*b1-22156322500/108531333*b2+55390806250000/51769445841*b3-138477015625000000/24694025666157*b4+346192539062500000000/11779050242756889*b5-1 b1j,b2j,b3j,b4j,b5j=solve(f7,f8,f9,f10,f11) b1j =1 b2j = 27824558815573855050308402/17509966474003480259529201 b3j =1958599880988746129641979/2501423782000497179932743b4j = -19421775618281782311371837/93386487861351894717489072 b5j = -7276583151737082934420395/124515317148469192956652096 B=double(b1j,b2j,b3j,b4j,b5j)B = 1.0000 1.5891 0.7830 -0.2080 -0.0584 Gezw=subs(Gezw,b1 b2 b3 b4 b5,B) Gezw =(1-1/z)2*(1+3578267065776849/2251799813685248/z+3526291609466701/4503599627370496/z2-7492981345843367/36028797018963968/z3-2105496282460953/36028797018963968/z4) 13、求针对单位速度信号输入的最少拍无波纹控制器Dw(z)并说明Dw(z)的可实现性。 Gz=0.68436*z-1*(1+2.904*z-1)*(1-0.6703*z-1)*(1-0.3011*z-1)*(1+0.1908*z-1)/(1-z-1)2/(1-0.8187*z-1)/(1-0.3679*z-1)/(1-0.2019*z-1) Gz =17109/25000/z*(1+363/125/z)*(1-6703/10000/z)*(1-3011/10000/z)*(1+477/2500/z)/(1-1/z)2/(1-8187/10000/z)/(1-3679/10000/z)/(1-2019/10000/z) Dwz=Gczw/Gezw/Gz Dwz =25000/17109*(925332561593647/2251799813685248+2353474703694299/4503599627370496/z)/(1+3578267065776849/2251799813685248/z+3526291609466701/4503599627370496/z2-7492981345843367/36028797018963968/z3-2105496282460953/36028797018963968/z4)*(1-8187/10000/z)*(1-3679/10000/z)*(1-2019/10000/z) 14、用程序仿真方法分析单位速度信号输入时闭环系统动态性能和稳态性能。 N,D=numden(simplify(Gczw); numc=sym2poly(N)numc =1.0e+029 * 0.5783 1.9635 0.2602 -1.6051 0.1282 0.0822 denc=sym2poly(D)denc =1.0e+029 * 1.4074 0 0 0 0 0 0 t=0:0.2:2t = 0 0.2000 0.4000 0.6000 0.8000 1.0000 u=tu = 0 0.2000 0.4000 0.6000 0.8000 1.0000 dlsim(numc,denc,u)15、用图形仿真方法(Simulink)分析单位速度信号输入时闭环系统动态性能和稳态性能。 N,D=numden(simplify(Dwz); numdw=sym2poly(N)numdw =1e+027 * 1.8507 -0.2162 -2.2670 1.1602 -0.1431.0 dendw=sym2poly(D)dendw = 1.0e+027 *3.0821 4.8976 2.4133 -0.6410 -0.180116、根据8、9、14、15、的分析,说明有波纹和无波纹的差别和物理意义。 有波纹和无波纹的差别在于有波纹控制器经过一定的采样周期之后,其输入跟踪输出,但 其数字控制器的曲线会出现小的波动,该波动为波纹,而无波纹数字控制器 D(z)的输出曲 线经过一定采样周期之后会为一常数,不会产生波纹。 最少拍控制系统的系统结构简单,设计完成后容易在计算机上实现。但它对系统的适应性 差,只能保证在采样点上的输出可以跟踪输入,有可能在采样点之间呈现波纹,会引起系统的振荡。产生波纹的原因是 U(z)不能在优先个采样周期内变为 0,即 u(KT)不等于零, 使系统的输出 y(t)产生波动。实现最少拍无波纹的控制,必须要系统在典型信号的作用 下,经有限个采样周期后,系统的稳态误差保持恒值或为 0,系统的数字控制器 D(z)的 输出 u(KT)也必须保持为0。 16、根据8、9、14、15、的分析,说明有波纹和无波纹的差别和物理意义。 有波纹和无波纹的差别在于有波纹控制器经过一定的采样周期之后,其输入跟踪输出,但其数字控制器的曲线会出现小的波动,该波动为波纹,而无波纹数字控制器 D(z)的输出曲线经过一定采样周期之后会为一常数,不会产生波纹。最少拍控制系统的系统结构简单,设计完成后容易在计算机上实现。但它对系统的适应性差,只能保证在采样点上的输出可以跟踪输入,有可能在采样点之间呈现波纹,会引起系统的振荡。产生波纹的原因是U(z)不能在优先个采样周期内变为0,即u(KT)不等于零,使系统的输出 y(t)产生波动。实现最少拍无波纹的控制,必须要系统在典型信号的作用下,经有限个采样周期后,系统的稳态误差保持恒值或为0,系统的数字控制器 D(z)的输出 u(KT)也必须保持为0。 心得体会短短两个周末的实习就这样结束了。从这次实习中,我认识到了亲身实践是我们大学生活很有用也很充实的一步,通过实习能学到在课堂上学不到的很多东西。在课堂上我们仅仅知道该怎么去做,但没有亲自实践,只能是靠想象,所以有很多东西都难以理解。我发现很多知识都是以前在课本上学过的,当时印象并不怎么深刻,但是,一经实习,似乎都能很容易理解其原理,并能对其有更深的记忆。总体来说,这次实习我受益匪浅。在摸索该如何运用MATLAB仿真课程实现所需功能的过程中,培养了我的逻辑思维能力,增加了实际编程的能力。在让我体会到了分析与设计系统的艰辛的同时,更让我体会到成功的喜悦和快乐。总之,这次实验我收获颇多。实践是获得知识的一种最好的手段!在这里,我非常地感谢老师的指导,没有您的指导,我们的实习也不会进行的如此顺利;也得感谢我的同学们,感谢他们给予我的帮助,才使得我们的实习很好的完成;同时也得感谢学校、学院给我们提供这么好的实验环境。我忠心的感谢候建华老师,请允许我由衷的送上一句,老师,您辛苦了! 参考文献:MATLAB 7.0从入门到精通(修订版)刘保柱 人民邮电MATLAB基础及应用 于润伟 机械工业MATLAB在自动控制中的应用 吴晓燕 西安电子科技大附录:(程序) z=-2 -6;p=0 0 -1 -5 -8;k=668; Gs=zpk(z,p,k) Gz=c2d(Gs,0.2,zoh) a,b,c=zpkdata(Gz) Gz=zpk(a,b,c,0.2,variable,z-1) syms z b0 b1 a0 a1 a2 a3 Gez=(1-z-1)3*(1+b0*z-1) Gcz=z-1*(1+2.904*z-1)*(a0+a1*z-1+a2*z-2) f1=subs(Gcz,z,1)-1 f2=subs(diff(Gcz,1),z,1) f3=subs(diff(Gcz,2),z,1) a0j a1j a2j=solve(f1,f2,f3
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 高三试卷:江浙皖高中(县中)发展共同体2024-2025学年高三上学期10月联考数学试题
- 二零二五年度物联网设备购销合同物联网应用
- 二零二五年度房屋买卖垫资借款及贷款利率调整合同
- 二零二五年度工人工资支付与环境保护合同
- 2025版离婚后共同财产分割及债务处理三方协议书
- 二零二五年度电动四轮车总代理服务合同样本
- 二零二五年度稻谷种植保险代理服务合同
- 二零二五年度创新型科技公司股权投资合同范本
- 二零二五年度鲜活海鲜冷链物流运输合同
- 2025版农业产业链无息贷款合作协议范本
- 超星尔雅学习通《脑洞大开背后的创新思维(大连理工大学)》2025章节测试答案
- 护理服务规范与礼仪标准:护理服务规范礼仪标准及考核评分标准
- T-CRHA 089-2024 成人床旁心电监测护理规程
- 黄豆苷元药理作用研究-深度研究
- 2025年电信人工智能学习考试题库(含答案)
- 监理实施细则模板(信息化、软件工程)
- 2025年无房产证二手房交易协议书样本
- CNAS-CL01:2018 检测和校准实验室能力认可准则
- 中考名著《唐诗三百首》习题集
- 人民卫生营养与食品卫生学第8版营养与食品卫生学
- 2025年安徽省交通控股集团有限公司招聘笔试参考题库含答案解析
评论
0/150
提交评论