CN104184165A - 一种适用于无主控中心的光伏发电机组协同控制方法 - Google Patents

Abstract

本发明公开了一种适用于无主控中心的光伏发电机组协同控制方法。引入辅助变量迭代计算，得到光伏发电机组的输出功率；通过高层控制器计算其储备率指令值；对每台光伏发电机输出的最大功率进行近似计算，再计算储备率的实测值；建立通信矩阵，通过底层控制器计算储备率的参考值以及并网逆变器直流母线电压的参考值；通过光伏并网逆变器电压电流双闭环PI控制方法，实现每台光伏发电机输出功率的稳定控制。本发明无需设置控制光伏发电机组的主控中心，节省通信设备和成本，提高运行的鲁棒性，实现“即插即用”功能；实现地理分散的大量光伏发电机输出功率的协同控制，为平滑电力网络的负载波动、提高电力网络中电压和频率的稳定性提供支持。

Description

一种适用于无主控中心的光伏发电机组协同控制方法

技术领域

本发明涉及一种光伏发电机组控制方法，尤其是涉及一种适用于无主控中心的光伏发电机组协同控制方法。

背景技术

太阳能具有无污染、可再生和储量大等优点，因此太阳能光伏发电必将得到广泛应用。然而，由于光伏发电机输出功率具有随机性和波动性，当大量的光伏发电机接入电力网络，会对电力网络的安全性和可靠性产生不利影响。因此，有必要研发一种协调光伏发电机输出功率的控制方法，统一控制有数量较多的光伏发电机组成的光伏发电机组的输出功率。

现有的光伏发电机组输出功率控制方法，包括集中式控制和分散式控制两种。集中式控制方法实现简单，并可以使光伏发电机组的输出功率快速收敛至给定值，然而该方法需要设置控制光伏发电机组的主控中心，而且需要所有光伏发电机都与主控中心实现可靠通信，这提高了通信成本，降低了光伏发电机组控制的可靠性和鲁棒性。

分散式控制方法无需主控中心，可以实现每台光伏发电机的自治控制，并可以实现光伏发电机的“即插即用”功能。然而，当电力网络中存在数量较多的光伏发电机时，该方法无法实现功率的有效分配与控制，因此无法为电力网络提供辅助服务，如平滑负载、频率和电压的控制等。

发明内容

为解决上述问题，本发明提出一种适用于无主控中心的光伏发电机组协同控制方法，以实现光伏发电机组输出功率的协调控制、提高光伏发电机组运行可靠性和鲁棒性、实现光伏发电机的“即插即用”功能、减小通信成本；具有综合集中式控制方法和分散式控制方法的优点，同时避免二者缺点的光伏发电机组输出功率的协同控制方法，实现光伏发电机组中每台发电机输出功率的协调控制，并为电力网络提供必要的辅助服务。

本发明的技术方案采用如下步骤：

1)通过对各个光伏发电机输出的有功功率引入辅助变量进行迭代计算，得到由多台光伏发电机组成的光伏发电机组的输出功率；

2)通过高层控制器计算得到光伏发电机组的储备率的指令值；

3)对光伏发电机组中每台光伏发电机输出的最大功率进行近似计算，然后计算光伏发电机组中每台光伏发电机的储备率的实测值；

4)建立通信矩阵，通过底层控制器计算光伏发电机组中每台光伏发电机的储备率的参考值，然后计算每台光伏发电机的并网逆变器直流母线电压u_dc的参考值；

5)通过光伏并网逆变器电压电流双闭环PI控制方法，实现光伏发电机组中每台光伏发电机输出功率的稳定控制。

所述的步骤1)中的多台光伏发电机组成的光伏发电机组的输出功率P_out采用以下公式1计算：

$\left\{\begin{array}{c}{\mathrm{\ω}}_{\mathrm{ij}}(k+1)={\mathrm{\ω}}_{i(j-1)}\left(k\right)+P_{\mathrm{ij}},j=\mathrm{2,3},...,N_i\\ {\mathrm{\ω}}_{i1}(k+1)=P_{i1}\\ P_{\mathrm{out}}=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}\underset{j=1}{\overset{N_i}{\mathrm{\Σ}}}{P}_j=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{\ω}}_{{\mathrm{iN}}_i}\end{array}\right.---\left(1\right)$

