Routh Example #2

MATLAB Code

%         ***** MATLAB Code Starts Here
%

p = linspace(-5,5,1001); % range of pole locations to be graphed

p_quad = p.^2 + p - 10; % quadradic function in the s^1 row

p1p2 = roots([1 1 -10]); % roots of the quadratic function

p_lin = p + 1; % first-order function in the s^2 row

p3 = roots([1 1]); % root of the first-order function

figure(1),clf,plot(p,p_quad,'b',p,p_lin,'r',[p1p2(1) p1p2(1)],[-15 20],'k--',...

[p1p2(2) p1p2(2)],[-15 20],'k--',[p3 p3],[-15 20],'k--'),plotax

xlabel('Pole Locations p'),ylabel('Values of Constraint Functions'),...

title('First-Order and Second-Order Constraints on Pole Location')

text(-4.5,13,'p < -3.7')

text(-0.9,13,'p > -1')

text(2.8,13,'p > 2.7')

%
%         ***** MATLAB Code Stops Here

Click the icon to return to the Dr. Beale's home page

Lastest revision on Wednesday, June 7, 2006 12:02 PM