model.src

一种求解目标函数最小化的MATLAB仿真程序.

model:=proc(X)
global n_gen,m_pv,s_pq,f1,f2;
n_gen:=54:
m_pv:=0:
s_pq:=64:
f1:=array(1..118,1..1):
f2:=array(1..64,1..1):
f1[1,1]:=((((((((e[1]*(e[1]*g[1,1]))+(((b[1,22]*e[1])*(e[22]*sin((0-th[21]))))+((cos((0-th[21]))*e[22])*
(e[1]*g[1,22]))))+(((b[1,31]*e[1])*(e[31]*sin((0-th[30]))))+((cos((0-th[30]))*e[31])*(e[1]*
g[1,31]))))+(((b[1,36]*e[1])*(e[36]*sin((0-th[35]))))+((cos((0-th[35]))*e[36])*(e[1]*g[1,36]))))+
(((b[1,81]*e[1])*(e[81]*sin((0-th[80]))))+((cos((0-th[80]))*e[81])*(e[1]*g[1,81]))))+(((b[1,93]*
e[1])*(e[93]*sin((0-th[92]))))+((cos((0-th[92]))*e[93])*(e[1]*g[1,93]))))+(((b[1,95]*e[1])*
(e[95]*sin((0-th[94]))))+((cos((0-th[94]))*e[95])*(e[1]*g[1,95]))))-p[1]):

f1[2,1]:=((((e[2]*(e[2]*g[2,2]))+(((b[2,55]*e[2])*(e[55]*sin((th[1]-th[54]))))+((cos((th[1]-th[54]))*e[55])*
(e[2]*g[2,55]))))+(((b[2,56]*e[2])*(e[56]*sin((th[1]-th[55]))))+((cos((th[1]-th[55]))*e[56])*(e[2]*
g[2,56]))))-p[2]):

f1[3,1]:=((((e[3]*(e[3]*g[3,3]))+(((b[3,57]*e[3])*(e[57]*sin((th[2]-th[56]))))+((cos((th[2]-th[56]))*e[57])*
(e[3]*g[3,57]))))+(((b[3,60]*e[3])*(e[60]*sin((th[2]-th[59]))))+((cos((th[2]-th[59]))*e[60])*(e[3]*
g[3,60]))))-p[3]):

f1[4,1]:=((((e[4]*(e[4]*g[4,4]))+(((b[4,57]*e[4])*(e[57]*sin((th[3]-th[56]))))+((cos((th[3]-th[56]))*e[57])*
(e[4]*g[4,57]))))+(((b[4,58]*e[4])*(e[58]*sin((th[3]-th[57]))))+((cos((th[3]-th[57]))*e[58])*(e[4]*
g[4,58]))))-p[4]):

f1[5,1]:=(((((e[5]*(e[5]*g[5,5]))+(((b[5,57]*e[5])*(e[57]*sin((th[4]-th[56]))))+((cos((th[4]-th[56]))*e[57])*
(e[5]*g[5,57]))))+(((b[5,59]*e[5])*(e[59]*sin((th[4]-th[58]))))+((cos((th[4]-th[58]))*e[59])*(e[5]*
g[5,59]))))+(((b[5,71]*e[5])*(e[71]*sin((th[4]-th[70]))))+((cos((th[4]-th[70]))*e[71])*(e[5]*g[5,71]))))-
p[5]):

f1[6,1]:=(((e[6]*(e[6]*g[6,6]))+(((b[6,59]*e[6])*(e[59]*sin((th[5]-th[58]))))+((cos((th[5]-th[58]))*e[59])*
(e[6]*g[6,59]))))-p[6]):

f1[7,1]:=(((((((((e[7]*(e[7]*g[7,7]))+(((b[7,55]*e[7])*(e[55]*sin((th[6]-th[54]))))+((cos((th[6]-th[54]))*e[55])*
(e[7]*g[7,55]))))+(((b[7,56]*e[7])*(e[56]*sin((th[6]-th[55]))))+((cos((th[6]-th[55]))*e[56])*(e[7]*
g[7,56]))))+(((b[7,58]*e[7])*(e[58]*sin((th[6]-th[57]))))+((cos((th[6]-th[57]))*e[58])*(e[7]*g[7,58]))))+
(((b[7,60]*e[7])*(e[60]*sin((th[6]-th[59]))))+((cos((th[6]-th[59]))*e[60])*(e[7]*g[7,60]))))+(((b[7,62]*
e[7])*(e[62]*sin((th[6]-th[61]))))+((cos((th[6]-th[61]))*e[62])*(e[7]*g[7,62]))))+(((b[7,63]*e[7])*
(e[63]*sin((th[6]-th[62]))))+((cos((th[6]-th[62]))*e[63])*(e[7]*g[7,63]))))+(((b[7,117]*e[7])*(e[117]*
sin((th[6]-th[116]))))+((cos((th[6]-th[116]))*e[117])*(e[7]*g[7,117]))))-p[7]):

