daona.m

电力系统中导纳矩阵的形成,用于计算潮流,比较重要

n=input('\n请输入节点数:n=');
m=input('\n请输入支路数：m=');
ph=input('\n请输入平衡母线的节点号:ph=');
B1=input('\n请输入支路信息：B1=');
%它以矩阵形式存贮支路的情况，每行存贮一条支路
%第一列存贮支路的一个端点
%第二列存贮支路的另一个端点
%第三列存贮支路的阻抗
%第四列存贮支路的对地导纳
%第五列存贮变压器的变比，注意支路为1
B2=input('请输入节点信息：B2=');
%第一列为电源侧的功率
%第二列为负荷侧的功率
%第三列为该点的电压值
%第四列为该点的类型：1为PQ节点,2为PV节点,3为平衡节点
A=input('请输入节点号及对地阻抗：A=');
Y=zeros(n);
e=zeros(1,n);
f=zeros(1,n);
for i=1:n
    if A(i,2)~=0
        p=A(i,1);
        Y(p,p)=1./A(i,2);
    end
end
for i=1:m
    p=B1(i,1);
    q=B1(i,2);
    Y(p,p)=Y(p,p)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
    Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
    Y(q,p)=Y(p,q);
    Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2;
end