潮流计算的计算机算法.pdfVIP

饭疏食，饮水，曲肱而枕之，乐亦在其中矣。不义而富且贵，于我如浮云。——《论语》

高等电力系统分析

（潮流计算的计算机算法）

PQ分解法潮流计算（IEEE14）

饭疏食，饮水，曲肱而枕之，乐亦在其中矣。不义而富且贵，于我如浮云。——《论语》

目录

一、MATLAB源程序

二、对支路参数（B1）、节点参数（B2）的说明

三、带入数据，运行结果

饭疏食，饮水，曲肱而枕之，乐亦在其中矣。不义而富且贵，于我如浮云。——《论语》

MATLAB源程序

clear

closeall

n=input(请输入节点数:n=);

：

n仁input(请输入支路数n1=)

)

isb=input(请输入平衡节点号:isb=)

pr=input(请输入误差精度:pr=);

：

B1=input(请输入支路参数B1=);

B2=input(请输入节点参数:B2=);

:n

n2=input(请输入PQ节点个数：

2=);

Y=zeros(n);

fori=1:n1

P=B1(i,1)；

q=B1(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;end

disp(导纳矩阵Y=);

disp(Y)

%

%P,Q,V,0

下面是求矩阵

V=zeros(1,n);O=zeros(1,n);P=zeros(1,n);Q=zeros(1,n);G=real(Y);B=imag(Y);

fori=1:n

P(i)=B2(i,3);