【正文】 rgS(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)))。B2(i,6)==2B2(i,6)==3 for j=1:n Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。for i=1:n endendB2Times=Times+1。3 0。j++) y[i]=a[i][j]*x[j]。 float h[N][N]={}。 float temp。 float lr,lx,lb1,lg,lb。 int ku=0,kr=0,kp=1,kq=1。 fp=fopen(,r)。 lx=tX。 g[j][j]+=lg。 b[i][j]=lb。 } float x1[N1],x2 。 printf(dp[%d]=%8f\n,i+1,dp[i])。 for (i=1。 for(j=0。i++) else goto top。 pv1=pv1+u[1]*u[i]*(g[1][i]*cos(r[1]r[i])+b[1][i]*sin(r[1]r[i]))。i++) } +u[j]*u[j]*(h[j][i]+b[j][i])+u[j]*u[i]*(g[j][i]*sin(r[j]r[i])b[j][i]*cos(r[j]r[i]))。 printf(\n=====xian lu sun hao=====\n)。 getch()。 printf(ku=%d\n,ku)。jN。 } r[%d]=%8f\n,i+1,u[i],i+1,r[i])。iN。iN。 {j++) { for(j=i+1。i++) { for(j=0。 printf(\n=====xian lu gong lv=======\n)。 getch()。 }i++) float q[N][N]={0}。 if(kp==0) kq=0。 y2[i]=dqu[i+NM]。 for(i=0。 { top: for(i=1。 } for(i=1。 }iN。 printf(\n=====die dai qiu jie=====\n)。 printf(\n=====gei ding chu zhi=====\n)。 printf(\n)。 {line[i].b=H[i][5]。 line[i].Lindex=(int)H[i][0]。j++) float u[N]={,}。0 1 0 0 2。 B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。 end k=0。 if i==j k=k+1。 end Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。h=0。for i=1:nfor i=1:n h=h+1。 endendJacbi%求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量DetaU=zeros(2*n2,1)。 %對角元素的處理 %對PV節(jié)點(diǎn)的處理 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 B2(i,6)==2 I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。 DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。 h=h+1。 endendPVU%創(chuàng)建DetaS，用于存儲有功功率、無功功率和電壓幅值的不平衡量h=0。節(jié)點(diǎn)號和對地參數(shù):X=39。n1=input(39。%第五烈為含變壓器支路的變壓器的變比%第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，%“0”為不含有變壓器。B1=input(39。 endendfor i=1:nt=0。B2(i,6)==3 k=k+1。 if i==j endendk=0。 Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。 %對PQ節(jié)點(diǎn)處理B2(i,6)==2for i=1:n Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。 endim。 for(j=0。 }for(i=m1。tline+K。 lr=tR。 g[i][i]+=lg。 { }iN。 float y1[N1],y2 。jN。 {i++) u[i]=u[i]+x2[i(NM)]。i++) printf(u[%d]=%8f\n,i+1,u[i])。 qv2=qv2+u[2]*u[i]*(g[2][i]*sin(r[2]r[i])b[2][i]*cos(r[2]r[i]))。 for(i=0。j++) { i++) fprintf(fp1,sp=%8f+j(%8f)\n,ps,qs)。 fprintf(fp1,sv1=%8f+j(%8f)\n,pv1,qv1)。 fprintf(fp1,wang luo zong sun hao:\n)。 for(j=i+1。 fprintf(fp1,xian lu cao liu:\n)。 dSq+=dq[i][j]。 printf(\n)。 for(i=0。 printf(\n======die dai hou dian ya yu xiang jiao zhi ======\n)。 +u[j]*u[j]*(g[j][i])+u[j]*u[i]*(g[j][i]*cos(r[j]r[i])+b[j][i]*sin(r[j]r[i]))。i++) pv2=pv2+u[2]*u[i]*(g[2][i]*cos(r[2]r[i])+b[2][i]*sin(r[2]r[i]))。 for(i=0。 float s[N][N]。 else } for(i=NM。 getc