完整版matlabPID比例微分控制KP,KI,KD讨论

武汉工程大学

实^报告LabReport

WuhanInstituteofTechnology课程：自动控制原理实验名称：控制系统的根轨迹分析学院：电气信息学院实验日期：11月28日实验报告日期：12月1日班级：2015级测控技术与仪器1-2班姓名：余文广学号： 1504200327报告退发（订正、重做）同组人：无教师审批签字 评分:）电气信息学院 武汉工程大学SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology实验名称：控制系统的根轨迹分析实验目的：初步了解基于Simulation的控制系统数字仿真方法掌握基于Matlab的控制系统的根轨迹分析方法实验任务：1）已知单位反馈控制系统的开环传递函数为G0（s）=§（§+］；§+5），试设计PID控制器，并进行单位阶跃响应分析。2）被控对象的传递函数为G2）被控对象的传递函数为G（s）=400 ，s(s2+30s+200)用Simulink建模并分析其单位阶跃响应。用MATLAB命令绘出其根轨迹图。实验步骤:k*总述：第一部分为课本上1.s（总述：第一部分为课本上1.s（s+1）（s+2）、2.s(s+2)(sA2+2s+2)以及例B-3题的Matlab根轨迹图仿真。第二部分详细讨论实验报告上的（1）（2）第三部分为总结一.熟悉基本知识Continuous（连续环节）Discontinuities（非线性环节）Discrete（离散环节）MathOperations（数学运算环节）Sinks(输出方式)PageI

电气信息学院武汉工程大学电气信息学院SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnologySources（输入源）口0上极点damping（阻尼）第一部分:根轨迹绘图第一部分:根轨迹绘图例题B-3截图图例题B-3截图图1-1»El113dl»切zpg力»rLtureCE）.说明：如上图可见，左图为零极点分布，右图为关于K的根轨迹图，中间为程序图。有两个对称的复数极点，零点一个，即可判断根轨迹有两条分支，如右图所示。同时也可根据计算计算出相应分离点渐进线交点、渐进线与实轴交点等。程序语句中pzmap（G）表示绘制零极点分布图指令，而rlocus则表示绘制根轨迹指令。PageII

电气信息学院武汉工程大学电气信息学院SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnologyk*图1-2有关s(s+1)(s+2)的Matlab根轨迹图仿真i工洞由Unri曲止

UmM毋+.H工LL」d2qmPIMP的：iL[ur■⑵iLwus'G：!.“G-T：FI[l]rI,]32DJ>>Javari"!])”TlM-lfSffii.说明：如上图可见，左图为零极点分布，右图为关于K的根轨迹图，中间为程序图。有3个极点，分别是0、-1、-2无零点。即可判断根轨迹有三条分支。可以初步判断根轨迹分离点在-1、-2之间，如右图所示。同时也可根据计算计算出相应分离点渐进线交点、渐进线与实轴交点等。程序语句中pzmap(G)表示绘制零极点分布图指令，而rlocus则表示绘制根轨迹指令。图1-3如图为根轨迹与虚轴交点说明图02由上的r要入1D7Freuncy(rad'sA:1.盟02由上的r要入1D7Freuncy(rad'sA:1.盟』 七 行 』 4T 」0 1 2 3RaalAxisfCDCcnds11.说明：如上如所示当点击根轨迹与虚轴交点时，出现如图文字框图可以看到相关数据，比如增益、阻尼等。PageIII

3 电气信息学院 武汉工程大学的Matlab的Matlab根轨迹图仿真图1-4有关 图 s(s+2)(sA2+2s+2)M-WiK常次/MATL&B电£RPde-Z&oMap7f亶SoiA>.inUlnOpi■口网wa>HijgramRsJ■|M-WiK常次/MATL&B电£RPde-Z&oMap7f亶SoiA>.inUlnOpi■口网wa>HijgramRsJ■|M1<巧为H』『IR0用|■口卢明评口工;整历由LABf.制|无皿二的2Q社H同Aks;Kccrids')二「Fnuiat。»门《3诉，miei”帧 工wn皿5bum,蘸却Lef触t■untibi«frHgILexusiJmhw"'J」■>k，、、白Wk£•二n目・四£>tJ([]].'I44401>一fL[ut*(LIPeb^i阳.fLiuri[21lLiKU5-l.5l..说明：如上图可见，左图为零极点分布，右图为关于K的根轨迹图，中间为程序图。有4个极点。即可判断根轨迹有四条分支。根轨迹，如右图所示。同时也可根据计算计算出相应分离点渐进线交点、渐进线与实轴交点等。由于更轨始发于极点而终归于零点，既可以知道如图随着k值不断增大，四条始与极点的根轨迹在k=1时聚集并分离并终于无穷远处。程序语句中pzmap（G）表示绘制零极点分布图指令，而rlocus则表示绘制根轨迹指令。第二部分，Simulink建模和仿真.建模思路：Sources（输入源） >MathOperations（数学运算环节）一Sinks—>（输出方式）一>gain（增益），可能还有控制部分和延迟部分等，其他部分。图2-2simulink中PID所在位置PageIV电气信息学院武汉工程大学电气信息学院武汉工程大学SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnologyfir&t*-SiinuliakSimiilinkLibraryBrowserIntegratorLimitedSinulim/CoAti回f：F7SilTiTllLlJcCaninonlyViedSlucksCgSimiilinkLibraryBrowserIntegratorLimitedSinulim/CoAti回f：F7SilTiTllLlJcCaninonlyViedSlucksCg匕irE.口口snssKboirdLookijpTableelifatkOparitiq^TodelYerificationHodel-VideUtilitiesPurts也SubsystemsSignalAttriButesSig^ialRoutingSinisSmiTCC3IJser-DefinedFuitcticuLs>AlJJlti>:"nalHatk&Di>Aet_i：«spaceBloclsetAu!:>:>S2?rstemTolLokCojTimjnicatiudsSystem1门口ILdxCojtinnjtlIqationsS;^etsm7oc-ILokKDL>CojTiputerVieioilSystemToijlboxContrt'lSystamToaltowHSFSystarnToolbo]-:，DSPSj/etarri7flolLa:-!HDLSupportEmbedied匚口derFuziyLogieTo>皿LCoderHTlF.ViayiPiATM^Ax+BuViriiblaTineIIelay说明：如上图，不仅仅包含了PID还有Continuous（连续环节）Discontinuities（非线性环节）Discrete（离散环节）MathOperations（数学运算环节）Sinks（输出方式）Sources（输入源）等环节.仿真流程打开simulink新建blankmodel，按设计流程建立如下图框图图2-3G图2-3G0(s)=s(s+1)(s+5)PIDsimulink仿真PageV电气信息学院武汉工程大学.SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology-Simulrnk\first*说明电气信息学院武汉工程大学.SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology-Simulrnk\first*说明:如图可见，输入为阶跃信号，函数为G0(s)= ,s(s+1)(s+5)输出示波器。图2-4图2-4无PID控制器时，示波器波形图如上图可见上升时间t约为20s时系统趋于稳定，同时纵坐标为1，纵坐标没有超出1的部分。接下来加上PID出现如下图.有关kp的测试PageVIe

a-电气信息学院e

a-电气信息学院武汉工程大学SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology图2-5KP=1、Ki=1,图2-5KP=1、Ki=1,Kd=1情况下示波器波形siI■工工晅由mn)iFaLtirrr-BfEi^nrriit1-SnmJIiniLdrt¥kw如图所示：系统为振幅逐渐变大（变大得比较小，但是如果增大t=200,就会比较明显）但是可以近似看做是系统阻尼等于一时的正弦波，波峰靠近2,波谷靠近0。可以判断此时误差大于无PID时，同时误差逐渐增大（可以认为是系统对于误差的响应过快出现了超调量）。接下来准备单独改变kp的值也即为比例微分部分，观察kp对系统性能的影响。图2-6kp=2图2-6kp=2,ki=1,kd=1图2-7kp=4,ki=1,kd=1说明：如图可见当kp=2,4时，相较与kp=1，系统逐渐稳定，并且可以发现kp=4时系统稳定的程度比kp=2优。PageVII

