Sunday, March 30, 2014

z domain root locus

data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAX8AAABXCAIAAABjrtLfAAAGaklEQVR4nO2cS5LbMAwFdf9LK4upcTyyCOFHAZS7V7EFPIIg+KxUqrLtAAAVbNUFAMCXgvsAQA24DwDUgPsAQA24DwDUgPsAQA24DwDUgPsAQA24DwDUgPsAQA24DwDUgPsAQA24Ty+2jROBFtwwisx6F7Ztw3qgFbNnknFvAdYDPZk6mUx8PVgPdGbefDL0xWA90J9JU8rcV4L1wCrMmFVGvwysB9YifWKZ/hqwHliR3LnlAhSw/VJdCICN3NHlAhSA9cC64D4Lg/XA6mTNMNfgbnAfWB3cZ1VwH1gd3GdJsB54BimTzE24Fdxnf/t3E1MrNFnC021MvLYvBPdZD8b6vQP6blxmXbqG1Xo4KRncZzEY6NGFz8pyfH/65/3NieTavpl4f2jufTDNo9eWrCyNywTXhRe4z0p8+TRbX0wcWaYO4z5BcJ+V+PJpXtF9RuseUK74MHCfZfjmMf2hlfvIYfKiWM+LYAfm/JdlAw4xM5ZuC5O6ivvoh/bLD3Tv5j7CD4LjR8OxsbYDgfus4j7vUnJhmoWeTRf3UXqKKcxdiS9xKrhPH/fRO5Swukbh8bRwH1MRl8Hxi9rwqjcs6WaWcx9hFWX646l3H2sFvmmbXdVsutVzP89wny8/xAPF7uNY/vLFJ1bRFKk4uM8+uMyXbdFnKdV0xQ6Xdqs9j0r32X6xZsmCkZLmqQVpVUwVK7qP/PHLz7TMfXzWcymYpTZP00efSmp574PPVoRORtxH85rzibzW48F97tb00aeScoTbe+kssj1d+oJwBLLCqfVwoDXuM8l6ZMHRBGiGMqtON03KAEgE98F9AGoocJ8ZV1ojaLUevfINNCkDIJHu7qOxDI2g/HZzT7VW4zPVCbAc3d3nM+szV+k+I83Ealu5j68YgAiOETWl/En35AQKPU2cunOfeDrBjgHcxtSp/pPuybFXKWdN3bav2nQiHQO4k6lT/Sfdk2OvUs6yCjqCIz1KoUMNALl8nfv47Ln85neoASCXAvfxrZriPrnrCsEm9MqmygGaU+k++oXleKWab6u4D8AkatzHunDcfeR099K30aEGgFyK3cf6TiEHOJ5GbOtOmpQBkEiZ++wWA1K6z2nA8C88utWbXPusMmbsRd9Ma5ZJWT8Ah6fKajVhJhx900+vvrGnkYO78j8yaxojOmkVjIpQnpCjiRHZEuKVzNjIYSKVS7yHjbKsysIAXH6jrzYLR980e9kVV0bZk8vLEp/GqEIk+bOUyw1fKuhlNfrxBiUSLGbGRpT34bIS/X2QK1GmmC5Sk77p9xL/Xoi5/EZP/HJ1uZn7HKd4jPtM2sVpSe5bFFEOus9IfIm+zXCf04BRYULZAo9ynz3bLFpZzx6oZ9JGrPOtz3LcnNOnp5Y00tGbV5BI35R7ufxer/MeIBQmJAqCj3KfPfXHqpX1/OA7sMtBfCdLNpJlUpazPvelF2/Vt92yF80jpc57jHUhWc2RdRQJ5qeTsqtEnVzcJz2Scl8hoRhZSpNlUlZmbb/vR6PCPp+26tth9VF5SjW9zive/dRRnlYkmD+D+MZ6Ws/uKkyYdeFjUFk596MsvfL7R2FdpVMoaxA+aoi4z+7aS0TnFRxZKyXlRCSYPw/33tpazw/W8jTxvv26b9Hn0B++meo+Qm2aGkar60l0H+VeIjqv4EspWSESP9SJS8zDscPOvvND+kl/3udLZGXlNB8EN4WPHL4/rfwyRShPszt59al9c+8lorMrLoXmxCPxQ524BJjIPenIELhv0WWKRvnSaOQwjX63vrn3EtE5jdevlRI/1IlLgBXT4UWmyqfscx+r8ibiq1DpPiV9i+wluHqu+2RZz477lCD/Up0Gjx4JH33iPuv5THEo61Ny3eeevlW5j2Z31mnEfdYm7j5x1zjNsuroL0MT96nq25PcRxOpUssSAiumIz/NPRAv41THOu56ZY2g6Z5/Rp6qVfXNuhf9I31P3DHyKj5wn0qUx3kIO71CkbGQRUarx5XlhU5FIneytm/6vWTF7KkvPsFenWvmyoEVzaHOOPjncek+cErhBHI29SiP/55i1mX0jnB/JWtR+OPH2bQAAwoiNIe+CdS+d3MwXcCA3NA3H+V/5edUesFFsqJsCH070GHSOBIAqAH3AYAacB8AqAH3AYAacB8AqAH3AYAacB8AqAH3AYAacB8AqAH3AYAa/gGWxnxIsyEOvAAAAABJRU5ErkJggg==

z = tf('z');
sys = (0.0173*z)/((z - 0.9048)*(z - 0.8187));
rlocus(sys)
axis([-1 1 -1 1])
zeta = 0.5;
Wn = 0.2;
zgrid(zeta,Wn)


 http://ctms.engin.umich.edu/CTMS/index.php?aux=Extras_Dsserror

Thursday, March 27, 2014

Design a compensator

HW #6

s = tf('s');
k=19.678
sys = (100*k)/(s*(s + 16.8)^2);
sys_cl = feedback(sys,1)
step(sys_cl)


Saturday, March 22, 2014

Root Locus for Lightly Damped Systems

 HW # 7
s = tf('s');
sys = (100)/((s + j*2)*(s - j*2));
rlocus(sys)
axis([-10 4 -6 10])
Zeta = 0.456;
Wn = 1.8;
sgrid(Zeta,Wn)
 


G(s)K(s)/(1+G(s)K(s)


r = roots(p)
p = [1 20 114.5 190.5]



    1.0000   20.0000  114.5000  190.5000


r =

  -11.4580
   -5.5420
   -3.0000

 

 s = tf('s');
k=3.02
sys = (110.5*k*(s + 1))/(s*(s + 1.53)*(s + 3)*(s + 5.542)*(s +11.45));
sys_cl = feedback(sys,1)
step(sys_cl)


Sunday, March 2, 2014

root locus -lead compensation

s = tf('s');
sys = (20*16.8)/(s*(s + 5)*(s + 20));
rlocus(sys)
axis([-11 0 -1 4])
Zeta = 0.5;
Wn = 4.8;
sgrid(Zeta,Wn)