【正文】 for k3=k1:N0 for k2=k1:N1 J(k3,k2)=J(k3,k2)J(k3,k)*J(k,k2)。 end J(k3,k)=0。 end end end%用高斯消去法解w=J*V for k=3:2:N01%修正各節(jié)點電壓 L=(k+1)./2。 e(L)=e(L)J(k,N)。 k1=k+1。 f(L)=f(L)J(k1,N)。 end for k=3:N0%判斷是否收斂 DET=abs(J(k,N))。 if DET=pr IT2=IT2+1。 endendToc附錄C 遺傳算法無功優(yōu)化matlab程序遺傳算法無功優(yōu)化matlab程序分為6個主子程序，每個子程序又包含其他在程序。6個主子程序分別為讀取數(shù)據(jù)read子程序；初始解種群initial子程序；選擇selection子程序；交叉crossover子程序；變異mutition子程序；計算適應(yīng)度adapting子程序。，若想運行其他節(jié)點系統(tǒng)，只需要改變這兩個數(shù)據(jù)文件，而不必修改程序。，不再增加說明；。read子程序：function [NGAmax,popsize,NV,Gnumber,NLC,NT,NCD,NLD,VGmax,VGmin,Vjmax,Vjmin,LC,T,CD,LD,rmin]=read(node) if node==3 fid = fopen(39。39。,39。r39。)。 elseif node==4 fid = fopen(39。39。,39。r39。)。 elseif node==5 fid = fopen(39。39。,39。r39。)。 elseif node==9 fid = fopen(39。39。,39。r39。)。 elseif node==10 fid = fopen(39。39。,39。r39。)。 elseif node==11 fid = fopen(39。39。,39。r39。)。 elseif node==13 fid = fopen(39。39。,39。r39。)。 elseif node==14 fid = fopen(39。39。,39。r39。)。 elseif node==30 fid = fopen(39。39。,39。r39。)。 elseif node==39 fid = fopen(39。39。,39。r39。)。 elseif node==43 fid = fopen(39。39。,39。r39。)。 elseif node==57 fid = fopen(39。39。,39。r39。)。 elseif node==118 fid = fopen(39。39。,39。r39。)。 elseif node==145 fid = fopen(39。39。,39。r39。)。 elseif node==162 fid = fopen(39。39。,39。r39。)。 elseif node==300 fid = fopen(39。39。,39。r39。)。 else disp(39。非標準IEEE數(shù)據(jù)！39。)。 beak。 end NGAmax=fscanf(fid,39。%d39。,[1,1])。 %最大進化代數(shù) popsize=fscanf(fid,39。%d39。,[1,1])。 %種群規(guī)模 NV=fscanf(fid,39。%d39。,[1,1])。 %可調(diào)發(fā)電機數(shù)量 NLC=fscanf(fid,39。%d39。,[1,1])。 %連續(xù)型無功補償裝置數(shù)量 NT=fscanf(fid,39。%d39。,[1,1])。 %可調(diào)變壓器數(shù)量 NCD=fscanf(fid,39。%d39。,[1,1])。 %離散電容補償裝置數(shù)量 NLD=fscanf(fid,39。%d39。,[1,1])。 %離散電抗補償裝置數(shù)量 VGmax=fscanf(fid,39。%g39。,[1,1])。 %發(fā)電機機端電壓上限 VGmin=fscanf(fid,39。%g39。,[1,1])。 %發(fā)電機機端電壓下限 Vjmax=fscanf(fid,39。%g39。,[1,1])。 %非發(fā)電機節(jié)點電壓上限 Vjmin=fscanf(fid,39。%g39。,[1,1])。 %非發(fā)電機節(jié)點電壓下限 Gnumber=fscanf(fid,39。%g39。,[3,NV])。 %發(fā)電機所在節(jié)點號、無功上下限 LC=fscanf(fid,39。%g39。,[3,NLC])。 %連續(xù)型無功補償所在節(jié)點號、補償容量上下限 T=fscanf(fid,39。%g39。,[5,NT])。 %可調(diào)變壓器首末節(jié)點號、上下限、單檔容量 CD=fscanf(fid,39。%g39。,[4,NCD])。 %離散電容補償裝置位置、上下限組數(shù)、單組容量 LD=fscanf(fid,39。%g39。,[4,NLD])。 %離散電抗補償裝置位置、上下限檔數(shù)、單檔容量 rmin=fscanf(fid,39。%d39。,[1,1])。 %懲罰因子的最小取值 fclose(fid)。%關(guān)閉數(shù)據(jù)文件 Gnumber=Gnumber39。 LC=LC39。 T=T39。 CD=CD39。 LD=LD39。initial子程序：function [x,fitness,n,n1,isb,pr,iternu,X,B1,B2]=initial(popsize,NV,Gnumber,NLC,NT,NCD,NLD,node,VGmax,VGmin,Vjmax,Vjmin,LC,T,CD,LD,rmin) x=zeros(popsize,NV+NLC+NT+NCD+NLD)。 DV=zeros(1,NV)。 DLC=zeros(1,NLC)。 DNT=zeros(1,NT)。 DCD=zeros(1,NCD)。 DLD=zeros(1,NLD)。 fitness=zeros(1,popsize)。 [B1,B2,n,n1,isb,pr,iternu,X]=readnet(node)。%讀取網(wǎng)絡(luò)參數(shù) k=1。%disp(LC) while k=popsize for i=1:NV DV(i)=rand*(VGmaxVGmin)+VGmin。 end for i=1:NLC DLC(i)=rand*(LC(i,2)LC(i,3))+LC(i,3)。 end for i=1:NT DNT(i)=round(rand*(T(i,3)T(i,4))+T(i,4))。 end for i=1:NCD DCD(i)=round(rand*(CD(i,2)CD(i,3))+CD(i,3))。 end for i=1:NLD DLD(i)=round(rand*(LD(i,2)LD(i,3))+LD(i,3))。 end [BB1,BB2]=renewb(NV,Gnumber,NLC,NT,NCD,NLD,DV,DLC,DNT,DCD,DLD,LC,T,CD,LD,B1,B2,n,n1)。%更新控制變量數(shù)據(jù)[PL,Vf,Vq,overflag,V,Gnum,GS]=flow(BB1,BB2,n,n1,isb,pr,iternu,X,NV,Gnumber,Vjmax,Vj min)。%計算潮流 if overflag==0%潮流收斂輸出 x(k,:)=[DV,DLC,DNT,DCD,DLD]。 sumLCD=0。 for i=1:NCD sumLCD=sumLCD+DCD(i)*CD(i,4)。%離散電容補償裝置 end for i=1:NLD sumLCD=sumLCD+DLD(i)*LD(i,4)。%離散電抗補償裝置 end sumLC=0。 for i=1:NLC sumLC=sumLC+DLC(i)。%連續(xù)型無功補償裝置(靜止補償器) end fitness(k)=1e6(sumLCD*5*1e4+PL*1e3*2000*+rmin*(Vf+Vq))。 k=k+1。 end endselection子程序：function [x,fitness,fitness_avg,fitness_max,fitness_newavg,fitness_newmax]=selection(x,fitness,popsize)[x,fitness,fitness_avg,fitness_max]=GAsort(x,fitness,popsize)。X_new=x。fitness_new=fitness。sel=popsize/4。for i=1:sel X_new(i,:)=x(i,:)。 X_new(i+sel,:)=x(i,:)。 fitness_new(i)=fitness(i)。 fitness_new(i+sel)=fitness(i)。endfor i=sel+1:popsizesel X_new((i+sel),:)=x(i,:)。 fitness_new(i+sel)=fitness(i)。endfitness_newsum=0。fitness_newmax=fitness_new(1)。for i=1:popsize fitness_newsum=fitness_newsum+fitness_new(i)。 if fitness_new(i)fitness_newmax fitness_newmax=fitness_new(i)。 endendfitness_newavg=fitness_newsum/popsize。x=X_new。fitness=fitness_new。crossover子程序：function [x,fitness]=crossover(x,fitness,popsize,Fitness_newavg,Fitness_newmax,NV,NLC,NT,NCD,NLD) pc1=。 pc2=。 for i=2:popsize/2 if fitness(i)fitness(popsize/2+i) f=fitness(i)。 else f=fitness(popsize/2+i)。 end if fFitness_newavg pc=pc1(pc1pc2)*(fFitness_newavg)/(Fitness_newmaxFitness_newavg)。 else pc=pc1。 end crossr=rand。 if crossrpc crossp=round(rand(1,NV+NLC+NT+NCD+NLD))。 for m=1:NV+NLC+NT+NCD+NLD a=rand。 temp=x(i,m)。 if crossp(1,m)==1 if mNV+NLC x(i,m)=round(a*x(i,m)+(1a)*x(popsize/2+i,m))。