其中，n为光伏发电机的总数量，n为正整数；n台光伏发电机分成m个小组，m表示小组的数量，且1≤m≤n，i为第i个小组，1≤i≤m；N_i表示第i个小组由N_i台光伏发电机组成；j表示每个小组中光伏发电机的数量，且j＝1,2,3,…,N_i；P_ij为第i个小组中第j台光伏发电机输出的有功功率的实测值；P_i1为第i个小组中第1台光伏发电机输出的有功功率的实测值；ω_ij为第i个小组中第j台光伏发电机有功功率的辅助变量；ω_i1为第i个小组中第1台光伏发电机有功功率的辅助变量；k为迭代运算的指数，且k＝0,1,2,…。

所述的步骤2)中的光伏发电机组的储备率的指令值β₀，采用以下公式2的高层控制器计算得到：

β₀＝∫K_P(P₀-P_out)dt    (2)

其中，K_P为高层控制器增益，P₀为光伏发电机组输出功率的指令值，t为积分计算的时间。

所述的步骤3)中的光伏发电机组中每台光伏发电机输出的最大功率P_imax，采用以下公式3进行近似计算：

$P_{i\max }=\left\{\begin{array}{c}0.{727U}_i,0<U_i\&le;22\\ 26U_i-\mathrm{556,22}<U_i\&le;25\\ {98U}_i-\mathrm{2356,25}<U_i\&le;26\\ {137U}_i-\mathrm{3370,26}<U_i\&le;28\\ 466,U_i>28\end{array}\right.---\left(3\right)$

其中，U_i为第i台光伏发电机的输出电压的实测值。

所述的步骤3)中的光伏发电机组中每台光伏发电机的储备率的实测值β_i，采用以下公式4计算得到：

${\mathrm{\&beta;}}_i=1-\frac{P_i}{P_{i\max }}---\left(4\right)$

其中，P_i为第i台光伏发电机的输出功率的实测值。

所述的步骤4)中的通信矩阵S采用以下公式5：

其中，i、j分别为通信矩阵S的行数和列数，s_ij(t)为通信矩阵S中的第i行、第j列元素；当i与j相同时，均有s_ij(t)＝1；当第j台光伏发电机的信息在t时刻可被第i台光伏发电机获知时，s_ij(t)＝1，否则s_ij(t)＝0；当第i台光伏发电机在t时刻可从高层控制中接受信息时，s_i0(t)＝1，否则s_i0(t)＝0。

所述的步骤4)中的光伏发电机组中每台光伏发电机的储备率的参考值采用以下公式6的底层控制器计算：

${\mathrm{\&beta;}}_i^{\mathrm{ref}}=K_0[-{\mathrm{\&beta;}}_i+D_{i0}{\mathrm{\&beta;}}_0+\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{D}_{\mathrm{ij}}{\mathrm{\&beta;}}_j]+{\mathrm{\&beta;}}_i---\left(6\right)$

其中，K₀为底层控制器增益，上式6中的通信矩阵S的算子D_ij采用以下公式7进行计算：

$D_{\mathrm{ij}}=s_{\mathrm{ij}}/\underset{j=0}{\overset{n}{\mathrm{\&Sigma;}}}{s}_{\mathrm{ij}},i=\mathrm{1,2},...,n---\left(7\right)$

其中，s_ij为通信矩阵S中的第i行、第j列元素。

所述的步骤4)中的每台光伏发电机的并网逆变器直流母线电压的参考值采用以下公式8计算：

$P_i=F(u_{\mathrm{dc}}^{\mathrm{ref}}/(1-{\mathrm{\&beta;}}_i^{\mathrm{ref}}))---\left(8\right)$

其中，F(·)为光伏发电机的最大功率曲线。

本发明具有的有益的效果是：

本发明实现地理分散、数量较多的光伏发电机输出功率快速收敛至指令值，提高光伏发电机组中熟练较多的光伏发电机的协同工作能力；无需主控中心，降低了通信成本，提高了光伏发电机组运行控制的鲁棒性；实现光伏发电机的“即插即用”功能，便于光伏发电机组模块化功能的实现；为电力网络提供辅助服务，如平滑负载、频率和电压的控制等。

附图说明

图1为本发明方法的控制结构图。

图2为本发明方法的控制框图。

图3为本发明方法的高层控制器框图。

图4为本发明方法的P_out计算过程结构图。

图5为实施例仿真例的结构连接示意图。

图6为实施例光伏发电机组的功率指令突变时的仿真图。

