一种基于BFOA的PMSM预测控制系统参数补偿方法 |
|||||||
申请号 | CN202311515629.5 | 申请日 | 2023-11-14 | 公开(公告)号 | CN117544033B | 公开(公告)日 | 2024-05-14 |
申请人 | 江南大学; | 发明人 | 沈艳霞; 杨佳丽; 赵芝璞; | ||||
摘要 | 本 发明 公开了一种基于BFOA的PMSM预测控制系统参数补偿方法,涉及 电机 控制领域,解决了传统PMSM有限控制集模型预测 电流 控制(FCS‑MPCC)存在控制性能依赖于电机参数准确性的问题。该方法基于参数畸变对于预测电流影响的分析,设计了BFOA的多适应度函数,并结合BFOA对电机各参数补偿值进行 迭代 优化,优化后的参数补偿值可降低预测电流偏差,提高PMSM预测控制系统鲁棒性,即使在电机参数发生失配情况下仍能一定程度上保证较好的动、静态控制性能,最终结合优化补偿值计算逆变器 电路 不同 开关 状态下的代价函数,输出代价函数取值最小时的开关状态。 | ||||||
权利要求 | 1.一种基于BFOA的PMSM预测控制系统参数补偿方法,其特征在于,所述方法包括: |
||||||
说明书全文 | 一种基于BFOA的PMSM预测控制系统参数补偿方法技术领域背景技术[0002] 由于永磁同步电机(Permanent Magnet Synchronous Motor,PMSM)是一个复杂的非线性、强耦合、多变量系统,其控制难度高、成本大。近年来,预测控制、滑模控制、智能控制等多种先进的控制策略被提出并运用到PMSM控制系统中。其中,有限控制集模型预测电流控制(Finite Control Set Model Predictive Current Control,FCS‑MPCC)因其对模型要求低、控制综合质量好、在线计算方便等特点脱颖而出。但FCS‑MPCC在实际应用中仍存在明显的缺点:PMSM物理结构存在齿槽效应、边端效应和饱和效应,在运行过程中会出现温升参数变化,而预测结果依赖于系统模型参数,参数失配会导致控制性能下降。针对这一问题,主要的解决方法可归纳为两种:一种即在线参数识别,另一种为预测误差补偿方法。参数识别需要精确的电压和电流测量,在实际实施中较为复杂。预测误差补偿方法通过设计观测器、调节器等去校正预测模型的参数,但补偿限于单项参数:d、q轴电感或磁体磁链,无法同时对多个参数进行补偿,且有时会削弱系统动态响应。因此,为进一步提升PMSM的发展和应用,PMSM预测控制系统参数鲁棒性的研究已成为当前热点。 发明内容[0003] 本发明人针对上述问题及技术需求,提出了一种基于BFOA的PMSM预测控制系统参数补偿方法,来维持电机在各类参数畸变情况下控制性能的稳定性,克服传统参数补偿策略难以实现参数全覆盖、计算量大的问题。本发明的技术方案如下: [0004] 一种基于BFOA的PMSM预测控制系统参数补偿方法,包括如下步骤: [0006] 根据电机各参数补偿前后的预测电流偏差确定BFOA的多适应度函数,多适应度函数为电机各参数的补偿值的约束条件; [0009] 其进一步的技术方案为,FCS‑MPCC模型为: [0010] ; [0011] 式中, 、 表示当前采样时刻k的d、q轴电流; 、 表示下一采样时刻k+1的d、q轴预测电流值, 、 表示当前采样时刻k的d、q轴电压; 为预测控制系统的采样时长; 分别为d、q轴电感; 为定子电阻; 为转子电角速度; 为永磁体磁链。 [0012] 其进一步的技术方案为,根据模型计算电机各参数补偿前后的预测电流偏差,表达式为: [0013] ; [0014] [0015] 式中, 、 表示下一采样时刻k+1、电机各参数补偿后的d、q轴预测电流值, 、 表示下一采样时刻k+1、电机各参数补偿前发生畸变时的d、q轴预测电流值, 、 表示当前采样时刻k的d、q轴电流, 、 表示当前采样时刻k的d、q轴电压; 、 为d轴电流偏差分量, 、 为q轴电流偏差分量; 为预测控制系统的采样时长; 、 分别为d、q轴电感标定值; 是定子电阻标定值; 为转子电角速度; 为永磁体磁链标定值; 、 、 、 为电机各参数的畸变值,、 、 、 为电机各参数的补偿值。 [0016] 其进一步的技术方案为,根据电机各参数补偿前后的预测电流偏差确定BFOA的多适应度函数,包括: [0017] 将电机各参数补偿前后的预测电流偏差与未补偿时模型理论和实际的预测电流偏差进行对比,为缩小预测电流偏差从而令 、 取值趋于零, 、 部分项趋于零,得到BFOA的多适应度函数,分别记作 、 、 、 ,表达式为: [0018] 。 [0019] 其进一步的技术方案为,初始细菌种群基于BFOA进行迭代优化,包括: [0020] 初始化BFOA各参数,包括种群个数、游动步数、趋化次数、繁殖次数、迁移概率和分散次数; [0021] 所有细菌进入趋化阶段,并根据多适应度函数值选择游动或翻滚,在达到趋化次数后所有细菌进入繁殖阶段; [0022] 在繁殖阶段,根据多适应度函数值选取趋化运动更合理的细菌进行自我复制,并淘汰其余细菌,在达到繁殖次数后所有细菌进入分散阶段; [0023] 在分散阶段,所有细菌以迁移概率执行分散动作,在达到分散次数后将输出的第一个细菌的位置信息作为电机各参数的优化补偿值。 [0024] 其进一步的技术方案为,根据多适应度函数值选择游动或翻滚,包括: [0025] 当细菌个体i完成一次趋化动作后更新其位置信息,并计算其多适应度函数值; [0026] 若细菌个体i在第j+1次趋化后的每个适应度函数值均小于细菌个体i在第j次趋化后的相应适应度函数值,则细菌个体i沿当前方向游动,直到达到游动步数;否则,细菌个体i沿随机方向进行翻滚,并重新执行当细菌个体i完成一次趋化动作后更新其位置信息的步骤。 [0027] 其进一步的技术方案为,根据多适应度函数值选取趋化运动更合理的细菌进行自我复制,包括: [0028] 计算每个细菌个体所处位置对应的多适应度函数值,将属于同一个适应度函数的计算值进行排序,并选取有一半以上计算值均趋于零的预定数量的细菌个体进行自我复制。 [0029] 其进一步的技术方案为,所有细菌以迁移概率执行分散动作,包括: [0030] 当某个细菌个体满足发生迁移的概率,则该细菌个体消失的同时在空间任意位置随机生成一个新的细菌个体,且新旧细菌个体具有不同的位置信息。 [0031] 其进一步的技术方案为,结合电机各参数的优化补偿值计算逆变器电路不同开关状态下的代价函数,包括: [0032] 计算逆变器电路不同开关状态z下的d、q轴电压, ; [0033] 将电机各参数的优化补偿值和d、q轴电压输入至模型中,得到不同开关状态下的电机各参数补偿后的预测电流值,记为 、 ; [0034] 基于电机各参数补偿后的预测电流值与d、q轴参考电流值 之间的差值计算代价函数,表达式为: [0035] 。 [0036] 本发明的有益技术效果是: [0037] 本方法实现了群智能优化算法与电机参数补偿的结合,可对定子电阻、d电感、q轴电感、永磁体磁链在内的四种参数实现优化补偿,其中群智能优化算法优选细菌觅食优化算法(Bacterial Foraging Optimization Algorithm,BFOA);其次,根据FCS‑MPCC模型生成电机参数补偿值的约束条件,作为BFOA的多适应度函数,有效提高了PMSM预测控制系统的鲁棒性,提升PMSM在高温等复杂环境中的控制性能;此外,此方法可实现对预测模型补偿参数的实时优化,PMSM调速系统的动、静态响应更为快速、稳定。 [0038] 上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂,以下特举较佳实施例,并配合附图,详细说明如下。 附图说明[0040] 图2是本申请提供的基于BFOA的PMSM预测控制系统参数补偿方法流程图。 [0041] 图3是本申请提供的细菌种群基于BFOA进行迭代优化的流程图。 [0042] 图4是无参数补偿的传统预测控制系统框图。 [0043] 图5是本申请提供的定子电阻发生畸变下系统响应仿真图。 [0044] 图6是本申请提供的d轴电感发生畸变下系统响应仿真图。 [0045] 图7是本申请提供的q轴电感发生畸变下系统响应仿真图。 [0046] 图8是本申请提供的永磁体磁链发生畸变下系统响应仿真图。 [0047] 图9是本申请提供的定子电阻、d轴电感、q轴电感、永磁体磁链均发生畸变下系统响应仿真图。 [0048] 其中图5 图9中的(a) 为三种情况下的电机转速响应图,(b)为参数畸变有无补偿~下的d轴电流跟踪图,(c)为参数畸变有无补偿下的q轴电流跟踪图。 具体实施方式[0049] 下面结合附图对本发明的具体实施方式做进一步说明。 [0050] 在本实施例中,控制对象为PMSM,控制目标为PMSM达到给定转速,预测控制系统框图如图1所示。由图1可知,PMSM的FCS‑MPCC模型采用双闭环控制,外环为转速环,内环为电流环,三相逆变器为中点钳位型逆变器。外环依据转速反馈进行PI控制得到电机q轴参考电流值 ;内环预测模型应基于电机标定参数(即定子电阻、dq轴电感、永磁体磁链的标定值)得到dq轴预测电流值。不同的预测电流值会导致代价函数的取值不同,当代价函数取值最小时,输出对应逆变器开关状态,从而实现电机对于给定转速的跟踪。实际应用时,当在高温等复杂实验环境中电机实际参数与标定参数发生偏差时,细菌觅食优化算法会根据设计的多个适应度函数来输出电机参数补偿值。对预测模型进行参数补偿后,代价函数会重新判定逆变器开关状态实现对于扰动的调节。 [0051] 基于上述原理,本实施例提供了一种基于BFOA的PMSM预测控制系统参数补偿方法,如图2所示,具体包括如下步骤: [0052] 步骤1:构建PMSM的FCS‑MPCC模型。 [0053] 构建FCS‑MPCC模型是预测控制系统实现功能的第一步骤,本实施例中的FCS‑MPCC模型是由PMSM电压方程和磁链方程推导而来。 [0054] 步骤1.1:PMSM电压方程坐标变换。 [0055] PMSM三相静止坐标系下的电压方程为: [0056] (1) [0057] 式中: 为定子电阻; 为三相定子电压; 为三相定子电流;为三相定子绕组磁链,磁链随着输入电流以及转子角度 的变化而变化。 将式(1)进行clark、park变换,得到在d(直轴)‑q(交轴)旋转坐标系中的电压方程为: [0058] (2) [0059] 式中: 分别为d、q轴电压; 分别为d、q轴电流; 分别为d、q轴绕组磁链; 为转子电角速度。 [0060] 步骤1.2:PMSM磁链方程坐标变换。 [0061] PMSM三相静止坐标系下的磁链方程为: [0062] (3) [0063] 式中: 、 、 为三相绕组自感系数, 、 、 、 、 、 为两两绕组间的互感系数, 、 、 为永磁体磁场交链三相定子绕组的永磁磁链分量。将式(3)进行clark、park变换,得到在d‑q旋转坐标系中的磁链方程为: [0064] (4) [0065] 式中: 分别为d、q轴电感; 为永磁体磁链。 [0066] 步骤1.3:PMSM的FCS‑MPCC模型推导。 [0067] 将式(4)带入式(2)得到电压方程为: [0068] (5) [0069] 将式(5)中记: [0070] (6) [0071] 式中: 为预测控制系统的采样时长; 、 表示当前采样时刻k的 d、q轴电流、 、 表示下一采样时刻k+1的d、q轴预测电流值。那么将式(6)带入式(5)经过推导得到FCS‑MPCC模型为: [0072] (7) [0073] 上述模型中,定子电阻 、dq轴电感 、永磁体磁链 均为电机参数,由于温度升高、外部干扰的存在常会发生畸变。 [0074] 步骤2:根据模型计算电机各参数补偿前后的预测电流偏差。 [0075] 步骤2.1:分析参数失配影响。 [0076] 当电机参数发生畸变时,由式(7)可知,预测电流会发生变化,理论预测电流与实际预测电流之间存在偏差,即为预测电流偏差。由PMSM的参数标定值可得FCS‑MPCC理论模型为: [0077] (8) [0078] 式中: 、 是理论上d、q轴预测电流值, 、 、 、 为电机各参数标定值。存在参数失配的FCS‑MPCC实际模型为: [0079] (9) [0080] 式中: 、 为实际工况中参数发生畸变时的d、q轴预测电流值,、 、 、 为电机各参数发生畸变后的值,且: [0081] (10) [0082] 式中: 、 、 、 为电机各参数的畸变值,取值有正有负。则由于参数畸变引起的理论和实际的d、q轴预测电流偏差为: [0083](11) [0084] 步骤2.2:对FCS‑MPCC模型进行参数补偿。 [0085] 为了减小 轴预测电流偏差,需对模型进行参数补偿,补偿后的预测模型为: [0086](12) [0087] 式中: 、 表示下一采样时刻k+1、电机各参数补偿后的d、q轴预测电流值; 、 、 、 为电机各参数的补偿值。 、 、 、 并非随机 生成,应满足一定条件从而起到降低预测电流偏差的作用。 [0088] 步骤2.3:计算电机各参数补偿前后的预测电流偏差,表达式为: [0089] (13) [0090] [0091] 式中, 、 为d轴电流偏差分量, 、 为q轴电流偏差分量,上标1和2表示偏差分量标号,无实际意义。 [0092] 步骤3:根据电机各参数补偿前后的预测电流偏差确定BFOA的多适应度函数。 [0093] 将电机各参数补偿前后的预测电流偏差式(13)与未补偿时模型理论和实际的预测电流偏差式(11)进行对比,为缩小预测电流偏差从而令 、 取值趋于零, 、部分项趋于零,得到BFOA的多适应度函数即为电机各参数的补偿值的约束条件,BFOA的适应度函数由4元1次的四组方程组成,分别记作 、 、 、 ,表达式为: [0094] (14) [0095] 在生成 、 、 、 时,需满足 、 、 、 均趋于零,使得式(13)的取值小于式(11)成为可能,有效降低了d、q轴预测电流偏差。 [0096] 步骤4:随机生成初始细菌种群,并基于BFOA进行迭代优化,得到满足多适应度函数均趋于零的电机各参数的优化补偿值。 [0097] BFOA对参数补偿值进行迭代优化过程中,具体可分为趋化、繁殖、分散三个阶段。设细菌的初始种群个数为N,种群中单个细菌所处的位置即对应现实问题的优选解,每个细菌i的位置信息包含电机各参数的补偿值组成的四维向量,表示为: [0098] (15) [0099] 式中: ,细菌 在第 次趋化步骤、第v次繁殖步骤和第 次分散步骤后的位置信息可用 表示。即: [0100] (16) [0101] 此时细菌 所处位置对应的优良度可用适应度函数 、 、、 体现。如图3所示,该步骤具体包括如下分步骤: [0102] 步骤4.1:初始化BFOA各参数,包括种群个数N、游动步数NS、趋化次数Nc、繁殖次数Nre、迁移概率Ped和分散次数Ned,具体数值均可以根据实际问题情况人为设置。 [0103] 步骤4.2:所有细菌进入趋化阶段,并根据多适应度函数值选择游动或翻滚,在达到趋化次数Nc后所有细菌进入繁殖阶段。 [0104] 步骤4.2.1:细菌的趋化阶段包含游动和翻滚两种行为,当细菌个体i完成一次趋化动作后更新其位置信息 ,并计算其多适应度函数值 、、 、 。 [0105] 其中 (17) [0106] 式中: 表示在第 次趋化步骤、第v次繁殖步骤和第 次分散步骤后的细菌位置信息; 为细菌游动时的步长单位,且 ; 为翻滚选择的随机方向,具体为: [0107] (18) [0108] 式中: 为随机生成的向量,且 。 [0109] 步骤4.2.2:若细菌个体i在第j+1次趋化后的每个适应度函数值均小于细菌个体i在第j次趋化后的相应适应度函数值,即 && & ,则细菌个 体i沿当前方向游动,直到达到游动步数NS,进入步骤4.2.3;否则,细菌个体i沿随机方向进行翻滚,并重新执行步骤4.2.1。 [0110] 步骤4.2.3:判断当前趋化次数j是否大于预设趋化次数Nc,若是则进入繁殖阶段,否则令j=j+1,并重新执行步骤4.2.1。 [0111] 步骤4.3:在繁殖阶段,根据多适应度函数值选取趋化运动更合理的细菌进行自我复制,并淘汰其余细菌,在达到繁殖次数Nre后所有细菌进入分散阶段。 [0112] 具体的,计算每个细菌个体i所处位置对应的多适应度函数值 、、 、 ,将属于同一个适应度函数的计算值进行排序,并选取有 一半以上计算值均趋于零的一半细菌个体进行自我复制,剩余一半细菌个体被淘汰,则繁殖后的细菌种群便包含了双倍较好适应度的细菌个体。判断当前繁殖次数v是否大于预设繁殖次数Nre,若是则进入分散阶段,否则令v=v+1,并重新执行步骤4.2.1。 [0113] 步骤4.4:在分散阶段,所有细菌以迁移概率Ped执行分散动作,在达到分散次数Ned后将输出的第一个细菌的位置信息作为电机各参数的优化补偿值。 [0114] 具体的,当某个细菌个体满足发生迁移的概率,则该细菌个体消失的同时在空间任意位置随机生成一个新的细菌个体,且新旧细菌个体具有不同的位置信息,即不同的适应度,这使得种群部分个体能够跳出原本的位置进行游动翻转,具有更强的全局搜索能力。判断当前分散次数l是否大于预设分散次数Ned,若是则算法结束,并输出第一个细菌的四维位置信息,否则令l=l+1,并重新执行步骤4.2.1。 [0115] 步骤5:结合电机各参数的优化补偿值计算逆变器电路不同开关状态下的代价函数,输出代价函数取值最小时的开关状态。 [0116] 步骤5.1:计算逆变器电路不同开关状态z下的d、q轴电压。 [0117] 三相逆变电路具有8种开关状态: ,不同开关状态下,三相定子电压也会随之发生变化: [0118] (19) [0119] 式中: 为逆变电路直流电压源; 为第z个三相逆变电路的开关状态,且。将 进行 变换得到 : [0120] (20) [0121] 由式(20)可知 的变化会反应在d、q轴电压 中,即不同开关状态下,会有不同的d、q轴电压。 [0122] 步骤5.2:将步骤4得到的电机各参数的优化补偿值和式(20)的d、q轴电压输入至式(12)的模型中,得到不同开关状态下的电机各参数补偿后的预测电流值,记为 、。 [0123] 步骤5.3:基于电机各参数补偿后的预测电流值与d、q轴参考电流值之间的差值计算代价函数,表达式为: [0124] (21) [0125] 式中, 是d、q轴参考电流值, , 为根据转速外环PI控制所得。不同的开关状态下,代价函数取值也不同,最终输出使代价函数取得最小值的开关状态。 [0126] 为验证上述所提基于BFOA的PMSM预测控制系统参数补偿方法的有效性,将PMSM传统预测控制模型与本申请提出的结合BFOA的预测控制模型进行仿真对比。具体包括: [0128] PMSM传统预测控制系统框图如图4所示。图4中无BFOA参数补偿模块,预测模型中均为标定值,其余同图1。按照框图搭建对比仿真实验,实验中PMSM的主要参数如表1所示,BFOA的主要参数如表2所示。 [0129] 表1 永磁同步电机仿真参数 [0130] [0131] 表2 BFOA主要参数 [0132] [0133] 2)设计对照组实验。 [0134] 分别设置5组参数畸变对照实验: [0135] (1) ; [0136] (2) ; [0137] (3) ; [0138] (4) ; [0139] (5) ; [0140] 在仿真电机中更改标定参数模拟参数畸变,前4组实验将单个参数进行更改,模拟畸变,最后一组实验使4组参数同时发生畸变。 [0141] 3)分析实验结果。 [0142] 5组实验的响应仿真图分别如图5 图9所示,其中,图5‑(a)为电机转速响应图,图~5‑(b)为d轴电流跟踪图,图5‑(c)为q轴电流跟踪图,转速响应图中包含了参数未发生畸变、参数发生畸变未补偿、参数发生畸变后进行BFOA补偿这三种情况。电流跟踪图分为参数发生畸变未补偿、参数发生畸变后进行BFOA补偿两种情况,且每幅图中包含了参考电流值、预测电流值和实际电流值三种电流跟踪情况,图6 图9的情况同图5一致,在此不再赘述。 ~ [0143] 分析仿真结果,定子电阻畸变会使得电机转速发生稳态偏差,q轴预测电流值与实际电流值发生偏差;d轴电感发生畸变会使得电机转速发生脉动,d轴电流出现明显脉动,q轴预测电流值与实际电流值存在偏差;q轴电感发生畸变会使得电机转速出现超调现象,q轴电流出现畸变,实际电流跟踪预测电流时出现偏差;永磁体磁链出现畸变会影响电机转速响应速度,且转速出现稳态偏差,q轴预测电流值与实际电流值发生偏差。在进行本申请提出的参数补偿方法后,5组实验中均有补偿效果,首先电机转速在稳定性、准确性等方面均有提升,超调幅度降低;其次q轴电流脉动幅度减小,预测电流跟踪实际电流更为准确,消除了静态偏差,d轴电流脉动幅度减小,更为平滑稳定;预测控制系统动、静态性能均出现提升。 [0144] 以上所述的仅是本申请的优选实施方式,本发明不限于以上实施例。可以理解,本领域技术人员在不脱离本发明的精神和构思的前提下直接导出或联想到的其他改进和变化,均应认为包含在本发明的保护范围之内。 |