石油大學(xué)電力系統(tǒng)分析matlab編程.doc

電 力 系 統(tǒng) 分 析大作業(yè) m=2; %m=PQ節(jié)點(diǎn)個(gè)數(shù)v=1;%v=PV節(jié)點(diǎn)個(gè)數(shù)P=-0.3 -0.55; %P=PQ節(jié)點(diǎn)的P值Q=-0.18 -0.13; %Q=PQ節(jié)點(diǎn)的Q值PP=0.5;%PP=PV節(jié)點(diǎn)的P值V=1.1;%V=PV節(jié)點(diǎn)的U值E=1 1 1.1 1.05; %E=PQ,PV,V節(jié)點(diǎn)e的初值 F=0 0 0 0 ; %F=PQ,PV,V節(jié)點(diǎn)f的初值G= 1.042093 -0.588235 0 -0.453858; -0.588235 1.069005 0 -0.480769; 0 0 0 0; -0.453858 -0.480769 0 0.934619 ; B= -8.24834 2.352941 3.666667 1.891074; 2.352941 -4.727377 0 2.403846; 3.666667 0 -3.333333 0; 1.891074 2.403846 0 -4.26142 ;Y=G+j*B;X=; %X=Xn=m+v+1;%總的節(jié)點(diǎn)數(shù)FX=ones(2*n-2,1);%F(x)矩陣F1=zeros(n-1,n-1);%F(x)導(dǎo)數(shù)矩陣a=0;%記錄迭代次數(shù)EF=zeros(n-1,n-1);%最后的節(jié)點(diǎn)電壓矩陣while max(FX)=10(-5)for i=1:m %PQ節(jié)點(diǎn) FX(2*i-1)=P(i);%P FX(2*i)=Q(i);%Q for w=1:n FX(2*i-1)= FX(2*i-1)-E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)-F(i)*G(i,w)*F(w)-F(i)*B(i,w)*E(w); %P FX(2*i)=FX(2*i)-F(i)*G(i,w)*E(w)+F(i)*B(i,w)*F(w)+E(i)*G(i,w)*F(w)+E(i)*B(i,w)*E(w); %Q endendfor i=m+1:n-1 %PV節(jié)點(diǎn) FX(2*i-1)=PP(i-m);%P FX(2*i)=V(i-m)2-E(i)2-F(i)2; %Q for w=1:n FX(2*i-1)= FX(2*i-1)-E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)-F(i)*G(i,w)*F(w)-F(i)*B(i,w)*E(w);%P endendfor i=1:m %PQ節(jié)點(diǎn) for w=1:n-1 if i=w F1(2*i-1,2*w-1)=-(G(i,w)*E(i)+B(i,w)*F(i); F1(2*i-1,2*w)=B(i,w)*E(i)-G(i,w)*F(i); F1(2*i,2*w-1)=B(i,w)*E(i)-G(i,w)*F(i); F1(2*i,2*w)=G(i,w)*E(i)+B(i,w)*F(i); else F1(2*i-1,2*w-1)=-G(i,i)*E(i)-B(i,i)*F(i); F1(2*i-1,2*w)=B(i,i)*E(i)-G(i,i)*F(i); F1(2*i,2*w-1)=B(i,i)*E(i)-G(i,i)*F(i); F1(2*i,2*w)=G(i,i)*E(i)+B(i,i)*F(i); for k=1:n F1(2*i-1,2*w-1)=F1(2*i-1,2*w-1)-G(i,k)*E(k)+B(i,k)*F(k); F1(2*i-1,2*w)= F1(2*i-1,2*w)-G(i,k)*F(k)-B(i,k)*E(k); F1(2*i,2*w-1)=F1(2*i,2*w-1)+G(i,k)*F(k)+B(i,k)*E(k); F1(2*i,2*w)=F1(2*i,2*w)-G(i,k)*E(k)+B(i,k)*F(k); end end endendfor i=m+1:n-1 %PV節(jié)點(diǎn) for w=1:n-1 if i=w F1(2*i-1,2*w-1)=-(G(i,w)*E(i)+B(i,w)*F(i); F1(2*i-1,2*w)=B(i,w)*E(i)-G(i,w)*F(i); F1(2*i,2*w-1)=0; F1(2*i,2*w)=0; else F1(2*i-1,2*w-1)=-G(i,i)*E(i)-B(i,i)*F(i); F1(2*i-1,2*w)=B(i,i)*E(i)-G(i,i)*F(i); F1(2*i,2*w-1)=-2*E(i); F1(2*i,2*w)=-2*F(i); for k=1:n F1(2*i-1,2*w-1)=F1(2*i-1,2*w-1)-G(i,k)*E(k)+B(i,k)*F(k); F1(2*i-1,2*w)= F1(2*i-1,2*w)-G(i,k)*F(k)-B(i,k)*E(k); end end endendX=inv(F1)*(-FX) ;for i=1:n-1 E(i)=E(i)+X(2*i-1); F(i)=F(i)+X(2*i);enda=a+1;fprintf(第%d次迭代后的節(jié)點(diǎn)電壓分別為:n,a);disp(E+j*F);fprintf(第%d次迭代后功率偏差P Q電壓偏差V的平方分別為:n,a);disp(FX);end disp(收斂后的節(jié)點(diǎn)電壓用極坐標(biāo)表示為:);EF=E+j*F;for i=1:n-1 fprintf(%d號(hào)節(jié)點(diǎn)電壓的幅值為:,i) disp(abs(EF(i); fprintf(%d號(hào)節(jié)點(diǎn)電壓的相角度數(shù)為,i) disp(angle(EF(i)*180/pi); endPPH=0;for i=1:n PPH=PPH+EF(n)*conj(Y(n,i)*conj(EF(i);endfprintf(平衡節(jié)點(diǎn)的功率);disp(PPH);運(yùn)行結(jié)果：第1次迭代后的節(jié)點(diǎn)電壓分別為: 0.9928 - 0.0087i 0.9759 - 0.1077i 1.1000 + 0.1268i 1.0500 第1次迭代后功率偏差P Q電壓偏差V的平方分別為: -0.2773 -0.0564 -0.5260 0.0196 0.5000 0第1次迭代后雅克比矩陣為: -1.0194 -8.3719 0.5882 2.3529 0 3.6667 -8.1248 1.0648 2.3529 -0.5882 3.6667 0 0.5882 2.3529 -1.0450 -4.8770 0 0 2.3529 -0.5882 -4.5778 1.0930 0 0 0 4.0333 0 0 0 -3.6667 0 0 0 0 -2.2000 0第2次迭代后的節(jié)點(diǎn)電壓分別為: 0.9840 - 0.0084i 0.9586 - 0.1083i 1.0924 + 0.1293i 1.0500 第2次迭代后功率偏差P Q電壓偏差V的平方分別為: -0.0015 -0.0028 -0.0137 -0.0547 0.0034 -0.0161第2次迭代后雅克比矩陣為: -0.8069 -8.3606 0.6044 2.3308 0.0318 3.6401 -7.9984 1.4051 2.3308 -0.6044 3.6401 -0.0318 0.8276 2.2330 -1.0182 -4.6347 0 0 2.2330 -0.8276 -4.3623 2.0871 0 0 -0.4651 4.0333 0 0 -0.0318 -3.6401 0 0 0 0 -2.2000 -0.2537第3次迭代后的節(jié)點(diǎn)電壓分別為: 0.9839 - 0.0084i 0.9583 - 0.1083i 1.0924 + 0.1293i 1.0500 第3次迭代后功率偏差P Q電壓偏差V的平方分別為: -0.0000 -0.0000 -0.0003 -0.0011 0.0001 -0.0001第3次迭代后雅克比矩陣為: -0.7913 -8.2929 0.5985 2.3103 0.0308 3.6079 -7.9219 1.3979 2.3103 -0.5985 3.6079 -0.0308 0.8188 2.1918 -0.9856 -4.6126 0 0 2.1918 -0.8188 -4.2190 2.0881 0 0 -0.4739 4.0055 0 0 -0.0308 -3.6079 0 0 0 0 -2.1848 -0.2585第4次迭代后的節(jié)點(diǎn)電壓分別為: 0.9839 - 0.0084i 0.9583 - 0.1083i 1.0924 + 0.1293i 1.0500 第4次迭代后功率偏差P Q電壓偏差V的平方分別為: 1.0e-006 * -0.0072 -0.0078 -0.0808 -0.3859 0.0106 -0.0015第4次迭代后雅克比矩陣為: -0.7912 -8.2920 0.5985 2.3100 0.0308 3.6075 -7.9209 1.3979 2.3100 -0.5985 3.6075 -0.0308 0.8186 2.1910 -0.9850 -4.6123 0 0 2.1910 -0.8186 -4.2163 2.0881 0 0 -0.4740 4.