fm_y.m

这是一个很适合研究和学习用的电力系统仿真软件

function fm_y
%FM_Y build admittance matrix
%
%FM_Y
%
%see also FM_LF and DAE, Bus and Line structures
%
%Author:    Federico Milano
%Date:      11-Nov-2002
%Version:   1.0.0
%
%E-mail:    [e-mail]
%Web-site:  http://thunderbox.uwaterloo.ca/~fmilano
%
% Copyright (C) 2002-2006 Federico Milano

global Line Bus Shunt Settings

if ~Line.n, return, end

% process line data and build admittance matrix Line.Y
chrg = Line.con(:,10);
y = 1./(Line.con(:,8) + j*Line.con(:,9));
ts = Line.con(:,11).*exp(j*Line.con(:,12)*pi/180);
ts2= ts.*conj(ts);

Line.Y = sparse(Line.from, Line.to, -y./conj(ts), Bus.n, Bus.n) + ...
         sparse(Line.to, Line.from, -y./ts, Bus.n, Bus.n) + ...
         sparse(Line.from, Line.from, (y+j*chrg/2)./ts2, Bus.n, Bus.n) + ...
         sparse(Line.to, Line.to, y+j*chrg/2, Bus.n, Bus.n);

% add diagonal shunt admittances
b = [1:Bus.n];
Line.Y = Line.Y + sparse(b,b,Shunt.g+j*Shunt.b,Bus.n,Bus.n);

% check for missing connections (0 diagonal elements)
b = find(diag(Line.Y)==0);
%c = b-1;
%if c(1) < 0, c(1) = b(1)+1; end
Line.Y = Line.Y - sparse(b,b,j*1e-6,Bus.n,Bus.n);

%         + sparse(c,b,j*1e-6,Bus.n,Bus.n) ...
%         + sparse(b,c,j*1e-6,Bus.n,Bus.n) ...
%         - sparse(c,c,j*1e-6,Bus.n,Bus.n);