双闭环直流调速系统的分数阶PI~λ控制
来源:五一七教育网
第15卷第8期2015年3月 科学技术与工程 Vo1.15 No.8 Mar.2015 1671—1815(2015)08—0085—06 Science Technology and Engineering @2015 Sci.Tech.Engrg. 机电技术 双闭环直流调速系统的分数阶PIA控制 张海明 缪仲翠 赵静琼 (兰州交通大学机电工程学院 ,自动化与电气学院 ,兰州730070) 摘要分数阶PI 控制算法有助于比较新的、控制性能比较高的控制器。将分数阶PI 算法应用到双闭环直流调速系统 的速度控制器中,利用粒子群优化算法确定分数阶PI 控制器的最优的 、K。和 参数;并和整数阶PI控制器进行比较。通 过Matlab/Simulink仿真验证了分数阶PI 控制器的稳定性、跟随性能、抗扰动性能都优于普通整数阶PI算法,同时验证了粒子 群优化算法参数整定的有效性。 关键词分数阶P1 粒子群 双闭环调速 速度控制 跟随性 中图法分类号TM571.61; 文献标志码A 目前,高性能的调速系统在不断地发展和完 一般都用整数PID。在传统整数阶PID基础上用分 善,而工业控制领域中对调速系统的稳定性、鲁棒 数阶PID来改进和提高双闭环直流调速系统的性 性、跟随性等性能的要求在不断地提高。实际应用 能,具有一定的工程应用意义。 中,各类型的调速系统都各有优势。直流电动机是 以他励直流电动机为控制对象实现直流电机的 单输入单输出系统,它的机械特性适合于不同负载 双闭环控制,其中电流环采用整数阶PI控制,而速 的要求。因此直流调速系统由于结构简单、调速性 度调节器采用分数阶PI 控制器,并和整数阶PI控 能优良一直广泛的应用于工业现场的不同领域,例 制器的控制性能进行了比较。对于速度控制器PI 如电力牵引、采矿业等。同时直流调速系统也是其 和PI 的参数整定都采用了粒子群优化算法,并采用 它类型调速系统研究的基础,应用于直流调速系统 相同的目标函数。通过Matlab/Simulink仿真实验 的各种算法也可以扩展应用于其他类型的调速系 验证分数阶控制系统的稳定性、跟随性和鲁棒性优 统。因此进一步的研究直流调速系统时必要的。 于整数阶PI控制系统。 在直流调速系统中,整数阶PID控制在工业生 产中应用最为广泛。因为整数阶PID控制器具有结 1分数阶控制器简介 构简单、容易实现、参数调整方便、鲁棒性强等特点。 PID类型的分数阶控制器,即PI D ,其控制器 但是当被控对象参数发生变化时由于PID参数不能 算法如式(1)所示 实时修改对控制性能有很大的影响,并且在工业生 产现场中不可预测的干扰很多,要满足不同负荷变 U(s)=K。+ s +Kps ;A, >0 (1) 化或受到干扰时的控制性能要求是很困难的。以分 式(1)中 。、 、 分别为比例、积分、微分系数,A 数阶微积分理论为基础的分数阶PID控制器继承并 为分数阶微积分的阶次,(A, )∈[0,1]。当(A, 发扬了整数阶PID的优点,并弥补了不足。 )={(1,1),(1,0),(0,1),(0,0)}时可以分别得 近几年,分数阶控制器在不同的领域取得了一 到经典的{PID,PI,PD,P}控制器,可以看出整数阶 定的成果¨ 。在运动控制领域,分数阶控器大部 PID是分数阶PI 的特例。而分数阶PI 控制 分对永磁同步电动机、直流电动机的速度调节,在对 器对控制系统的性能的调节自由度更大,为系统获 直流电动机调节中直接的对速度调节 ’ 。双闭环 得更优性能提供了新的可能性 J。 调速系统在工程实际应用已非常广发,速度调节器 时域内分数阶控制器算法为 “(t)=K。e(t)+ iD一 e(t)+K。D e(t) (2) 2014年4月28日收到 国家自然科学基金(61 162016)、 式(2)中,D 代表了微积分,。=(一A, )。分数阶 甘肃省自然基金(1310RJZA037)、 微积分定义形式很多,比较常用的有以下二种: 甘肃省高等学校基本科研业务费专项资金资助 Grtinwald—Letnikov(GL)定义和Riemann.Liouville 第一作者简介:张海明(1963一),男,江苏徐州人,高级工程师。研 究方向:电动机智能控制。E—mail:gslztyzhm@163.tom。 (RL)定义是两种最常用的分数阶微积分的 定义 ] 科学技术与工程 15卷 RL定义如式(3)所示。 2分数阶数字控制器的数字实现 3 根据文献[12]可以用不同的方法和技术用IIR 。 =币 式(3)中,n是正数,厂( )是 的Gamma函数,即 F( )=I e-tt ~dt,厂( +1)= ,( ); >0。 GL定义如式(4)所示。 和FIR滤波器来近似分数阶控制器。近几年的研究 表明,用有限带宽的高阶滤波器来近似逼近分数阶 控制器可达到满意效果。Oustaloup’S在频域范围 内的近似滤波器可以在一定的频段内近似表示分数 D ) 。 (_1) ( —jh)(4) 阶控制器。假设需要逼近的频段为(O ,O ),则可 \ 厂\一 ! ± ) r(j+1)厂(0一 +1)。 式(4)中h是步长,[.]是取整符号。 从控制和信号处理的角度GL定义更实用而且 直观,尤其是在离散化的时候 ]。用式(4)离散化 处理时可直接取t ( 为采样周期),式(2)可离 散为式(5)。 “(k)=K。e(k)+KiD一 e(k)+K。D e(k)(5) 式(5)中 肌( ) 』 =i 0 (一 ㈠、 ,1 (6) 因此,离散的PI 控制算法式(6)可表示为 ) )+每毫(_1) A + ( ) (7) 式(7)表示当前控制量u(k)与过去的所有偏差值 有关,即随着时间的推移计算量、存储量大,这在实 际应用中不实用。但是式(7)可以看出分数阶微积 分的特征,例如不仅分数阶积分有记忆性而且分数 阶微分也具有记忆性。分数阶微积分具有渐变特 性,能更精确的描述系统的动态过程。在实际应用 中,可以利用短记忆原理表示如式(8)¨ 。 u( ): e( )+ ∑c e( 一 )+ ∑ e( 一 ) (8) 式(8)中当 <L/T时 :0,当k> 时v=k一 ,J/ 。L是记忆时间长度, 。’=(一1) f、 1是二项 式系数,可以用递归方式计算,即 c =l; =(1一半) ; =1,2,3,… (9) 可以看出式(8)是FIR滤波器。同样也可以用 IIR滤波器近似离散形式的分数阶PID控制器 。 以构造出连续滤波器的传递函数为 : N . Gf(s)=Kn 1 十叫^ (10) O)k O :∞b∞ ∞b∞ ∞b∞u 卜 ,, ∞b∞u ¨ ,, KO h∞ = ∞u 。 式(10)中Ot为分数阶的阶次,Ⅳ为滤波器的阶次, 一般取2~5;∞ 、 为用户选定的拟合频率的上限 和下限,一般在该区域内能较好的拟合分数阶微分 算子。根据当前的很多文献,递推滤波器Oustaloup 一般在频段(10~,10 )范围内,滤波器的阶次取 5阶 3双闭环直流调速系统的数学模型 以双闭环直流调速系统为控制对象,采用他励 直流电动机,电动机参数如表1所示。双闭环中电 流调节器采用普通PI,速度控制器用分数阶PI 代 替普通PI,主要研究系统的速度控制性能。双闭环 直流调速系统的数学模型如图1所示。为了正确评 价分数阶PI 的控制性能,与普通PI作为速度控制 器的控制效果进行比较。 表1 直流电机及控制系统的参数 Table 1 DC Motor Drive and control system[ 参数名称 参数值 参数名称 参数值 4速度控制器PI/PI 设计 PI/PI 控制器设计的关键是对参数的确定。整 数阶PI参数只要确定K 、Ki参数,而分数阶PI 控 制器不仅要确定 。、Ki参数,而且要确定分数阶次 A,比普通PI多了一个自由度,为系统获得更优性能 8期 张海明,等:双闭环直流调速系统的分数阶PI 控制 87 图1 双闭环直流电机调速系统 Fig.1 DC double closed—loop control system 提供了新的可能性。采用粒子群优化算法对PI和 PI 控制器参数确定,为了公平的评价两个控制器, 在参数优化过程中采用相同的目标函数。 4.1粒子群优化算法 作为评价指标。 PSO确定PI/PI 控制器参数流程如下¨ : (1)初始化粒子群,随机产生所有粒子(即PI/ PI 控制器参数)的位置和速度,并确定粒子的最优 位置。位置和速度公式如式(13)所示。 『 £+1 = +Clr1(P 一 f)+c2r2(G£一 £) 粒子群优化算法是基于群体智能理论的一种新 兴进化计算技术。粒子群算法通过群体微粒问的合 作与竞争而产生群体智能来指导优化搜索,具有较 强的收敛性、全局优化的特性,同时计算效率也比较 高¨ 。利用Matlab/Simulink仿真平台对PI//PI 控 制器参数用粒子群算法优化确定。双闭环的数学用 Simulink搭建,粒子群算法和双闭环数学模型之间 连接的桥梁是粒子和粒子对应的适应度,而粒子就 是要优化的PI/PI 控制器参数,适应度就是评价控 制系统的性能指标。选择合适的评价性能指标很关 键,对应实际应用的双闭环系统,即要考虑控制精 I …: + … (13) 式(13)中, 表示粒子的位置; 表示粒子的速度,W 为惯性因子,W=0.6;c。、C:为加速常数,r。、r 为[0, 1]区间的随机数;P 为粒子搜索到的最优位置,G 是整个粒子群搜索到的最优位置; (2)对于每个粒子,将其适应值(即控制系统 的性能指标)与整个粒子所经历过的最优位置Pt的 适应值进行比较,如较好则作为当前的最优位置; (3)对每个粒子,将其适应值与整个粒子群所 经历过的最优位置G,的适应值比较,如较好则作为 当前的最优位置; (4)按速度和位置公式更新粒子的速度和 位置; 度、跟随性能,又要考虑鲁棒性,同时还要考虑控制 量不能超过执行机构的运行范围。为得到满意的过 渡过程动态特性,采用误差绝对时间积分性能指标 为参数选择的最小目标函数。为了防止控制能量过 大,在目标函数中加入控制输入的平方项,因此选 用式(1 1)作为参数优化的评价指标 。 r∞ . J=J (W1 le(t)I+W2 I ( ) I dt+W3t (11) (5)若没有满足终止条件,则返回步骤(2),否 则退出得到最优解。 经确定PI和PI 控制器为式(14)所示。 为了尽可能减小超调,采用了惩罚功能,即一旦 f ( ):5.18+ , 控制 产生超调,将超调量作为最优指标的一项,此时最优 指标为: Jl ( ):7.25+ 5仿真验证分析 ,PI 控制 (14) 。。 t,=fo w I e( )l+ :l u( ) I+ W4 l e(t)I出+w3t (12) 为了分析比较整数阶PI速度控制器和分数阶 式(12)中W 、W 、加 、W 分别为权重值。通过在 MATLAB/SIMULINK仿真平台用PSO算法对PI和 PI 控制器的控制系统的控制性能,在Matlab/Simu— link环境下,按照实际的直流电机模型建模与仿真, 采用相同的电流内环,速度控制器分别采用了分数 PI 控制器的参数优化,并都利用式(11)和式(12) 88 科学技术与工程 15卷 阶PI 控制器和整数阶PI控制器。为了验证两种控 制器的跟随性和稳定性,首先以阶跃信号作为给定 输入信号,其输出响应曲线如图2所示。由图2可 知,分数阶PI 速度控制系统比普通PI的速度响应 更快速的跟随上了给定输入,而且控制过程很平滑。 为了更加直观的比较两种控制器的性能,表2列出 了主要性能指标。对于双闭环直流调速系统说,速 度控制器就是电流控制器的给定输入,而直流电动 机的电枢电流与电磁转矩成正比,因此电枢电流的 变化决定了电动机运行的动态性能。PI 控制器和 PI控制器中电动机电枢电流的动态变化如图3所 示。由图3可以看出PI控制系统中的电枢电流在 启动过程有一定的波动,而PI 控制系统的电枢电 流比较平滑。 l 800 l塑 堕堕塑婆墨篓壁堕 仿真时间/s 图2阶跃响应曲线 Fig.2 Step responses 表2跟随性能指标 Table 2 Tracking performance index 仿真时『日j/s 图3 PI/PI 控制时电枢电流响应 Fig.3 Armature current response of PI/PI system 为了进一步的验证系统的跟随性能,在两种控 制系统中分别输入方波信号。方波输入表示系统启 动到额定转速,在2 S时突然停机,在2.5 S时又启 动,枷姗㈣ 启动后又分别运行在高速、低速段,最后额定转 速运行。图4、图5分别是PI㈣ 0 控制系统和PI控制系 统速度输出响应,可以看出在每次速度变化时PI控 制系统总会有比较大的波动和振荡,而PI 控制系 统能比较平稳的跟随上给定输入的变化。图6是两 种控制系统的电枢电流,由图6可明显看出PI 控制 系统的电枢电流波形平滑动态性能较好。 2 000 l 500 旨1 000 500 0 l 2 3 4 5 6 7 8 9 l【1 f/s 图4方波输入时PI控制系统输出 Fig.4 Responses of PI system when square input 宴 g 图5方波输入时PI 控制系统输出 Fig.5 Responses of PI system when square input 为了比较两种控制系统的鲁棒性,在上述方波 输入时,4 S时给系统的电枢电压加一扰动,图7和 图8分别是突加扰动后系统速度输出响应。由图7 和图8可以看出,PI 控制系统受到扰动后速度输出 虽然有一定的超调,但是和PI控制系统相比没有振 荡,而且恢复时间短。可以看出分数阶PI 控制器 的抗扰动性比普通Pl控制器好。 综上所述,使用分数阶PI 控制器双闭环直流 调速系统的速度控制与使用普通PI控制器的系统 相比,在跟随性以及抗扰动性方面均有很大提高;采 用粒子群优化算法对于分数阶PI 控制器的各参数 的整定是非常有效的。 8期 张海明,等:双闭环直流调速系统的分数阶PI 控制 89 图6 方波输入时系统的电枢电流 Fig.6 Armature current response of PI/PI system when square input 2 000 0 500 O 图7 PI控制系统4 S扰动时速度响应 Fig.7 Responses of PI system when 4 S disturbance 图8 PI 控制系统4 s扰动时速度响应 Fig.8 Responses of PI system when 4 S disturbance 6 结论 分数阶PID是对整数阶PID的扩展。以直流电 动机控制对象采用双闭环调速,电流环采用普通整 数阶PI控制器,为了提高系统控制性能,速度环采 用了分数阶PI 控制器,并和普通整数阶PI控制器 进行了比较。对分数阶PI 和整数PI控制器的参数 整定方法都采用粒子群算法对其进行整定,并采用 了相同的优化目标函数。为了研究比较两种控制器 c【山/J 口I圳,J 稳定性、跟随性,在进行仿真实验时分别给系统输入 了阶跃、方波信号,为更进一步的比较系统的抗扰动 性能在系统运行过程中突加扰动信号。经仿真分析 分数阶PI 作为速度控制器时系统的稳定性、跟随 性以及抗扰动性能都优于整数阶PI控制器,因此对 控制性能要求高的场合可以使用分数阶控制器。 参考文献 1 Podlubny I.Fractional—order systems and PI D .controllers.IEEE Transactions on Automatic Control,1999;44(1):2O8—_214 2 Xue D,Zhao C,Chen Y Q.Fractional order PID control of a DC nlo- tor with elastic shaft:a case study.Proc of the 2006 American Con— trolconferenee,Minneapolis,Minnesota,USA,2006:3182—3187 3 Chen Yangquan.Ubiquitous fractional order controllers.Proceedings of the 2nd IFAC Symposium on Fractional Derivatives and Applica— tions(IFAC FDA06,Plenary Paper),2006:19_21 4张瑞萍,皮佑国.基于分数阶PI速度控制器的永磁同步电动机控 制.电子技术学报,2012;(27)11:69——75 Zhang Ruiping,Pi Youguo.Fractional-order PI speed controller for permanent magnet synchronous motor.Transactiona of China Electro- technical Society,2012;27(1):69—_75 5 Huang Jiacai,Li Hongsheng,Teng Fulin.Fractional order sliding mode controller for the speed control of a permanent magnet synchro— Dons motor. 24th Chinese Control and Decision Conference (CCDC),2012:1203--1208 6朱呈祥,邹云.分数阶控制研究综述.控制与决策,2009;24 (2):161—169 Zhu Chengxiang,Zou Yun.Summary of research on fractiona.1 order contro1.Control and Decision,2009;24(2):161—169 7 Barbosa R S,Jesus I S,Silva M F.Fuzzy reasoning in fractional-or- der PD controllers.Proceedings of AIC10,10th WSEAS Internation— al Conference on Applied Informatics and Communications..New As.. pects of Applied Informatics,Biomedical Electronics Informatics and Informatics and Communications,Taipei,Taiwan,2010:252—_257 8 Barbosa R S.Machado J A T.Silva M F.Time domain design of fractional differintegrators using least—squares.Signal Processing, 2006:86:2567--2581 9 Machado J AT,Analysis and design of fractional-order digital control systems.SAMS Journal of Systems Analysis,Modelling and Simula- tion,1997;27:107—122 10 Levron O A,Mathieu F.Frequency—band complex noninteger differ— entiator:Characterization and synthesis.IEEE Trans on Circuits and Systems I:Fundamental Theory and Applications2000:25—39 ,11薛定宇.控制系统计算机辅助设计.第3版.北京:清华大学出 版社,2012 Xue Dingyu.Computer aided design of control system.Beijing:Tsin- ghua University Press,2012 12陈伯时.电力拖动自动控制系统:运动控制系统.北京:机械工 业出版社,2003 Chen Boshi.Automatic control system of electirc drive:motion con- trol system.Beijing:Mechanical Industrial Publishing House,2003 (下转第98页) 98 科学技术与工程 15卷 An Improved 2W-ary Fast Scalar Multiplication Algorithm ZHANG Jie (College of Electronic and Information Engineering,Lanzhou Jiaotong University,Lanzhou 730070,P.R.China) [Abstract]The scalar multiplication is the basic operations of the elliptic CHIVe cryptography(ECC),and is a1 SO the most time consuming and vulnerable operation.According to the flaw OH using recoding and R.L for fast corn. puting scalar multiplication algorithm,the orignal algorithm is improved.Under the premise of the resisting SPAthe improved algorithm has some advantages,such as to avoid the negative number to use in an array subscript,.. to reduce the amount of 50%storage,to delete some operations about elimination after assignmentThe theoretical analysis indicates that the operation of the improved algorithm is less 2 than original algorithm.The digital veriica.f tion shows that the operation complexity of improved algorithm reduce to 17%of original algorithm. [Key words] elliptic curve cryptosystem tp p ) ≯ ) ≯ ) ≯ ) ) scalar muhiplication scalar addition side channel attack p p ≯ ≯ p (上接第89页) l3王凌,刘波.微粒群优化与调度算法.北京:清华大学出版 ler using hybrid particle swarm optimization.Journal of Computer 社,2008:16一l9 Wang Ling,Liu Bo.Particle swalT ̄optimization and scheduling[go— Applications,2009;(29)9:2483—_2486 l5史峰,王辉,胡斐,等.MATLAB智能算法3O个案例分析. rithms.Beijing:Tsinghua University Press.2008:l6—19 北京:北京航空航天大学出版社,2013:138—139 Shi Feng,Wang Hui,Hu Fei.Intelligent algorithm MATLAB 30 case 14胡海波,黄友锐.混合粒子群算法优化分数阶PID控制参数研 究.计算机应用,2009;29(9)9:2483--2486 Hu Haibo,Huang Yourui.Research of fractional order PID control— analysis.Beijing:Beijing University of Aeronautics and Astronautics Press.2013:l38—139 Fractional Order PIA Controller of Double Closed Loop DC Motor Speed Control ZHANG Hai—ming ,MIAO Zhong.cui ,ZHAO Jing.qiong (Lanzhou Jiaotong University School of Automation&Electircal Engineering . School of Mechanical and Electrical Engineering2,Lanzhou 730070,P.R.China) [Abstract] Fractional order PI D is a relatively new and high control performance controller.The Fractional or— der PI was applied to the double closed—loop DC speed regulation system of speed controller,using particles swarm optimization algorithm to determine the parameters of the controller.The stability,follow performance and disturb— ance resistance of the fractional order PI is superior to integer order P1 which was demonstrated by Matlab/Simu. 1ink.The effectiveness of particle swarm optimization algorithm parameters setting was verified either. [Key words] frational order PI formance particle swarm double closed—loop speed control following per.