基于matlab的電力系統(tǒng)潮流計算(參考版)
2025-01-19 10:52本頁面
【正文】 getch()。 printf(\n===========THE END==============\n)。i++) for(i=0。 fprintf(fp1,jie dian dian ya yu xiang jiao:\n)。 fprintf(fp1,\n)。 fprintf(fp1,sp=%8f+j(%8f)\n,ps,qs)。 fprintf(fp1,ping heng jie dian gong lv:\n)。 fprintf(fp1,\n)。 fprintf(fp1,sv2=%8f+j(%8f)\n,pv2,qv2)。 fprintf(fp1,sv1=%8f+j(%8f)\n,pv1,qv1)。 fprintf(fp1,PV jie dian gong lv:\n)。 fprintf(fp1,\n)。 fprintf(fp1,dS=%8f+j(%8f)\n,dSp,dSq)。 fprintf(fp1,wang luo zong sun hao:\n)。 fprintf(fp1,\n)。 fprintf(fp1,ds[%d][%d]=%8f+j%8f\n,i+1,j+1,dp[i][j],dq[i][j])。j++) for(j=i+1。 { i++) for(i=0。 fprintf(fp1,zhi lu sun hao:\n)。 fprintf(fp1,\n)。 fprintf(fp1,s[%d][%d]=%8f+j%8f\n,i+1,j+1,p[i][j],q[i][j])。j++) for(j=0。 { i++) for(i=0。 fprintf(fp1,xian lu cao liu:\n)。 } dSq+=dq[i][j]。jN。iN。 }jN。 { iN。 } printf(\n)。 }jN。 getch()。 { i++) for(i=0。 r[%d]=%8f\n,i+1,u[i],i+1,r[i])。 printf(u[%d]=%8fi++) 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]))。 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])) 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]))。 {j++) for(j=0。i++) 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]))。 pv2=pv2+u[2]*u[i]*(g[2][i]*cos(r[2]r[i])+b[2][i]*sin(r[2]r[i]))。iN。 printf(sv1=%8f+j(%8f)\n,pv1,qv1)。 } {i++) for(i=0。 printf(\n=====PV jie dian gong lv=====\n)。 qs=qs+u[0]*u[i]*(g[0][i]*sin(r[0]r[i])b[0][i]*cos(r[0]r[i]))。 ps=ps+u[0]*u[i]*(g[0][i]*cos(r[0]r[i])+b[0][i]*sin(r[0]r[i]))。iN。 getch()。 FILE *fp1。 float dSp=0,dSq=0。 float dq[N][N]={0}。 float s[N][N]。 float p[N][N]={0}。 float ps=0,pv1=0,pv2=0。 { }}void val(float u[N],float g[N][N],float b[N][N],float r[N],int ku, int kr,float h[N][N]) } val(u,g,b,r,ku,kr,h)。 if(kq==0) else } val(u,g,b,r,ku,kr,h)。 { else } ku=ku+1。 printf(u[%d]=%8f\n,i+1,u[i])。i++) for(i=NM。 u[i]=u[i]+x2[i(NM)]。i++) for(i=NM。iM。 if(max=eps)iN。 } printf(\n)。 printf(dq[%d]=%8f\n,i+1,dq[i])。 dqu[i]=dq[i]/u[i]。 dq[i]=q[i]u[i]*iq。 iq=iq+u[j]*(g[i][j]*sin(r[i]r[j])b[i][j]*cos(r[i]r[j]))。jN。 {i++) /*算dq對應(yīng)B,僅M維,除去平衡結(jié)點和PV結(jié)點*/ for(i=NM。 kr=kr+1。 getch()。iN。 xiang jiao chu zhi =====\n,kr+1)。 r[i]=r[i]+dr[i]。 dr[i]=x1[i1]/u[i]。 {iN。 guass(1,N1,y1,B,x1)。 /*起值不同,為了對應(yīng),故加一*/ y1[i]=dpu[i+1]。i++) for(i=0。 { if (max=eps) } if(fabs(dpu[i])max)iN。 max=fabs(dpu[1])。jN。 {i++) /*算dp對應(yīng)B39。 for(i=1。 float y1[N1],y2 。 float dp[N],dq[N],dr[N]。 q[%d]=%8f\n,i+1,u[i],i+1,p[i],i+1,q[i])。 printf(u[%d]=%8fiN。 } B[i][j]=b[i][j]。 {j++) for(j=0。i++) for(i=0。 printf(\n=====jie dian dao na ju zhen=====\n)。 getch()。 b[j][i]=lb。 g[j][i]=lg。 b[j][j]+=lb+lb1。 g[i][i]+=lg。 lb=lx/(lr*lr+lx*lx)。 l
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1