图7为实施例光伏发电机组的光照强度突变时的仿真图。

具体实施方式

下面结合附图及具体实施例对本发明作进一步详细说明。

如图1、图2所示，本发明包括以下步骤：

1)通过对各个光伏发电机输出的有功功率引入辅助变量进行迭代计算，得到由多台光伏发电机组成的光伏发电机组的输出功率；

2)通过基于积分反馈控制的高层控制器计算得到光伏发电机组的储备率的指令值；

3)对光伏发电机组中每台光伏发电机输出的最大功率进行近似计算，然后计算光伏发电机组中每台光伏发电机的储备率的实测值；

4)建立通信矩阵，通过底层控制器计算光伏发电机组中每台光伏发电机的储备率的参考值，然后计算每台光伏发电机的并网逆变器直流母线电压udc的参考值；

5)通过光伏并网逆变器电压电流双闭环PI控制方法，实现光伏发电机组中每台光伏发电机输出功率的稳定控制，进而实现无主控中心的光伏发电机组的协同控制。

如图4所示，步骤1)中的多台光伏发电机组成的光伏发电机组的输出功率P_out采用以下公式1计算：

$\left\{\begin{array}{c}{\mathrm{\&omega;}}_{\mathrm{ij}}(k+1)={\mathrm{\&omega;}}_{i(j-1)}\left(k\right)+P_{\mathrm{ij}},j=\mathrm{2,3},...,N_i\\ {\mathrm{\&omega;}}_{i1}(k+1)=P_{i1}\\ P_{\mathrm{out}}=\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}\underset{j=1}{\overset{N_i}{\mathrm{\&Sigma;}}}{P}_j=\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{\mathrm{\&omega;}}_{{\mathrm{iN}}_i}\end{array}\right.---\left(1\right)$

其中，n为光伏发电机的总数量，n为正整数；n台光伏发电机分成m个小组，m表示小组的数量，且1≤m≤n，i为第i个小组，1≤i≤m；N_i表示第i个小组由N_i台光伏发电机组成；j表示每个小组中光伏发电机的数量，且j＝1,2,3,…,N_i；P_ij为第i个小组中第j台光伏发电机输出的有功功率的实测值；P_i1为第i个小组中第1台光伏发电机输出的有功功率的实测值；ω_ij为第i个小组中第j台光伏发电机有功功率的辅助变量；ω_i1为第i个小组中第1台光伏发电机有功功率的辅助变量；k为迭代运算的指数，且k＝0,1,2,…。

步骤2)中的光伏发电机组的储备率的指令值β₀，采用以下公式2的基于积分反馈控制的高层控制器计算得到，如图3所示，图中s为拉普拉斯算子：

β₀＝∫K_P(P₀-P_out)dt    (2)

其中，K_P为高层控制器增益，P₀为光伏发电机组输出功率的指令值，t为积分计算的时间。

步骤3)中的光伏发电机组中每台光伏发电机输出的最大功率P_imax，采用以下公式3进行近似计算：

$P_{i\max }=\left\{\begin{array}{c}0.{727U}_i,0<U_i\&le;22\\ 26U_i-\mathrm{556,22}<U_i\&le;25\\ {98U}_i-\mathrm{2356,25}<U_i\&le;26\\ {137U}_i-\mathrm{3370,26}<U_i\&le;28\\ 466,U_i>28\end{array}\right.---\left(3\right)$

其中，U_i为第i台光伏发电机的输出电压的实测值。

步骤3)中的光伏发电机组中每台光伏发电机的储备率的实测值β_i，采用以下公式4计算得到：

${\mathrm{\&beta;}}_i=1-\frac{P_i}{P_{i\max }}---\left(4\right)$

其中，P_i为第i台光伏发电机的输出功率的实测值。

步骤4)中的通信矩阵S采用以下公式5：

其中，i、j分别为通信矩阵S的行数和列数，分别对应于小组数量和每个小组中光伏发电机的数量，s_ij(t)为通信矩阵S中的第i行、第j列元素；当i与j相同时，均有s_ij(t)＝1；当s_ij(t)表示的第j台光伏发电机的信息在t时刻可被第i台光伏发电机获知时，s_ij(t)＝1，否则s_ij(t)＝0；当s_i0(t)表示的第i台光伏发电机在t时刻可从高层控制中接受信息时，s_i0(t)＝1，否则s_i0(t)＝0。

