31 lines
2.2 KiB
Matlab
31 lines
2.2 KiB
Matlab
function [G,B,GB,Y,r,c,Angle] = admmatrix(Busnum,Linei,Linej,Liner,Linex,Lineb,Transfori...
|
||
,Transforj,Transforr,Transforx,Transfork0,Branchi,Branchb)
|
||
%**************************************************************************
|
||
% 程序功能 : 子函数——形成节点导纳矩阵Y
|
||
% 编 者：
|
||
% 编制时间：2010.12
|
||
%**************************************************************************
|
||
%% 支路导纳计算
|
||
G = -sparse(Linei,Linej,Liner./(Liner.^2+Linex.^2),Busnum,Busnum) - sparse(Linej,Linei,Liner./(Liner.^2+Linex.^2),Busnum,Busnum);
|
||
G = G - sparse(1:Busnum,1:Busnum,sum(G,2)'); % 计算各线路支路电导
|
||
B = sparse(Linei,Linej,Linex./(Liner.^2+Linex.^2),Busnum,Busnum)+sparse(Linej,Linei,Linex./(Liner.^2+Linex.^2),Busnum,Busnum);
|
||
B = B - sparse(1:Busnum,1:Busnum,sum(B,2)')+sparse(Linei,Linei,Lineb,Busnum,Busnum)+sparse(Linej,Linej,Lineb,Busnum,Busnum);
|
||
%% 变压器支路计算
|
||
if Transfori>0
|
||
mr = Transforr./(Transforr.^2+Transforx.^2); % 计算变压器支路电导
|
||
mx = -Transforx./(Transforr.^2+Transforx.^2); % 计算变压器支路电纳
|
||
G = G-sparse(Transfori,Transforj,mr./Transfork0,Busnum,Busnum)-sparse(Transforj,Transfori,mr./Transfork0,Busnum,Busnum)...
|
||
+sparse(Transfori,Transfori,mr./Transfork0./Transfork0,Busnum,Busnum)+sparse(Transforj,Transforj,mr,Busnum,Busnum);
|
||
B = B-sparse(Transfori,Transforj,mx./Transfork0,Busnum,Busnum)-sparse(Transforj,Transfori,mx./Transfork0,Busnum,Busnum)...
|
||
+sparse(Transfori,Transfori,mx./Transfork0./Transfork0,Busnum,Busnum)+sparse(Transforj,Transforj,mx,Busnum,Busnum);
|
||
end
|
||
%% 接地支路计算
|
||
if Branchi>0 % 判断有无接地支路
|
||
B = B+sparse(Branchi,Branchi,Branchb,Busnum,Busnum);
|
||
end
|
||
%% 化作极坐标形式
|
||
GB = G+B.*1i; %将电导，电纳合并，写成复数形式
|
||
Y = abs(GB); %求节点导纳幅值
|
||
[r,c] = find(Y);
|
||
Angle = angle(GB(GB~=0)); %求节点导纳角度
|
||
%Angle=angle(GB); |