直流电机PI控制器设计与性能分析

1、课程设计题目：直流电机PI控制器设计与性能分析初始条件：WR一直流电机控制系统的方框图如图所示，其中Y为电机转速，Va为电枢电压，W为负载转矩。令电枢电压由PI控制律给出，PI表达式为：tVa=(kpe+k|Jedt),其中e=r-y。0要求完成的主要任务:(包括课程设计工作量及其技术要求，以及说明书撰写等具体要求)(1)写出以Va和W为输入的直流电机系统微分方程；(2)试求kp和ki的值，使闭环系统的特征方程的根包括-30±30j；计算系统性能指标。(3)计算在单位阶跃参考输入、单位斜坡参考输入、单位阶跃扰动输入、单位斜坡扰动输入时系统的稳态误差；(4)用Matlab验证上述设计，

2、并给出系统响应曲线，分析其跟踪性能；(5)对上述任务写出完整的课程设计说明书，说明书中必须写清楚分析的过程，附Matlab源程序或Simulink仿真模型，说明书的格式按照教务处标准书写。时间安排:任务时间(天)指导老师下达任务书，审题、查阅相关资料2分析、计算2编写程序1撰写报告2论文答辩1指导教师签名:系主任（或责任教师）签名:目录摘要11设计意义及要求11.1 设计意义11.2 设计要求12设计原理32.1 PI校正控制原理32.2 叠加原理33PI控制器的具体设计及分析53.1 写控制系统微分方程53.2 求Kp和KI的值及动态性能指标53.2.1 求Kp和KI的值53.2.2 动态性

3、能指标63.3 计算稳态误差83.3.1 单位阶跃参考输入时系统稳态误差83.3.2 单位斜坡参考输入时系统稳态误差83.3.3 单位阶跃扰动输入时系统稳态误差93.3.4 单位斜坡扰动输入时系统稳态误差103.4 验证程序和图形103.4.1 单位阶跃参考输入时系统响应103.4.2 单位斜坡参考输入时系统响应123.4.3 单位阶跃扰动输入时系统响应133.4.4 单位斜坡扰动输入时系统响应错误！未定义书签。4体会和总结错误！未定义书签。参考文献错误！未定义书签。摘要扰动是一个现实系统无法避免的因素，扰动能够影响一个的反馈系统的正常运行，这时就必须要利用系统的闭环控制自行实现调节，使系统从

4、新归为稳态。为了实现系统的自动控制，只有处理好扰动对系统的影响，系统才能实现所希望的功能。本设计是一个直流电机控制系统，电机转速是系统的输出信号，负载转矩是系统的扰动信号，目的是为了探讨参考输入和扰动输入分别对系统的影响。本次课程设计针对PI闭环反馈，做了原理分析，推导了整个系统的传递函数,通过给予典型信号作为输入及扰动，分析了系统的相关动态性能、稳态误差及稳定性。并用MATLAB软件验证自己的计算，从而进行校正、检验。最终检验结果与数值计算一样，结果令人满意，可见掌握自动控制原理的一些基本知识的重要性。关键字：自动控制扰动PI闭环反馈MATLAB1设计意义及要求1.1 设计意义本设计是一个直

