Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- %definicje macierzy Jordana dla
- J_lista(:,:,1)=[ 0 0;0 0];
- J_lista(:,:,2)=[ 0 1;0 0];
- J_lista(:,:,3)=[ 0 0;0 -1];
- J_lista(:,:,4)=[-1 0;0 -2];
- J_lista(:,:,5)=[-1 0;0 -1];
- J_lista(:,:,6)=[-1 1;0 -1];
- J_lista(:,:,7)=[-1 0;0 1];
- J_lista(:,:,8)=[ 0 -1;1 0];
- J_lista(:,:,9)=[-1 -1;1 -1];
- %definicje nazw dla interfejsu greficznego
- typeNames = {
- '1. Punkt'
- '2. Poziome Proste'
- '3. Pionowe Proste'
- '4. Węzeł'
- '5. Węzeł 2'
- '6. Zdegenerowany Węzeł'
- '7. Siodło'
- '8. Okrąg'
- '9. Ognisko'
- };
- systemNames = {
- '1. \lambda_{1,2}= 0 - dzielnik liniowy'
- '2. \lambda_{1,2} = 0 -- dzielnik nieliniowy'
- '3. \lambda_1 = -1, \lambda_2 = 0'
- '4. \lambda_1 = -2, \lambda_2 = -1'
- '5. \lambda_{1,2} = -1 - dzielnik liniowy'
- '6. \lambda_{1,2} = -1 - dzielnik nieliniowy'
- '7. \lambda_1 = -1, \lambda_2 = 1'
- '8. \lambda_{1,2} = 0 \pm 1i'
- '9. \lambda_{1,2} = -1 \pm 1i'
- };
- %interfejs graficzny do wyboru typu wykresu
- %mozliwy wielokrotny wybór
- [systems,~]=listdlg('ListString',typeNames);
- %otwarcie modelu -celem prezentacji
- open('model2');
- %2x2 macierz przejścia P
- P= [1 0;
- 0 1];%maci
- %wartości początkowe
- i=1;
- for x=-9.5:1:9.5
- for y=-9.5:1:9.5
- WPC_lista(:,:,i)=[x;y];
- i=i+1;
- end
- end
- for system = systems
- figure(system)
- clf
- J=J_lista(:,:,system);
- A=P*J*P^-1
- hold on;
- for k=1:i-1
- WPC=WPC_lista(:,:,k);
- sim('model2')
- x1 = x(:, 1); % wybierz x1
- x2 = x(:, 2); % wybierz x2
- plot(x1,x2,'k')
- end
- % plotv(P(:,1),'-r')
- % plotv(P(:,2),'-r')
- points=P(:,1).*linspace(-10,10);
- line(points(1,:),points(2,:),'Color','black','LineStyle','--')
- points=P(:,2).*linspace(-10,10);
- line(points(1,:),points(2,:),'Color','black','LineStyle','--')
- quiver([0 0],[0 0],P(1,:),P(2,:),1)
- % quiver([0;0],P(:,2),P(:,2),0)
- hold off;
- % formatujemy wykres
- title(systemNames(system)); % tytuł
- xlabel('x1'); % osie: x1 i x2
- ylabel('x2');
- xlim([-10 10]); % zakres: od -10.0 do 9.0
- ylim([-10 10]);
- grid on; % włącz siatkę
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement