采用乘子罚函数的牛顿法最优无功潮流(二)

<正> 四、采用乘子罚函数法的牛顿型无功潮流算法 1 乘子罚函数法原理 考虑等式约束问题 Minf(x) S.T.gi(x)=0 i=1,2,……m 设其极小点x~*,考虑增广拉格朗日函数:M(x,u)=f(x)-sum from i=1 to m(u_ig_i(x))+C/2 sum from i=1 to m[g_i(x)]~2 存在u~*使(x~*,u~*)为L(x,u)的稳定点,即▽_xL(x~*,u~*)=▽f(x~*)-sum from i=1 to m(u_i~*▽g_i(x~*)=0 而附加项1/2 sum from i=1 to m[g_i(x)]~2在x~*处的梯度为零,因此,▽_xM(x~*,u~*)