Advertisement
Guest User

Untitled

a guest
May 29th, 2017
259
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 0.61 KB | None | 0 0
  1. pkg load symbolic
  2. syms x1 x2 x3 x4 x5 x6
  3.  
  4. k1 = 30;
  5. k2 = 50;
  6. k3 = 70;
  7. m1 = 5;
  8. m2 = 20;
  9. m3 = 40;
  10. b1 = 17.14;
  11. b2 = 57.83;
  12. b3 = 116.42;
  13.  
  14. dx1 = x2
  15. dx2 = - b1*(x2 - x4)/m1 - k1*(x1 - x3)/m1
  16. dx3 = x4
  17. dx4 = b1*(x2 - x4)/m2 + k1*(x1 - x3)/m2 - k2*x3/m2 - b2*x4/m2
  18. dx5 = x6
  19. dx6 = b2*(x4 - x6)/m3 + k2*(x3 - x5)/m3 - k3*x5/m3 - b3*x6/m3
  20.  
  21. # State space matrix - A
  22. A = [0 1 0 0 0 0 ; -6 -857/250 6 857/250 0 0; 0 0 0 1 0 0; 3/2 857/1000 -4 -7497/2000 0 0; 0 0 0 0 0 1; 0 05/4 5783/4000 -3 -697/160]
  23.  
  24. # Eigenvalues
  25. e = eig(A)
  26.  
  27. # I want the same eigenvalues for all position so alla position are equal in overshoot.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement