基于matlab的電力系統(tǒng)潮流計算-文庫吧資料
2025-01-22 10:52本頁面
【正文】 g=lr/(lr*lr+lx*lx)。 lb1=tb。 lr=tR。 {i=tHeadnode1。tline+K。}line[i].X=H[i][4]。line[i].Headnode=(int)H[i][1]。 {iK。 } }temp)。j6。 { iK。 /*函數(shù)申明*/ float g[N][N]={},b[N][N]={}。 float q[N]={1,1,1,}。 float p[N]={1,}。 float r[N]={}。struct line *t。 float b。 float X。 float R。 int Endnode。 int Headnode。}}struct line{int Lindex。jm。i){ }for(i=m1。 y[i]=y[i]y[k]*a[i][k]/a[k][k]。jm。i++) for(i=k+1。km1。j++) a[i][j]=b[i+n][j+n]。 for(j=0。 y[i]=c[i]。im。int i,j,k。4 0。2 0。0 5+0i 0 3。0 21i 1 0 0 2。4 2 0 1。1 3 + 0 1 0。 if B2(i,6)==3 if B2(i,6)==2DetaU%修正節(jié)點電壓j=0。 endendJacbiDetaU=zeros(2*n2,1)。 end if k==(n1) else Jacbi(2*h,2*k)=2*real(B2(i,3))。 Jacbi(2*h,2*k1)=2*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))+real(I(h,1))。 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))。 if j~=isb if B2(i,6)==3 endendk=0。 end if k==(n1) else Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。 Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(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*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。 if j~=isb if B2(i,6)==2k=0。 endendI%創(chuàng)建JacbiJacbi=zeros(2*n2)。 I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。 h=h+1。for i=1:n endendDetaS%創(chuàng)建Ii=zeros(n1,1)。 DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。 DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。 t=t+1。 h=h+1。 if i~=isbamp。 endendt=0。 DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。 DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。 h=h+1。 if i~=isbamp。 endendOrgS%創(chuàng)建DetaSh=0。 for j=1:n h=h+1。 if i~=isbamp。 end 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)))。 OrgS(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)))。B2(i,6)==2for i=1:n %!!!初始功率參數(shù)在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣h=0。 endendB2%開始循環(huán)**********************************************************************while abs(max(DetaU))prOrgS=zeros(2*n2,1)。 B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。 j=j+1。 %對PV節(jié)點的處理 endendfor i=1:n B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。 j=j+1。 %對PQ節(jié)點處理for i=1:nDetaU=inv(Jacbi)*DetaS。 end end k=0。 end Jacbi(2*h,2*k)=0。 Jacbi(2*h,2*k1)=0。 Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。 Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。 else Jacbi(2*h,2*k)=2*real(B2(i,3))。 Jacbi(2*h,2*k1)=2*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
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1