【正文】 j,3)))real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。j=0。 if B2(i,6)==2 if k==(n1) if i==j k=k+1。 end k=0。 else Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*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))。 endendI%創(chuàng)建Jacbi(雅可比矩陣)Jacbi=zeros(2*n2)。 endendDetaS%創(chuàng)建I，用于存儲節(jié)點電流參數(shù)i=zeros(n1,1)。 PVU(t,1)=B2(i,3)。 t=t+1。 %不含變壓器的支路B2=input(39。請輸入平衡節(jié)點號:isb=39。其中第一列為節(jié)點編號，第二列為節(jié)點對地%參數(shù)?；贛ATLAB的電力系統(tǒng)潮流計算%簡單潮流計算的小程序，相關(guān)的原始數(shù)據(jù)數(shù)據(jù)數(shù)據(jù)輸入格式如下：%B1是支路參數(shù)矩陣，第一列和第二列是節(jié)點編號。n=input(39。)。請輸入節(jié)點參數(shù):B2=39。 p=B1(i,1)。 q=B1(i,2)。 Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。 Y(q,p)=Y(p,q)。 Y(p,p)=Y(p,p)+1/B1(i,3)。 Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3))。 %將OrgS、DetaS初始化%創(chuàng)建OrgS，用于存儲初始功率參數(shù)h=0。 if i~=isbamp。 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)))。 end if B2(i,6)==3 %對PV節(jié)點的處理，注意這時不可再將h初始化為0h=0。h=0。 %對PQ節(jié)點的處理 %對角元素的處理 if j~=isb 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*k)=0。 end end %對PV節(jié)點的處理for i=1:n if i~=isbamp。 h=h+1。 for j=1:n if B2(i,6)==2 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*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。 else if k==(n1) endendJacbiDetaU=zeros(2*n2,1)。4 2 0 1。4 0。km1。 for(i=k+1。}}struct line{int Lindex。 int Headnode。 int Endnode。 float R。 float X。 float b。 float g[N][N]={},b[N][N]={}。 /*函數(shù)申明*/iK。 { { lg=lr/(lr*lr+lx*lx)。 b[j][j]+=lb+lb1。 g[j][i]=lg。 getch()。 printf(\n=====jie dian dao na ju zhen=====\n)。 for(i=0。j++) { max=fabs(dpu[1])。 } if (max=eps) {iN。 kr=kr+1。i++) /*算dq對應(yīng)B,僅M維,除去平衡結(jié)點和PV結(jié)點*/jN。 iq=iq+u[j]*(g[i][j]*sin(r[i]r[j])b[i][j]*cos(r[i]r[j]))。 printf(\n)。 if(max=eps) if(kq==0) val(u,g,b,r,ku,kr,h)。 }}void val(float u[N],float g[N][N],float b[N][N],float r[N],int ku, int kr,float h[N][N]) { float p[N][N]={0}。 FILE *fp1。 qs=qs+u[0]*u[i]*(g[0][i]*sin(r[0]r[i])b[0][i]*cos(r[0]r[i]))。 printf(\n=====PV jie dian gong lv=====\n)。i++) } { printf(u[%d]=%8f r[%d]=%8f\n,i+1,u[i],i+1,r[i])。i++) } } { j++) fprintf(fp1,s[%d][%d]=%8f+j%8f\n,i+1,j+1,p[i][j],q[i][j])。 { j++) fprintf(fp1,ds[%d][%d]=%8f+j%8f\n,i+1,j+1,dp[i][j],dq[i][j])。 printf(\n===========THE END==============\n)。 for(i=0。 fprintf(fp1,\n)。 fprintf(fp1,ping heng jie dian gong lv:\n)。 fprintf(fp1,sv2=%8f+j(%8f)\n,pv2,qv2)。 fprintf(fp1,PV jie dian gong lv:\n)。 fprintf(fp1,dS=%8f+j(%8f)\n,dSp,dSq)。 fprintf(fp1,\n)。 for(i=0。 fprintf(fp1,\n)。 for(i=0。jN。iN。 }iN。 { for(j=0。 printf(sv1=%8f+j(%8f)\n,pv1,qv1)。 { ps=ps+u[0]*u[i]*(g[0][i]*cos(r[0]r[i])+b[0][i]*sin(r[0]r[i]))。iN。 float dSp=0,dSq=0。 float ps=0,pv1=0,pv2=0。 } val(u,g,b,r,ku,kr,h)