《潮流计算的计算机算法》由会员分享，可在线阅读，更多相关《潮流计算的计算机算法（17页珍藏版）》请在金锄头文库上搜索。
1、-高等电力系统分析 潮流计算的计算机算法PQ分解法潮流计算IEEE14)目录一、 MATLAB源程序二、 对支路参数B1、节点参数B2的说明三、 带入数据，运行结果一、MATLAB源程序clearclose alln=input(请输入节点数:n=);n1=input(请输入支路数:n1=);isb=input(请输入平衡节点号:isb=);pr=input(请输入误差精度:pr=);B1=input(请输入支路参数:B1=);B2=input(请输入节点参数:B2=);n2=input(请输入PQ节点个数:n2=);Y=zeros(n);for i=1:n1 p=B1(i,1); q=B1(
2、i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)+B1(i,4)*1j); %非对角元 Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j; %对角元 Y(q,q)=Y(q,q)+1/(B1(i,3)+B1(i,4)*1j)+B1(i,6)*1j;enddisp(导纳矩阵Y=);disp(Y)%-%-下面是求P,Q,V,O矩阵-V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,n);G=real(Y);B=imag(Y);for i=1:n P(i)=B2(i,
3、3); Q(i)=B2(i,4); V(i)=B2(i,5); O(i)=B2(i,6); end B3=B(1:n-1,1:n-1); %不含平衡节点，由节点导纳虚部构成B4=B(1:n2,1:n2); %所有PQ节点%-%-下面是求P,Q矩阵-D*=0;ICT=1;Mp=1;Mq=1; while ICT=0 m1=1;m2=1; for i=1:n if i=isb C(i)=0; D(i)=0; for j1=1:n C(i)=C(i)+V(i)*V(j1)*(G(i,j1)*cos(O(i)-O(j1)+B(i,j1)*sin(O(i)-O(j1); D(i)=D(i)+V(i)*V
4、(j1)*(G(i,j1)*sin(O(i)-O(j1)-B(i,j1)*cos(O(i)-O(j1); end DP(m1)=P(i)-C(i); m1=m1+1; if B2(i,2)=1 DQ(m2)=Q(i)-D(i); m2=m2+1; end end end m1=m1-1; %所有节点数 m2=m2-1; %PQ节点数 DPQ=DP;DQ; %求DP,DQ V1=V(:,1:m1); V2=diag(V1); V3=inv(V2); %对V矩阵求逆 H=V3*DP; %P/V K=-inv(B3)*H; %-P/V/B3 deltO=V3*K; %角=-P/V/V/B3ma*1=
5、ma*(abs(DP); for i=1:m1 if ma*1pr Mp=0; else O(i)=O(i)+deltO(i); Mq=1; end end V4=V(:,1:m2); V5=diag(V4); V6=inv(V5); L=V6*DQ; N=-inv(B4)*L; deltV=N; %V=-Q/V/B ma*2=ma*(abs(DQ); for i=1:m2 if ma*2pr Mq=0; else if B2(i,2)=1; V(i)=V(i)+deltV(i); Mp=1; end end end if Mp=0&Mq=0 ICT=0; else ICT=1; end D*
6、=D*+1;end%-%-迭代完毕，开场输出结果-disp(-);disp(迭代次数为：);disp(D*);for i=1:n E(i)=V(i)*cos(O(i)+1j*V(i)*sin(O(i); o(i)= 180*angle(E(i)/pi;enddisp(-);disp(修正后各节点电压标么值为节点号从小到大排列:);disp(V);disp(-);disp(修正后各节点电压相角为节点号从小到大排列:);disp(o);%-计算各个节点的功率-disp(-);disp(各节点的功率为:); for p=1:n C(p)=0; for q=1:n C(p)=C(p)+conj(Y(p
7、,q)*conj(E(q); end S(p)=E(p)*C(p);enddisp(S);%-计算各支路的功率-for i=1:n1 p=B1(i,1);q=B1(i,2); Si(p,q)=E(p)*(conj(E(p)*conj(Y(p,p)-Y(p,q)+(conj(E(p)-conj(E(q)*conj(Y(p,q);disp(-);disp(各条支路的首端功率为:); disp(Si(p,q);Si(q,p)=E(q)*(conj(E(q)*conj(Y(q,q)-Y(p,q)+(conj(E(q)-conj(E(p)*conj(Y(p,q);disp(-);disp(各条支路的末端
8、功率为:); disp(Si(q,p);DS(i)=Si(p,q)+Si(q,p);disp(-);disp(各条支路的功率损耗为:); disp(DS(i);end%-计算平衡节点功率-Sp=0;for i=1:n Sp=Sp+V(n)*conj(Y(n,i)*conj(V(i);enddisp(-);disp(平衡节点功率为:);disp(Sp);二、 对支路参数B1、节点参数B2的说明1.节点数：142.支路数：203.支路矩阵B1的各支路参数：起点编号，终点编号，电阻，电抗，电导，电纳120.013350.0421100;1300.2091200;1400.5561800;1100.0
9、58110.1763200.034;1110.067010.1710300.0128;2100.056950.1738800.0346;21200.2520200;2140.054030.2230400.0492;3400.1100100;31300.1761500;450.031810.084500;490.127110.2703800;560.082050.1920700;6120.094980.198900;780.220920.1998800;7120.122910.2558100;890.170930.3480200;8120.066150.1302700;10110.046990.
10、1979700.0438;10140.019380.0591700.0528;4. 节点参数矩阵B2的各节点参数：对应的每一列为节点编号，类型，注入有功，注入无功，电压幅值，电压相位其中节点类型：1=PQ节点，2=PV节点，0=平衡节点11-0.4780.03910;21-0.076-0.01610;310010;41-0.295-0.16610;51-0.09-0.05810;61-0.035-0.01810;71-0.061-0.01610;81-0.135-0.05810;91-0.149-0.0510;1020.18301.0450;112-0.94201.010;122-0.1120.0471.70;13200.1741.90;1400
