2/16/14 File Exchange - MATLAB Central

www.mathworks.in/matlabcentral/fileexchange/28954-y-bus/content/ybus_final.m 1/2

Search: File Exchange

Site Help Patents Trademarks Privacy Policy Preventing Piracy Terms of Use

Featured MathWorks.com Topics: New Products Support Documentation Training Webinars Newsletters MATLAB Trials Careers

Create Account Log In

File Exchange Answers Newsgroup Link Exchange Blogs Trendy Cody Contest MathWorks.com

by arvind singh

Y-bus

06 Oct 2010

admittance matrix

ybus_final.m

Dow nload Submission

Code covered by the BSD License

Highlights from Y-bus

ybus_final.mProgram to form Admittance Bus Formation with Transformer Tap setting..

View all files

% Program to form Admittance Bus Formation with Transformer Tap setting..

clc

clear all

nbranch=input('enter total no. of lines in network');

disp('enter line data')

for p=1:1:nbranch

fb=input('from bus no. :');

tb=input('to bus no. :');

r=input('resistance of line =');

x=input('reactance of line =');

B=input('line charging admittance"b/2"=');

T=input('tap=');

z = r + i*x; % Z matrix...

y = 1./z; % To get inverse of each element...

Ldata(p,:)=[fb tb r x B T y];

end

fb = Ldata(:,1); % From bus number...

tb = Ldata(:,2); % To bus number...

r = Ldata(:,3); % Resistance, R...

x = Ldata(:,4); % Reactance, X...

b = Ldata(:,5); % Ground Admittance, B/2...

a = Ldata(:,6); % Tap setting value..

y = Ldata(:,7);

b = i*b; % Make B imaginary...

nbus = max(max(fb),max(tb)); % no. of buses...

Y = zeros(nbus,nbus); % Initialise YBus...

% Formation of the Off Diagonal Elements...

for k=1:nbranch

Y(fb(k),tb(k)) = Y(fb(k),tb(k))-y(k);

Y(tb(k),fb(k)) = Y(fb(k),tb(k));

end

% Formation of Diagonal Elements....

for m =1:nbus

for n =1:nbranch

if fb(n) == m

Y(m,m) = Y(m,m) + y(n) + b(n);

elseif tb(n) == m

Y(m,m) = Y(m,m) + y(n) + b(n);

end

end

end

Y % Bus Admittance Matrix..

Contact us

© 1994-2014 The MathWorks, Inc.

2/16/14 File Exchange - MATLAB Central

www.mathworks.in/matlabcentral/fileexchange/28954-y-bus/content/ybus_final.m 2/2