基于matlab的電力系統(tǒng)潮流計(jì)算-免費(fèi)閱讀
【正文】 } fp1=fopen(,w+)。 printf(kr=%d\n,kr)。 printf(dS=%8f+j(%8f)\n,dSp,dSq)。 for(j=i+1。 for(i=0。 } for(i=0。 getch()。 printf(sv2=%8f+j(%8f)\n,pv2,qv2)。 { for(i=0。 float dp[N][N]={0}。 kp=0。 if(fabs(dqu[i])max) max=fabs(dqu[NM])。j++) iq=0。kq=1。 printf(r[%d]=%8f\n,i+1,r[i])。i++) { dpu[i]=dp[i]/u[i]。 float dpu[N],dqu[N]。 { while(kp==1) printf(\n)。 printf(%8f,B[i][j])。iN。 } j=tEndnode1。 line[i].R=H[i][3]。 for(j=0。i++)}line[K]。 for(j=k+1。 {im。k++) /*高斯函數(shù)*/{float a[N][N],y[N]。 /*線路數(shù)*/define eps 1e4void guass(int n,int m,float c[],float b[][N],float x[]) /*PQ結(jié)點(diǎn)數(shù)*/define K 9 /*總結(jié)點(diǎn)數(shù)*/define M 65 0]電力系統(tǒng)潮流計(jì)算——9結(jié)點(diǎn)算例PQ法:標(biāo)號,起始結(jié)點(diǎn),終止結(jié)點(diǎn),支路電阻參數(shù),支路電抗參數(shù),支路對地導(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 95 3 0 1]B2(節(jié)點(diǎn)參數(shù)矩陣)[0 1 0 0 2。 endendfor i=1:nDetaU=inv(Jacbi)*DetaS。 end Jacbi(2*h,2*k)=0。 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。 end k=0。 if i==j k=k+1。B2(i,6)==3 end %非對角元素的處理 end Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。 Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。k=0。for i=1:n DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。 DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。 h=h+1。 if i~=isbamp。 for j=1:n h=h+1。 if i~=isbamp。 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(j,3)))real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。B2(i,6)==2j=0。 %含有變壓器的支路 if B1(i,6)==0)。pr=input(39。請輸入節(jié)點(diǎn)數(shù):n=39。節(jié)點(diǎn)編號由小到大編寫%對于含有變壓器的支路，第一列為低壓側(cè)節(jié)點(diǎn)編號，第二列為高壓側(cè)節(jié)點(diǎn)%編號，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理。%X為節(jié)點(diǎn)號和對地參數(shù)矩陣。isb=input(39。)。Times=1。 %置迭代次數(shù)為初始值%創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣for i=1:n1 p=B1(i,1)。 q=B1(i,2)。 Y(p,q)=Y(p,q)1/B1(i,3)。 Y(q,p)=Y(p,q)。 Y(p,p)=Y(p,p)+1/B1(i,3)+*B1(i,4)。 Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。 %對PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0 %對PQ節(jié)點(diǎn)的處理for i=1:n %非對角元素的處理 if B2(i,6)==3 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)=2*imag(B2(i,3))。 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)。 endendB2%開始循環(huán)**********************************************************************while abs(max(DetaU))prOrgS=zeros(2*n2,1)。 %!!!初始功率參數(shù)在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣h=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)))。 if i~=isbamp。 h=h+1。 DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。 DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。 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,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。 if k==(n1) endendk=0。 Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3)
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1