[发明专利]一种基于直角坐标解法的极坐标牛顿法潮流算法

权利要求书

1.一种基于直角坐标解法的极坐标牛顿法潮流算法，其特征包括以下步骤：

步骤1：打开数据文件，读取Y阵数据文件到Y(n,2n)数组；

步骤2：根据Y(n,2n)数组，用对极坐标牛顿法的计算式进行数学变换得到基于直角坐标解法的计算式计算ΔP_i、ΔQ_i和J阵元素；

数学变换后的基于直角坐标解法的计算式如下：

ΔPi=Pi-Σj=1j=n(Gijeiej+Gijfifj+Bijejfi-Bijeifj)]]>

ΔQi=Qi-Σj=1j=n(Gijejfi-Gijeifj-Bijeiej-Bijfifj)]]>

H_ij＝-G_ije_jf_i+G_ije_if_j+B_ije_ie_j+B_ijf_if_j

N_ij＝-G_ije_ie_j-G_ijf_if_j-B_ije_jf_i+B_ije_if_j

M_ij＝G_ije_ie_j+G_ijf_if_j+B_ije_jf_i-B_ije_if_j＝-N_ij

L_ij＝-G_ije_jf_i+G_ije_if_j+B_ije_ie_j+B_ijf_if_j＝H_ij

Hii=Σj=1j=n(Gijejfi-Gijeifj-Bijeiej-Bijfifj)+Biiei2+Biifi2]]>

Nii=-Σj=1j=n(Gijeiej+Gijfifj+Bijejfi-Bijeifj)+Gijeiej+Gijfifj-2(ei2+fi2)Gii]]>

Mii=-Σj=1j=n(Gijeiej+Gijfifj+Bijejfi-Bijeifj)+Gijeiej+Gijfifj=Nii+2(ei2+fi2)Gii]]>

Lii=-Σj=1j=n(Gijejfi-Gijeifj-Bijeiej-Bijfifj)-Biiei2-Biifi2+2Vi2Bii=-Hii+2(ei2+fi2)Bii]]>

步骤3：对J阵进行消元和回代求取ΔV_i、Δδ_i；

通过ΔV_i、Δδ_i求出电压幅值和相角的值再通过极坐标与直角坐标转换得到电压的实部和虚部然后计算ΔP_i、ΔQ_i；

步骤4：判断是否满足收敛条件；

如果不满足收敛条件，则跳转到步骤2；如果满足收敛条件，则执行步骤5；

步骤5：结束迭代并输出结果。