基于matlab的電力系統(tǒng)潮流計(jì)算-全文預(yù)覽
【正文】 ode1。line[i].X=H[i][4]。iK。 } }temp)。j6。 float q[N]={1,1,1,}。 float p[N]={1,}。 float r[N]={}。struct line *t。jm。i){jm。i++)j++) a[i][j]=b[i+n][j+n]。int i,j,k。2 0。0 21i 1 0 0 2。1 3 + 0 1 0。 if B2(i,6)==2DetaU%修正節(jié)點(diǎn)電壓j=0。 end else Jacbi(2*h,2*k)=2*real(B2(i,3))。 Jacbi(2*h,2*k1)=2*imag(B2(i,3))。 Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。 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))。 if B2(i,6)==3 endendk=0。 if k==(n1) if j~=isb endendI%創(chuàng)建JacbiJacbi=zeros(2*n2)。 I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。 h=h+1。 endendDetaS%創(chuàng)建Ii=zeros(n1,1)。 DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。 DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。 t=t+1。 h=h+1。 if i~=isbamp。 endendt=0。 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。 endendOrgS%創(chuàng)建DetaSh=0。 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)))。 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)))。 %!!!初始功率參數(shù)在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣h=0。 endendB2%開始循環(huán)**********************************************************************while abs(max(DetaU))prOrgS=zeros(2*n2,1)。 B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。 j=j+1。 endendfor i=1:n B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。 j=j+1。for i=1:n end k=0。 else Jacbi(2*h,2*k)=2*real(B2(i,3))。 Jacbi(2*h,2*k1)=2*imag(B2(i,3))。 Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。 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))。 if B2(i,6)==3 end %非對角元素的處理 for j=1:n h=h+1。for i=1:nfor i=1:nB2(i,6)==2 %對PQ節(jié)點(diǎn)的處理B2(i,6)==3 %對PV節(jié)點(diǎn)的處理，注意這時(shí)不可再將h初始化為0 for j=1:n h=h+1。for i=1:n else Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。 Y(p,p)=Y(p,p)+1/B1(i,3)+*B1(i,4)。 Y(q,p)=Y(p,q)。 Y(p,q)=Y(p,q)1/B1(i,3)。 q=B1(i,2)。 p=B1(i,1)。 %置迭代次數(shù)為初始值%創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣for i=1:n1Times=1。X=input(39。)。請輸入誤差精度:pr=39。isb=input(39。)。%X為節(jié)點(diǎn)號和對地參數(shù)矩陣。%第三列為支路的串列阻抗參數(shù)。節(jié)點(diǎn)編號由小到大編寫%對于含有變壓器的支路，第一列為低壓側(cè)節(jié)點(diǎn)編號，第二列為高壓側(cè)節(jié)點(diǎn)%編號，將變壓器的串聯(lián)阻抗置于低壓側(cè)處理。%B2為節(jié)點(diǎn)參數(shù)矩陣，其中第一列為節(jié)點(diǎn)注入發(fā)電功率參數(shù)；第二列為節(jié)點(diǎn)%負(fù)荷功率參數(shù)；第三列為節(jié)點(diǎn)電壓參數(shù)；第六列為節(jié)點(diǎn)類型參數(shù)，其中%“1”為平衡節(jié)點(diǎn)，“2”為PQ節(jié)點(diǎn)，“3”為PV節(jié)點(diǎn)參數(shù)。請輸入節(jié)點(diǎn)數(shù):n=39。)。pr=input(39。請輸入支路參數(shù):B1=39。)。Y=zeros(n)。 if B1(i,6)==0 %含有變壓器的支路DetaS=zeros(2*n2,1)。j=0。 %對PQ節(jié)點(diǎn)的處理B2(i,6)==2 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)))。 OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1