基于matlab的電力系統(tǒng)潮流計(jì)算(存儲(chǔ)版)
【正文】 )+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。 Jacbi(2*h,2*k1)=2*imag(B2(i,3))。 end if B2(i,6)==21 3 + 0 1 0。2 0。j++) a[i][j]=b[i+n][j+n]。i){jm。 float p[N]={1,}。temp)。 } }iK。 {i=tHeadnode1。 lb=lx/(lr*lr+lx*lx)。 for(j=0。 printf(u[%d]=%8f float dp[N],dq[N],dr[N]。 if(fabs(dpu[i])max)i++) { xiang jiao chu zhi =====\n,kr+1)。 getch()。 for(i=NM。 dq[i]=q[i]u[i]*iq。 }iN。 else else } float s[N][N]。 getch()。 for(i=0。iN。 pv2=pv2+u[2]*u[i]*(g[2][i]*cos(r[2]r[i])+b[2][i]*sin(r[2]r[i]))。i++) p[i][j]=u[i]*u[i]*(g[i][j])+u[i]*u[j]*(g[i][j]*cos(r[i]r[j])+b[i][j]*sin(r[i]r[j]))。 +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]))。 printf(\n======die dai hou dian ya yu xiang jiao zhi ======\n)。i++) for(i=0。 getch()。 printf(\n)。 dSq+=dq[i][j]。 fprintf(fp1,xian lu cao liu:\n)。 for(j=0。 fprintf(fp1,zhi lu sun hao:\n)。 for(j=i+1。 fprintf(fp1,wang luo zong sun hao:\n)。 fprintf(fp1,\n)。 fprintf(fp1,sv1=%8f+j(%8f)\n,pv1,qv1)。 fprintf(fp1,\n)。 fprintf(fp1,sp=%8f+j(%8f)\n,ps,qs)。 fprintf(fp1,jie dian dian ya yu xiang jiao:\n)。 getch()。i++)i++)i++)jN。 { }jN。 for(i=0。 } dp[i][j]=u[i]*u[i]*(g[i][j])+u[i]*u[j]*(g[i][j]*cos(r[i]r[j])+b[i][j]*sin(r[i]r[j]))j++) for(i=0。 qv2=qv2+u[2]*u[i]*(g[2][i]*sin(r[2]r[i])b[2][i]*cos(r[2]r[i]))。 float dq[N][N]={0}。 { printf(u[%d]=%8f\n,i+1,u[i])。i++) u[i]=u[i]+x2[i(NM)]。i++)iM。 dqu[i]=dq[i]/u[i]。 {iN。 r[i]=r[i]+dr[i]。 dr[i]=x1[i1]/u[i]。 for(i=0。iN。jN。i++) /*算dp對(duì)應(yīng)B39。 float y1[N1],y2 。 q[%d]=%8f\n,i+1,u[i],i+1,p[i],i+1,q[i])。iN。 } { g[i][i]+=lg。 lr=tR。tline+K。}line[i].Headnode=(int)H[i][1]。 }for(i=m1。 y[i]=y[i]y[k]*a[i][k]/a[k][k]。 for(j=0。 y[i]=c[i]。im。0 5+0i 0 3。 if B2(i,6)==3 if j~=isb end Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*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))。k=0。for i=1:nB2(i,6)==2 %對(duì)PQ節(jié)點(diǎn)處理DetaU=inv(Jacbi)*DetaS。 Jacbi(2*h,2*k1)=0。 Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。 endendk=0。 if k==(n1) %將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行 if i==j k=k+1。 if i~=isbB2(i,6)==3for i=1:nt=0。 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)))。 endendfor i=1:n endendYOrgS=zeros(2*n2,1)。)。B1=input(39。請(qǐng)輸入支路數(shù):n1=39。%第五烈為含變壓器支路的變壓器的變比%第六列為變壓器是否是否含有變壓器的參數(shù)，其中“1”為含有變壓器，%“0”為不含有變壓器。%第四列為支路的對(duì)地導(dǎo)納參數(shù)。n1=input(39。)。節(jié)點(diǎn)號(hào)和對(duì)地參數(shù):X=39。 endendOrgS%創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓PVU=zeros(nh1,1)。 endendPVU%創(chuàng)建DetaS，用于存儲(chǔ)有功功率、無(wú)功功率和電壓幅值的不平衡量h=0。 if i~=isbamp。 h=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。 h=h+1。 I(h,1)=(OrgS(2*h1,
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1