基于matlab的電力系統(tǒng)潮流計算-文庫吧在線文庫
【正文】 1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。 if B2(i,6)==2 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*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。 %對PV節(jié)點的處理 %對角元素的處理 %將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行 endendJacbi%求解修正方程，獲取節(jié)點電壓的不平衡量DetaU=zeros(2*n2,1)。 if B2(i,6)==3 if i~=isbamp。 h=h+1。 for j=1:n 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)))。for i=1:nfor i=1:nh=0。h=0。 Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。 Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。 end k=k+1。 if i==j k=0。 endfor i=1:n j=j+1。 B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。0 1 0 0 2。for(i=0。j++) a[i][j]=a[i][j]a[k][j]*a[i][k]/a[k][k]。main(){ float u[N]={,}。j++) H[i][j]=temp。 line[i].Lindex=(int)H[i][0]。line[i].b=H[i][5]。 for(t=line。 { printf(\n)。 getch()。 printf(\n=====gei ding chu zhi=====\n)。 for(i=0。 printf(\n=====die dai qiu jie=====\n)。 float ip,iq,max。iN。 for(j=0。 dp[i]=p[i]u[i]*ip。 } for(i=1。 max=fabs(dpu[i])。 } for(i=1。 printf(\n\n)。 top: { for(i=0。 y2[i]=dqu[i+NM]。 /*同上,對應(yīng)加NM*/ guass(NM,M,y2,B,x2)。iN。 printf(\n=====di %d ci die dai dian ya chu zhi=====\n,ku+1)。iN。 printf(\n\n)。kp=1。 kq=0。 if(kp==0) { float q[N][N]={0}。 printf(\n=====ping heng jie dian gong lv =====\n)。i++) } getch()。jN。 printf(\n=====xian lu gong lv=======\n)。 { for(j=0。 }i++) { for(j=i+1。 printf(\n)。 printf(\n)。i++)j++) {iN。iN。iN。 r[%d]=%8f\n,i+1,u[i],i+1,r[i])。 fprintf(fp1,u[%d]=%8f }jN。 }jN。 printf(\n=======shu ju bao cun=====\n)。 printf(ku=%d\n,ku)。 { dSp+=dp[i][j]。 getch()。 getch()。 printf(\n=====xian lu sun hao=====\n)。 printf(\n)。iN。 +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]))。 }i++) pv1=pv1+u[1]*u[i]*(g[1][i]*cos(r[1]r[i])+b[1][i]*sin(r[1]r[i]))。iN。 float ds[N][N]。 goto top。 else max=fabs(dqu[i])。i++) for(j=0。iN。 printf(\n===== di %d ci die dai hou dian ya for (i=1。 printf(dp[%d]=%8f\n,i+1,dp[i])。 float x1[N1],x2 。 }jN。 b[i][j]=lb。 h[i][j]=h[j][i]=lb1。 g[j][j]+=lg。 lx=tX。i++) for(i=0。 fp=fopen(,r)。 FILE *fp。 void val(float u[N],float g[N][N],float b[N][N],float r[N],int ku, int kr,float h[N][N])。 int ku=0,kr=0,kp=1,kq=1。 int i,j。 float lr,lx,lb1,lg,lb。 float H[K][6]。 float temp。 float B[N][N]={}。 float h[N][N]={}。 x[i]=y[i]/a[i][i]。j++) y[i]=a[i][j]*x[j]。}for(k=0。3 0。2 3 + 1 0。 endendB2Times=Times+1。for i=1:n end Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。 Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。 for j=1:n h=h+1。 if i~=isbB2(i,6)==3B2(i,6)==2 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(
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1