model.src

电力系统电压稳定研究的图形化软件

p[80]):

f2[26,1]:=(((((((g[80,21]*e[80])*(e[21]*sin((th[79]-th[20]))))-((cos((th[79]-th[20]))*e[21])*(e[80]*b[80,21])))+(((g[80,22]*
e[80])*(e[22]*sin((th[79]-th[21]))))-((cos((th[79]-th[21]))*e[22])*(e[80]*b[80,22]))))+(((g[80,79]*e[80])*
(e[79]*sin((th[79]-th[78]))))-((cos((th[79]-th[78]))*e[79])*(e[80]*b[80,79]))))+(0-(e[80]*(e[80]*
b[80,80]))))-q[80]):

f1[81,1]:=(((((((b[81,1]*e[81])*(e[1]*sin(th[80])))+((cos(th[80])*e[1])*(e[81]*g[81,1])))+(((b[81,21]*e[81])*(e[21]*
sin((th[80]-th[20]))))+((cos((th[80]-th[20]))*e[21])*(e[81]*g[81,21]))))+(((b[81,22]*e[81])*(e[22]*sin((th[80]-
th[21]))))+((cos((th[80]-th[21]))*e[22])*(e[81]*g[81,22]))))+(e[81]*(e[81]*g[81,81])))-p[81]):

f2[27,1]:=(((((((g[81,1]*e[81])*(e[1]*sin(th[80])))-((cos(th[80])*e[1])*(e[81]*b[81,1])))+(((g[81,21]*e[81])*(e[21]*
sin((th[80]-th[20]))))-((cos((th[80]-th[20]))*e[21])*(e[81]*b[81,21]))))+(((g[81,22]*e[81])*(e[22]*sin((th[80]-
th[21]))))-((cos((th[80]-th[21]))*e[22])*(e[81]*b[81,22]))))+(0-(e[81]*(e[81]*b[81,81]))))-q[81]):

f1[82,1]:=((((((b[82,21]*e[82])*(e[21]*sin((th[81]-th[20]))))+((cos((th[81]-th[20]))*e[21])*(e[82]*g[82,21])))+(((b[82,22]*
e[82])*(e[22]*sin((th[81]-th[21]))))+((cos((th[81]-th[21]))*e[22])*(e[82]*g[82,22]))))+(e[82]*(e[82]*
g[82,82])))-p[82]):

f2[28,1]:=((((((g[82,21]*e[82])*(e[21]*sin((th[81]-th[20]))))-((cos((th[81]-th[20]))*e[21])*(e[82]*b[82,21])))+(((g[82,22]*
e[82])*(e[22]*sin((th[81]-th[21]))))-((cos((th[81]-th[21]))*e[22])*(e[82]*b[82,22]))))+(0-(e[82]*
(e[82]*b[82,82]))))-q[82]):

f1[83,1]:=((((((b[83,22]*e[83])*(e[22]*sin((th[82]-th[21]))))+((cos((th[82]-th[21]))*e[22])*(e[83]*g[83,22])))+(e[83]*
(e[83]*g[83,83])))+(((b[83,87]*e[83])*(e[87]*sin((th[82]-th[86]))))+((cos((th[82]-th[86]))*e[87])*(e[83]*
g[83,87]))))-p[83]):

f2[29,1]:=((((((g[83,22]*e[83])*(e[22]*sin((th[82]-th[21]))))-((cos((th[82]-th[21]))*e[22])*(e[83]*b[83,22])))+(0-
(e[83]*(e[83]*b[83,83]))))+(((g[83,87]*e[83])*(e[87]*sin((th[82]-th[86]))))-((cos((th[82]-th[86]))*e[87])*
(e[83]*b[83,87]))))-q[83]):

f1[84,1]:=(((((((b[84,22]*e[84])*(e[22]*sin((th[83]-th[21]))))+((cos((th[83]-th[21]))*e[22])*(e[84]*g[84,22])))+(e[84]*
(e[84]*g[84,84])))+(((b[84,85]*e[84])*(e[85]*sin((th[83]-th[84]))))+((cos((th[83]-th[84]))*e[85])*(e[84]*
g[84,85]))))+(((b[84,88]*e[84])*(e[88]*sin((th[83]-th[87]))))+((cos((th[83]-th[87]))*e[88])*(e[84]*g[84,88]))))-
p[84]):