f1[8,1]:=(((((((e[8]*(e[8]*g[8,8]))+(((b[8,10]*e[8])*(e[10]*sin((th[7]-th[9]))))+((cos((th[7]-th[9]))*e[10])*
(e[8]*g[8,10]))))+(((b[8,61]*e[8])*(e[61]*sin((th[7]-th[60]))))+((cos((th[7]-th[60]))*e[61])*(e[8]*
g[8,61]))))+(((b[8,62]*e[8])*(e[62]*sin((th[7]-th[61]))))+((cos((th[7]-th[61]))*e[62])*(e[8]*g[8,62]))))+
(((b[8,64]*e[8])*(e[64]*sin((th[7]-th[63]))))+((cos((th[7]-th[63]))*e[64])*(e[8]*g[8,64]))))+(((b[8,72]*
e[8])*(e[72]*sin((th[7]-th[71]))))+((cos((th[7]-th[71]))*e[72])*(e[8]*g[8,72]))))-p[8]):

f1[9,1]:=((((e[9]*(e[9]*g[9,9]))+(((b[9,10]*e[9])*(e[10]*sin((th[8]-th[9]))))+((cos((th[8]-th[9]))*e[10])*
(e[9]*g[9,10]))))+(((b[9,64]*e[9])*(e[64]*sin((th[8]-th[63]))))+((cos((th[8]-th[63]))*e[64])*(e[9]*
g[9,64]))))-p[9]):

f1[10,1]:=((((((((b[10,8]*e[10])*(e[8]*sin((th[9]-th[7]))))+((cos((th[9]-th[7]))*e[8])*(e[10]*g[10,8])))+(((b[10,9]*
e[10])*(e[9]*sin((th[9]-th[8]))))+((cos((th[9]-th[8]))*e[9])*(e[10]*g[10,9]))))+(e[10]*(e[10]*
g[10,10])))+(((b[10,17]*e[10])*(e[17]*sin((th[9]-th[16]))))+((cos((th[9]-th[16]))*e[17])*(e[10]*g[10,17]))))+
(((b[10,65]*e[10])*(e[65]*sin((th[9]-th[64]))))+((cos((th[9]-th[64]))*e[65])*(e[10]*g[10,65]))))-p[10]):

f1[11,1]:=(((((e[11]*(e[11]*g[11,11]))+(((b[11,31]*e[11])*(e[31]*sin((th[10]-th[30]))))+((cos((th[10]-th[30]))*e[31])*
(e[11]*g[11,31]))))+(((b[11,32]*e[11])*(e[32]*sin((th[10]-th[31]))))+((cos((th[10]-th[31]))*e[32])*(e[11]*
g[11,32]))))+(((b[11,68]*e[11])*(e[68]*sin((th[10]-th[67]))))+((cos((th[10]-th[67]))*e[68])*(e[11]*g[11,68]))))-
p[11]):

f1[12,1]:=(((((e[12]*(e[12]*g[12,12]))+(((b[12,13]*e[12])*(e[13]*sin((th[11]-th[12]))))+((cos((th[11]-th[12]))*e[13])*
(e[12]*g[12,13]))))+(((b[12,14]*e[12])*(e[14]*sin((th[11]-th[13]))))+((cos((th[11]-th[13]))*e[14])*(e[12]*
g[12,14]))))+(((b[12,68]*e[12])*(e[68]*sin((th[11]-th[67]))))+((cos((th[11]-th[67]))*e[68])*(e[12]*g[12,68]))))-
p[12]):