步骤4)中的光伏发电机组中每台光伏发电机的储备率的参考值采用以下公式6的底层控制器计算：

${\mathrm{\&beta;}}_i^{\mathrm{ref}}=K_0[-{\mathrm{\&beta;}}_i+D_{i0}{\mathrm{\&beta;}}_0+\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{D}_{\mathrm{ij}}{\mathrm{\&beta;}}_j]+{\mathrm{\&beta;}}_i---\left(6\right)$

其中，K₀为底层控制器增益。

以上公式6中的通信矩阵S的算子D_ij采用以下公式7进行计算：

$D_{\mathrm{ij}}=s_{\mathrm{ij}}/\underset{j=0}{\overset{n}{\mathrm{\&Sigma;}}}{s}_{\mathrm{ij}},i=\mathrm{1,2},...,n---\left(7\right)$

其中，s_ij为通信矩阵S中的第i行、第j列元素。

步骤4)中的每台光伏发电机的并网逆变器直流母线电压的参考值采用以下公式8计算：

$P_i=F(u_{\mathrm{dc}}^{\mathrm{ref}}/(1-{\mathrm{\&beta;}}_i^{\mathrm{ref}}))---\left(8\right)$

其中，F(·)为光伏发电机的最大功率曲线。

步骤5)中的光伏并网逆变器电压电流双闭环PI控制方法，采用文献“张国月,曲轶龙,齐冬莲,等.基于重复控制的三电平光伏逆变技术[J].浙江大学学报(工学版),2012,46(7):1339—1344.”中的电压电流双闭环PI控制方法，对光伏发电机组中每台光伏发电机输出功率通过控制光伏并网逆变器来实现稳定控制。

本发明尤其适用于无主控中心的光伏(Photovoltaic，PV)发电机组协同控制方法。如图1所示，计算由多台光伏发电机组成的光伏发电机组的输出功率；设计高层控制器，计算光伏发电机组的储备率的指令值；估计光伏发电机组中每台光伏发电机输出的最大功率，并计算光伏发电机组中每台光伏发电机的储备率的实测值；设计通信矩阵和底层控制器，计算光伏发电机组中每台光伏发电机的储备率的参考值，并计算每台光伏发电机的并网逆变器直流母线电压的参考值；通过光伏并网逆变器电压电流双闭环PI控制方法，实现光伏发电机组中每台光伏发电机输出功率的稳定控制。

本发明的步骤1)通过辅助变量和迭代方法计算光伏发电机组的输出功率，可以消除分布式控制方法中公共检测点的限制，实现光伏发电机组的输出功率的准确计算。

如图2所示，本发明的步骤2)设计高层控制器计算光伏发电机组的储备率的指令值，可以实现光伏发电机组中各台光伏发电机在电力网络中的公平竞争。

本发明的步骤3)通过估计光伏发电机组中每台光伏发电机输出的最大功率来计算光伏发电机组中每台光伏发电机的储备率的实测值，不仅可以避免光伏电池板参数实时采集较为复杂的情况，而且可以较为准确的计算每台光伏发电机的储备率，为分布式控制方法的设计提供基础条件。

如图2所示，本发明的步骤4)通过设计通信矩阵和底层控制器，并计算光伏发电机组中每台光伏发电机的储备率的参考值，可以实现光伏发电机组中每台光伏发电机的储备率快速趋于同一数值，协调控制光伏发电机组输出的功率，并可以准确计算出每台光伏发电机的并网逆变器直流母线电压的参考值，为光伏并网逆变器的控制提供依据。

本发明的步骤5)采用工程中常用的光伏并网逆变器电压电流双闭环PI控制方法，可以实现光伏发电机组中每台光伏发电机输出功率的稳定控制。

在由数量较多的光伏发电机组成的光伏发电机组的平衡点处，每台光伏发电机的储备率如下：

1-P₁/P_1max＝...＝1-P_n/P_nmax＝β₀

其中，P_imax为第i台光伏发电机的最大功率。

在分析和计算每台光伏发电机的储备率的参考值时，将每台光伏发电机均视为一阶惯性环节，如下：

$T_i{\stackrel{\&CenterDot;}{\mathrm{\&beta;}}}_i={-\mathrm{\&beta;}}_i+{\mathrm{\&beta;}}_i^{\mathrm{ref}}$

其中，T_i为近似惯性时间常数，表示β_i的微分量。

