发明内容
本申请通过提供一种基于NSGA-Ⅱ参数寻优的ECPT***输出稳压控制方法,采用PID控制器对半桥式ECPT***进行输出稳压控制,建立了ECPT***的广义状态空间平均模型,并基于此模型提出了一种利用NSGA-Ⅱ多目标多约束遗传算法对PID控制器的控制参数进行自动寻优的方法。该优化方法在多目标多约束条件下自动搜寻出一组最优PID控制参数,保证闭环***在该PID控制器作用下,当ECPT***负载和传输距离在一定范围内发生变化时,输出电压仍然能保持上升时间短、稳态误差小、鲁棒性能好的效果,从而提高ECPT***的鲁棒性与传输性能。
为解决上述技术问题,本申请采用以下技术方案予以实现:
一种基于NSGA-Ⅱ参数寻优的ECPT***输出稳压控制方法,其关键在于,采用广义状态空间平均法建立ECPT***的广义状态空间平均模型(简称GSSA模型);基于GSSA模型利用PID控制器构建闭环网络对***进行输出稳压控制;采用NSGA-Ⅱ多目标多约束遗传算法对所述PID控制器的控制参数进行全局寻优,得到最优控制参数;在采用NSGA-Ⅱ多目标多约束遗传算法指标计算过程中,将种群中的每个个体参数分别带入PID控制器中,PID控制器根据GSSA模型的输出值与电压设定值比较后得到的误差量生成相应控制量,将该控制量作为GSSA模型的输入量,由此构成闭环控制网络,得到***在不同个体作用下的指标值;所述NSGA-Ⅱ多目标多约束遗传算法根据这些指标值对个体进行分层以及遗传操作,通过不断迭代最终在最后一代的种群个体中选择一组控制指标最优的个体作为PID控制器的控制参数对ECPT***进行输出稳压控制。
ECPT***因其阶数高、开关非线性等特性使得其对参数变化较为敏感,而且因为结构复杂,也使其理论分析变得困难。因此,有必要建立一个精确的数学模型对其内在动力学行为进行描述与分析,并在模型的基础上设计一定的控制律来保证***输出的稳定性。
采用广义状态空间平均法对ECPT***进行建模的具体步骤为:
A1:搭建半桥式ECPT***,确定***状态变量并建立时域微分模型;
A2:利用傅里叶变换将ECPT***的时域微分模型转换为频域线性微分模型;
A3:将所述频域线性微分模型中各状态变量的傅里叶系数的实虚部依次展开并定义为ECPT***的广义状态变量x;
A4:建立以ECPT***中各变量的傅里叶系数实虚部为广义状态变量的广义状态空间平均模型:其中,A为***矩阵,B为输入矩阵,x为***的广义状态变量,u为***的直流输入,为***的广义状态变量的一阶导数。
进一步地,步骤A1中搭建的半桥式ECPT***由半桥逆变电路、阻抗匹配电路、单级谐振网络电路、LCL复合谐振网络电路、整流滤波电路以及负载依次连接而成,其中所述半桥逆变电路包括两个功率开关器件Q1和Q2,阻抗匹配电路包括电感L1、电容C1a以及电容C1b,单级谐振网络电路包括电感L2a、第一耦合极板等效电容Cs1以及第二耦合极板等效电容Cs2,LCL复合谐振网络电路包括电感L2b、电容C2以及电感L3,***直流输入Edc的正极连接开关管Q1的漏极,开关管Q1的源极连接开关管Q2的漏极,开关管Q2的源极连接***直流输入Edc的负极,开关管Q1和开关管Q2的连接点作为输出连接电感L1的一端,电感L1的另一端依次连接电感L2a、第一耦合极板等效电容Cs1、电感L2b、以及电感L3,电感L3的另一端连接整流滤波电路的一个输入端,在电感L1与电感L2a之间连接并联的电容C1a和电容C1b,并联的电容C1a和电容C1b的另一端连接开关管Q2的源极,第二耦合电容Cs2的发射极板连接开关管Q2的源极,第二耦合电容Cs2的接收极板连接整流滤波电路的另一输入端,电容C2的一端连接在电感L2b与电感L3之间,电容C2的另一端连接第二耦合电容Cs2的接收极板与整流滤波电路之间的连线上;
构建的ECPT***由发射侧和拾取侧组成,其中Cs1、Cs2表示ECPT***中耦合极板的等效电容,该ECPT***的直流输入经过半桥逆变电路和发射侧的单级谐振网络电路后得到高频交流电,在其作用下,耦合机构的两极板之间形成交互电场,在交互电场的作用下产生位移电流“流过”极板,实现电能的无线传输,接收极板接收到的电能经过拾取侧的LCL复合谐振网络电路与整流滤波电路变换后提供给用电设备。若要对***输出电压进行闭环控制,可以通过能量信号并行传输或射频传输等方式将输出电压采样信号从拾取侧传输到发射侧,采样获得的实际值与期望设定值比较后得到偏差量,将此偏差量作为控制器输入,再根据具体控制策略生成对应的控制量来控制逆变电路开关管Q1和Q2的导通与关断,由此构成闭环控制***保证输出电压稳定性。
如图2所示,从拓扑结构上来看,该***主要由半桥逆变电路、阻抗匹配电路、单级谐振网络电路、LCL复合谐振网络电路、整流滤波电路以及负载组成。其中,单级谐振网络中的电感L2a用于对第一耦合极板等效电容Cs1和第二耦合极板等效电容Cs2进行补偿;L2b、C2、L3在***拾取侧构成LCL复合谐振网络,其作用主要有以下两点:一是提高A、B端右侧网络的功率因数,二是在耦合极板电流较小的条件下,提升了流经负载RL的输出电流且由于耦合极板的电压差较小,***的运行安全性得到提高;电容C1b用于补偿其右侧网络的阻抗使其呈现纯阻性;由于LCL复合谐振网络电路的引入使***拾取侧的总阻抗增加,因此在逆变电路输出端利用电感L1及电容C1a对其右侧所有电路的等效阻抗进行匹配,提升***输出功率。
将电容C1a和电容C1b简化为电容C1,即C1=C1a+C1b,将电感L2a和电感L2b简化为电感L2,即L2=L2a+L2b,将第一耦合极板等效电容Cs1和第二耦合极板等效电容Cs2简化为电容Cs,即Cs=Cs1Cs2/(Cs1+Cs2),并分别定义RL1、RL2、RL3、Rcs为简化***中电感L1、L2、L3以及耦合极板的内阻值,将简化***中电感电流和电容电压作为状态变量,采用零次谐波分量表示ECPT***中的直流状态变量ucf,采用基波分量表示ECPT***中的交流状态变量iL1,uc1,ucs,iL2,uc2,iL3,得到ECPT***时域微分模型为:
以及各时域状态变量的傅里叶级数表达形式:
式中,ω0为基波角频率,iL1为流经电感L1的电流,uc1为电容C1两端的电压,ucs为耦合极板两端的电压,iL2为流经电感L2的电流,uc2为电容C2两端的电压,iL3为流经电感L3的电流,ucf为滤波电容Cf两端的电压,交流状态变量iL1,uc1,ucs,iL2,uc2,iL3分别用其基波分量表示,<iL1>1和<iL1>-1为电流iL1的一阶傅里叶变换系数,<uC1>1和<uC1>-1为电压uc1的一阶傅里叶变换系数,<uCs>1和<uCs>-1为电压ucs的一阶傅里叶变换系数,<iL2>1和<iL2>-1为电流iL2的一阶傅里叶变换系数,<uC2>1和<uC2>-1为电压uc2的一阶傅里叶变换系数,<iL3>1和<iL3>-1为电流iL3的一阶傅里叶变换系数,直流状态变量ucf用其零次谐波分量表示,<ucf>0为电压ucf的零阶傅里叶变换系数;
步骤A4中广义状态空间平均模型:中***的广义状态变量x(t)为:
x(t)=[Re<iL1>1,Im<iL1>1,Re<uc1>1,Im<uc1>1,
Re<ucs>1,Im<ucs>1,Re<iL2>1,Im<iL2>1,
Re<uc2>1,Im<uc2>1,Re<iL3>1,Im<iL3>1,
<ucf>0]T
式中,Re表示实部,Im表示虚部;
***矩阵A具体为:
输入矩阵B具体为:
为了使ECPT***输出电压稳定,本发明通过采用PID控制器对逆变电路开关管的通断状态进行控制来实现对ECPT***的输出稳压控制,PID控制器的传递函数一般为:其中,Kp表示比例系数,Ti表示积分时间常数,Td表示微分时间常数。在PID控制器中,Kp、Ti、Td三个参数的选择直接影响PID控制器的性能,因此PID控制器设计的关键问题是如何选择这三个参数。到目前为止,有许多PID参数设计方法,其中比较常见的是ZN方法,然而这些参数设计方法通常只适用于低阶对象,对于ECPT***这样高阶复杂的***,想要设计出一组能满足多个控制目标的PID参数,显然这些方法是不太适用的。针对PID参数整定问题的多目标性,本发明选择多目标多约束遗传算法NSGA-Ⅱ对控制器参数进行全局优化,自动搜寻出一系列非支配最优解,并从中找出一组合适的控制参数,使***在该组控制参数的作用下当***可变参数发生变化时其输出依旧能达到预期控制目标:响应快、误差小、超调小、对变化参数鲁棒性好。
基于GSSA模型利用NSGA-Ⅱ多目标多约束遗传算法对PID控制器的控制参数进行全局寻优的具体步骤包括:
S1:根据经验预估PID控制器传递函数中三个控制参数:比例系数Kp,积分时间常数Ti,微分时间常数Td的取值范围,并在此范围内随机生成大小为N的初始种群;
S2:指标计算:将初始种群中的每个个体带入由GSSA模型与PID控制器构成的闭环控制***中,根据目标函数和约束函数分别计算出每个个体所对应的目标函数值与约束函数值,即***在该PID控制器作用下的对应指标值;
S3:根据个体的指标值对种群个体进行约束性非支配排序分层,并分层计算每个个体的拥挤度距离;
S4:设置计数器Gen=1;
S5:通过选择、交叉、变异操作后得到子代种群;
S6:将子代种群与父辈种群合并形成大小为2N的临时种群;
S7:对临时种群中的每个个体进行指标计算;
S8:根据临时种群中的每个个体的指标值对临时种群个体进行约束性非支配排序分层,并分层计算其拥挤度距离;
S9:根据每个个体与第一非支配前沿的接近程度及其拥挤度距离大小,从临时种群中选择前N个个体组成新的种群,作为下一代遗传操作的父辈;
S10:判断算法是否满足结束条件,如果满足,则进入步骤S11,否则将计数器Gen加1后,跳转到步骤S5;
S11:得到PID控制器参数最优解。
所谓PID参数优化,实际就是利用算法来优化Kp、Ti、Td三个参数,其本质是基于一定的目标函数和约束条件的参数寻优问题。因此,在遗传算法里关键的一步就是确定目标函数与约束条件,它决定了一个染色体(一组PID参数)的好坏。步骤S2中选取三个目标函数作为性能评价指标,分别是:
***输出从稳态值的10%到稳态值的90%所需的上升时间t_rise:J1=t90%-t10%;
时间乘平方误差积分ess:其中e(t)为误差;
鲁棒性能指标err:在ECPT***实际工程应用中,耦合机构两极板间的相对位置经常发生变化,且***常常工作在变负载或多负载应用条件下,因此,本发明选取***的可变参数为负载电阻RL和耦合机构的等效电容Cs,设定可变参数RL和Cs的标称值分别为和可变参数RL和Cs的相对变化范围分别为ΔRL和ΔCs,设定ΔRL=±20%,ΔCs=±10%,鲁棒性能指标err为J3=J3,1+J3,2+J3,3+J3,4,其中:
式中,J3,1,J3,2,J3,3和J3,4分别表示当***负载阻值和耦合机构距离发生变化,即当可变参数(RL,Cs)从标称值到极端最坏值变化时,***输出从切换时刻起始的时间乘平方误差积分,t表示切换时刻;所设定的切换场景为:
(1)RL切换:标称值→切换值(1+20%)
Cs切换:标称值→切换值(1+10%)
(2)RL切换:标称值→切换值(1-20%)
Cs切换:标称值→切换值(1+10%)
(3)RL切换:标称值→切换值(1+20%)
Cs切换:标称值→切换值(1-10%)
(4)RL切换:标称值→切换值(1-20%)
Cs切换:标称值→切换值(1-10%)
为了提高NSGA-Ⅱ算法的计算效率与非支配最优解的实用性,本发明引入约束条件对NSGA-Ⅱ算法进行优化处理,旨在将每一代的非支配最优解都限制到其所有目标函数值均接近各自平均值的那些非支配优化解。
步骤S2中选取四个经验条件对非支配最优解的目标函数值作约束处理:
Cons1:t_rise≤6e-3;
Cons2:ess≤3e-4;
Cons3:err≤1e-1;
式中,Cons1、Cons2、Cons3分别对三个目标函数值上升时间t_rise、时间乘平方误差积分ess以及鲁棒性能指标err进行边界条件约束,Cons4对***输出在响应过程中的超调量σP进行约束,当***输出超过稳态值h(∞)时,输出的最大值h(tp)减去稳态值h(∞)的差除以稳态值h(∞)乘以百分之百得到的一个输出最大偏差比小于等于5%。
与现有技术相比,本申请提供的技术方案,具有的技术效果或优点是:***输出响应快、稳态误差小,并且当***参数在一定范围内变化时其输出仍然能有效跟踪期望值。
实施例
一种基于NSGA-Ⅱ参数寻优的ECPT***输出稳压控制方法,如图1所示,采用广义状态空间平均法建立ECPT***的广义状态空间平均(GSSA)模型;基于GSSA模型利用PID控制器构建闭环网络对***进行输出稳压控制;采用NSGA-Ⅱ多目标多约束遗传算法对所述PID控制器的控制参数进行全局寻优,得到最优控制参数;在采用NSGA-Ⅱ多目标多约束遗传算法指标计算过程中,将种群中的每个个体参数分别带入PID控制器中,PID控制器根据GSSA模型的输出值与电压设定值比较后得到的误差量生成相应控制量,将该控制量作为GSSA模型的输入量,由此构成闭环控制网络,得到***在不同个体作用下的指标值;所述NSGA-Ⅱ多目标多约束遗传算法根据这些指标值对个体进行分层以及遗传操作,通过不断迭代最终在最后一代的种群个体中选择一组控制指标最优的个体作为PID控制器的控制参数对ECPT***进行输出稳压控制。
ECPT***因其阶数高、开关非线性等特性使得其对参数变化较为敏感,而且因为结构复杂,也使其理论分析变得困难。因此,有必要建立一个精确的数学模型对其内在动力学行为进行描述与分析,并在模型的基础上设计一定的控制律来保证***输出的稳定性。
采用广义状态空间平均法对ECPT***进行建模的具体步骤为:
A1:搭建半桥式ECPT***,确定***状态变量并建立时域微分模型;
A2:利用傅里叶变换将ECPT***的时域微分模型转换为频域线性微分模型;
A3:将所述频域线性微分模型中各状态变量的傅里叶系数的实虚部依次展开并定义为ECPT***的广义状态变量x;
A4:建立以ECPT***中各变量的傅里叶系数实虚部为广义状态变量的广义状态空间平均模型:其中,A为***矩阵,B为输入矩阵,x为***的广义状态变量,u为***的直流输入。
进一步地,步骤A1中搭建的半桥式ECPT***由半桥逆变电路、阻抗匹配电路、单级谐振网络电路、LCL复合谐振网络电路、整流滤波电路以及负载依次连接而成,其中所述半桥逆变电路包括两个功率开关器件Q1和Q2,阻抗匹配电路包括电感L1、电容C1a以及电容C1b,单级谐振网络电路包括电感L2a、第一耦合极板等效电容Cs1以及第二耦合极板等效电容Cs2,LCL复合谐振网络电路包括电感L2b、电容C2以及电感L3,***直流输入Edc的正极连接开关管Q1的漏极,开关管Q1的源极连接开关管Q2的漏极,开关管Q2的源极连接***直流输入Edc的负极,开关管Q1和开关管Q2的连接点作为输出连接电感L1的一端,电感L1的另一端依次连接电感L2a、第一耦合极板等效电容Cs1、电感L2b、以及电感L3,电感L3的另一端连接整流滤波电路的一个输入端,在电感L1与电感L2a之间连接并联的电容C1a和电容C1b,并联的电容C1a和电容C1b的另一端连接开关管Q2的源极,第二耦合电容Cs2的发射极板连接开关管Q2的源极,第二耦合电容Cs2的接收极板连接整流滤波电路的另一输入端,电容C2的一端连接在电感L2b与电感L3之间,电容C2的另一端连接第二耦合电容Cs2的接收极板与整流滤波电路之间的连线上;
构建的ECPT***由发射侧和拾取侧组成,其中Cs1、Cs2表示ECPT***中耦合极板的等效电容,该ECPT***的直流输入经过半桥逆变电路和发射侧的单级谐振网络电路后得到高频交流电,在其作用下,耦合机构的两极板之间形成交互电场,在交互电场的作用下产生位移电流“流过”极板,实现电能的无线传输,接收极板接收到的电能经过拾取侧的LCL复合谐振网络电路与整流滤波电路变换后提供给用电设备。若要对***输出电压进行闭环控制,可以通过能量信号并行传输或射频传输等方式将输出电压采样信号从拾取侧传输到发射侧,采样获得的实际值与期望设定值比较后得到偏差量,将此偏差量作为控制器输入,再根据具体控制策略生成对应的控制量来控制逆变电路开关管Q1和Q2的导通与关断,由此构成闭环控制***保证输出电压稳定性。
如图2所示,从拓扑结构上来看,该***主要由半桥逆变电路、阻抗匹配电路、单级谐振网络电路、LCL复合谐振网络电路、整流滤波电路以及负载组成。其中,单级谐振网络中的电感L2a用于对第一耦合极板等效电容Cs1和第二耦合极板等效电容Cs2进行补偿;L2b、C2、L3在***拾取侧构成LCL复合谐振网络,其作用主要有以下两点:一是提高A、B端右侧网络的功率因数,二是在耦合极板电流较小的条件下,提升了流经负载RL的输出电流且由于耦合极板的电压差较小,***的运行安全性得到提高;电容C1b用于补偿其右侧网络的阻抗使其呈现纯阻性;由于LCL复合谐振网络电路的引入使***拾取侧的总阻抗增加,因此在逆变电路输出端利用电感L1及电容C1a对其右侧所有电路的等效阻抗进行匹配,提升***输出功率。
将电容C1a和电容C1b简化为电容C1,即C1=C1a+C1b,将电感L2a和电感L2b简化为电感L2,即L2=L2a+L2b,将第一耦合极板等效电容Cs1和第二耦合极板等效电容Cs2简化为电容Cs,即Cs=Cs1Cs2/(Cs1+Cs2),并分别定义RL1、RL2、RL3、Rcs为简化***中电感L1、L2、L3以及耦合极板的内阻值,将简化***中电感电流和电容电压作为状态变量,采用零次谐波分量表示ECPT***中的直流状态变量ucf,采用基波分量表示ECPT***中的交流状态变量iL1,uc1,ucs,iL2,uc2,iL3,得到ECPT***时域微分模型为:
以及各时域状态变量的傅里叶级数表达形式:
式中,ω0为基波角频率,iL1为流经电感L1的电流,uc1为电容C1两端的电压,ucs为耦合极板两端的电压,iL2为流经电感L2的电流,uc2为电容C2两端的电压,iL3为流经电感L3的电流,ucf为滤波电容Cf两端的电压,交流状态变量iL1,uc1,ucs,iL2,uc2,iL3分别用其基波分量表示,<iL1>1和<iL1>-1为电流iL1的一阶傅里叶变换系数,<uC1>1和<uC1>-1为电压uc1的一阶傅里叶变换系数,<uCs>1和<uCs>-1为电压ucs的一阶傅里叶变换系数,<iL2>1和<iL2>-1为电流iL2的一阶傅里叶变换系数,<uC2>1和<uC2>-1为电压uc2的一阶傅里叶变换系数,<iL3>1和<iL3>-1为电流iL3的一阶傅里叶变换系数,直流状态变量ucf用其零次谐波分量表示,<ucf>0为电压ucf的零阶傅里叶变换系数;
利用傅里叶变换将ECPT***时域微分模型转换为频域线性微分模型如下:
式中,Edc为***的直流输入,<Edc>0为直流输入Edc的零阶傅里叶变换系数,g(t)、p(t)为非线性转换函数,分别表示发射侧半桥式逆变环节和拾取侧整流环节的转换函数,转换函数的定义如下:
式中,T为逆变器工作周期,θ1、θ2为函数的初相角,n为周期的整数倍。
步骤A4中广义状态空间平均模型:中***的广义状态变量x(t)为:
x(t)=[Re<iL1>1,Im<iL1>1,Re<uc1>1,Im<uc1>1,
Re<ucs>1,Im<ucs>1,Re<iL2>1,Im<iL2>1,
Re<uc2>1,Im<uc2>1,Re<iL3>1,Im<iL3>1,
<ucf>0]T
式中,Re表示实部,Im表示虚部;
***矩阵A具体为:
输入矩阵B具体为:
为了使ECPT***输出电压稳定,本发明通过采用PID控制器对逆变电路开关管的通断状态进行控制来实现对ECPT***的输出稳压控制,PID控制器的传递函数一般为:其中,Kp表示比例系数,Ti表示积分时间常数,Td表示微分时间常数。在PID控制器中,Kp、Ti、Td三个参数的选择直接影响PID控制器的性能,因此PID控制器设计的关键问题是如何选择这三个参数。到目前为止,有许多PID参数设计方法,其中比较常见的是ZN方法,然而这些参数设计方法通常只适用于低阶对象,对于ECPT***这样高阶复杂的***,想要设计出一组能满足多个控制目标的PID参数,显然这些方法是不太适用的。针对PID参数整定问题的多目标性,本发明选择多目标多约束遗传算法NSGA-Ⅱ对控制器参数进行全局优化,自动搜寻出一系列非支配最优解,并从中找出一组合适的控制参数,使***在该组控制参数的作用下当***可变参数发生变化时其输出依旧能达到预期控制目标:响应快、误差小、超调小、对变化参数鲁棒性好。
如图3所示,基于GSSA模型利用NSGA-Ⅱ多目标多约束遗传算法对PID控制器的控制参数进行全局寻优的具体步骤包括:
S1:根据经验预估PID控制器传递函数中三个控制参数:比例系数Kp,积分时间常数Ti,微分时间常数Td的取值范围,并在此范围内随机生成大小为N的初始种群;
S2:指标计算:将初始种群中的每个个体带入由GSSA模型与PID控制器构成的闭环控制***中,根据目标函数和约束函数分别计算出每个个体所对应的目标函数值与约束函数值,即***在该PID控制器作用下的对应指标值;
S3:根据个体的指标值对种群个体进行约束性非支配排序分层,并分层计算每个个体的拥挤度距离;
S4:设置计数器Gen=1;
S5:通过选择、交叉、变异操作后得到子代种群;
S6:将子代种群与父辈种群合并形成大小为2N的临时种群;
S7:对临时种群中的每个个体进行指标计算;
S8:根据临时种群中的每个个体的指标值对临时种群个体进行约束性非支配排序分层,并分层计算其拥挤度距离;
S9:根据每个个体与第一非支配前沿的接近程度及其拥挤度距离大小,从临时种群中选择前N个个体组成新的种群,作为下一代遗传操作的父辈;
S10:判断算法是否满足结束条件,如果满足,则进入步骤S11,否则将计数器Gen加1后,跳转到步骤S5;
S11:得到PID控制器参数最优解。
所谓PID参数优化,实际就是利用算法来优化Kp、Ti、Td三个参数,其本质是基于一定的目标函数和约束条件的参数寻优问题。因此,在遗传算法里关键的一步就是确定目标函数与约束条件,它决定了一个染色体(一组PID参数)的好坏。步骤S2中选取三个目标函数作为性能评价指标,分别是:
***输出从稳态值的10%到稳态值的90%所需的上升时间t_rise:J1=t90%-t10%;
时间乘平方误差积分ess:其中e(t)为误差;
鲁棒性能指标err:在ECPT***实际工程应用中,耦合机构两极板间的相对位置经常发生变化,且***常常工作在变负载或多负载应用条件下,因此,本发明选取***的可变参数为负载电阻RL和耦合机构的等效电容Cs,设定可变参数RL和Cs的标称值分别为和可变参数RL和Cs的相对变化范围分别为ΔRL和ΔCs,设定ΔRL=±20%,ΔCs=±10%,鲁棒性能指标err为J3=J3,1+J3,2+J3,3+J3,4,其中:
式中,J3,1,J3,2,J3,3和J3,4分别表示当***负载阻值和耦合机构距离发生变化,即当可变参数(RL,Cs)从标称值到极端最坏值变化时,***输出从切换时刻起始的时间乘平方误差积分,t表示切换时刻;所设定的切换场景为:
(1)RL切换:标称值→切换值(1+20%)
Cs切换:标称值→切换值(1+10%)
(2)RL切换:标称值→切换值(1-20%)
Cs切换:标称值→切换值(1+10%)
(3)RL切换:标称值→切换值(1+20%)
Cs切换:标称值→切换值(1-10%)
(4)RL切换:标称值→切换值(1-20%)
Cs切换:标称值→切换值(1-10%)
为了提高NSGA-Ⅱ算法的计算效率与非支配最优解的实用性,本发明引入约束条件对NSGA-Ⅱ算法进行优化处理,旨在将每一代的非支配最优解都限制到其所有目标函数值均接近各自平均值的那些非支配优化解。
步骤S2中选取四个经验条件对非支配最优解的目标函数值作约束处理:
Cons1:t_rise≤6e-3;
Cons2:ess≤3e-4;
Cons3:err≤1e-1;
式中,Cons1、Cons2、Cons3分别对三个目标函数值上升时间t_rise、时间乘平方误差积分ess以及鲁棒性能指标err进行边界条件约束,Cons4对***输出在响应过程中的超调量σP进行约束,当***输出超过稳态值h(∞)时,输出的最大值h(tp)减去稳态值h(∞)的差除以稳态值h(∞)乘以百分之百得到的一个输出最大偏差比小于等于5%。
为了进一步了解本发明的显著效果,现将基于NSGA-Ⅱ算法优化的PID控制器带入ECPT闭环***进行仿真及实验分析,具体如下:
NSGA-Ⅱ算法的参数具体设置为:种群大小100,遗传代数30,交叉概率0.9,变异概率0.1。图4为当遗传代数为30时,100个种群个体的目标值分布情况的仿真图,由仿真结果图可以看出,当种群进化到30代时,其非支配最优解均匀分布在最优解集中,种群中每个个体的三个目标值均达到较小值,控制性能效果较好。在100个种群个体中选择等级为1,拥挤度距离较大的一个非支配最优解作为最终选定的最优解Sopt:
Kp=8.8654471e-1,Ti=0.002,Td=0
即选用(0.887,0.002,0)作为此次PID控制参数优化的最终解。
为了验证基于NSGA-Ⅱ算法优化的PID控制器对ECPT***输出控制效果,本实施例基于上述所提出的PID控制器参数寻优方法及目标函数和约束条件设计出一组最优PID控制参数(0.887,0.002,0),将其带入闭环控制***进行仿真分析和实验验证。在图2所示的***电路拓扑基础上,加入PID控制回路,根据***期望输出电压值与采样得到的实际值比较后获得偏差及偏差变化率,经过PID控制器模块运算后得到相应的控制量,占空比计算模块将其转换为逆变电路开关管的控制信号,从而构成闭环控制***,保证输出电压稳压控制。
令输出电压设定值为10V,在Matlab/Simulink平台下对闭环***进行仿真分析,***电路参数取值如表1所示。
在仿真运行中利用开关器件同时切换可变参数RL、Cs,使负载RL变化±20%,Cs变化±10%,分别模拟四种参数变化的最坏情况。负载两端电压u0的仿真波形如图5、6、7、8所示。其中,图5为(20Ω,1nF)到(24Ω,1.1nF)跳变时负载电压u0的仿真波形图,图6为(20Ω,1nF)到(16Ω,1.1nF)跳变时负载电压u0的仿真波形图,图7为(20Ω,1nF)到(24Ω,0.9nF)跳变时负载电压u0的仿真波形图,图8为(20Ω,1nF)到(16Ω,0.9nF)跳变时负载电压u0的仿真波形图。从以上四幅图可以看出,当***耦合机构等效电容Cs以及负载电阻RL在允许变化范围内发生四种最坏情况变化时,闭环***输出经过最快8ms最长20ms左右的调节过程后均能稳定下来。而在***启动过程,输出电压经过8ms左右从稳态值的10%上升到稳态值的90%,与图4中的优化目标值一致。从而验证了经过NSGA-Ⅱ算法优化后的PID控制器的控制效果:***输出上升速度快,超调量小,且闭环***在可变参数扰动情况下鲁棒性能好。
基于图2的拓扑结构和表1的***参数,搭建了半桥式ECPT***实验平台。本实施例选用STM32F407单片机作为***的控制器实现PID控制器运算与占空比计算。单片机输出逆变电路开关管的控制信号,驱动电路将其转换为实际的驱动脉冲。***启动过程的负载输出电压u0、Q1管栅极触发脉冲uGS、逆变器输出电压u1实验波形如图9所示。从图中可以看出,Q1管栅极触发脉冲uGS作为闭环***的直接控制信号,当***开始运行时便根据PID控制器输出的控制量转换为占空比可变的触发脉冲信号。***在半桥逆变环节的作用下,逆变器输出电压u1为占空比可变的方波输出,其占空比与触发脉冲uGS的占空比一致。与仿真结果类似,PID控制器经过12ms左右的动态调节使***输出电压从1V上升到9V,且进入稳态后,其稳态误差基本为零。在整个动态过程中***超调很小,调节时间与仿真基本一致。由此说明,在实际应用中,该PID控制器使***有较好的动态特性和静态特性。
当***负载阻值和耦合机构等效电容发生变化时,***输出电压u0与流经电感L3的电流iL3的实验波形如图10、11、12、13所示。其中,图10为(20Ω,1nF)到(24Ω,1.1nF)跳变时***输出电压u0与流经电感L3的电流iL3的实验波形图,图11为(20Ω,1nF)到(16Ω,1.1nF)跳变时***输出电压u0与流经电感L3的电流iL3的实验波形图,图12为(20Ω,1nF)到(24Ω,0.9nF)跳变时***输出电压u0与流经电感L3的电流iL3的实验波形图,图13为(20Ω,1nF)到(16Ω,0.9nF)跳变时***输出电压u0与流经电感L3的电流iL3的实验波形图。图10-13的实验结果与图5-8的仿真结果基本一致,当***参数发生四种最坏情况跳变时,输出电压均能经过一个较短的调节过程后进入10V稳态阶段,调节时间最快为12ms,最长为30ms左右。当负载阻值发生变化时,流经L3的电流iL3会发生相应变化来保证***输出电压值恒定。当负载阻值增大时,iL3的幅值会相应减小,反之,当负载阻值减小时,iL3的幅值会相应增大以保证u0值不变。
本申请的上述实施例中,一种基于NSGA-Ⅱ参数寻优的ECPT***输出稳压控制方法,采用PID控制器对半桥式ECPT***进行输出稳压控制,建立了ECPT***的广义状态空间平均模型,并基于此模型提出了一种利用多目标多约束遗传算法对PID控制器的控制参数进行自动寻优的方法。该优化方法在多目标多约束条件下自动搜寻出一组最优PID控制参数,保证闭环***在该PID控制器作用下,当ECPT***负载和传输距离在一定范围内发生变化时,输出电压仍然能保持上升时间短、稳态误差小、鲁棒性能好的效果,从而提高ECPT***的鲁棒性与传输性能。
应当指出的是,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的普通技术人员在本发明的实质范围内所做出的变化、改性、添加或替换,也应属于本发明的保护范围。