f2[30,1]:=(((((((g[84,22]*e[84])*(e[22]*sin((th[83]-th[21]))))-((cos((th[83]-th[21]))*e[22])*(e[84]*b[84,22])))+(0-
(e[84]*(e[84]*b[84,84]))))+(((g[84,85]*e[84])*(e[85]*sin((th[83]-th[84]))))-((cos((th[83]-th[84]))*e[85])*
(e[84]*b[84,85]))))+(((g[84,88]*e[84])*(e[88]*sin((th[83]-th[87]))))-((cos((th[83]-th[87]))*e[88])*(e[84]*
b[84,88]))))-q[84]):

f1[85,1]:=((((((b[85,84]*e[85])*(e[84]*sin((th[84]-th[83]))))+((cos((th[84]-th[83]))*e[84])*(e[85]*g[85,84])))+(e[85]*
(e[85]*g[85,85])))+(((b[85,86]*e[85])*(e[86]*sin((th[84]-th[85]))))+((cos((th[84]-th[85]))*e[86])*(e[85]*
g[85,86]))))-p[85]):

f2[31,1]:=((((((g[85,84]*e[85])*(e[84]*sin((th[84]-th[83]))))-((cos((th[84]-th[83]))*e[84])*(e[85]*b[85,84])))+(0-
(e[85]*(e[85]*b[85,85]))))+(((g[85,86]*e[85])*(e[86]*sin((th[84]-th[85]))))-((cos((th[84]-th[85]))*e[86])*
(e[85]*b[85,86]))))-q[85]):

f1[86,1]:=((((((b[86,23]*e[86])*(e[23]*sin((th[85]-th[22]))))+((cos((th[85]-th[22]))*e[23])*(e[86]*g[86,23])))+(((b[86,85]*
e[86])*(e[85]*sin((th[85]-th[84]))))+((cos((th[85]-th[84]))*e[85])*(e[86]*g[86,85]))))+(e[86]*(e[86]*
g[86,86])))-p[86]):

f2[32,1]:=((((((g[86,23]*e[86])*(e[23]*sin((th[85]-th[22]))))-((cos((th[85]-th[22]))*e[23])*(e[86]*b[86,23])))+(((g[86,85]*
e[86])*(e[85]*sin((th[85]-th[84]))))-((cos((th[85]-th[84]))*e[85])*(e[86]*b[86,85]))))+(0-(e[86]*
(e[86]*b[86,86]))))-q[86]):

f1[87,1]:=((((((b[87,25]*e[87])*(e[25]*sin((th[86]-th[24]))))+((cos((th[86]-th[24]))*e[25])*(e[87]*g[87,25])))+(((b[87,83]*
e[87])*(e[83]*sin((th[86]-th[82]))))+((cos((th[86]-th[82]))*e[83])*(e[87]*g[87,83]))))+(e[87]*(e[87]*
g[87,87])))-p[87]):

f2[33,1]:=((((((g[87,25]*e[87])*(e[25]*sin((th[86]-th[24]))))-((cos((th[86]-th[24]))*e[25])*(e[87]*b[87,25])))+(((g[87,83]*
e[87])*(e[83]*sin((th[86]-th[82]))))-((cos((th[86]-th[82]))*e[83])*(e[87]*b[87,83]))))+(0-(e[87]*
(e[87]*b[87,87]))))-q[87]):

f1[88,1]:=((((((b[88,25]*e[88])*(e[25]*sin((th[87]-th[24]))))+((cos((th[87]-th[24]))*e[25])*(e[88]*g[88,25])))+(((b[88,84]*
e[88])*(e[84]*sin((th[87]-th[83]))))+((cos((th[87]-th[83]))*e[84])*(e[88]*g[88,84]))))+(e[88]*(e[88]*
g[88,88])))-p[88]):

f2[34,1]:=((((((g[88,25]*e[88])*(e[25]*sin((th[87]-th[24]))))-((cos((th[87]-th[24]))*e[25])*(e[88]*b[88,25])))+(((g[88,84]*
e[88])*(e[84]*sin((th[87]-th[83]))))-((cos((th[87]-th[83]))*e[84])*(e[88]*b[88,84]))))+(0-(e[88]*
(e[88]*b[88,88]))))-q[88]):

