【文章內(nèi)容簡(jiǎn)介】 rgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。endend end OrgS。%四：創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓PVU=zeros(nh1,1)。t=0。for i=1:nif B2(i,6)==3t=t+1。PVU(t,1)=B2(i,3)。end end PVU。%五：創(chuàng)建DetaS,用于存儲(chǔ)有功功率、無(wú)功功率和電壓幅值的不平衡量h=0。for i=1:n%對(duì)PQ節(jié)點(diǎn)的處理if i~=isbamp。B2(i,6)==2h=h+1。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。end end t=0。for i=1:n%六：對(duì)PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0if i~=isbamp。B2(i,6)==3h=h+1。t=t+1。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。end end DetaS。%七：創(chuàng)建I，用于存儲(chǔ)節(jié)點(diǎn)電流參數(shù)i=zeros(n1,1)。h=0。for i=1:nif i~=isbh=h+1。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。end end I。%八：創(chuàng)建Jacbi(雅可比矩陣)Jacbi=zeros(2*n2)。h=0。k=0。for i=1:n%對(duì)PQ節(jié)點(diǎn)的處理if B2(i,6)==2h=h+1。for j=1:nif j~=isbk=k+1。if i==j%對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。else%非對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。endif k==(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。endendendend end k=0。for i=1:n%對(duì)PV節(jié)點(diǎn)的處理if B2(i,6)==3h=h+1。for j=1:nif j~=isbk=k+1。if i==j%對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。Jacbi(2*h,2*k)=2*real(B2(i,3))。else%非對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=0。Jacbi(2*h,2*k)=0。endif k==(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。endendendend end Jacbi。%九：求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量DetaU=zeros(2*n2,1)。DetaU=inv(Jacbi)*DetaS。DetaU。%修正節(jié)點(diǎn)電壓 j=0。for i=1:n%對(duì)PQ節(jié)點(diǎn)處理if B2(i,6)==2j=j+1。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。end end for i=1:n%對(duì)PV節(jié)點(diǎn)的處理if B2(i,6)==3j=j+1。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。end end B2。%十：開始循環(huán)********************************************************************** while abs(max(DetaU))pr OrgS=zeros(2*n2,1)。%！!初始功率參數(shù)在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣h=0。j=0。for i=1:nif i~=isbamp。B2(i,6)==2h=h+1。for j=1:nOrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。endend end for i=1:nif i~=isbamp。B2(i,6)==3h=h+1。for j=1:nOrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。endend end OrgS。%創(chuàng)建DetaS h=0。for i=1:nif i~=isbamp。B2(i,6)==2h=h+1。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。end end t=0。for i=1:nif i~=isbamp。B2(i,6)==3h=h+1。t=t+1。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。end end DetaS。%創(chuàng)建I i=zeros(n1,1)。h=0。for i=1:nif i~=isbh=h+1。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。end end I。%創(chuàng)建Jacbi Jacbi=zeros(2*n2)。h=0。k=0。for i=1:nif B2(i,6)==2h=h+1。for j=1:nif j~=isbk=k+1。if i==jJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。elseJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。Jacbi(2*h,2*k)=Jacbi(2*h