5、流电机控制系统，电机转速是系统的输出信号，负载转矩是系统的扰动信号，目的是为了探讨参考输入和扰动输入分别对系统的影响，当今社会自动控制技术发展很快，而通过一学期对自动控制原理的学习，这一次的课程设计即可以很好的检验学习成果，也可以把理论知识运用到实际，为以后工作学习打下基础。1.2 设计要求图1-1直流电机控制系统的方框图一直流电机控制系统的方框图如图1-1所示，其中Y为电机转速，va为电枢电压，W为负载转矩。令电枢电压由PI控制定律求取，PI表达式为tva=(kPe+K(edt):,其中e=r-y要求完成的主要任务：(包括课程设计工作量及其技术要求，以及说明书撰写等具体要求)(1)写出以Va

6、和W为输入的直流电机系统微分方程；(2)试求kp和ki的值，使闭环系统的特征方程的根包括30±30j；计算系统性能指标。(3)计算在单位阶跃参考输入、单位斜坡参考输入、单位阶跃扰动输入、单位斜坡扰动输入时系统的稳态误差；(4)用Matlab验证上述设计，并给出系统响应曲线，分析其跟踪性能；(5)对上述任务写出完整的课程设计说明书，说明书中必须写清楚分析的过程，附Matlab源程序或Simulink仿真模型，说明书的格式按照教务处标准书写。2设计原理2.1 PI校正控制原理比例一积分控制规律的控制器，称PI控制器，输出信号m(t)同时成比例的反应输入信号e(t)及其积分，即m(t)=K

7、p&t)Kp：e(t)dt(1.1)Ti式中Kp为可调比例系数；Ti为可调积分时间常数。在串联校正时，PI控制器相当于在系统中增加了一个位于原点的开环极点，同时也增加了一个位于S左半平面的开环零点。位于原点的极点可以提高系统的型别，可以消除或减小系统的稳态误差，改善系统的稳态性能；而增加的负实数零点则用来减小系统的阻尼程度，缓和PI控制器极点对系统稳定性及动态过程产生的不利影响。只要积分时间常数Ti足够大，PI控制器对系统稳定性的不利影响可大为减弱。在控制工程实践中，PI控制器主要用来改善控制系统的稳态性能1。PI控制系统不改变系统自然频率，但可以增大系统阻尼比。适当选择开环增益和积分

8、时间常数，既可以减小系统在斜坡输入时的稳态误差，又可以使系统在阶跃输入时有满意的动态性能。2.2 叠加原理图1-1所给的系统2构图，R(S)是有用输入，称为输入信号；W(S)是扰动信号；Y(S)是系统的输出信号。为了研究有用输入作用对系统输出信号Y(S)的影响，需要求有用输入作用下的闭环传递函数Y(S)/R(S)o同样为了研究扰动作用W(S)对系统输出信号Y(S)的影响，也需要求取扰动作用下的闭环传递函数Y(S)/W(S)o此外，在控制系统的分析和设计中，还常用到在输入信号R(S)或扰动W(S)作用下，以误差信号E(S)作为输出量的闭环误差传递函数或E(S)/W(S),以此用于研究系统的稳态误

9、差。其中，以误差信号E(S)作为输出量的传递函数称为误差传递函数。叠加原理即在W(S)=0的条件下，求出输入信号作用下闭环传递函数Y(S)/R(S),在R(S)=0的条件下，求扰动输入信号作用下闭环传递函数Y(S)/W(S);随后将二者线性叠加便得到最终传递函数。3PI控制器的具体设计及分析3.1 写控制系统微分方程令图1-1中的W(S)输入为零，可以得到方程：一130CVa(S)Y(S)(3.1)S30令图1-1中的V(S)输入为零，可以得到方程：1000_W(S)=Y(S)(3.2)S30由叠加定理得：3001000Y(S)=Va(S)kW(S)(3.3)S30S30两边同时乘以S+30可

10、得：300Va(S)-1000W(S)=(S30)Y(S)(3.4)最后进行拉式反变换得：(300Va(t)-1000W(t)=y30y(t)(3.5)3.2 求Kp和KI的值及动态性能指标3.2.1 求Kp和KI的值先求PI控制器的传递函数t已知Va=kpe-kIoedt,两边同时进行拉式变换可得：Va(S)=KpE(S)?E(S)S得其相应的传递函数：D(S)Va(S)E(S)(2)求系统开环传递函数和闭环传递函数当扰动输入W(s)=0时系统的开环传递函数为:G(S)=300KPS300KlS(S30)(3.6)由开环传递函数及梅森增益公式可求得闭环传递函数为:中(S)=GO)1G(S)3

11、00KPS300KI2""S(30300KP)S300Kl(3.7)(3)求闭环特征方程及其解闭环特征方程为：S2(30300KP)S300KI-0特征根为：S=-(15150Kp)_(15150Kp)2-300K|(4)求解Kp和K|因为闭环系统特征根包括-30.30j,可得如下方程组。15+150Kp=30J15+150KP)2-300KI=900解方程得：KP=0.1,KI=63.2.2 动态性能指标当输入是阶跃函数时，本系统传递函数为:对于有实零点的二阶系统：.：,(S)2n(Sz_z(S22nSn)(3.8)其中Wn-302z=60_22根据上升时间t.，峰值时间

12、tp,最大超调量6%,调节时间ts的定义,由式分别求得二阶系统具有闭环实零极点时的欠阻尼单位阶跃响应的各项指标。冗-Q中tr二Wn.1-1n一中PWnJ-1、%=1(112r)-二(二-)212c,2*e100%(3.9),I、1(4叱)一(.:=2%)(.=-5%)l1ts=(3ln)z-'n,1,21-arctanWnr二z+Wn.1-2:arctanzWnl=，z2-2zwnw2可得:tr=0.0580s(3.10)tp=0.0842sp、.%=5.65%ts=0.1218s(=&%)ts=0.0884s(.-：'5%)I上升时间tr=0.0580s,峰值时间tP

13、=0.0842s,最大超调量6%=5.65%,在=&%调节时间ts=0.1218s3.3 计算稳态误差3.3.1 单位阶跃参考输入时系统稳态误差/1令W(S)=0,R(S)=，并将Kp=0.1、KI=6带入开环传递函数得:S60(1)G(S)=60SS(1)30于是v=1、开环增益K=60,为I型系统。静态位置误差系数为：kp=Si产6?稳态误差为:，、1Cess(二)=01kp因此当输入为单位阶跃信号时系统的稳态误差为03.3.2 单位斜坡参考输入时系统稳态误差令w=0,r(t)=t,并将Kp=0.1、K|=弗入开环传递函数得:60号1)G(S)二一S0S(1)30于是Y=1、开环增

14、益K=60,为I型系统。静态速度误差系数为:kv=g曰SG(S)=60160稳态误差为:,、1ess(二)=kv则系统对单位斜坡输入信号的稳态误差约为0.017。3.3.3 单位阶跃扰动输入时系统稳态误差一1当R=°,W=时系统在扰动下输出量的实际值为:-1000SYn(S)W(S)S260S1800而输出量的希望值为零，因此误差信号为：1000SEn(S)W(S)S260S1800一.、，一1系统在R(S)=0,W(s)二一时的稳态误差:Sc11000S1Cessn(二)二limSEn(S)=lim1S=0ssnss0S260S1800S因此系统在单位阶跃扰动转矩作用下的稳态误差为

15、_1R(S)=0,W(s)=力系统在S时的稳态误差:essng)=雪5日(5)=雪1000S-2-2S60S1800SS=0.56因此系统斜坡扰动转矩作用下的稳态误差为0.5603.4验证程序和图形3.4.1 单位阶跃参考输入时系统响应扰动W(S)=0时系统的闭环传递函数为G(s)30S1800''(S)=-21G(s)S60S1800单位阶跃参考输入时的响应，即负载转矩W(s)=0,输入为R(s)=1/s的直流电机控制系统得输出响应，MATLAB的Editor/Debugger输入程序为:num=30,1800;%分子多项式den=1,60,1800%分母多项式sys=tf(

16、num,den)ltiview('step',sys)System.那RtwHrrt(imrdf)-3D5W_：富y赋士”灯gFonkirrplhjdfi:"07CkftFDd(知16T制He闺kce)口可否3SttpRttp4nMSeeingli?Q1押Hr“，皿不丁.Fhalduf!1岩nl-3.E<Trre；seccndsj得图3所示的系统仿真响应曲线图3-1单位阶跃参考输入时系统响应曲线在仿真图3-1中可以直观看出：t=0.0524stp=0.0783sp、%=6.7%ts=0.124s(A=±2%)ts=0.0994s(:=5%)所得结果与公

17、式(3.10)的数值计算结果几乎一样；当经过一段调节时间后系统输出趋近于1,达到稳态，与输入相等，因此在单位阶跃参考输入时系统的稳态误差为零，与计算得出的数值相同求单位斜坡参考输入时响应,机控制系统得输出响应，num=30,1800;den=1,60,1800;sys=tf(num,den);t=0:0.0001:0.5;u=t;lsim(sys,u,t,0);gridxlabel('t');ylabel('y(t)');title('rampresponse')即负载转矩MATLAB的Editor/Debugger得图3-2所示的系统仿真响应曲线W(s)=0,输入为R(s)=1/S2的直流电输入程序为：%分子多项式%分母多项式%系统表示%仿真时间