>> k=1 k = 1 >> den= poly([1 -6 -10]) den = 1 15 44 -60 >> Gs = tf(k , den) Gs = 1 ------------------------ s^3 + 15 s^2 + 44 s - 60 Continuous-time transfer function. >> den= poly([0 -6 -10]) den = 1 16 60 0 >> Gs = tf(k , den) Gs = 1 ------------------- s^3 + 16 s^2 + 60 s Continuous-time transfer function. >> rlocus(Gs) >> 60/193 ans = 0.3109 >> -1.79*2 ans = -3.5800 >> 3.51*2 ans = 7.0200 >> oldsys= feedback(Gs , 1) oldsys = 1 ----------------------- s^3 + 16 s^2 + 60 s + 1 Continuous-time transfer function. >> step(oldsys) >> 4/1.79 ans = 2.2346 >> k=193 k = 193 >> Gs = tf(k , den) Gs = 193 ------------------- s^3 + 16 s^2 + 60 s Continuous-time transfer function. >> oldsys= feedback(Gs , 1) oldsys = 193 ------------------------- s^3 + 16 s^2 + 60 s + 193 Continuous-time transfer function. >> step(oldsys) >> k=1 k = 1 >> Gs = tf(k , den) Gs = 1 ------------------- s^3 + 16 s^2 + 60 s Continuous-time transfer function. >> num1 = [1 6] num1 = 1 6 >> den2 = poly([-29.03]) den2 = 1.0000 29.0300 >> Lead = tf(num1 , den2) Lead = s + 6 --------- s + 29.03 Continuous-time transfer function. >> sys = series(Lead , Gs) sys = s + 6 ------------------------------------ s^4 + 45.03 s^3 + 524.5 s^2 + 1742 s Continuous-time transfer function. >> rlocus(sys) >> Gs = tf(k , den) Gs = 1 ------------------- s^3 + 16 s^2 + 60 s Continuous-time transfer function. >> num2 =[1 0.0479] num2 = 1.0000 0.0479 >> den3 = poly([-0.01]) den3 = 1.0000 0.0100 >> lag = tf(num2 , den3) lag = s + 0.0479 ---------- s + 0.01 Continuous-time transfer function. >> lead-lag = series(lag , sys) lead-lag = series(lag , sys) | Error: The expression to the left of the equals sign is not a valid target for an assignment. >> leadlag = series(lag , sys) leadlag = s^2 + 6.048 s + 0.2874 ------------------------------------------------ s^5 + 45.04 s^4 + 524.9 s^3 + 1747 s^2 + 17.42 s Continuous-time transfer function. >> newsys = feedback(leadlag , 1) newsys = s^2 + 6.048 s + 0.2874 --------------------------------------------------------- s^5 + 45.04 s^4 + 524.9 s^3 + 1748 s^2 + 23.47 s + 0.2874 Continuous-time transfer function. >> leadlag leadlag = s^2 + 6.048 s + 0.2874 ------------------------------------------------ s^5 + 45.04 s^4 + 524.9 s^3 + 1747 s^2 + 17.42 s Continuous-time transfer function. >> rlocus(leadlag) >> k2 =1850 k2 = 1850 >> iii = tf(k2 ,1) iii = 1850 Static gain. >> leadlag2 = series(iii , leadlag) leadlag2 = 1850 s^2 + 1.119e04 s + 531.7 ------------------------------------------------ s^5 + 45.04 s^4 + 524.9 s^3 + 1747 s^2 + 17.42 s Continuous-time transfer function. >> newsys = feedback(leadlag2 , 1) newsys = 1850 s^2 + 1.119e04 s + 531.7 ----------------------------------------------------------- s^5 + 45.04 s^4 + 524.9 s^3 + 3597 s^2 + 1.121e04 s + 531.7 Continuous-time transfer function. >> step(oldsys , newsys)
Posted on: Sat, 21 Sep 2013 04:25:08 +0000
Trending Topics
Recently Viewed Topics
© 2015