资源描述：
《潮流计算源程序及运行结果》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、»%木程序的功能是用牛顿——拉夫逊法进行潮流计算%B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳%5、支路的变比；6、支路首端处于K侧为1，1侧为0%B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电乐初始位%4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容fi%6、甘点分类标号：1为平衡节点（应为1号节点）；2为PQ节点;%3为PV节点；clear;n=10;%input（’请输入节点数:n-）;nl=10;%inputf请输入支路数:nl-）;isb=l;%inp＜请输入平衡母线节点号:isb=’）;pr=0.00001;%i叩ut
2、（•请输入误差精度:pr='）;Bl=[l2230.03512+0.08306i0.13455i0.0068+0.18375i011.023810；1;140.05620+0.13289i0.05382i10；450.008ll+0.24549i01.023811;160.05620+0.13289i0.05382i10；460.04215+0.09967i0.04037i10；670.0068+0.18375i01.023811;680.02810+0.06645i0.10764i10；8100.0081l+0.24549i011;890.03512+0.08306i0.134
3、55i10]B2=[001.11.101;00100200.343+0.21256i1002；00100200.204+0.12638i1002；00100200.306+0.18962i1002；0010020.501.11.103；00.343+0.21256i1002];%hiputC请输入各节点参数形成的矩阵：B2=');Y=zcros(n);c=zcros(1，n);f=zeros(l,n);V=zcros(l，n);sida=zcros(1，n);Sl=zcros(nl);%%%fori=l:nl%支路数ifBl(i，6)==0%左节点处于1侧p=Bl(i,l);q=
4、Bl(i,2);else%左节点处于K侧p=Bl(i,2);q=Bl(i,l);Y(p,q)=Y(p,q)-l./(Bl(i,3)*Bl(i,5));%非对角元Y(q，p)=Y(p，q);%非对角元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)A2)+B1(i,4)./2;%对角元K侧Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;%对角元1侧end%求导纳矩阵dispC导纳矩阵Y=.);disp(Y)%分解出导纳阵的实部和虚部%给定各节点初始电压的实部和虚部%PV节点电压给定模位%给定各节点注入功率%i节点注入功率SG-SL%i节点无
5、功补偿量%G=real(Y);B=imag(Y);fori=l:ne(i)=real(B2(i,3))；f(i)=imag(B2(i,3))；V(i)=B2(i,4);endfori=l:nS(i)=B2(i,l)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);end%===========—=—================—===========P=rcal(S);Q=imag(S);%分解出各节点注入的冇功和无功功率ICTl=0;IT2=l;N0=2*n;N=N0+l;a=0;%迭代次数ICT1、a;不满足收敛要求的•点数IT2whileIT2〜=0%N0=
6、2*n雅可比矩阵的阶数；N=N0+l扩展列IT2=0;a=a+l;fori=l:nif一isb%非平衡节点C(i)=0;D(i)=0;forjl=l:nD(i)=D(i)+G(ij1)*f(jl)+B(i,j1)*e(jl);%E(Gij*fj+Bij*ej)endPl=C(i)*e(i)+f(i)*D(i);%节点功率P计算ei2(Gij*ej-Bij*fj)+fiE(Gij*fj+Bij*ej)Ql=C(i)*f(i)-e(i)*D(i);%节点功率Q计算fi2(Gij*ej-Bij*fj)-eiE(Gij*fj+Bij*ej)%求i节点有功和无功功率P’，Q
7、的计算值V2=
8、e(i)A2+f(i)A2;%电压模平方%=====以下针对非PV节点来求取功率差及Jacobi矩阵元素=====ifB2(i,6)^30/dh:PV节点DP=P(i)-Pl;%节点有功功率差DQ=Q(i)-Ql;%节点无功功率差%=========以上为除平衡节点外K它节点的功率计算==========%==========求取Jacobi矩阵============forjl=l:nifjl〜=isb&jl^=i%非T衡节点&非对角元Xl=-G(i,jl)*e(i)-B(i,
