lab4dynamiczne

%Jan Zasadny

open('model1');
sim('model1');

figure;
plot(tout,simout);
title('Częstotliwość 0.1 Hz');

figure;
plot(tout,simout1);
title('Częstotliwość 0.5 Hz');

figure;
plot(tout,simout2);
title('Częstotliwość 1 Hz');



G1= tf([1],[1 1 2]);
G2= tf([1 0],[1 -0.1 1]);

figure;
hold on;
nyquist(G1,'r');
nyquist(G2,'g');
plot(-1,0);
hold off;

%Jan Zasadny


tabGG=[[1],[0.4 0],[1 0],[1]];
tabGD=[[4 1],[0.4 1],[1 2 1],[1 0.1 1]];
tabA=[2,3,1,2,4];
tabOmega=[0.001,1,30,0.01,5];
tabFi=[1,5,2,1,0];

for i=1:1:4
    figure;
    G=tf(tabGG(i),tabGD(i));
    bode(G,{1e1,1e4});
end;

%for i=1:1:4


%Jan Zasadny



%Jan Zasadny

freq1=0;
freq2=freq1/64;

RC=1/(2*pi*freq2);

G=tf([1],[RC 1]);

figure;
hold on;
bode(G,freq1,'ro');
bode(G,freq2,'go');
bode(G,{1e1,1e4});
hold off;



Jan Zasadny

Zadanie 1
Na początku wykresy wszystkich systemów oscylują, po pewnym czasie przyjmują postać sinusoidy.
Wraz ze wzrostem częstotliwości wymuszeń sinusoidalnych zwiększa się czas jaki system potrzebuje do stabiliacji.

Zadanie 2
Pierwszy układ jest stabilny ponieważ wykres charakterystyki nie obejmuje punktu [-1,0]

Drugi układ nie jest stabilny ponieważ wykres obejmuje ten punkt.

Zadanie 3