电力系统分析潮流计算问题 内附详细题目 10
这道题应该如何处理应该先将环网拆开算出拆开后算两侧的自然功率分布在从功率分点开拆算电压降落吗？求指导。。。。。这个并联支路真心不知道怎么处理了。。。这是电力系统稳态的题。...
这道题应该如何处理 应该先将环网拆开 算出 拆开后算两侧的自然功率分布 在从功率分点开拆算电压降落吗？ 求指导。。。。。这个并联支路真心不知道怎么处理了。。。
这是电力系统稳态的题。。。。请用手算法 分析一下 感谢 展开
这是电力系统稳态的题。。。。请用手算法 分析一下 感谢 展开
展开全部
＃包括“stdio.h中”
的＃include“math.h的”
＃定义M 20 / *节点数目，分行数目限制* /
＃定义N 5 / *迭代次数的限制* /
整数N，M，DD = 0，PQ = 0，PV = 0，TT，QY;
浮动EPS / *节点数目，地支的，PQ节点，PV节点，控制打印，互联网一些精密* /
结构jiedian / *节点的数据结构* /
{int类型; / *节点的分支数类型:1-PQ节点; 2-PV节点; 3 - 平衡节点* /
持股量P，Q，E，F，V; / *节点主动和无功功率，电压，实部，虚部发言电压值* /
} jiedian [M] / *作为一个参考节点的无功功率，相比jiedian [1]。 Q * /
结构植绿/ *分支的数据结构* / {P1，P2，S / *分支结束节点号，分支类型：1 - 正常滑中国 - 变压器分支3 - 地支鲁* /
浮动R，X，B，克拉; / *防滑性，抗，导纳，非标准变压器匝数比* / }植绿[M]; /> 结构hulianwang / *网络结构* /
{数，PV / * Internet区域的号码，指定PV节点数量* /
持股量P，EPS2 / *指定激活时，其允许误差* /
诠释计数; / *每个区域包括相互连接的节点* /
INT A [M] / *包括节点* /
} hulianwang [M];
静态浮动G [M] [M]，B [M] [M]，G1 [M] [M]，B1 [M] [M] / *节点导纳矩阵矩阵* /
静浮YKB [2 * M] [2 * M] / *节点导纳矩阵矩阵* /
静态浮动饮子[2 * M] [2 * M];系数表* / * /
静态持股量P [M] [M]，Q [M] [M] / *流* /
FILE * FP1，FP2 * / *文件的输入和输出指针* / / a>
无效输入（）/ *输入文件in.txt基本参数，节点数据，分支机构数据线* /
{
INT I，J，H / *循环变量节点（分公司）类型* /
/ **********************************开放电网参数文件********************** /
FP1 =的FOPEN（英寸TXT“，”R“）;
如果（FP1 == NULL）
{
的printf（“无法打开文件in.txt！\ N”）;
出口（0）;
}
/ ***********阅读的节点，多家分支机构，互联网号码，PQ节点，节点和准确性的光伏*********** /
> fscanf（FP1，“％D，％D，％D，％F \ N”，＆N，＆M＆QY＆EPS）;
/ ******* *****************读取节点信息****************************** ******* /
（1 <= N;我+ +）
{
fscanf（FP1，“％D”，＆H）; BR />（H == 1）/ *类型H = 1的PQ节点* /
{
PQ + +;
fscanf（FP1，“％F％F ％F，％F \ N“，＆jiedian的[I] P，＆jiedian的[I] Q，＆jiedian的[I]。＆jiedian的[I] F）;
jiedian [ ]。 S = 1;，
jiedian [I]。 = SQRT（jiedian [I]，E * jiedian [I] E + jiedian [I]，F * jiedian [I] F）;
}
其他（H == 2） / * H型= 2是PV节点* /
{
PV + +;
fscanf（FP1，“％F，％F，％F，％F \”＆ jiedian [] P，＆jiedian的[I] V，＆jiedian的[I]，E，及jiedian [I]，F）;
jiedian [I]。 S = 2;，
jiedian [I]。 Q = 0;
}
其他/ *型h = 3是平衡的节点* /
{ jiedian [I]。 P = 0;，
jiedian [I]。 Q = 0;，
jiedian [I]。 E = 1;，
jiedian [I]。 F = 0;，
jiedian [I]。 V = 1;，
jiedian [I]。 S = 3;
}
}
/ *************************** ****阅读分支信息************************************* /
（I = 1; <= M; + +）
{
fscanf（FP1，“为％d”，＆H）;
（H == 1）/ *类型H = 1是一个普通的分支* /
{
fscanf（FP1，“％D％D％F％F％F \ N”，植绿[I] P1，植绿[I] P2，＆植绿[I]，R＆植绿[I]，＆植绿[I]）;
植绿[I]。 KT = 1;，
植绿[I]。 = 1;
}
（H == 2）/ * H型= 2是变压器支路* /
{
fscanf（FP1，“％D， ％％F，％F，％D，F \ N“，植绿[I]。P1，植绿[I]。P2，植绿[I]，R＆植绿[I]，X，＆植绿I]。克拉）;，
植绿[I]。 S = 2;
}
（H == 3）/ * H型= 3接地分支* /
{
fscanf（FP1，“％D，％ D，％F，％F，％F \ N“，植绿[I]。P1，植绿[I]。P2，植绿[I]，R＆植绿[I]，X，＆植绿[ ，B）;
植绿[I]。 KT = 1;，
植绿[I]。 S = 3;
DD + +;
}
}
/ ******************* ************读取互联网信息*********************************** / **
（qy! = 0）
（I = 1; <= QY; + +）/ *输入互联网状态* /
{
fscanf（FP1，“％D，％D，％F，％F，为％d”，＆hulianwang [i]于民，＆hulianwang的[I]光伏hulianwang [I]，P＆hulianwang [I] 。EPS2，hulianwang的[I]计数）;
（J = 1; <= hulianwang [I]。计数; J + +）
fscanf（FP1，“％D”， （hulianwang [I]一个[J]））;
}
FCLOSE（FP1）;
/ ********* *打开输出文件********************************************* *********** /
FP2 =打开（“out.txt的”，“W”）;
（FP2 == NULL）
{
>的printf（“无法打开文件！\ N”）;
出口（0）;
}
fprintf（FP2，“\ N ****原始数据如下************************************************ *********** \ n \ n已“）;
fprintf（FP2，”节点：％2D科数：％2D正确的地支的：％2D PQ节点编号：％2D PV节点：％2D精度1％F \ N“，N，M，DD，PQ，PV，EPS）;
fprintf（FP2，”\ N ---------- -------------------------------------------------- ----------------- \ N \ N“）;
（1 <= PQ + +）
fprintf（FP2 ，“％2D PQ节点节点P [％d〕=％F Q [％] =％F \”，我，我，jiedian [I]，P，I，jiedian []。Q）; p>就（i = PQ +1 <= PQ + PV + +）
fprintf（FP2，“节点％2D PV节点P [％] =％F V [％ e] =％f \ N“，我，我，jiedian [I]，P，I，jiedian [I]，V）;
fprintf（FP2，”节点％2D平衡节点\ N“，我）
fprintf（FP2，“，\ N -------------------------------------- ----------------------------------------- \ N \ N“）; />（I = 1; <= M; + +）
{
如果（植绿[I]。== 1）
fprintf（FP2，“支部％2D ：普通分支相关节点：％2D，2D R％=％F X =％F B =％F \“，我，植绿[I] P1，植绿[I]。P2，植绿[I]，R植绿[I]，植绿[I]）;
其他（植绿[I]。== 2）
fprintf（FP2，“分支％2D：变压器分会相关节点： ％2D，2D R％=％F X =％f KT =％F \“，我植绿[I] P1，植绿[I]。P2，植绿[I]，R，植绿[I]，X植绿[I]。克拉）;
其他
fprintf（FP2，“支部％2D：相关节点树枝在地上：％2D=％F X =％F B =％F \ “，我植绿[I]。P1，植绿[I]，R，植绿[I]，植绿[] b）条;
}
（I = 1 <= QY; + +）
{
fprintf（FP2，“Internet区域中：％2D指定PV节点：％2D提供有效：％8.5F允许误差：包括节点％8.5F ：“，hulianwang [I]数，hulianwang [I]。光伏，hulianwang [I]，P，hulianwang []。EPS2）;
（J = 1; <= hulianwang [I]。计数; J + +）
fprintf（FP2，“2D”，hulianwang [I]一个[J]）; fprintf（FP2，“\ N”）;
} <br的的 BR />}
无效油花（）/ *使用节点数据，分支机构的数据，形成一个新的号码* / {
INT I，K，J，JD [N / * JD [N] [记录]节点* / struct jiedian的温度; / *中间临时党支部连接的节点数量* /
（I = 1; <= N + +）/ *分支的节点上连接数赋0 * /
{
JD [0] = 0; }
（我= 1 <=我+ +）/ *考虑每个节点* /
{
（J = 1，J <= M，J + +）/ *考虑每个分支* /
{
（（植绿[J]。P1 = | |植绿[J]。P2 ==）&&（植绿[J]。！= 3））/ *如果这个分支的节点，不接地分支* /
JD [0] + +;
}
}
（I = 1 <PQ; + +）/ *正确的PQ节点根据分支机构连接到重新编号的数量* /
{
（J = i +1的J <= PQ研究+ +）/ *升序排序算法* /
{
（JD [I]> JD [J]）
{
TEM = jiedian [I];
jiedian [I] = jiedian [J];
jiedian [J] = TEM
（K = 1; K <= M，K + +）/ *更新分支信息* /
{
（植绿[K]。P1 == J）
植绿[K]。 P1 = I 否则，如果（植绿[K]。P2 == J）
植绿[K]。 P2 = I;
其他（植绿[K]。P1 ==我）
植绿[K]。 P1 = J;
否则，如果（植绿[K]。P2 ==）
植绿[K]。 P2 = J;
}
}
}
}
就（i = PQ +1; <PQ + PV; + +）/ * PV节点连接上重新编排的分支机构数量* /
{
（J = i +1的; J <= PQ + PV; J + +）/ *小大排序算法* /
{
如果（JD [我]> JD [J]）
{
温度= jiedian [我];
jiedian [I] = jiedian [J];
jiedian [J] = TEM
（K = 1; K <= M，K + +）/ *更新分支信息* /
> {
（植绿[K]。P1 == J）
植绿[K]。 P1 = I;
否则，如果（植绿[K]。P2 == J）
植绿[K]。 P2 = I;
其他（植绿[K]。P1 ==我）
植绿[K]。 P1 = J;
否则，如果（植绿[K]。P2 ==）
植绿[K]。 P2 = J;
}
}
}
}
fprintf（FP2，“\ N ****节点优化结果如下*********************************************** ********** \ n \ n已“）;
（I = 1; <=我+ +）
{
（jiedian [I]。== 1）
fprintf（FP2，“％2D PQ节点，节点P [％d〕=％F Q [％] =％F \”，我，我，jiedian [I]，P，I，jiedian [I]。Q）;
如果（jiedian [I]。== 2）
fprintf（FP2，“节点％2D PV节点P [％d的=％F V [％d〕=％F \“，我，我，jiedian [I]，P，我jiedian [I]，V）;
（jiedian [I]。= = 3）
fprintf（FP2，“％2D平衡节点的节点\ N”）;
}
}
的无效form_y（）/ *使用分支数据形成Y，要注意地面腿* /
{
I，J，K;
float类型的
（I = 1 <= N;我+ +）
（J = 1; <= N，J + +）
G [我] [J] = B [I] [J] = 0;
（1 <= M; + +）/ *节点导纳矩阵主对角线准入* /
（J = 1; <= N; J + +）
（植绿[I] P1 ==| |植绿[I]。P2 == J）
{
S =植绿[I]。 R *植绿[I]。 R +植绿[I]。 X *植绿[I]。 X;
（S == 0）;
G [J] [J] + =植绿[I]。 R / S;
B [J] [J] + =植绿[我]。 X / S;
（植绿[I]。== 1）/ *如果是普通的分支* / B [J] [J] + =植绿[I]。 B / 2;
}
其他（植绿[I]。== 2）/ *如果是普通变压器支路* /
{
（植绿[ P1 == J）
{
G [J] [J] + =（植绿[I]，R / S *（1-植绿[]克拉））/（植绿我克拉*植绿[I]克拉）;
B [J] [J] + =（-植绿[I] X / S *（1-植绿[]。克拉））/（植绿[I] KT *植绿[I]克拉）;
}
其他（植绿[I] P2 == j）条
{
[J] [J] + =（植绿[I]，R / S（植绿[I]千吨-1））/植绿[I]。克拉;
B [J] [J] + =（-植绿[I]，X / S *（植绿[I] KT-1））/植绿[I]。克拉;
}
}
其他（植绿[I]。== 3）/ *如果是在地面上分支* /
{
乙[J] [J] + =植绿[I]。 b;
}
}
（K = 1; K <= M，K + +）/ *节点导纳矩阵非准入的主对角线上* /
{
=植绿[K]。 P1;
=植绿[K]。 P2;
S =植绿[K]。 R *植绿[K]。 R +植绿[K]。 X *植绿[K]。 X;
（S == 0）;
G [我] [J] + =植绿[K]。 R / S;
B [我] [J] + =植绿[K]。 X / S;
（植绿[K]。克拉= 1.0）
{
G [我] [J] / =植绿[K]。克拉;
B [我] [J] / =植绿[K]。克拉;
}
G [J] [I] = G [I] [J];
B [J] [] = B [I] [J];
>}
fprintf（FP2，“\ N \ N ******************************节点导纳矩阵********************************** \ N“）;
（I = 1 <= N;我+ +）
{
fprintf（FP2，“\”）;
（J = 1; <= N研究+） BR /> fprintf（FP2，“％+ J％8.5F 8.5F，G [I] [J]，B [I] [J]）;
}
}
BR />无效form_j，（）/ *使用节点数据和Y形式的J * /
{
浮法EI，FI，A = 0，B = 0;
I1，J1，K1，I，J，K;
（I = 1; <= 2 *（PQ + PV）+1; + +）
（J = 1; J <= 2 *（PQ + PV）+1; J + +）
YKB [I] [J] = 0;
（I = 1; < = PQ + +）
（J = 1; J <; J + +）
{I1 =;
J1 = J;
EI = jiedian I]。 E;
Fi = jiedian的[I]。 F;
（i! = J）/ *找我！ = j的H，N，J，L * /
{YKB *-1 [2] [2 * J-1] = G [期1] [J1] * EI + B [期1] [J1] * FI / * H * /
YKB [2 *的i-1] [2 * J] = B [I1] [J1] * EI + G [I1] [J1] * FI / * N * /
YKB [2] [2 * J-1] = B [I1]？[J1] * EI + G [I1] [J1] * FI / * J * / YKB [2] [2 * J] =-G [I1] [J1] * EI-B [I1] [J1] * FI / * L * /
}
BR /> / *其他寻求当我= J H，N，J，K * /
{A = 0，B = 0;
（K = 1; K <= N，K + +）
（k! = I）
{K1 = K;
= A + G [I1] [K1] * jiedian [K]。 EB [I1] [K1] * jiedian [K]。 F;
B = B + G [I1] [K1] * jiedian [K]。 F + B [期1] [K1] * jiedian [●]。 E;
}
YKB [2 *的i-1] [2 * J-1] = 2 * G [I1] [I1] * jiedian [I]。 E + A / * H * /
YKB [2] [2 * J] = -2 * B [I1] [I1] * jiedian [I]。 F +一个; / * L * /
YKB [2 *的i-1] [2 * j]的= 2 * G [期1] [1期1] * jiedian []。 F + B; / * N * /
YKB [2] [2 * J-1] = -2 * B [I1] [I1] * jiedian [I]。 EB / * J * /
} }
就（i = PQ +1; <= PQ + PV; + +）/ *子矩阵形成光伏节点* /
（十= 1; J <; J + +）
{I1 =;
J1 = J;
EI = jiedian [我]。 E;
Fi = jiedian的[I]。 F;
（i! = J）/ *找我！ =当H，N * /
{
YKB [2 *的i-1] [2 * J-1] = G [期1] [J1] * EI + B [期1] [J1] * FI / * H * /
YKB [2 *的i-1] [2 * J] = B [I1] [J1] * EI + G [I1] [J1] * FI / * N * /
}
其他/ *求时，我= J H，N，R，S * /
{A = 0，B = 0;
（K = 1，K <= N，K +）
（k! = I）
{K1 = K;
+ = G [I1] [K1] * jiedian [K]。 EB [I1] [K1] * jiedian [K]。 F;
B + = G [I1] [K1] * jiedian [K]。 F + B [期1] [K1] * jiedian [●]。 E;
}
YKB [2 *的i-1] [2 * J-1] = 2 * G [I1] [I1] * jiedian [I]。 E + A / * H * /
YKB [2 *的i-1] [2] = 2 * G [I1] [I1] * jiedian [I]。 F + B; / * N * /
YKB [2] [2 * J-1] = 2 * EI / * R * /
YKB [2] [2 * J ] = 2 * FI / * S * /
}
}
fprintf（FP2，“\ N \ N \ N ********** *********************雅可比************** *********** \ N“）;
（I = 1; <= 2 *（PQ + PV）; + +） {
fprintf（FP2，“\”）;
（J = 1; <= 2 *（PQ + PV）; J + +）
fprintf（FP2，“％8.5 F“，YKB [I] [J]）;
}
}
无效form_yz（）/ *使用的J的外形表* / { BR /> INT I，J，K，
浮法A [2 * M] [2 * M]，L [2 * M] [2 * M]，R [2 * M] [2 * M]。
浮法X，Y，Z;
（I = 1; <= 2 *（PQ + PV）; + +）
（J = 1;研究<= 2 *（PQ + PV）; J + +）
{[I] [J] = 0，R [I] [J] = 0，L [I] [J] = 0 ;}
（I = 1; <= 2 *（PQ + PV）; + +）
（J = 1; <= 2 *（PQ + PV）; J + +） A [I] [J] = YKB [I] [J];
（I = 1; <= 2 *（PQ + PV）+ + ）
{
（J = 1; <= 2 *（PQ + PV），J + +）
{
如果（> J）
{
R [我] [J] = 0;
X = 0;
（K = 1; K <= J-1，K +）
{ BR /> X = X + L [I] [K] * R [K] [J];
}
L [我] [J] = A [I] [J] X;
}
（I == J）
{
R [我] [J] = 1;
Y = 0;
（K = 1; K <= I-1，K +）
{
Y = Y + L [我] [K] * R [K] [I];
} L [I] [I] = [I] [J];
}
（I <j）的
{
L [我] [J ] = 0;
Z = 0;
（K = 1; K <= I-1，K +）
Z = Z + L [我] [K] * R [K] [J]; }
R [我] [J] =（[I] [J] Z）/ L [] []; BR />}
}
}
（I = 1; <= 2 *（PQ + PV）; + +）
{（J = 1; <= 2 *（PQ + PV）; J + +）
{
（I <J）
饮子[I] [J] = R [我] [J];
（I == J）
饮子[I] [J] = 1 / L [I] [J];
（> j的）
饮子[我] [J] = L [I] [J];
}
}
fprintf（FP2，“\ N \ N \ N ********************************因子表**************** ************************* \ N“）;
（I = 1; <= 2 *（ PQ + PV）; + +）{
fprintf（FP2，“\”）;
（J = 1; <= 2 *（PQ + PV）; + +）
fprintf（FP2，“％8.5F”，银子岩[I] [J]）;
}
}
的＃include“math.h的”
＃定义M 20 / *节点数目，分行数目限制* /
＃定义N 5 / *迭代次数的限制* /
整数N，M，DD = 0，PQ = 0，PV = 0，TT，QY;
浮动EPS / *节点数目，地支的，PQ节点，PV节点，控制打印，互联网一些精密* /
结构jiedian / *节点的数据结构* /
{int类型; / *节点的分支数类型:1-PQ节点; 2-PV节点; 3 - 平衡节点* /
持股量P，Q，E，F，V; / *节点主动和无功功率，电压，实部，虚部发言电压值* /
} jiedian [M] / *作为一个参考节点的无功功率，相比jiedian [1]。 Q * /
结构植绿/ *分支的数据结构* / {P1，P2，S / *分支结束节点号，分支类型：1 - 正常滑中国 - 变压器分支3 - 地支鲁* /
浮动R，X，B，克拉; / *防滑性，抗，导纳，非标准变压器匝数比* / }植绿[M]; /> 结构hulianwang / *网络结构* /
{数，PV / * Internet区域的号码，指定PV节点数量* /
持股量P，EPS2 / *指定激活时，其允许误差* /
诠释计数; / *每个区域包括相互连接的节点* /
INT A [M] / *包括节点* /
} hulianwang [M];
静态浮动G [M] [M]，B [M] [M]，G1 [M] [M]，B1 [M] [M] / *节点导纳矩阵矩阵* /
静浮YKB [2 * M] [2 * M] / *节点导纳矩阵矩阵* /
静态浮动饮子[2 * M] [2 * M];系数表* / * /
静态持股量P [M] [M]，Q [M] [M] / *流* /
FILE * FP1，FP2 * / *文件的输入和输出指针* / / a>
无效输入（）/ *输入文件in.txt基本参数，节点数据，分支机构数据线* /
{
INT I，J，H / *循环变量节点（分公司）类型* /
/ **********************************开放电网参数文件********************** /
FP1 =的FOPEN（英寸TXT“，”R“）;
如果（FP1 == NULL）
{
的printf（“无法打开文件in.txt！\ N”）;
出口（0）;
}
/ ***********阅读的节点，多家分支机构，互联网号码，PQ节点，节点和准确性的光伏*********** /
> fscanf（FP1，“％D，％D，％D，％F \ N”，＆N，＆M＆QY＆EPS）;
/ ******* *****************读取节点信息****************************** ******* /
（1 <= N;我+ +）
{
fscanf（FP1，“％D”，＆H）; BR />（H == 1）/ *类型H = 1的PQ节点* /
{
PQ + +;
fscanf（FP1，“％F％F ％F，％F \ N“，＆jiedian的[I] P，＆jiedian的[I] Q，＆jiedian的[I]。＆jiedian的[I] F）;
jiedian [ ]。 S = 1;，
jiedian [I]。 = SQRT（jiedian [I]，E * jiedian [I] E + jiedian [I]，F * jiedian [I] F）;
}
其他（H == 2） / * H型= 2是PV节点* /
{
PV + +;
fscanf（FP1，“％F，％F，％F，％F \”＆ jiedian [] P，＆jiedian的[I] V，＆jiedian的[I]，E，及jiedian [I]，F）;
jiedian [I]。 S = 2;，
jiedian [I]。 Q = 0;
}
其他/ *型h = 3是平衡的节点* /
{ jiedian [I]。 P = 0;，
jiedian [I]。 Q = 0;，
jiedian [I]。 E = 1;，
jiedian [I]。 F = 0;，
jiedian [I]。 V = 1;，
jiedian [I]。 S = 3;
}
}
/ *************************** ****阅读分支信息************************************* /
（I = 1; <= M; + +）
{
fscanf（FP1，“为％d”，＆H）;
（H == 1）/ *类型H = 1是一个普通的分支* /
{
fscanf（FP1，“％D％D％F％F％F \ N”，植绿[I] P1，植绿[I] P2，＆植绿[I]，R＆植绿[I]，＆植绿[I]）;
植绿[I]。 KT = 1;，
植绿[I]。 = 1;
}
（H == 2）/ * H型= 2是变压器支路* /
{
fscanf（FP1，“％D， ％％F，％F，％D，F \ N“，植绿[I]。P1，植绿[I]。P2，植绿[I]，R＆植绿[I]，X，＆植绿I]。克拉）;，
植绿[I]。 S = 2;
}
（H == 3）/ * H型= 3接地分支* /
{
fscanf（FP1，“％D，％ D，％F，％F，％F \ N“，植绿[I]。P1，植绿[I]。P2，植绿[I]，R＆植绿[I]，X，＆植绿[ ，B）;
植绿[I]。 KT = 1;，
植绿[I]。 S = 3;
DD + +;
}
}
/ ******************* ************读取互联网信息*********************************** / **
（qy! = 0）
（I = 1; <= QY; + +）/ *输入互联网状态* /
{
fscanf（FP1，“％D，％D，％F，％F，为％d”，＆hulianwang [i]于民，＆hulianwang的[I]光伏hulianwang [I]，P＆hulianwang [I] 。EPS2，hulianwang的[I]计数）;
（J = 1; <= hulianwang [I]。计数; J + +）
fscanf（FP1，“％D”， （hulianwang [I]一个[J]））;
}
FCLOSE（FP1）;
/ ********* *打开输出文件********************************************* *********** /
FP2 =打开（“out.txt的”，“W”）;
（FP2 == NULL）
{
>的printf（“无法打开文件！\ N”）;
出口（0）;
}
fprintf（FP2，“\ N ****原始数据如下************************************************ *********** \ n \ n已“）;
fprintf（FP2，”节点：％2D科数：％2D正确的地支的：％2D PQ节点编号：％2D PV节点：％2D精度1％F \ N“，N，M，DD，PQ，PV，EPS）;
fprintf（FP2，”\ N ---------- -------------------------------------------------- ----------------- \ N \ N“）;
（1 <= PQ + +）
fprintf（FP2 ，“％2D PQ节点节点P [％d〕=％F Q [％] =％F \”，我，我，jiedian [I]，P，I，jiedian []。Q）; p>就（i = PQ +1 <= PQ + PV + +）
fprintf（FP2，“节点％2D PV节点P [％] =％F V [％ e] =％f \ N“，我，我，jiedian [I]，P，I，jiedian [I]，V）;
fprintf（FP2，”节点％2D平衡节点\ N“，我）
fprintf（FP2，“，\ N -------------------------------------- ----------------------------------------- \ N \ N“）; />（I = 1; <= M; + +）
{
如果（植绿[I]。== 1）
fprintf（FP2，“支部％2D ：普通分支相关节点：％2D，2D R％=％F X =％F B =％F \“，我，植绿[I] P1，植绿[I]。P2，植绿[I]，R植绿[I]，植绿[I]）;
其他（植绿[I]。== 2）
fprintf（FP2，“分支％2D：变压器分会相关节点： ％2D，2D R％=％F X =％f KT =％F \“，我植绿[I] P1，植绿[I]。P2，植绿[I]，R，植绿[I]，X植绿[I]。克拉）;
其他
fprintf（FP2，“支部％2D：相关节点树枝在地上：％2D=％F X =％F B =％F \ “，我植绿[I]。P1，植绿[I]，R，植绿[I]，植绿[] b）条;
}
（I = 1 <= QY; + +）
{
fprintf（FP2，“Internet区域中：％2D指定PV节点：％2D提供有效：％8.5F允许误差：包括节点％8.5F ：“，hulianwang [I]数，hulianwang [I]。光伏，hulianwang [I]，P，hulianwang []。EPS2）;
（J = 1; <= hulianwang [I]。计数; J + +）
fprintf（FP2，“2D”，hulianwang [I]一个[J]）; fprintf（FP2，“\ N”）;
} <br的的 BR />}
无效油花（）/ *使用节点数据，分支机构的数据，形成一个新的号码* / {
INT I，K，J，JD [N / * JD [N] [记录]节点* / struct jiedian的温度; / *中间临时党支部连接的节点数量* /
（I = 1; <= N + +）/ *分支的节点上连接数赋0 * /
{
JD [0] = 0; }
（我= 1 <=我+ +）/ *考虑每个节点* /
{
（J = 1，J <= M，J + +）/ *考虑每个分支* /
{
（（植绿[J]。P1 = | |植绿[J]。P2 ==）&&（植绿[J]。！= 3））/ *如果这个分支的节点，不接地分支* /
JD [0] + +;
}
}
（I = 1 <PQ; + +）/ *正确的PQ节点根据分支机构连接到重新编号的数量* /
{
（J = i +1的J <= PQ研究+ +）/ *升序排序算法* /
{
（JD [I]> JD [J]）
{
TEM = jiedian [I];
jiedian [I] = jiedian [J];
jiedian [J] = TEM
（K = 1; K <= M，K + +）/ *更新分支信息* /
{
（植绿[K]。P1 == J）
植绿[K]。 P1 = I 否则，如果（植绿[K]。P2 == J）
植绿[K]。 P2 = I;
其他（植绿[K]。P1 ==我）
植绿[K]。 P1 = J;
否则，如果（植绿[K]。P2 ==）
植绿[K]。 P2 = J;
}
}
}
}
就（i = PQ +1; <PQ + PV; + +）/ * PV节点连接上重新编排的分支机构数量* /
{
（J = i +1的; J <= PQ + PV; J + +）/ *小大排序算法* /
{
如果（JD [我]> JD [J]）
{
温度= jiedian [我];
jiedian [I] = jiedian [J];
jiedian [J] = TEM
（K = 1; K <= M，K + +）/ *更新分支信息* /
> {
（植绿[K]。P1 == J）
植绿[K]。 P1 = I;
否则，如果（植绿[K]。P2 == J）
植绿[K]。 P2 = I;
其他（植绿[K]。P1 ==我）
植绿[K]。 P1 = J;
否则，如果（植绿[K]。P2 ==）
植绿[K]。 P2 = J;
}
}
}
}
fprintf（FP2，“\ N ****节点优化结果如下*********************************************** ********** \ n \ n已“）;
（I = 1; <=我+ +）
{
（jiedian [I]。== 1）
fprintf（FP2，“％2D PQ节点，节点P [％d〕=％F Q [％] =％F \”，我，我，jiedian [I]，P，I，jiedian [I]。Q）;
如果（jiedian [I]。== 2）
fprintf（FP2，“节点％2D PV节点P [％d的=％F V [％d〕=％F \“，我，我，jiedian [I]，P，我jiedian [I]，V）;
（jiedian [I]。= = 3）
fprintf（FP2，“％2D平衡节点的节点\ N”）;
}
}
的无效form_y（）/ *使用分支数据形成Y，要注意地面腿* /
{
I，J，K;
float类型的
（I = 1 <= N;我+ +）
（J = 1; <= N，J + +）
G [我] [J] = B [I] [J] = 0;
（1 <= M; + +）/ *节点导纳矩阵主对角线准入* /
（J = 1; <= N; J + +）
（植绿[I] P1 ==| |植绿[I]。P2 == J）
{
S =植绿[I]。 R *植绿[I]。 R +植绿[I]。 X *植绿[I]。 X;
（S == 0）;
G [J] [J] + =植绿[I]。 R / S;
B [J] [J] + =植绿[我]。 X / S;
（植绿[I]。== 1）/ *如果是普通的分支* / B [J] [J] + =植绿[I]。 B / 2;
}
其他（植绿[I]。== 2）/ *如果是普通变压器支路* /
{
（植绿[ P1 == J）
{
G [J] [J] + =（植绿[I]，R / S *（1-植绿[]克拉））/（植绿我克拉*植绿[I]克拉）;
B [J] [J] + =（-植绿[I] X / S *（1-植绿[]。克拉））/（植绿[I] KT *植绿[I]克拉）;
}
其他（植绿[I] P2 == j）条
{
[J] [J] + =（植绿[I]，R / S（植绿[I]千吨-1））/植绿[I]。克拉;
B [J] [J] + =（-植绿[I]，X / S *（植绿[I] KT-1））/植绿[I]。克拉;
}
}
其他（植绿[I]。== 3）/ *如果是在地面上分支* /
{
乙[J] [J] + =植绿[I]。 b;
}
}
（K = 1; K <= M，K + +）/ *节点导纳矩阵非准入的主对角线上* /
{
=植绿[K]。 P1;
=植绿[K]。 P2;
S =植绿[K]。 R *植绿[K]。 R +植绿[K]。 X *植绿[K]。 X;
（S == 0）;
G [我] [J] + =植绿[K]。 R / S;
B [我] [J] + =植绿[K]。 X / S;
（植绿[K]。克拉= 1.0）
{
G [我] [J] / =植绿[K]。克拉;
B [我] [J] / =植绿[K]。克拉;
}
G [J] [I] = G [I] [J];
B [J] [] = B [I] [J];
>}
fprintf（FP2，“\ N \ N ******************************节点导纳矩阵********************************** \ N“）;
（I = 1 <= N;我+ +）
{
fprintf（FP2，“\”）;
（J = 1; <= N研究+） BR /> fprintf（FP2，“％+ J％8.5F 8.5F，G [I] [J]，B [I] [J]）;
}
}
BR />无效form_j，（）/ *使用节点数据和Y形式的J * /
{
浮法EI，FI，A = 0，B = 0;
I1，J1，K1，I，J，K;
（I = 1; <= 2 *（PQ + PV）+1; + +）
（J = 1; J <= 2 *（PQ + PV）+1; J + +）
YKB [I] [J] = 0;
（I = 1; < = PQ + +）
（J = 1; J <; J + +）
{I1 =;
J1 = J;
EI = jiedian I]。 E;
Fi = jiedian的[I]。 F;
（i! = J）/ *找我！ = j的H，N，J，L * /
{YKB *-1 [2] [2 * J-1] = G [期1] [J1] * EI + B [期1] [J1] * FI / * H * /
YKB [2 *的i-1] [2 * J] = B [I1] [J1] * EI + G [I1] [J1] * FI / * N * /
YKB [2] [2 * J-1] = B [I1]？[J1] * EI + G [I1] [J1] * FI / * J * / YKB [2] [2 * J] =-G [I1] [J1] * EI-B [I1] [J1] * FI / * L * /
}
BR /> / *其他寻求当我= J H，N，J，K * /
{A = 0，B = 0;
（K = 1; K <= N，K + +）
（k! = I）
{K1 = K;
= A + G [I1] [K1] * jiedian [K]。 EB [I1] [K1] * jiedian [K]。 F;
B = B + G [I1] [K1] * jiedian [K]。 F + B [期1] [K1] * jiedian [●]。 E;
}
YKB [2 *的i-1] [2 * J-1] = 2 * G [I1] [I1] * jiedian [I]。 E + A / * H * /
YKB [2] [2 * J] = -2 * B [I1] [I1] * jiedian [I]。 F +一个; / * L * /
YKB [2 *的i-1] [2 * j]的= 2 * G [期1] [1期1] * jiedian []。 F + B; / * N * /
YKB [2] [2 * J-1] = -2 * B [I1] [I1] * jiedian [I]。 EB / * J * /
} }
就（i = PQ +1; <= PQ + PV; + +）/ *子矩阵形成光伏节点* /
（十= 1; J <; J + +）
{I1 =;
J1 = J;
EI = jiedian [我]。 E;
Fi = jiedian的[I]。 F;
（i! = J）/ *找我！ =当H，N * /
{
YKB [2 *的i-1] [2 * J-1] = G [期1] [J1] * EI + B [期1] [J1] * FI / * H * /
YKB [2 *的i-1] [2 * J] = B [I1] [J1] * EI + G [I1] [J1] * FI / * N * /
}
其他/ *求时，我= J H，N，R，S * /
{A = 0，B = 0;
（K = 1，K <= N，K +）
（k! = I）
{K1 = K;
+ = G [I1] [K1] * jiedian [K]。 EB [I1] [K1] * jiedian [K]。 F;
B + = G [I1] [K1] * jiedian [K]。 F + B [期1] [K1] * jiedian [●]。 E;
}
YKB [2 *的i-1] [2 * J-1] = 2 * G [I1] [I1] * jiedian [I]。 E + A / * H * /
YKB [2 *的i-1] [2] = 2 * G [I1] [I1] * jiedian [I]。 F + B; / * N * /
YKB [2] [2 * J-1] = 2 * EI / * R * /
YKB [2] [2 * J ] = 2 * FI / * S * /
}
}
fprintf（FP2，“\ N \ N \ N ********** *********************雅可比************** *********** \ N“）;
（I = 1; <= 2 *（PQ + PV）; + +） {
fprintf（FP2，“\”）;
（J = 1; <= 2 *（PQ + PV）; J + +）
fprintf（FP2，“％8.5 F“，YKB [I] [J]）;
}
}
无效form_yz（）/ *使用的J的外形表* / { BR /> INT I，J，K，
浮法A [2 * M] [2 * M]，L [2 * M] [2 * M]，R [2 * M] [2 * M]。
浮法X，Y，Z;
（I = 1; <= 2 *（PQ + PV）; + +）
（J = 1;研究<= 2 *（PQ + PV）; J + +）
{[I] [J] = 0，R [I] [J] = 0，L [I] [J] = 0 ;}
（I = 1; <= 2 *（PQ + PV）; + +）
（J = 1; <= 2 *（PQ + PV）; J + +） A [I] [J] = YKB [I] [J];
（I = 1; <= 2 *（PQ + PV）+ + ）
{
（J = 1; <= 2 *（PQ + PV），J + +）
{
如果（> J）
{
R [我] [J] = 0;
X = 0;
（K = 1; K <= J-1，K +）
{ BR /> X = X + L [I] [K] * R [K] [J];
}
L [我] [J] = A [I] [J] X;
}
（I == J）
{
R [我] [J] = 1;
Y = 0;
（K = 1; K <= I-1，K +）
{
Y = Y + L [我] [K] * R [K] [I];
} L [I] [I] = [I] [J];
}
（I <j）的
{
L [我] [J ] = 0;
Z = 0;
（K = 1; K <= I-1，K +）
Z = Z + L [我] [K] * R [K] [J]; }
R [我] [J] =（[I] [J] Z）/ L [] []; BR />}
}
}
（I = 1; <= 2 *（PQ + PV）; + +）
{（J = 1; <= 2 *（PQ + PV）; J + +）
{
（I <J）
饮子[I] [J] = R [我] [J];
（I == J）
饮子[I] [J] = 1 / L [I] [J];
（> j的）
饮子[我] [J] = L [I] [J];
}
}
fprintf（FP2，“\ N \ N \ N ********************************因子表**************** ************************* \ N“）;
（I = 1; <= 2 *（ PQ + PV）; + +）{
fprintf（FP2，“\”）;
（J = 1; <= 2 *（PQ + PV）; + +）
fprintf（FP2，“％8.5F”，银子岩[I] [J]）;
}
}
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
