Advertisement
Guest User

Untitled

a guest
Dec 13th, 2018
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.18 KB | None | 0 0
  1.  
  2.  
  3. %definicje macierzy Jordana dla
  4. J_lista(:,:,1)=[ 0 0;0 0];
  5. J_lista(:,:,2)=[ 0 1;0 0];
  6. J_lista(:,:,3)=[ 0 0;0 -1];
  7. J_lista(:,:,4)=[-1 0;0 -2];
  8. J_lista(:,:,5)=[-1 0;0 -1];
  9. J_lista(:,:,6)=[-1 1;0 -1];
  10. J_lista(:,:,7)=[-1 0;0 1];
  11. J_lista(:,:,8)=[ 0 -1;1 0];
  12. J_lista(:,:,9)=[-1 -1;1 -1];
  13.  
  14. %definicje nazw dla interfejsu greficznego
  15. typeNames = {
  16. '1. Punkt'
  17. '2. Poziome Proste'
  18. '3. Pionowe Proste'
  19. '4. Węzeł'
  20. '5. Węzeł 2'
  21. '6. Zdegenerowany Węzeł'
  22. '7. Siodło'
  23. '8. Okrąg'
  24. '9. Ognisko'
  25.  
  26. };
  27. systemNames = {
  28. '1. \lambda_{1,2}= 0 - dzielnik liniowy'
  29. '2. \lambda_{1,2} = 0 -- dzielnik nieliniowy'
  30. '3. \lambda_1 = -1, \lambda_2 = 0'
  31. '4. \lambda_1 = -2, \lambda_2 = -1'
  32. '5. \lambda_{1,2} = -1 - dzielnik liniowy'
  33. '6. \lambda_{1,2} = -1 - dzielnik nieliniowy'
  34. '7. \lambda_1 = -1, \lambda_2 = 1'
  35. '8. \lambda_{1,2} = 0 \pm 1i'
  36. '9. \lambda_{1,2} = -1 \pm 1i'
  37. };
  38. %interfejs graficzny do wyboru typu wykresu
  39. %mozliwy wielokrotny wybór
  40. [systems,~]=listdlg('ListString',typeNames);
  41. %otwarcie modelu -celem prezentacji
  42. open('model2');
  43. %2x2 macierz przejścia P
  44. P= [1 0;
  45. 0 1];%maci
  46.  
  47. %wartości początkowe
  48. i=1;
  49. for x=-9.5:1:9.5
  50. for y=-9.5:1:9.5
  51. WPC_lista(:,:,i)=[x;y];
  52. i=i+1;
  53. end
  54. end
  55.  
  56. for system = systems
  57. figure(system)
  58. clf
  59. J=J_lista(:,:,system);
  60. A=P*J*P^-1
  61. hold on;
  62. for k=1:i-1
  63. WPC=WPC_lista(:,:,k);
  64. sim('model2')
  65. x1 = x(:, 1); % wybierz x1
  66. x2 = x(:, 2); % wybierz x2
  67. plot(x1,x2,'k')
  68. end
  69.  
  70. % plotv(P(:,1),'-r')
  71. % plotv(P(:,2),'-r')
  72.  
  73. points=P(:,1).*linspace(-10,10);
  74. line(points(1,:),points(2,:),'Color','black','LineStyle','--')
  75. points=P(:,2).*linspace(-10,10);
  76. line(points(1,:),points(2,:),'Color','black','LineStyle','--')
  77.  
  78. quiver([0 0],[0 0],P(1,:),P(2,:),1)
  79. % quiver([0;0],P(:,2),P(:,2),0)
  80.  
  81.  
  82. hold off;
  83.  
  84. % formatujemy wykres
  85. title(systemNames(system)); % tytuł
  86. xlabel('x1'); % osie: x1 i x2
  87. ylabel('x2');
  88. xlim([-10 10]); % zakres: od -10.0 do 9.0
  89. ylim([-10 10]);
  90. grid on; % włącz siatkę
  91. end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement