ex7_9cap_sec

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

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

% 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;

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

View Full Document
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
Ask a homework question - tutors are online