[0128] 下面结合附图及实施例对本发明做进一步说明。
[0129] 如图1、2所示,本实施例提供了一种光伏多水泵可达集估计与补偿协调控制方法,包括如下步骤:
[0130] 步骤S1:提供一光伏多水泵物理系统模型;其中,所述光伏多水泵物理系统模型包括光伏、多个DC/AC变换器、多个水泵和多个水管;所述每个DC/AC变换器均连接一个水泵;所述每个水泵均连接一根水管;所述光伏为所述多个DC/AC变换器及多个水泵提供电能;
[0131] 步骤S2:根据物理学原理以及中立II型T-S模糊模型,建立所述光伏多水泵物理系统的非线性动态模型;
[0132] 步骤S3:考虑到每台水泵受外界环境干扰因素不同;用统一的定性与定量方法去抑制不同的外界环境干扰显得非常困难。针对这个问题设计估计器去估计所述光伏多水泵物理系统模型的外界干扰信号;
[0133] 步骤S4:基于步骤S3建立的估计控制器所估计的外界干扰信号,设计基于补偿的反馈控制器,用以使所述外界干扰信号能够被抑制并实现光伏多水泵的稳定工作。
[0134] 在步骤S2中,根据物理学原理以及中立II型T-S模糊模型的表达方法,建立光伏多水泵物理系统的模糊动态模型。具体步骤如下:
[0135] 步骤S21:首先建立光伏单水泵非线性系统模型,如公式(1)所示:
[0136]
[0137] 式中,k6=ωe-ωr;
isd、isq表
示d轴、q轴电流; 分别表示定子d轴、q轴磁链; 分别表示转子d轴、q轴
磁链;ωr表示转子角速度;Udc是光伏发电输出电压;TL是负载转矩;ωe表示电磁场角速度;
Lm表示电枢互感; 表示转子时间常数;Lr和Rr分别是转子电感和电阻;Cdc是直流侧电容; 是漏感导数;Ls和Rs分别表示定子的电感和电阻;Cdc表示直流侧电容;
Te表示电磁转矩;p表示极对数。
[0138] 步骤S22:接着考虑光伏多水泵系统,每个水泵系统定义为角标i;根据基尔霍夫电流定理,得到:
[0139]
[0140] 将公式(2)代入(1)得到如下的光伏多水泵耦合非线性系统,
[0141]
[0142] 式中,
[0143]
[0144]
[0145]k6(i)=ωe(i)-ωr(i);
[0146] 步骤S23:将isd(i),isq(i),ωr(i),udc(i)选择作为光伏多水泵耦合非线性系统的输出,并且考虑输出测量通道卷入干扰,表达如下:
[0147]
[0148] 式中, ωi(t)是输出测量通道干扰。
[0149] 选择 作为模糊前件变量,并对其欧拉离散化得到如下光伏多水泵系统中立II型T-S模糊模型:
[0150]
[0151] 其中,是非线性Aii(t),Bi(t)Aij(t)函数 线性化后
的参数矩阵; 是中立II型的模糊集。
[0152] 在步骤S3中,考虑到每台水泵受外界环境干扰因素不同;用统一的定性与定量方法去抑制不同的外界环境干扰显得非常困难。针对这个问题设计估计器去估计这些外界的干扰信号。具体实施步骤如下:
[0153] 步骤S31:首先定义 并且引入一个变量那么随着公式(5)得到:
[0154]
[0155] 式中,
[0156]
[0157] 现在引入模糊观测器,如下:
[0158]
[0159] 式中, 并且 是一个辅助的状态矢量,是要设计的观测器增益。
[0160] 现在,我们进一步定义
[0161]
[0162]
[0163] 式中,Si是非奇异矩阵,那么得到:
[0164]
[0165] 随着公式(6)-(10)得到:
[0166]
[0167] 因为Si是非奇异矩阵,因此系统(11)可以再次表达为:
[0168]
[0169] 式中,
[0170] 步骤S32:接着定义如下的李雅普诺夫函数,
[0171]
[0172] 式中, 是正定对称矩阵。取以上李雅普诺夫函数的差分,得到:
[0173]
[0174] 由于
[0175]
[0176] 式中, 和标量k>0;
[0177] 定义正定对称矩阵 和矩阵 那么由公式(12)得到:
[0178]
[0179]
[0180] 式中,
[0181] 现在,定义如下的性能指标函数:
[0182]
[0183] 式中,α∈[0,1].
[0184] 混合(13)-(17),J(t)<0成立,假如以下的不等式成立,
[0185]
[0186] 式中,
[0187]
[0188]
[0189]
[0190]
[0191] 步骤S33:进一步为了将不等式(18)转化线性矩阵不等式,定义矩阵 为:
[0192]
[0193] 式中, 是非奇异矩阵。
[0194] 现在将(20)代入(18),并且定义 并且抽出模糊前进变量,得到:
[0195]
[0196] 式中,l∈Li, 是模糊规则,
[0197]
[0198]
[0199] 接着,不等式(21)成立可得到J(k)<0,那么
[0200] V(k+1)-1<α(V(k)-1). (23)
[0201] 从公式(23)得到:
[0202] V(k)<αk(V(0)-1)+1, (24)
[0203] 对于零初始的情况,得到:
[0204]
[0205] 式中,
[0206] 设计求解估计控制器的算法如下:
[0207] 和(21),其中 估计器的增益可以求解如下:
[0208]
[0209] 在步骤S4中,考虑到每台水泵受外界环境干扰因素不同;用统一的定性与定量方法去抑制不同的外界环境干扰显得非常困难。针对这个问题设计估计器去估计这些外界的干扰信号。具体实施步骤如下:
[0210] 步骤S41:首先考虑如下的补偿控制器:
[0211]
[0212] 式中, 和是要设计的补偿控制器增益。
[0213] 将公式(27)代入(5),得到如下的光伏多水泵闭环控制系统:
[0214]
[0215] 式中,
[0216] 步骤S42:接着考虑如下的李雅普诺夫函数:
[0217]
[0218] 式中, 是正定对称矩阵。通过取李雅普诺夫函数V(k)的差分得到:
[0219]
[0220] 定义对称正定矩阵 和矩阵 随着公式(28)得到:
[0221]
[0222] 式中,
[0223] 随着步骤S3得到 我们进一步考虑如下的性能指标函数:
[0224]
[0225] 式中,β∈[0,1].
[0226] 混合(29)-(31),以下的不等式成立,那么保证J(t)<0,
[0227]
[0228] 式中,是对称正定矩阵, 和 是适当维数矩阵,标量β
∈[0,1].
[0229] 定义,通过使用锥补引理并且抽出模糊前进变量后得到:
[0230]
[0231]
[0232] 式中,
[0233]
[0234]
[0235]
[0236]
[0237] 步骤S43:不等式(33)和(34)成立得到J(k)<0,
[0238]
[0239] 式中,
[0240] 对于零初始的情况,得到:
[0241]
[0242] 设计求解反馈补偿控制器的算法如下:
[0243] 和(33)和(34), 其中 和 是控制器的增益,那么这一系列设计步骤可以实现水泵的扰动被反向补偿,稳定光伏多水泵系统的稳定工作。
[0244] 以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。