f1[89,1]:=(((((((b[89,26]*e[89])*(e[26]*sin((th[88]-th[25]))))+((cos((th[88]-th[25]))*e[26])*(e[89]*g[89,26])))+(((b[89,27]*
e[89])*(e[27]*sin((th[88]-th[26]))))+((cos((th[88]-th[26]))*e[27])*(e[89]*g[89,27]))))+(((b[89,28]*e[89])*
(e[28]*sin((th[88]-th[27]))))+((cos((th[88]-th[27]))*e[28])*(e[89]*g[89,28]))))+(e[89]*(e[89]*g[89,89])))-
p[89]):

f2[35,1]:=(((((((g[89,26]*e[89])*(e[26]*sin((th[88]-th[25]))))-((cos((th[88]-th[25]))*e[26])*(e[89]*b[89,26])))+(((g[89,27]*
e[89])*(e[27]*sin((th[88]-th[26]))))-((cos((th[88]-th[26]))*e[27])*(e[89]*b[89,27]))))+(((g[89,28]*e[89])*
(e[28]*sin((th[88]-th[27]))))-((cos((th[88]-th[27]))*e[28])*(e[89]*b[89,28]))))+(0-(e[89]*(e[89]*
b[89,89]))))-q[89]):

f1[90,1]:=((((((b[90,26]*e[90])*(e[26]*sin((th[89]-th[25]))))+((cos((th[89]-th[25]))*e[26])*(e[90]*g[90,26])))+(e[90]*
(e[90]*g[90,90])))+(((b[90,91]*e[90])*(e[91]*sin((th[89]-th[90]))))+((cos((th[89]-th[90]))*e[91])*(e[90]*
g[90,91]))))-p[90]):

f2[36,1]:=((((((g[90,26]*e[90])*(e[26]*sin((th[89]-th[25]))))-((cos((th[89]-th[25]))*e[26])*(e[90]*b[90,26])))+(0-
(e[90]*(e[90]*b[90,90]))))+(((g[90,91]*e[90])*(e[91]*sin((th[89]-th[90]))))-((cos((th[89]-th[90]))*e[91])*
(e[90]*b[90,91]))))-q[90]):

f1[91,1]:=(((((((b[91,27]*e[91])*(e[27]*sin((th[90]-th[26]))))+((cos((th[90]-th[26]))*e[27])*(e[91]*g[91,27])))+(((b[91,29]*
e[91])*(e[29]*sin((th[90]-th[28]))))+((cos((th[90]-th[28]))*e[29])*(e[91]*g[91,29]))))+(((b[91,90]*e[91])*
(e[90]*sin((th[90]-th[89]))))+((cos((th[90]-th[89]))*e[90])*(e[91]*g[91,90]))))+(e[91]*(e[91]*g[91,91])))-
p[91]):

f2[37,1]:=(((((((g[91,27]*e[91])*(e[27]*sin((th[90]-th[26]))))-((cos((th[90]-th[26]))*e[27])*(e[91]*b[91,27])))+(((g[91,29]*
e[91])*(e[29]*sin((th[90]-th[28]))))-((cos((th[90]-th[28]))*e[29])*(e[91]*b[91,29]))))+(((g[91,90]*e[91])*
(e[90]*sin((th[90]-th[89]))))-((cos((th[90]-th[89]))*e[90])*(e[91]*b[91,90]))))+(0-(e[91]*(e[91]*
b[91,91]))))-q[91]):

f1[92,1]:=((((((b[92,28]*e[92])*(e[28]*sin((th[91]-th[27]))))+((cos((th[91]-th[27]))*e[28])*(e[92]*g[92,28])))+(((b[92,30]*
e[92])*(e[30]*sin((th[91]-th[29]))))+((cos((th[91]-th[29]))*e[30])*(e[92]*g[92,30]))))+(e[92]*(e[92]*
g[92,92])))-p[92]):

f2[38,1]:=((((((g[92,28]*e[92])*(e[28]*sin((th[91]-th[27]))))-((cos((th[91]-th[27]))*e[28])*(e[92]*b[92,28])))+(((g[92,30]*
e[92])*(e[30]*sin((th[91]-th[29]))))-((cos((th[91]-th[29]))*e[30])*(e[92]*b[92,30]))))+(0-(e[92]*
(e[92]*b[92,92]))))-q[92]):

