用MATLAB进行控制系统的滞后

.课程设计任务书学生姓名：李超专业班级：电气1001班指导教师：刘志立工作单位：自动化学院题目:用MATLAB进行控制系统的滞后－超前校正设计初始条件：已知一单位反馈系统的开环传递函数是G(s)Ks(s1)(s2)要求系统的静态速度误差系数Kv10S1，45。要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1、MATLAB作出满足初始条件的最小K值的系统伯德图，计算系统的幅值裕量和相位裕量。2、前向通路中插入一相位滞后－超前校正，确定校正网络的传递函数。3、用MATLAB画出未校正和已校正系统的根轨迹。4、用Matlab对校正前后的系统进行仿真分析，画出阶跃响应曲线，计算其时域性能指标。5、课程设计说明书中要求写清楚计算分析的过程，列出MATLAB程序和MATLAB输出。说明书的格式按照教务处 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 书写。时间安排:任务时间（天）审题、查阅相关资料1分析、计算3编写程序1撰写报告2 论文 政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载 答辩0.5指导教师签名:年月日系主任（或责任教师）签名:年月日..摘要MATLAB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用的算法都是科研和工程计算中的最新研究成果，而前经过了各种优化和容错处理。在通常情况下，可以用它来代替底层编程语言，如C和C++。在计算要求相同的情况下，使用MATLAB的编程工作量会大大减少。MATLAB的这些函数集包括从最简单最基本的函数到诸如矩阵，特征向量、快速傅立叶变换的复杂函数。函数所能解决的问题其大致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组操作以及建模动态仿真等。此次课程设计就是利用MATLAB对一单位反馈系统进行滞后-超前校正。通过运用MATLAB的相关功能，绘制系统校正前后的伯德图、根轨迹和阶跃响应曲线，并计算校正后系统的时域性能指标。关键字：超前-滞后校正MATLAB伯德图时域性能指标..目录1滞后-超前校正设计目的和原理..............................................................................11.1滞后-超前校正设计目的...............................................................................11.2滞后-超前校正设计原理...............................................................................12滞后-超前校正的设计过程......................................................................................22.1校正前系统的参数........................................................................................22.1.1用MATLAB绘制校正前系统的伯德图..............................................32.1.2用MATLAB求校正前系统的幅值裕量和相位裕量..........................32.1.3用MATLAB绘制校正前系统的根轨迹..............................................42.1.4对校正前系统进行仿真分析.............................................................52.2滞后-超前校正设计参数计算.......................................................................62.2.1选择校正后的截止频率c................................................................62.2.2确定校正参数、T2和T1.................................................................62.3滞后-超前校正后的验证...............................................................................72.3.1用MATLAB求校正后系统的幅值裕量和相位裕量..........................72.3.2用MATLAB绘制校正后系统的伯德图..............................................82.3.3用MATLAB绘制校正后系统的根轨迹..............................................92.3.4用MATLAB对校正前后的系统进行仿真分析................................103心得体会..................................................................................................................12参考文献......................................................................................................................13..用MATLAB进行控制系统的滞后－超前校正设计滞后-超前校正设计目的和原理1.1滞后-超前校正设计目的所谓校正就是在系统不可变部分的基础上，加入适当的校正元部件，使系统满足给定的性能指标。校正 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 主要有串联校正、并联校正、反馈校正和前馈校正。确定校正装置的结构和参数的 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 主要有两类：分析法和综合法。分析法是针对被校正系统的性能和给定的性能指标，首先选择合适的校正环节的结构，然后用校正方法确定校正环节的参数。在用分析法进行串联校正时，校正环节的结构通常采用超前校正、滞后校正和滞后-超前校正这三种类型。超前校正通常可以改善控制系统的快速性和超调量，但增加了带宽，而滞后校正可以改善超调量及相对稳定度，但往往会因带宽减小而使快速性下降。滞后-超前校正兼用两者优点，并在结构设计时设法限制它们的缺点。1.2滞后-超前校正设计原理滞后-超前校正RC网络电路图如图1-1所示。图1-1滞后-超前校正RC网络..下面推导它的传递函数：R21MssC2Gcs1E(s)R11sC11R2sC2R1sC11R1C1s1R2C2s1R1C1R2C2R1C2sR1C1R2C2s2令T1R1C1,T2R2C2,T1T2R1C1R2C2R1C2,1，则1T1s1T2sGcsT1s1Ts12其中T1为超前部分的参数，T2为滞后部分。滞后-超前校正的频域设计实际是超前校正和滞后校正频域法设计的综合，基本方法是利用滞后校正将系统校正后的穿越频率调整到超前部分的最大相角处的频率。具体方法是先合理地选择截止频率c，先设计滞后校正部分，再根据已经选定的设计超前部分。应用频率法确定滞后超前校正参数的步骤：1、根据稳态性能指标，绘制未校正系统的伯德图；2、选择校正后的截止频率c；3、确定校正参数；4、确定滞后部分的参数T2；5、确定超前部分的参数T1；6、将滞后部分和超前部分的传递函数组合在一起，即得滞后-超前校正的传递函数；7、绘制校正后的伯德图，检验性能指标。滞后-超前校正的设计过程2.1校正前系统的参数..根据初始条件，调整开环传递函数：Gs0.5Ks10.5ss1当系统的静态速度误差系数Kv10S1时，0.5KKv。则K20s1满足初始条件的最小K值时的开环传递函数为Gs10s1s10.5s2.1.1用MATLAB绘制校正前系统的伯德图绘制伯德图可用命令bode(num,den)程序：num=[10];den=[0.5,1.5,1,0];bode(num,den)grid得到的伯德图如图2-1所示。图2-1校正前系统的伯德图2.1.2用MATLAB求校正前系统的幅值裕量和相位裕量..用命令margin(G)可以绘制出G的伯德图，并标出幅值裕量、相位裕量和对应的频率。用函数[kg,r,wg,wc]=margin(G)可以求出G的幅值裕量、相位裕量和幅值穿越频率。程序：num=[10];den=[0.5,1.5,1,0];G=tf(num,den);margin(G)[kg,r,wg,wc]=margin(G)得到的幅值裕量和相位裕量如图2-2所示。图2-2校正前系统的幅值裕量和相位裕量运行结果：kg=0.3000r=-28.0814wg=1.4142wc=2.4253即幅值裕量h20lg0.310.5dB，相位裕量=-28.0814o。2.1.3用MATLAB绘制校正前系统的根轨迹..MATLAB中专门提供了绘制根轨迹的有关函数。[p,z]=pzmap(num,den)的功能是绘制连续系统的零、极点图。[r,k]=rlocus(num,den)的功能是绘制k0部分的根轨迹。程序：num=[10];den=[0.5,1.5,1,0];rlocus(num,den)得到校正前系统的根轨迹如图2-3所示。图2-3校正前系统的根轨迹2.1.4对校正前系统进行仿真分析Simulink是可以用于连续、离散以及混合的线性、非线性控制系统建模、仿真和分析的软件包，并为用户提供了用方框图进行建模的图形接口，很适合于控制系统的仿真。仿真后得到的结果如图2-4和图2-5所示。..图2-4校正前系统的仿真图图2-5校正前系统仿真的阶跃响应曲线2.2滞后-超前校正设计参数计算2.2.1选择校正后的截止频率c若性能指标中对系统的快速性未提明确要求时，一般对应作为c。Gj180的频率在本题中，从图3中得，c=1.5。这样，未校正系统的相位裕量为o，与要求值仅0差+45o，这样大小的超前相角通过简单的超前校正是很容易实现的。2.2.2确定校正参数、T2和T11sin。在本题中，由超前部分应产生超前相角而定，即sin145550，因此1sin507.551sin50取11c，以使滞后相角控制在-5o以内，在本题中10.1，因此滞后部分的T215T2传递函数为s0.1。s0.01..过c,20lgGjc，作20dB/dec直线，由该直线与0dB线交点坐标确定T1T1。未校正系统的伯德图在c=1.5处的增益是13dB。所以过点（1.5，-13）画一条20dB/dec的直线，与0dB线的交点确定转折频率。经计算得，转折频率10.89，另一转T1折频率为6.7。所以超前部分的传递函数为s0.89。T1s6.7将滞后校正部分和超前校正部分的传递函数组合在一起，得滞后-超前校正的传递函数为s0.89s0.1Gcss6.7s0.01系统校正后的传递函数为10s0.89s0.1GsGcsss10.5s1s6.7s0.012.3滞后-超前校正后的验证由于校正过程中，多处采用的是近似计算，可能会造成滞后-超前校正后得到的系统的传递函数不满足题目要求的性能指标。所以需要对滞后-超前校正后的系统进行验证。下面用MATLAB求已校正系统的相角裕量和幅值裕量。2.3.1用MATLAB求校正后系统的幅值裕量和相位裕量程序：num=[10,9.9,0.89];den=[0.5,4.855,11.0985,6.8055,0.067,0];G=tf(num,den);margin(G)[kg,r,wg,wc]=margin(G)得到的校正后系统的幅值裕量和相位裕量如图2-6所示。运行结果：kg=5.9195r=47.6239..wg=3.6762wc=1.2072即校正后系统的相位裕量47.6239，KvlimsGs10满足指标。s0BodeDiagramGm=15.4dB(at3.68rad/sec),Pm=47.6deg(at1.21rad/sec)150100)Bd(edutingaM)ged(esahP500-50-100-150-90-135-180-225-270-4-3-2-101231010101010101010Frequency(rad/sec)图2-6校正后系统的幅值裕量和相位裕量假设验证结果不满足指标，重新选择校正后的截止频率，重复上述过程，直到满足性能指标为止。2.3.2用MATLAB绘制校正后系统的伯德图程序：num=[10,9.9,0.89];den=[0.5,4.855,11.0985,6.8055,0.067,0];bode(num,den)grid得到的伯德图如图2-7所示。..图2-7校正后系统的伯德图2.3.3用MATLAB绘制校正后系统的根轨迹程序：num=[10,9.9,0.89];den=[0.5,4.855,11.0985,6.8055,0.067,0];rlocus(num,den)得到的校正后系统的根轨迹如图2-8所示。图2-8校正后系统的根轨迹..2.3.4用MATLAB对校正前后的系统进行仿真分析用Simulink对校正后的系统仿真。仿真后得到的结果如图2-9和图2-10所示。图2-9校正后系统的仿真图图2-10校正后系统仿真的阶跃响应曲线用MATLAB编程计算校正后系统的时域性能指标。程序：k=10;num=conv([1,0.89],[1,0.1]);den=conv(conv(conv(conv([1,0],[1,1]),[0.5,1]),[1,6.7]),[1,0.01]);sys=tf(k*num,den);Lsys=feedback(sys,1,-1);[y,t,x]=step(Lsys);..plot(t,y);ltiview得到的阶跃响应曲线如图2-11所示。图2-11校正后阶跃响应曲线由图12可知，超调量%23.8%上升时间tr1.35s峰值时间tp2.33s调节时间ts13s其中调节时间是取的2%误差范围。..心得体会随着科学技术发展的日新月异，MATLAB已成为当今应用软件中空前活跃的领域，在生活中的应用可以说是无处不在，因此掌握MATLAB这个软件基本的使用方法对我们是十分有益的。MATLAB可用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。当然，MATLAB也可以用对反馈系统进行校正。此次课程设计的内容对一个单位反馈系统进行滞后-超前校正。回顾此次实践的整个过程，虽然只有短短的几天，但是真的在这个自己独立学习的过程中学到了好多东西。课程设计开始阶段比较顺利，但是做到计算校正后系统的时域性能指标这里时，遇到了不小麻烦，不会用MATLAB编程得阶跃响应曲线。后来，在同学的帮助下，终于看懂了课件了的那段程序，于是就求出了滞后-超前校正后系统的时域性能指标。这次的课程设计，不仅让我们更好的更深一步的了解MATLAB这个十分有用的软件，也能运用他对某一电路图进行仿真，与理论上相结合，从而进一步验证理论的正确性，也是理论运用于实践的很好的证明。与此同时，通过此次课程设计，加深了系统进行滞后-超前设计过程的理解，还掌握了用MATLAB编程计算系统时域性能指标和系统幅值裕量、相位裕量的方法。总而言之，这次的课程设计的确让我受益匪浅，还让我把许多新知识尽收囊中。..参考文献陈杰主.MATLAB宝典.电子工业出版.2007.吴天明，谢小竹等.MATLAB电力系统设计与分析.国防工业出版社.2004.李国勇.智能控制与其MATLAB实现.电子工业出版社.2005.胡寿松.自动控制原理.科学出版社.2007.王万良.自动控制原理.高等教育出版社.2008.黄坚.自动控制原理及其应用[M].高等教育出版社.2004...本科生课程设计成绩评定表姓名李超性别男专业、班级电气工程及其自动化1001班课程设计题目：用MATLAB进行控制系统的滞后－超前校正设计课程设计答辩或质疑记录：1、滞后-超前校正设计目的和原理所谓校正就是在系统不可变部分的基础上，加入适当的校正元部件，使系统满足给定的性能指标。滞后-超前校正的频域设计实际是超前校正和滞后校正频域法设计的综合，基本方法是利用滞后校正将系统校正后的穿越频率调整到超前部分的最大相角处的频率。2、怎么进行滞后－超前校正设计应用频率法确定滞后超前校正参数的步骤：根据稳态性能指标，绘制未校正系统的伯德图。选择校正后的截止频率c。确定校正参数。确定滞后部分的参数T2。确定超前部分的参数T1。将滞后部分和超前部分的传递函数组合在一起，即得滞后-超前校正的传递函数。绘制校正后的伯德图，检验性能指标。3、MATLAB的功能使用用命令可以绘制出的伯德图，并标出幅值裕量、相位裕量和对应的频率。用函数可以求出G的幅值裕量、相位裕量和幅值穿越频率。用Simulink功能可以进行系统仿真。最后根据计算得出滞后-超前校正的传递函数，再用MATLAB绘制伯德图与进行系统仿真分析确认结果。成绩评定依据：评定项目评分成绩1．选题合理、目的明确（10分）2． 设计方案 关于薪酬设计方案通用技术作品设计方案停车场设计方案多媒体教室设计方案农贸市场设计方案 正确、具有可行性、创新性（20分）3．设计结果（例如：系统设计程序、仿真程序）（20分）4．态度认真、学习刻苦、遵守纪律（15分）5．设计报告的规范化、参考文献充分（不少于5篇）（10分）6．答辩（25分）总分最终评定成绩（以优、良、中、及格、不及格评定）指导教师签字：年月日.