lecture_23_exmp

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

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

View Full Document Right Arrow Icon
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;
Background image of page 1

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

View Full DocumentRight 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: 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

This document was uploaded on 03/23/2010.

Page1 / 2

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

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