这学期在上大四专业选修课时，一个课后作业是写一个电力系统的导纳矩阵计算程序，我找了找资料肝了半天的程序如下图所示，希望能对大家有用。
1.具有非标准变比的变压器处理 如果变压器阻抗归算到非标准变比的K侧，则在参数输入脚本文件data中变压器参数第五列记为0；如果归算到非标准变比的1侧，则在参数输入脚本文件data中变压器参数第五列记为1。
2.线路为双回线时，需要将两个节点之间线路阻抗减半，导纳加倍；同时，如果线路是双回线则在data线路参数中第五列记为1，否则记为0。
程序data脚本用于输入参数，buildYbus是计算导纳矩阵函数 data.m
%-----Input Data-------输入参数脚本 %线路参数：第一、二列为母线节点编号 第三列为线路阻抗 第四列为线路对地导纳(Y) %---------第五列（0——单回线，1——双回线） lineData=[1 2 0.04+0.25i 0.5i 0; 1 3 0.1+0.35i 0 0; 2 3 0.08+0.3i 0.5i 0;]; %变压器参数：第一、二列为母线节点编号 第三列变压器阻抗 第四列变压器非标准变比(1:K) %---------第五列（0——变压器阻抗归算到K侧即第二列母线节点编号，1——归算到1侧即第一列节点编号） transData=[4 2 0.015i 1.05 1; 5 3 0.03i 1.05 1;];主程序buildYbus.m
function Ybus=buildYus() %调用data输入参数 data; %变压器π型等效导纳参数 trans=zeros(size(transData ,1),5); trans(:,1:2)=transData(:,1:2); for k=1:size(transData,1) if transData(k,5) %变压器阻抗归算到1侧计算导纳矩阵 trans(k,3)=1/(transData(k,3)*transData(k,4)); trans(k,4)=(transData(k,4)-1)/(transData(k,3)*transData(k,4)); trans(k,5)=(1-transData(k,4))/(transData(k,4)^2*transData(k,3)); else %变压器阻抗归算到K侧计算导纳矩阵 trans(k,3)=transData(k,4)/transData(k,3); trans(k,4)=(transData(k,4).^2-transData(k,4))/transData(k,3); trans(k,5)=(1-transData(k,4))/transData(k,3); end end %节点数 nodeNum=max([lineData(:,1)' lineData(:,2)' transData(:,1)' transData(:,2)']); %计算节点导纳矩阵 Ybus=zeros(nodeNum); %根据线路参数修改节点导纳矩阵 n1=lineData(:,1); n2=lineData(:,2); for k=1:size(lineData,1) %判断线路参数若为双回线，则阻抗减半，导纳加倍 if lineData(k,5) lineData(k,3)=lineData(k,3)/2; lineData(k,4)=2*lineData(k,4); end Ybus(n1(k),n2(k))=Ybus(n1(k),n2(k))-1/lineData(k,3); Ybus(n2(k),n1(k))=Ybus(n2(k),n1(k))-1/lineData(k,3); Ybus(n1(k),n1(k))=Ybus(n1(k),n1(k))+1/lineData(k,3)+lineData(k,4)/2;%对地导纳Y/2 Ybus(n2(k),n2(k))=Ybus(n2(k),n2(k))+1/lineData(k,3)+lineData(k,4)/2; end %根据变压器参数修改节点导纳矩阵 n1=trans(:,1); n2=trans(:,2); for k=1:size(trans,1) Ybus(n1(k),n2(k))=Ybus(n1(k),n2(k))-trans(k,3); Ybus(n2(k),n1(k))=Ybus(n2(k),n1(k))-trans(k,3); Ybus(n1(k),n1(k))= Ybus(n1(k),n1(k))+trans(k,4)+trans(k,3); Ybus(n2(k),n2(k))= Ybus(n2(k),n2(k))+trans(k,5)+trans(k,3); end end