Great Deal! Get Instant $10 FREE in Account on First Order + 10% Cashback on Every Order Order Now

This assignment requires a MATLAB code using background from power system analysis (Y BUS Matrix). Read instructions 1 of 1 and 1 of 2. The rest of the file attachments are examples of problems to...

1 answer below »
This assignment requires a MATLAB code using background from power system analysis (Y BUS Matrix). Read instructions 1 of 1 and 1 of 2. The rest of the file attachments are examples of problems to test your code.
Answered Same Day Sep 02, 2021

Solution

Vishvajeet answered on Sep 02 2021
169 Votes
First version "Transmission line only "
% Bus Bus R X G B a
% No. No. p.u. p.u. p.u. p.u.
linedata=[1 2 0.1 0.2 0 0 1
1 4 0.05 0.2 0 0.22 1
2 3 0.05 0.25 0.05 0.11 1
2 4 0.05 0.1 0.1 0.25 1
2 6 0.07 0.2 0.2 0.28 1
3 5 0.12 0.26 0.01 0.3 1
3 6 0.02 0.1 0 0 1
4 5 0.02 0.4 0 0 1
5 6 0.1 0.3 0 0.22 1];
nl = linedata(:,1);
nr = linedata(:,2);
R = linedata(:,3);
X = linedata(:,4);
B = j*linedata(:,6);
G = j*linedata(:,5);
Bc = .5*( G + j*B);
a = linedata(:,7);
n
=length(nl); %linedata(:,1)
nbus = max(max(nl), max(nr));
Z = R + j*X;
y= ones(n
,1)./Z;
Ybus=zeros(nbus,nbus);%
anch admittance
% initialize Ybus to zero
% formation of the off diagonal elements
% formation of the diagonal elements
for n=1:nbus ;
for k=1:n
Ybus(nl(k),nr(k))=Ybus(nl(k),nr(k))-y(k)/a(k);
Ybus(nr(k),nl(k))=Ybus(nl(k),nr(k));
if nl(k)==n
Ybus(n,n) = Ybus(n,n)+y(k)/(a(k)^2) + Bc(k);
elseif nr(k)==n
Ybus(n,n) = Ybus(n,n)+y(k) +Bc(k);
else,...
SOLUTION.PDF

Answer To This Question Is Available To Download

Related Questions & Answers

More Questions »

Submit New Assignment

Copy and Paste Your Assignment Here