课程设计-用MATLAB进行控制系统的滞后-超前校正设计

武汉理工大学《自动控制原理》课程设计说明书课程设计任务书学生姓名：***专业班级：自动化0805指导教师：*****工作单位：自动化学院题目:用MATLAB进行控制系统的滞后－超前校正设计初始条件：已知一单位反馈系统的开环传递函数是)2)(1()(sssKsG要求系统的静态速度误差系数110SKv，45。要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1、MATLAB作出满足初始条件的最小K值的系统伯德图，计算系统的幅值裕量和相位裕量。2、前向通路中插入一相位滞后－超前校正，确定校正网络的传递函数。3、用MATLAB画出未校正和已校正系统的根轨迹。4、用Matlab对校正前后的系统进行仿真分析，画出阶跃响应曲线，计算其时域性能指标。5、课程设计说明书中要求写清楚计算分析的过程，列出MATLAB程序和MATLAB输出。说明书的格式按照教务处标准书写。时间安排:任务时间（天）审题、查阅相关资料2分析、计算3编写程序2撰写报告2论文答辩1指导教师签名:年月日系主任（或责任教师）签名:年月日武汉理工大学《自动控制原理》课程设计说明书摘要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确定校正参数、2T和1T.................................................................62.3滞后-超前校正后的验证...............................................................................72.3.1用MATLAB求校正后系统的幅值裕量和相位裕量..........................72.3.2用MATLAB绘制校正后系统的伯德图..............................................82.3.3用MATLAB绘制校正后系统的根轨迹..............................................92.3.4用MATLAB对校正前后的系统进行仿真分析................................103心得体会..................................................................................................................12参考文献......................................................................................................................13武汉理工大学《自动控制原理》课程设计说明书1用MATLAB进行控制系统的滞后－超前校正设计1滞后-超前校正设计目的和原理1.1滞后-超前校正设计目的所谓校正就是在系统不可变部分的基础上，加入适当的校正元部件，使系统满足给定的性能指标。校正方案主要有串联校正、并联校正、反馈校正和前馈校正。确定校正装置的结构和参数的方法主要有两类：分析法和综合法。分析法是针对被校正系统的性能和给定的性能指标，首先选择合适的校正环节的结构，然后用校正方法确定校正环节的参数。在用分析法进行串联校正时，校正环节的结构通常采用超前校正、滞后校正和滞后-超前校正这三种类型。超前校正通常可以改善控制系统的快速性和超调量，但增加了带宽，而滞后校正可以改善超调量及相对稳定度，但往往会因带宽减小而使快速性下降。滞后-超前校正兼用两者优点，并在结构设计时设法限制它们的缺点。1.2滞后-超前校正设计原理滞后-超前校正RC网络电路图如图1所示。图1滞后-超前校正RC网络武汉理工大学《自动控制原理》课程设计说明书2下面推导它的传递函数：222112122112211221111221111111)(sCRCRsCRCRCRsCRsCRsCRsCRsCRsCRsEsMsGc令1,,,21221121222111CRCRCRTTCRTCRT，则sTsTsTsTsGc21211111其中1T为超前部分的参数，2T为滞后部分的参数。滞后-超前校正的频域设计实际是超前校正和滞后校正频域法设计的综合，基本方法是利用滞后校正将系统校正后的穿越频率调整到超前部分的最大相角处的频率。具体方法是先合理地选择截止频率c，先设计滞后校正部分，再根据已经选定的设计超前部分。应用频率法确定滞后超前校正参数的步骤：1、根据稳态性能指标，绘制未校正系统的伯德图；2、选择校正后的截止频率c；3、确定校正参数；4、确定滞后部分的参数2T；5、确定超前部分的参数1T；6、将滞后部分和超前部分的传递函数组合在一起，即得滞后-超前校正的传递函数；7、绘制校正后的伯德图，检验性能指标。2滞后-超前校正的设计过程2.1校正前系统的参数武汉理工大学《自动控制原理》课程设计说明书3根据初始条件，调整开环传递函数：sssKsG5.0115.0当系统的静态速度误差系数110SKv时，vKK5.0。则120sK满足初始条件的最小K值时的开环传递函数为ssssG5.011102.1.1用MATLAB绘制校正前系统的伯德图程序：num=[10];den=[0.5,1.5,1,0];bode(num,den)grid得到的伯德图如图2所示。图2校正前系统的伯德图2.1.2用MATLAB求校正前系统的幅值裕量和相位裕量武汉理工大学《自动控制原理》课程设计说明书4用命令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)得到的幅值裕量和相位裕量如图3所示。图3校正前系统的幅值裕量和相位裕量运行结果：kg=0.3000r=-28.0814wg=1.4142wc=2.4253即幅值裕量dBh5.103.0lg20，相位裕量=-28.0814o。2.1.3用MATLAB绘制校正前系统的根轨迹武汉理工大学《自动控制原理》课程设计说明书5MATLAB中专门提供了绘制根轨迹的有关函数。[p,z]=pzmap(num,den)的功能是绘制连续系统的零、极点图。[r,k]=rlocus(num,den)的功能是绘制0k部分的根轨迹。程序：num=[10];den=[0.5,1.5,1,0];rlocus(num,den)得到校正前系统的根轨迹如图4所示。图4校正前系统的根轨迹2.1.4对校正前系统进行仿真分析Simulink是可以用于连续、离散以及混合的线性、非线性控制系统建模、仿真和分析的软件包，并为用户提供了用方框图进行建模的图形接口，很适合于控制系统的仿真。仿真后得到的结果如图5和图6所示。武汉理工大学《自动控制原理》课程设计说明书6图5校正前系统的仿真图图6校正前系统仿真的阶跃响应曲线2.2滞后-超前校正设计参数计算2.2.1选择校正后的截止频率c若性能指标中对系统的快速性未提明确要求时，一般对应180jG的频率作为c。从图3中得，c=1.5。这样，未校正系统的相位裕量为0o，与要求值仅差+45o，这样大小的超前相角通过简单的超前校正是很容易实现的。2.2.2确定校正参数、2T和1T由超前部分应产生超前相角而定，即sin1sin1。在本题中，50545，因此55.750sin150sin1取cT15112，以使滞后相角控制在-5o以内，因此1.012T，滞后部分的传递函数为01.01.0ss。武汉理工大学《自动控制原理》课程设计说明书7过ccjGlg20,，作20dB/dec直线，由该直线与0dB线交点坐标1T确定1T。未校正系统的伯德图在c=1.5处的增益是13dB。所以过点（1.5，-13）画一条20dB/dec的直线，与0dB线的交点确定转折频率。经计算得，转折频率89.011T，另一转折频率为7.61T。所以超前部分的传递函数为7.689.0ss。将滞后校正部分和超前校正部分的传递函数组合在一起，得滞后-超前校正的传递函数为01.01.