1.摘要
光伏模型未知参数提取的准确性直接关系到光伏系统建模、仿真和控制的有效性。本文提出了一种新的元启发式算法,称为广义正态分布优化算法,该方法受到广义正态分布模型的启发;每个个体使用广义正态分布曲线来更新其位置。与大多数元启发式算法不同,所提出的方法只需要基本的种群大小和终止条件就可以解决优化问题。
2.算法原理
GNDO的灵感来自于正态分布理论,其概率密度函数可表示为:
f
(
x
)
=
1
2
π
δ
e
x
p
(
−
(
x
−
μ
)
2
2
δ
2
)
f(x)=\frac{1}{\sqrt{2\pi}\delta}exp\Bigg(-\frac{(x-\mu)^2}{2\delta^2}\Bigg)
f(x)=2πδ1exp(−2δ2(x−μ)2)
GNDO中设计的信息共享策略包括局部开发和全局探索。局部开发基于建立的广义正态分布模型,该模型以当前平均位置和当前最优位置为指导,全局探索与三个随机选择的个体有关。
局部开发
局部开发是指在由所有个体当前位置组成的搜索空间周围寻找更好解决方案的过程。基于种群中个体的分布与正态分布之间的关系,可以通过以下方式构建一个用于优化的广义正态分布模型:
v
i
t
=
μ
i
+
δ
i
×
η
,
i
=
1
,
2
,
3
,
.
.
.
,
N
v_i^t=\mu_i + \delta_i\times\eta,i=1,2,3,...,N
vit=μi+δi×η,i=1,2,3,...,N
其中vi为第i个个体在t时刻的轨迹向量,μi为第i个个体的广义平均位置,δi为广义标准差,η为惩罚因子。μi、δi和η可以定义为:
μ
i
=
1
3
(
x
i
t
+
x
B
e
s
t
t
+
M
)
δ
i
=
1
3
[
(
x
i
t
−
μ
)
2
+
(
x
B
e
s
t
t
−
μ
)
2
+
(
M
−
μ
)
2
]
η
=
{
−
l
o
g
(
λ
1
)
×
c
o
s
(
2
π
λ
2
)
,
i
f
a
<
=
b
−
l
o
g
(
λ
1
)
×
c
o
s
(
2
π
λ
2
+
π
)
,
o
t
h
e
r
w
i
s
e
\left.\begin{aligned}&\mu_{i} = \frac{1}{3}\left(\boldsymbol{x}_{i}^{t}+\boldsymbol{x}_{\mathrm{Best}}^{t}+\boldsymbol{M}\right)\\&\delta_{i} = \sqrt{\frac{1}{3}\left[\left(\boldsymbol{x}_{i}^{t}-\mu\right)^{2}+\left(\boldsymbol{x}_{\mathrm{Best}}^{t}-\mu\right)^{2}+\left(\boldsymbol{M}-\mu\right)^{2}\right]}\\&\eta = \left\{\begin{array}{ll}\sqrt{-log(\lambda_1)}&\times cos(2\pi\lambda_2),\mathrm{if~}a<=b\\\sqrt{-log(\lambda_1)}&\times cos(2\pi\lambda_2+\pi),\mathrm{otherwise}\end{array}\right.\end{aligned}\right.
μi=31(xit+xBestt+M)δi=31[(xit−μ)2+(xBestt−μ)2+(M−μ)2]η={−log(λ1)−log(λ1)×cos(2πλ2),if a<=b×cos(2πλ2+π),otherwise
全局探索
全局探索是指全面搜索一个空间以找到有前景的区域,GNDO中的全局探索是基于三个随机选择的个体,可以表达为:
v
i
′
=
x
i
′
+
β
×
(
∣
λ
3
∣
×
v
1
)
⏟
Local information sharing
+
(
1
−
β
)
×
(
∣
λ
4
∣
×
v
2
)
⏟
Global information sharing
\begin{array}{cccc}\boldsymbol{v}_i'=\boldsymbol{x}_i'&+&\underbrace{\beta\times(|\lambda_3|\times\boldsymbol{v}_1 )}_{\text{Local information sharing}}&+\underbrace{(1-\beta)\times(|\lambda_4|\times\boldsymbol{v}_2 )}_{\text{Global information sharing}}\end{array}
vi′=xi′+Local information sharing
β×(∣λ3∣×v1)+Global information sharing
(1−β)×(∣λ4∣×v2)
v1和v2是两个跟踪向量:
v
1
=
{
x
i
t
−
x
p
1
t
,
if
f
(
x
i
t
)
<
f
(
x
p
1
t
)
x
p
1
t
−
x
i
t
,
otherwise
\boldsymbol{v}_{1} = \left\{\begin{matrix}\boldsymbol{x}_{i}^{t}-\boldsymbol{x}_{\text{p}1}^{t},\text{if}f(\boldsymbol{x}_{i}^{t})<f(\boldsymbol{x}_{\text{p}1}^{t})\\\boldsymbol{x}_{\text{p}1}^{t}-\boldsymbol{x}_{i}^{t},\text{otherwise}\end{matrix}\right.
v1={xit−xp1t,iff(xit)<f(xp1t)xp1t−xit,otherwise
v
2
=
{
x
p
2
t
−
x
p
3
t
,
if
f
(
x
p
2
t
)
<
f
(
x
p
3
t
)
x
p
3
t
−
x
p
2
t
,
otherwise
\boldsymbol{v}_{2} = \begin{cases}\boldsymbol{x}_{\text{p}2}^{t}-\boldsymbol{x}_{\text{p}3}^{t},\text{if} f(\boldsymbol{x}_{\text{p}2}^{t})<f(\boldsymbol{x}_{\text{p}3}^{t})\\\boldsymbol{x}_{\text{p}3}^{t}-\boldsymbol{x}_{\text{p}2}^{t},\text{otherwise}\end{cases}
v2={xp2t−xp3t,iff(xp2t)<f(xp3t)xp3t−xp2t,otherwise
伪代码
3.结果展示
4.参考文献
[1] Zhang Y, Jin Z, Mirjalili S. Generalized normal distribution optimization and its applications in parameter extraction of photovoltaic models[J]. Energy Conversion and Management, 2020, 224: 113301.