f1[13,1]:=((((((b[13,12]*e[13])*(e[12]*sin((th[12]-th[11]))))+((cos((th[12]-th[11]))*e[12])*(e[13]*g[13,12])))+(e[13]*
(e[13]*g[13,13])))+(((b[13,71]*e[13])*(e[71]*sin((th[12]-th[70]))))+((cos((th[12]-th[70]))*e[71])*(e[13]*
g[13,71]))))-p[13]):

f1[14,1]:=((((((((b[14,12]*e[14])*(e[12]*sin((th[13]-th[11]))))+((cos((th[13]-th[11]))*e[12])*(e[14]*g[14,12])))+(e[14]*
(e[14]*g[14,14])))+(((b[14,16]*e[14])*(e[16]*sin((th[13]-th[15]))))+((cos((th[13]-th[15]))*e[16])*(e[14]*
g[14,16]))))+(((b[14,69]*e[14])*(e[69]*sin((th[13]-th[68]))))+((cos((th[13]-th[68]))*e[69])*(e[14]*g[14,69]))))+
(((b[14,116]*e[14])*(e[116]*sin((th[13]-th[115]))))+((cos((th[13]-th[115]))*e[116])*(e[14]*g[14,116]))))-p[14]):

f1[15,1]:=(((((e[15]*(e[15]*g[15,15]))+(((b[15,16]*e[15])*(e[16]*sin((th[14]-th[15]))))+((cos((th[14]-th[15]))*e[16])*
(e[15]*g[15,16]))))+(((b[15,64]*e[15])*(e[64]*sin((th[14]-th[63]))))+((cos((th[14]-th[63]))*e[64])*(e[15]*
g[15,64]))))+(((b[15,70]*e[15])*(e[70]*sin((th[14]-th[69]))))+((cos((th[14]-th[69]))*e[70])*(e[15]*g[15,70]))))-
p[15]):

f1[16,1]:=(((((((((b[16,14]*e[16])*(e[14]*sin((th[15]-th[13]))))+((cos((th[15]-th[13]))*e[14])*(e[16]*g[16,14])))+(((b[16,15]*
e[16])*(e[15]*sin((th[15]-th[14]))))+((cos((th[15]-th[14]))*e[15])*(e[16]*g[16,15]))))+(e[16]*(e[16]*
g[16,16])))+(((b[16,53]*e[16])*(e[53]*sin((th[15]-th[52]))))+((cos((th[15]-th[52]))*e[53])*(e[16]*g[16,53]))))+
(((b[16,68]*e[16])*(e[68]*sin((th[15]-th[67]))))+((cos((th[15]-th[67]))*e[68])*(e[16]*g[16,68]))))+(((b[16,115]*
e[16])*(e[115]*sin((th[15]-th[114]))))+((cos((th[15]-th[114]))*e[115])*(e[16]*g[16,115]))))-p[16]):

f1[17,1]:=((((((((b[17,10]*e[17])*(e[10]*sin((th[16]-th[9]))))+((cos((th[16]-th[9]))*e[10])*(e[17]*g[17,10])))+(e[17]*
(e[17]*g[17,17])))+(((b[17,18]*e[17])*(e[18]*sin((th[16]-th[17]))))+((cos((th[16]-th[17]))*e[18])*(e[17]*
g[17,18]))))+(((b[17,74]*e[17])*(e[74]*sin((th[16]-th[73]))))+((cos((th[16]-th[73]))*e[74])*(e[17]*g[17,74]))))+
(((b[17,78]*e[17])*(e[78]*sin((th[16]-th[77]))))+((cos((th[16]-th[77]))*e[78])*(e[17]*g[17,78]))))-p[17]):

f1[18,1]:=((((((b[18,17]*e[18])*(e[17]*sin((th[17]-th[16]))))+((cos((th[17]-th[16]))*e[17])*(e[18]*g[18,17])))+(e[18]*
(e[18]*g[18,18])))+(((b[18,73]*e[18])*(e[73]*sin((th[17]-th[72]))))+((cos((th[17]-th[72]))*e[73])*(e[18]*
g[18,73]))))-p[18]):

f1[19,1]:=((((((e[19]*(e[19]*g[19,19]))+(((b[19,20]*e[19])*(e[20]*sin((th[18]-th[19]))))+((cos((th[18]-th[19]))*e[20])*
(e[19]*g[19,20]))))+(((b[19,74]*e[19])*(e[74]*sin((th[18]-th[73]))))+((cos((th[18]-th[73]))*e[74])*(e[19]*
g[19,74]))))+(((b[19,76]*e[19])*(e[76]*sin((th[18]-th[75]))))+((cos((th[18]-th[75]))*e[76])*(e[19]*g[19,76]))))+
(((b[19,77]*e[19])*(e[77]*sin((th[18]-th[76]))))+((cos((th[18]-th[76]))*e[77])*(e[19]*g[19,77]))))-p[19]):