f1[93,1]:=((((((((b[93,1]*e[93])*(e[1]*sin(th[92])))+((cos(th[92])*e[1])*(e[93]*g[93,1])))+(((b[93,29]*e[93])*(e[29]*
sin((th[92]-th[28]))))+((cos((th[92]-th[28]))*e[29])*(e[93]*g[93,29]))))+(((b[93,54]*e[93])*(e[54]*sin((th[92]-
th[53]))))+((cos((th[92]-th[53]))*e[54])*(e[93]*g[93,54]))))+(e[93]*(e[93]*g[93,93])))+(((b[93,98]*e[93])*
(e[98]*sin((th[92]-th[97]))))+((cos((th[92]-th[97]))*e[98])*(e[93]*g[93,98]))))-p[93]):

f2[39,1]:=((((((((g[93,1]*e[93])*(e[1]*sin(th[92])))-((cos(th[92])*e[1])*(e[93]*b[93,1])))+(((g[93,29]*e[93])*(e[29]*
sin((th[92]-th[28]))))-((cos((th[92]-th[28]))*e[29])*(e[93]*b[93,29]))))+(((g[93,54]*e[93])*(e[54]*sin((th[92]-
th[53]))))-((cos((th[92]-th[53]))*e[54])*(e[93]*b[93,54]))))+(0-(e[93]*(e[93]*b[93,93]))))+(((g[93,98]*
e[93])*(e[98]*sin((th[92]-th[97]))))-((cos((th[92]-th[97]))*e[98])*(e[93]*b[93,98]))))-q[93]):

f1[94,1]:=(((((((b[94,31]*e[94])*(e[31]*sin((th[93]-th[30]))))+((cos((th[93]-th[30]))*e[31])*(e[94]*g[94,31])))+(((b[94,32]*
e[94])*(e[32]*sin((th[93]-th[31]))))+((cos((th[93]-th[31]))*e[32])*(e[94]*g[94,32]))))+(((b[94,33]*e[94])*
(e[33]*sin((th[93]-th[32]))))+((cos((th[93]-th[32]))*e[33])*(e[94]*g[94,33]))))+(e[94]*(e[94]*g[94,94])))-
p[94]):

f2[40,1]:=(((((((g[94,31]*e[94])*(e[31]*sin((th[93]-th[30]))))-((cos((th[93]-th[30]))*e[31])*(e[94]*b[94,31])))+(((g[94,32]*
e[94])*(e[32]*sin((th[93]-th[31]))))-((cos((th[93]-th[31]))*e[32])*(e[94]*b[94,32]))))+(((g[94,33]*e[94])*
(e[33]*sin((th[93]-th[32]))))-((cos((th[93]-th[32]))*e[33])*(e[94]*b[94,33]))))+(0-(e[94]*(e[94]*
b[94,94]))))-q[94]):

f1[95,1]:=(((((((((b[95,1]*e[95])*(e[1]*sin(th[94])))+((cos(th[94])*e[1])*(e[95]*g[95,1])))+(((b[95,31]*e[95])*(e[31]*
sin((th[94]-th[30]))))+((cos((th[94]-th[30]))*e[31])*(e[95]*g[95,31]))))+(((b[95,34]*e[95])*(e[34]*sin((th[94]-
th[33]))))+((cos((th[94]-th[33]))*e[34])*(e[95]*g[95,34]))))+(((b[95,36]*e[95])*(e[36]*sin((th[94]-th[35]))))+
((cos((th[94]-th[35]))*e[36])*(e[95]*g[95,36]))))+(e[95]*(e[95]*g[95,95])))+(((b[95,118]*e[95])*(e[118]*
sin((th[94]-th[117]))))+((cos((th[94]-th[117]))*e[118])*(e[95]*g[95,118]))))-p[95]):