武汉工程大学匕 电气信息学院武汉工程大学SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology图2-8kp=0.5ki=1kp=1图图2-8kp=0.5ki=1kp=1图2-9kp=20ki=1kd=1说明：由上图以及多次试验可得当kP小于1和大于5时都会导致系统稳定性逐渐下降可见k并不能无限制增大或者减小，而是在某个适当范围内具有该善系统性能的作用。图2-10,当kp=5ki=1kd=1时系统稳曲线。±_Scope -nXFitTodeVim*SinnuhlKHiHdp *.对于kd同理：改变多组ki的值发现：ki（微分环节）大于1时系统稳定性下降，并且此时如果ki=2,kd=1再进行kp值调整仍然发现kp=5时系统稳定性最优，但是不如ki=1时。也就是改变ki并不会影响kp的作用。一下是对于ki的进一步测试。P容VIII电气信息学院武汉工程大学电气信息学院SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology图2-11kp=1ki=0.5kd=1图图2-11kp=1ki=0.5kd=1图2-14kp=1ki=0.1kd=1注意：我缩短了时间图3kp=5ki=0.5kd=1图2-12kp=1ki=0.01kd=1图图2-12kp=1ki=0.01kd=1图2-13kp=5ki=0.01kd=1PageIXP电气信息学院 武汉工程大学SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology说明：如上图依次减小ki分别为0.5、0.1、。01发现ki（积分环节）具有消除y大于1部分的能力。并且最后加上当kp=5（比例环节）时发现此时仍然具有大于y的部分但若是继续减小ki仍然会继续减小y》0的部分，但是无法消除。3.对于kd的测试图2-14kp=1ki=1kd=0.5图2-15kp=1ki=1kd=0.1RkTods同EqmulfflfcmHelp□通过以上，发现kd的减小会导致系统逐渐不稳定，（振幅逐渐增大，并且随着kd的减小增幅增大的幅度越来越大）以下对于kd逐渐增大进行讨论。图2-16图图2-16图kp=1ki=1kd=2图2-17kp=1ki=1kd=4PageX电气信息学院武汉工程大学电气信息学院SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology说明，当kd逐渐增大的时候似乎出现了与kp逐渐增大相同的情形，但是明显振幅只是减小而且周期似乎也变大了。接下来继续扩大kd图2-18kp=1ki=1kd=8图图2-18kp=1ki=1kd=8图2-19kp=1ki=1kd=16综合上面kd分别为2、4、8、16。可以发现在系统初期，出现了一个新的信号，这个信号加快了系统的调整。图2-20kp=1ki=1kd=32PageXISchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology.⑥①心・鲁SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology.⑥①心・鲁，si,□T-5O.QO0FileToolsMewEimulad口曲Help电气信息学院武汉工程大学电气信息学院如上，可以初步判断kd的上升加入了一个新的调节量，这个量使得系统的调节时间更短。但是截图也可以发现kd的增大不仅仅加入了一个新的信号也使得振幅变小，对系统有一定的调节作用。5.kp。Ki，kd总结。通过网上查阅相关资料kpkikd之间的关系。PID通过以下公式对函数进行控制。Ns1+A-有了以下的几点发现（参考网上信息）：①.比例环节的作用：成比例地反映控制系统的偏差信号e（t），偏差一旦产生，控制器立即产生控制作用，以减小偏差，（通过增大或者减小总的系统阻尼来调节系统）②积分环节的作用：主要用于消除静差（余差），提高系统的无差度。积分作用的强弱取决于积分时间常数ki，ki越大，积分作用越弱，反之则越强。③彳微分环节：反映偏差信号的变化趋势，并能在偏差信号变得太大之前，在系统中引入一个有效的早期修正信号，从而加快系统的动作速度，减少调节时间。6.图2-21有关G图2-21有关G(s)=400s(s2+30s+200)的simulink仿真PageXII匕 电气信息学院 武汉工程大学SchoolofElectricalandInformationEngineeringWuhanInstituteofTechnology400图2-22G(s)400图2-22G(s)=s(s2+30s+200)示波器的输出波形p容Xm电气信息学院武汉工程大学电气信息学院武汉工程大学SchoolofElectricalandInformationEngi