f1[20,1]:=(((((((b[20,19]*e[20])*(e[19]*sin((th[19]-th[18]))))+((cos((th[19]-th[18]))*e[19])*(e[20]*g[20,19])))+(e[20]*
(e[20]*g[20,20])))+(((b[20,22]*e[20])*(e[22]*sin((th[19]-th[21]))))+((cos((th[19]-th[21]))*e[22])*(e[20]*
g[20,22]))))+(((b[20,77]*e[20])*(e[77]*sin((th[19]-th[76]))))+((cos((th[19]-th[76]))*e[77])*(e[20]*g[20,77]))))-
p[20]):

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

f1[22,1]:=(((((((((((((b[22,1]*e[22])*(e[1]*sin(th[21])))+((cos(th[21])*e[1])*(e[22]*g[22,1])))+(((b[22,20]*e[22])*(e[20]*
sin((th[21]-th[19]))))+((cos((th[21]-th[19]))*e[20])*(e[22]*g[22,20]))))+(e[22]*(e[22]*g[22,22])))+(((b[22,23]*
e[22])*(e[23]*sin((th[21]-th[22]))))+((cos((th[21]-th[22]))*e[23])*(e[22]*g[22,23]))))+(((b[22,30]*e[22])*
(e[30]*sin((th[21]-th[29]))))+((cos((th[21]-th[29]))*e[30])*(e[22]*g[22,30]))))+(((b[22,80]*e[22])*(e[80]*
sin((th[21]-th[79]))))+((cos((th[21]-th[79]))*e[80])*(e[22]*g[22,80]))))+(((b[22,81]*e[22])*(e[81]*sin((th[21]-
th[80]))))+((cos((th[21]-th[80]))*e[81])*(e[22]*g[22,81]))))+(((b[22,82]*e[22])*(e[82]*sin((th[21]-th[81]))))+
((cos((th[21]-th[81]))*e[82])*(e[22]*g[22,82]))))+(((b[22,83]*e[22])*(e[83]*sin((th[21]-th[82]))))+((cos((th[21]-
th[82]))*e[83])*(e[22]*g[22,83]))))+(((b[22,84]*e[22])*(e[84]*sin((th[21]-th[83]))))+((cos((th[21]-th[83]))*
e[84])*(e[22]*g[22,84]))))-p[22]):

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

f1[24,1]:=(((((((b[24,23]*e[24])*(e[23]*sin((th[23]-th[22]))))+((cos((th[23]-th[22]))*e[23])*(e[24]*g[24,23])))+(e[24]*
(e[24]*g[24,24])))+(((b[24,25]*e[24])*(e[25]*sin((th[23]-th[24]))))+((cos((th[23]-th[24]))*e[25])*(e[24]*
g[24,25]))))+(((b[24,26]*e[24])*(e[26]*sin((th[23]-th[25]))))+((cos((th[23]-th[25]))*e[26])*(e[24]*g[24,26]))))-
p[24]):

f1[25,1]:=(((((((((b[25,23]*e[25])*(e[23]*sin((th[24]-th[22]))))+((cos((th[24]-th[22]))*e[23])*(e[25]*g[25,23])))+(((b[25,24]*
e[25])*(e[24]*sin((th[24]-th[23]))))+((cos((th[24]-th[23]))*e[24])*(e[25]*g[25,24]))))+(e[25]*(e[25]*
g[25,25])))+(((b[25,26]*e[25])*(e[26]*sin((th[24]-th[25]))))+((cos((th[24]-th[25]))*e[26])*(e[25]*g[25,26]))))+
(((b[25,87]*e[25])*(e[87]*sin((th[24]-th[86]))))+((cos((th[24]-th[86]))*e[87])*(e[25]*g[25,87]))))+(((b[25,88]*
e[25])*(e[88]*sin((th[24]-th[87]))))+((cos((th[24]-th[87]))*e[88])*(e[25]*g[25,88]))))-p[25]):

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

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

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

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

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

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