基于matlab的電力系統(tǒng)潮流計(jì)算(編輯修改稿)
2025-02-12 10:52 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 i,3))。 Jacbi(2*h,2*k1)=0。 Jacbi(2*h,2*k)=0。 end if k==(n1) k=0。 end end end endendJacbiDetaU=zeros(2*n2,1)。DetaU=inv(Jacbi)*DetaS。DetaU%修正節(jié)點(diǎn)電壓j=0。for i=1:n if B2(i,6)==2 j=j+1。 B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。 endendfor i=1:n if B2(i,6)==3 j=j+1。 B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。 endendB2Times=Times+1。 %迭代次數(shù)加1endTimes一個(gè)原始數(shù)據(jù)的例子節(jié)點(diǎn)數(shù) 5支路數(shù) 5平衡節(jié)點(diǎn)編號(hào) 5精度pr B1(支路參數(shù)矩陣)[1 2 + 1 0。1 3 + 0 1 0。2 3 + 1 0。4 2 0 1。5 3 0 1]B2(節(jié)點(diǎn)參數(shù)矩陣)[0 1 0 0 2。0 21i 1 0 0 2。0 1 0 0 2。0 5+0i 0 3。0 0 0 1]X(節(jié)點(diǎn)號(hào)和對(duì)地參數(shù))[1 0。2 0。3 0。4 0。5 0]電力系統(tǒng)潮流計(jì)算——9結(jié)點(diǎn)算例PQ法:標(biāo)號(hào),起始結(jié)點(diǎn),終止結(jié)點(diǎn),支路電阻參數(shù),支路電抗參數(shù),支路對(duì)地導(dǎo)納參數(shù)1,2,5,2,5,9,3,6,9, 4,3,6,5,6,8,6,4,8, 7,5,7,8,4,7,9,1,4,:includeincludedefine N 9 /*總結(jié)點(diǎn)數(shù)*/define M 6 /*PQ結(jié)點(diǎn)數(shù)*/define K 9 /*線路數(shù)*/define eps 1e4void guass(int n,int m,float c[],float b[][N],float x[]) /*高斯函數(shù)*/{float a[N][N],y[N]。int i,j,k。for(i=0。im。i++){ y[i]=c[i]。 for(j=0。jm。j++) a[i][j]=b[i+n][j+n]。}for(k=0。km1。k++) for(i=k+1。im。i++) { for(j=k+1。jm。j++) a[i][j]=a[i][j]a[k][j]*a[i][k]/a[k][k]。 y[i]=y[i]y[k]*a[i][k]/a[k][k]。 }for(i=m1。i=0。i){ for(j=i+1。jm。j++) y[i]=a[i][j]*x[j]。 x[i]=y[i]/a[i][i]。}}struct line{int Lindex。 int Headnode。 int Endnode。 float R。 float X。 float b。}line[K]。struct line *t。main(){ float r[N]={}。 float u[N]={,,}。 float p[N]={1,,}。 float q[N]={1,1,1,,}。 float g[N][N]={},b[N][N]={}。 float h[N][N]={}。 float B[N][N]={}。 float temp。 float H[K][6]。 float 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])。 /*函數(shù)申明*/ FILE *fp。 fp=fopen(,r)。 for(i=0。iK。i++) { for(j=0。j6。j++) { fscanf(fp,%f,amp。temp)。 H[i][j]=temp。 } } fclose(fp)。 for(i=0。iK。i++) { line[i].Lindex=(int)H[i][0]。line[i].Headnode=(int)H[i][1]。line[i].Endnode=(int)H[i][2]。 line[i].R=H[i][3]。line[i].X=H[i][4]。line[i].b=H[i][5]。} for(t=line。tline+K。t++) {i=tHeadnode1。 j=tEndnode1。 lr=tR。 lx=tX。 lb1=tb。 lg=lr/(lr*lr+lx*lx)。 lb=lx/(lr*lr+lx*lx)。 g[i][i]+=lg。 g[j][j]+=lg。 b[i][i]+=lb+lb1。 b[j][j]+=lb+lb1。 h[i][j]=h[j][i]=lb1。 g[i][j]=lg。 g[j][i]=lg。 b[i][j]=lb。 b[j][i]=lb。 } getch()。 printf(\n=====jie dian dao na ju zhen=====\n)。 for(i=0。iN。i++) { for(j=0。jN。j++) { B[i][j]=b[i][j]。 printf(%8f,B[i][j])。 } printf(\n)。 } printf(\n)。 getch()。 printf(\n=====gei ding chu zhi=====\n)。 for(i=0。iN。i++) printf(u[%d]=%8f p[%d]=%8f q[%d]=%8f\n,i+1,u[i],i+1,p[i],i+1,q[i])。 printf(\n=====die dai qiu jie=====\n)。 while(kp==1) { float ip,iq,max。 float dp[N],dq[N],dr[N]。 float dpu[N],dqu[N]。 float y1[N1],y2 。 float x1[N1],x2 。 for(i=1。iN。i++) /*算dp對(duì)應(yīng)B39。,N1維,除去平衡結(jié)點(diǎn)*/1
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1