Modify this MATLAB code from previous Order ID: 90466 with Excel functions where the user inputs values onto an Excel sheet and the code outputsthe Y BUS.
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); nbr=length(nl); %linedata(:,1) nbus = max(max(nl), max(nr)); Z = R + j*X; y= ones(nbr,1)./Z; Ybus=zeros(nbus,nbus);%branch admittance % initialize Ybus to zero % formation of the off diagonal elements % formation of the diagonal elements for n=1:nbus ; for k=1:nbr 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, end end end Ybus Ybus = Columns 1 through 4 3.0665 - 8.7059i -12.0000 +24.0000i 0.0000 + 0.0000i -7.0588 +28.2353i 1 -12.0000 +24.0000i 8.0083 -20.1255i -4.6154 +23.0769i -24.0000 +48.0000i 0.0000 + 0.0000i -4.6154 +23.0769i 3.9507 -16.6023i 0.0000 + 0.0000i -7.0588 +28.2353i -24.0000 +48.0000i 0.0000 + 0.0000i 5.0662 -15.1496i 0.0000 + 0.0000i 0.0000 + 0.0000i -8.7805 +19.0244i -0.7481 +14.9626i 0.0000 + 0.0000i -9.3541 +26.7261i -11.5385 +57.6923i 0.0000 + 0.0000i Columns 5 through 6 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -9.3541 +26.7261i -8.7805 +19.0244i -11.5385 +57.6923i -0.7481 +14.9626i 0.0000 + 0.0000i 2.3281 - 8.6595i -6.0000 +18.0000i -6.0000 +18.0000i 4.2321 -16.9697i second version "Equipment impedance includ- ed " % Bus Bus R X G B a % No. No. p.u. p.u. p.u. p.u. linedata=[1 1 0 1.25 0 0 1 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 4 0 1.15 0 0 1 4 5 0.02 0.4 0 0 1 5 6 0.1 0.3 0 0.22 1 6 6 0 1.25 0 0 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); nbr=length(nl); %linedata(:,1) nbus = max(max(nl), max(nr)); Z = R + j*X; y= ones(nbr,1)./Z; Ybus=zeros(nbus,nbus);%branch admittance 2 % initialize Ybus to zero % formation of the off diagonal elements % formation of the diagonal elements for n=1:nbus ; for k=1:nbr 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, end end end Ybus Ybus = Columns 1 through 4 3.0665 - 4.7059i -12.0000 +24.0000i 0.0000 + 0.0000i -7.0588 +28.2353i -12.0000 +24.0000i 8.0083 -20.1255i -4.6154 +23.0769i -24.0000 +48.0000i 0.0000 + 0.0000i -4.6154 +23.0769i 3.9507 -16.6023i 0.0000 + 0.0000i -7.0588 +28.2353i -24.0000 +48.0000i 0.0000 + 0.0000i 5.0662 -10.8018i 0.0000 + 0.0000i 0.0000 + 0.0000i -8.7805 +19.0244i -0.7481 +14.9626i 0.0000 + 0.0000i -9.3541 +26.7261i -11.5385 +57.6923i 0.0000 + 0.0000i Columns 5 through 6