以上两式为光伏发电机组中每台光伏发电机储备率及其参考值的分配和计算规则，该规则是通过本发明实现光伏发电机组输出功率协同功率控制的公知依据。

本发明的具体实施例如下：

在图5所示的仿真例图上对本发明提出的控制方法进行了仿真实验。图中的每个节点表示电力网络中的节点，例如图中数字1～114、149、150、195、250、251、300、35、450、451，PV1～PV10所在的节点处设有光伏发电机。

仿真试验参数如表1所示。

表1

在Matlab/Simulink中建立仿真模型，采用本发明提出的控制方法，所得仿真实验数据：

光伏发电机组功率指令突变时：每台光伏发电机的储备率β_i跟踪至指令值β₀的时间小于2.5s，光伏发电机组输出功率的计算值P_out跟踪至指令值P₀的时间小于1.6s，每台光伏发电机输出最大功率的估计值P_estimate跟踪至实测值P_m的时间小于2s；

光伏发电机组光照强度突变时：第1光伏电厂与第2光伏电厂的储备率恢复至相同值的时间小于0.9s，第1光伏电厂与第2光伏电厂的输出功率恢复稳定的时间小于2.3s，光伏发电机组输出功率的计算值P_out跟踪至指令值P₀的时间小于1.8s，每台光伏发电机输出最大功率的估计值P_estimate跟踪至实测值P_m的时间小于1.2s。

实验截图如下：

(1)图6显示了光伏发电机组的功率指令突变时，每台光伏发电机的储备率与指令值、光伏发电机组输出功率的计算值与指令值以及每台光伏发电机输出最大功率的估计值与实测值的仿真图。可见，以上三组测量参数都可以在较短的时间内趋于相同。仿真实验表明，本发明提出的控制方法可以提高光伏发电机组中每台光伏发电机运行的稳定性，实现光伏发电机输出功率的协同控制。

(2)图7显示了光伏发电机组的光照强度突变时，第1光伏电厂与第2光伏电厂的储备率、第1光伏电厂与第2光伏电厂的输出功率、光伏发电机组输出功率的计算值与指令值以及每台光伏发电机输出最大功率的估计值与实测值的仿真图。可见，以上四组测量参数都可以在较短的时间内趋于稳定。仿真实验表明，本发明提出的控制方法可以提高光伏发电机组中每台光伏发电机运行的鲁棒性，实现光伏发电机组的安全运行。

上述具体实施方式用来解释说明本发明，而不是对本发明进行限制，在本发明的精神和权利要求的保护范围内，对本发明做出的任何修改和改变，都落入本发明的保护范围。

Claims (8)

1.一种适用于无主控中心的光伏发电机组协同控制方法，其特征在于：包括以下步骤：

1)通过对各个光伏发电机输出的有功功率引入辅助变量进行迭代计算，得到由多台光伏发电机组成的光伏发电机组的输出功率；

2)通过高层控制器计算得到光伏发电机组的储备率的指令值；

3)对光伏发电机组中每台光伏发电机输出的最大功率进行近似计算，然后计算光伏发电机组中每台光伏发电机的储备率的实测值；

4)建立通信矩阵，通过底层控制器计算光伏发电机组中每台光伏发电机的储备率的参考值，然后计算每台光伏发电机的并网逆变器直流母线电压u_dc的参考值；

5)通过光伏并网逆变器电压电流双闭环PI控制方法，实现光伏发电机组中每台光伏发电机输出功率的稳定控制。

2.根据权利要求1所述的一种适用于无主控中心的光伏发电机组协同控制方法，其特征在于：所述的步骤1)中的多台光伏发电机组成的光伏发电机组的输出功率P_out采用以下公式1计算：

$\left\{\begin{array}{c}{\mathrm{\&omega;}}_{\mathrm{ij}}(k+1)={\mathrm{\&omega;}}_{i(j-1)}\left(k\right)+P_{\mathrm{ij}},j=\mathrm{2,3},...,N_i\\ {\mathrm{\&omega;}}_{i1}(k+1)=P_{i1}\\ P_{\mathrm{out}}=\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}\underset{j=1}{\overset{N_i}{\mathrm{\&Sigma;}}}{P}_j=\underset{i=1}{\overset{m}{\mathrm{\&Sigma;}}}{\mathrm{\&omega;}}_{{\mathrm{iN}}_i}\end{array}\right.---\left(1\right)$

其中，n为光伏发电机的总数量，n为正整数；n台光伏发电机分成m个小组，m表示小组的数量，且1≤m≤n，i为第i个小组，1≤i≤m；N_i表示第i个小组由N_i台光伏发电机组成；j表示每个小组中光伏发电机的数量，且j＝1,2,3,…,N_i；P_ij为第i个小组中第j台光伏发电机输出的有功功率的实测值；P_i1为第i个小组中第1台光伏发电机输出的有功功率的实测值；ω_ij为第i个小组中第j台光伏发电机有功功率的辅助变量；ω_i1为第i个小组中第1台光伏发电机有功功率的辅助变量；k为迭代运算的指数，且k＝0,1,2,…。

3.根据权利要求1所述的一种适用于无主控中心的光伏发电机组协同控制方法，其特征在于：所述的步骤2)中的光伏发电机组的储备率的指令值β₀，采用以下公式2的高层控制器计算得到：

β₀＝∫K_P(P₀-P_out)dt    (2)

其中，K_P为高层控制器增益，P₀为光伏发电机组输出功率的指令值，t为积分计算的时间。

4.根据权利要求1所述的一种适用于无主控中心的光伏发电机组协同控制方法，其特征在于：所述的步骤3)中的光伏发电机组中每台光伏发电机输出的最大功率P_imax，采用以下公式3进行近似计算：

$P_{i\max }=\left\{\begin{array}{c}0.{727U}_i,0<U_i\&le;22\\ 26U_i-\mathrm{556,22}<U_i\&le;25\\ {98U}_i-\mathrm{2356,25}<U_i\&le;26\\ {137U}_i-\mathrm{3370,26}<U_i\&le;28\\ 466,U_i>28\end{array}\right.---\left(3\right)$

其中，U_i为第i台光伏发电机的输出电压的实测值。

5.根据权利要求1所述的一种适用于无主控中心的光伏发电机组协同控制方法，其特征在于：所述的步骤3)中的光伏发电机组中每台光伏发电机的储备率的实测值β_i，采用以下公式4计算得到：

${\mathrm{\&beta;}}_i=1-\frac{P_i}{P_{i\max }}---\left(4\right)$

其中，P_i为第i台光伏发电机的输出功率的实测值。

6.根据权利要求1所述的一种适用于无主控中心的光伏发电机组协同控制方法，其特征在于：所述的步骤4)中的通信矩阵S采用以下公式5：

其中，i、j分别为通信矩阵S的行数和列数，s_ij(t)为通信矩阵S中的第i行、第j列元素；当i与j相同时，均有s_ij(t)＝1；当第j台光伏发电机的信息在t时刻可被第i台光伏发电机获知时，s_ij(t)＝1，否则s_ij(t)＝0；当第i台光伏发电机在t时刻可从高层控制中接受信息时，s_i0(t)＝1，否则s_i0(t)＝0。

7.根据权利要求1所述的一种适用于无主控中心的光伏发电机组协同控制方法，其特征在于：所述的步骤4)中的光伏发电机组中每台光伏发电机的储备率的参考值采用以下公式6的底层控制器计算：

${\mathrm{\&beta;}}_i^{\mathrm{ref}}=K_0[-{\mathrm{\&beta;}}_i+D_{i0}{\mathrm{\&beta;}}_0+\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{D}_{\mathrm{ij}}{\mathrm{\&beta;}}_j]+{\mathrm{\&beta;}}_i---\left(6\right)$

其中，K₀为底层控制器增益，上式6中的通信矩阵S的算子D_ij采用以下公式7进行计算：

$D_{\mathrm{ij}}=s_{\mathrm{ij}}/\underset{j=0}{\overset{n}{\mathrm{\&Sigma;}}}{s}_{\mathrm{ij}},i=\mathrm{1,2},...,n---\left(7\right)$

其中，s_ij为通信矩阵S中的第i行、第j列元素。

8.根据权利要求1所述的一种适用于无主控中心的光伏发电机组协同控制方法，其特征在于：所述的步骤4)中的每台光伏发电机的并网逆变器直流母线电压的参考值采用以下公式8计算：

$P_i=F(u_{\mathrm{dc}}^{\mathrm{ref}}/(1-{\mathrm{\&beta;}}_i^{\mathrm{ref}}))---\left(8\right)$

其中，F(·)为光伏发电机的最大功率曲线。