{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

lecture_23_exmp

# lecture_23_exmp - My_Sys_3 = tf[1 2*xi_z*omega_z omega_z^2...

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

C:\Documents and Settings\me\Desktop\Lecture_23 (1).m Page 1 November 7, 2006 10:19:42 AM % Root locus for system 1; xi_z = 0.2; omega_z = 5; p_1 = 3; p_2 = 10; xi_p = 0.5; omega_p = 20; % Construct the open loop transfer function. NOTE how the command % ``conv'' is used here to find the coefficients of the denominator. My_Sys_1 = tf([1 2*xi_z*omega_z omega_z^2], conv(conv([1 p_1], [1 ... p_2]), [1 2*xi_p*omega_p omega_p^2])); figure(4), rlocus(My_Sys_1); % Root locus for system 2; xi_z = -0.2; % The zeros are on the % right half plan. omega_z = 5; p_1 = 3; p_2 = 10; xi_p = 0.5; omega_p = 20; My_Sys_2 = tf([1 2*xi_z*omega_z omega_z^2], conv(conv([1 p_1], [1 ... p_2]), [1 2*xi_p*omega_p omega_p^2])); figure(5), rlocus(My_Sys_2); % Root locus for system 3: The open loop system is unstable; xi_z = 0.2; omega_z = 5; p_1 = -3; p_2 = 10; xi_p = 0.5; omega_p = 20;

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: My_Sys_3 = tf([1 2*xi_z*omega_z omega_z^2], conv(conv([1 p_1], [1 . .. p_2]), [1 2*xi_p*omega_p omega_p^2])); figure(6), rlocus(My_Sys_3); % Integral control; % The open loop system: Ki = 1; My_Go = tf([5 10], [1 2 36]); figure(7), step(My_Go,'r--'), hold on; My_Gc = tf([Ki], [1 0]); My_Gcl = feedback(My_Go*My_Gc, 1); figure(7), step(My_Gcl); % test the ramp response; My_time = 0:1e-02:25; C:\Documents and Settings\me\Desktop\Lecture_23 (1).m Page 2 November 7, 2006 10:19:42 AM figure(8), subplot(211), lsim(My_Go, My_time, My_time); figure(8), subplot(212), lsim(My_Gcl, My_time, My_time); % Now add another integrator; My_Gc = tf([Ki], [1 0 0]); My_Gcl = feedback(My_Go*My_Gc, 1); figure(9), subplot(211), step(My_Gcl); figure(9), subplot(212), lsim(My_Gcl, My_time, My_time);...
View Full Document

{[ snackBarMessage ]}

### Page1 / 2

lecture_23_exmp - My_Sys_3 = tf[1 2*xi_z*omega_z omega_z^2...

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

View Full Document
Ask a homework question - tutors are online