f2[41,1]:=(((((((((g[95,1]*e[95])*(e[1]*sin(th[94])))-((cos(th[94])*e[1])*(e[95]*b[95,1])))+(((g[95,31]*e[95])*(e[31]*
sin((th[94]-th[30]))))-((cos((th[94]-th[30]))*e[31])*(e[95]*b[95,31]))))+(((g[95,34]*e[95])*(e[34]*sin((th[94]-
th[33]))))-((cos((th[94]-th[33]))*e[34])*(e[95]*b[95,34]))))+(((g[95,36]*e[95])*(e[36]*sin((th[94]-th[35]))))-
((cos((th[94]-th[35]))*e[36])*(e[95]*b[95,36]))))+(0-(e[95]*(e[95]*b[95,95]))))+(((g[95,118]*e[95])*
(e[118]*sin((th[94]-th[117]))))-((cos((th[94]-th[117]))*e[118])*(e[95]*b[95,118]))))-q[95]):

f1[96,1]:=((((((b[96,36]*e[96])*(e[36]*sin((th[95]-th[35]))))+((cos((th[95]-th[35]))*e[36])*(e[96]*g[96,36])))+(e[96]*
(e[96]*g[96,96])))+(((b[96,97]*e[96])*(e[97]*sin((th[95]-th[96]))))+((cos((th[95]-th[96]))*e[97])*(e[96]*
g[96,97]))))-p[96]):

f2[42,1]:=((((((g[96,36]*e[96])*(e[36]*sin((th[95]-th[35]))))-((cos((th[95]-th[35]))*e[36])*(e[96]*b[96,36])))+(0-
(e[96]*(e[96]*b[96,96]))))+(((g[96,97]*e[96])*(e[97]*sin((th[95]-th[96]))))-((cos((th[95]-th[96]))*e[97])*
(e[96]*b[96,97]))))-q[96]):

f1[97,1]:=((((((b[97,37]*e[97])*(e[37]*sin((th[96]-th[36]))))+((cos((th[96]-th[36]))*e[37])*(e[97]*g[97,37])))+(((b[97,96]*
e[97])*(e[96]*sin((th[96]-th[95]))))+((cos((th[96]-th[95]))*e[96])*(e[97]*g[97,96]))))+(e[97]*(e[97]*
g[97,97])))-p[97]):

f2[43,1]:=((((((g[97,37]*e[97])*(e[37]*sin((th[96]-th[36]))))-((cos((th[96]-th[36]))*e[37])*(e[97]*b[97,37])))+(((g[97,96]*
e[97])*(e[96]*sin((th[96]-th[95]))))-((cos((th[96]-th[95]))*e[96])*(e[97]*b[97,96]))))+(0-(e[97]*
(e[97]*b[97,97]))))-q[97]):

f1[98,1]:=((((((b[98,37]*e[98])*(e[37]*sin((th[97]-th[36]))))+((cos((th[97]-th[36]))*e[37])*(e[98]*g[98,37])))+(((b[98,93]*
e[98])*(e[93]*sin((th[97]-th[92]))))+((cos((th[97]-th[92]))*e[93])*(e[98]*g[98,93]))))+(e[98]*(e[98]*
g[98,98])))-p[98]):

f2[44,1]:=((((((g[98,37]*e[98])*(e[37]*sin((th[97]-th[36]))))-((cos((th[97]-th[36]))*e[37])*(e[98]*b[98,37])))+(((g[98,93]*
e[98])*(e[93]*sin((th[97]-th[92]))))-((cos((th[97]-th[92]))*e[93])*(e[98]*b[98,93]))))+(0-(e[98]*
(e[98]*b[98,98]))))-q[98]):

f1[99,1]:=(((((((b[99,36]*e[99])*(e[36]*sin((th[98]-th[35]))))+((cos((th[98]-th[35]))*e[36])*(e[99]*g[99,36])))+(e[99]*
(e[99]*g[99,99])))+(((b[99,100]*e[99])*(e[100]*sin((th[98]-th[99]))))+((cos((th[98]-th[99]))*e[100])*(e[99]*
g[99,100]))))+(((b[99,107]*e[99])*(e[107]*sin((th[98]-th[106]))))+((cos((th[98]-th[106]))*e[107])*(e[99]*g[99,107]))))-
p[99]):

f2[45,1]:=(((((((g[99,36]*e[99])*(e[36]*sin((th[98]-th[35]))))-((cos((th[98]-th[35]))*e[36])*(e[99]*b[99,36])))+(0-