基于matlab的電力系統(tǒng)潮流計算-閱讀頁
2025-01-31 10:52本頁面
【正文】 at lr,lx,lb1,lg,lb。 int i,j。 int ku=0,kr=0,kp=1,kq=1。 void val(float u[N],float g[N][N],float b[N][N],float r[N],int ku, int kr,float h[N][N])。 FILE *fp。 fp=fopen(,r)。 for(i=0。i++) for(j=0。j++) { fscanf(fp,%f,amp。 H[i][j]=temp。 fclose(fp)。 for(i=0。i++) line[i].Lindex=(int)H[i][0]。line[i].Endnode=(int)H[i][2]。 line[i].R=H[i][3]。line[i].b=H[i][5]。 for(t=line。t++) j=tEndnode1。 lx=tX。 g[j][j]+=lg。 b[i][i]+=lb+lb1。 h[i][j]=h[j][i]=lb1。 g[i][j]=lg。 b[i][j]=lb。 }iN。 {jN。 printf(%8f,B[i][j])。 } printf(\n)。 printf(\n)。 getch()。 printf(\n=====gei ding chu zhi=====\n)。 for(i=0。i++) p[%d]=%8f printf(\n=====die dai qiu jie=====\n)。 while(kp==1) { float ip,iq,max。 float dpu[N],dqu[N]。 float x1[N1],x2 。iN。,N1維,除去平衡結(jié)點*/ ip=0。 for(j=0。j++) ip=ip+u[j]*(g[i][j]*cos(r[i]r[j])+b[i][j]*sin(r[i]r[j]))。 dp[i]=p[i]u[i]*ip。 dpu[i]=dp[i]/u[i]。 printf(dp[%d]=%8f\n,i+1,dp[i])。 printf(\n)。 } getch()。 for(i=1。i++) { max=fabs(dpu[i])。iN1。 for (i=1。i++) } printf(\n===== di %d ci die dai hou dian ya for(i=1。i++) printf(r[%d]=%8f\n,i+1,r[i])。 printf(\n\n)。kq=1。 top:iN。 iq=0。 for(j=0。j++) max=fabs(dqu[NM])。 for (i=NM。i++) { if(fabs(dqu[i])max) max=fabs(dqu[i])。 } { for(i=0。i++) 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) } { kp=0。 else goto top。 float qs=0,qv1=0,qv2=0。 float q[N][N]={0}。 float dp[N][N]={0}。 float ds[N][N]。 int i,j。 printf(\n=====ping heng jie dian gong lv =====\n)。 for(i=0。i++) { } printf(sp=%8f+j(%8f)\n,ps,qs)。 getch()。iN。 pv1=pv1+u[1]*u[i]*(g[1][i]*cos(r[1]r[i])+b[1][i]*sin(r[1]r[i]))。 qv1=qv1+u[1]*u[i]*(g[1][i]*sin(r[1]r[i])b[1][i]*cos(r[1]r[i]))。 for(i=0。i++) { } printf(sv2=%8f+j(%8f)\n,pv2,qv2)。iN。jN。 q[i][j]=u[i]*u[i]*(h[i][j]+b[i][j])+u[i]*u[j]*(g[i][j]*sin(r[i]r[j])b[i][j]*cos(r[i]r[j]))。 dq[i][j]=u[i]*u[i]*(h[i][j]+b[i][j])+u[i]*u[j]*(g[i][j]*sin(r[i]r[j])b[i][j]*cos(r[i]r[j])) +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]))。 getch()。iN。 printf(\n=====xian lu gong lv=======\n)。iN。 { for(j=0。j++) { printf(s[%d][%d]=%8f+j(%8f)\n,i+1,j+1,p[i][j],q[i][j])。 printf(\n)。 } printf(\n)。 printf(\n=====xian lu sun hao=====\n)。 for(i=0。i++) getch()。 { for(j=i+1。j++) { printf(ds[%d][%d]=%8f+j%8f\n,i+1,j+1,dp[i][j],dq[i][j])。 printf(\n)。 } printf(\n)。 } printf(\n)。 printf(\n=====wang luo zong sun hao=====\n)。 getch()。 for(i=0。i++) { for(j=i+1。j++) { dSp+=dp[i][j]。 } printf(dS=%8f+j(%8f)\n,dSp,dSq)。 printf(\n======die dai ci shu========\n)。 printf(ku=%d\n,ku)。 printf(kr=%d\n,kr)。 printf(\n=======shu ju bao cun=====\n)。 fp1=fopen(,w+)。 {iN。jN。 }iN。jN。 }iN。 fprintf(fp1,u[%d]=%8f r[%d]=%8f\n,i+1,u[i],i+1,r[i])。 }
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1