ex7_9cap_sec - A(6,2)=1 y40=ys24 ys45 A(7,4)=1 y50=ys25...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
% Example 7.9 bus system data % % Number of generator buses ng = 2; % Number of load buses nl = 3; % Number of buses n = ng+nl; n % Ybus (from the last homework) A=zeros(8,5); y15=.9323-j*12.4301; A(1,1)=1; A(1,5)=-1; y25=.4464-j*4.9598; A(2,2)=1; A(2,5)=-1; y45=.8928-j*9.9197; A(3,4)=1; A(3,5)=-1; y24=.2232-j*2.4799; A(4,2)=1; A(4,4)=-1; y34=1.8645-j*24.8602; A(5,3)=1; A(5,4)=-1; ys25=j*.22/2; ys24=j*.43/2; ys45=j*.11/2; % Take out transformer from 1-5 % y15=0.5*y15; % - Exceeds transformer rating % Take out transformer from 3-4 % y34=0.5*y34; % - Voltage slightly low, say secure % Take out line from 2-5 % y25=0.5*y25; % - Load flow doesn't converge % ys25=0.5*ys25; % Take out line from 2-4 % y24=0.5*y24; % - Load flow converges but voltage too low % ys24=0.5*ys24; % Take out line from 4-5 % y45=0.5*y45; % - Voltage slightly low, say secure % ys45=0.5*ys45; % % Form Y y20=ys24+ys25;
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full Document Right Arrow Icon
Background image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: A(6,2)=1; y40=ys24+ys45; A(7,4)=1; y50=ys25+ys45; A(8,5)=1; Yl=diag([y15; y25; y45; y24; y34; y20; y40; y50]); Y=A'*Yl*A; Y % Power injections (unknown values set to 0) P = [0; 1.1; -2.00; 0; 0]; Q = [0; 0; -0.70; 0; 0]; Q % Specified voltages (unknown magnitudes set to 1 unknown angles set to 0) V=[1.0; 1.05; 1.0; 1.0; 1.0]; delta=zeros(5,1); d % Relabel buses so bus 1 is slack and bus 2 is the generator bus Ybus=[Y(1,1) Y(1,3) Y(1,2) Y(1,4:5); Y(3,1) Y(3,3) Y(3,2) Y(3,4:5); Y(2,1) Y(2,3) Y(2,2) Y(2,4:5); Y(4:5,1) Y(4:5,3) Y(4:5,2) Y(4:5,4:5)]; % Install capacitor at bus 3 Ybus(3,3)=Ybus(3,3)+j*0.5; Y % B1 and B2 for calculating angle and voltage updates respectively Bp=-imag(Ybus-diag(sum(Ybus))); %B' with neglecting any shunt capacitance Bpp = -imag(Ybus); % Neglect phase-shifting transformers but none here B...
View Full Document

{[ snackBarMessage ]}

Page1 / 2

ex7_9cap_sec - A(6,2)=1 y40=ys24 ys45 A(7,4)=1 y50=ys25...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online