原创力文档%本程序的功能是用P-Q分解法进行潮流计算
n=input(请输入节点数:n=);
nl=input(请输入支路数:nl=);
isb=input(请输入平衡母线节点号:isb=);
pr=input(请输入误差精度:pr=);
B1=input(请输入由支路参数形成的矩阵:B1=);
B2=input(请输入由节点参数形成的矩阵:B2=);
na=input(请输入PQ节点数na=);
Y=zeros(n);YI=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);
for i=1:nl
if B1(i,6)==0
p=B1(i,1); q=B1(i,2);
else p=B1(i,2); q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
YI(p,q)=YI(p,q)-1./B1(i,3);
Y(q,p)=Y(p,q);
YI(q,p)=YI(p,q);
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
YI(q,q)=YI(q,q)+1./B1(i,3);
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
Y
文档评论（0）