3545336548

Laboratorium Modelowania i symulacji 2008 r.

%odpowiedź skokowa figurę(XI)

axis([0,10,0,2])    % dopasować zakresy osi

hołd on grid on

[wy,czas]=step(system,max(t)); plot(czas,wy,opis (i, :)) xlabel('czas [s] '); ylabel(['T=',num2str(T)]) title('odpowiedź skokowa')

%odpowiedź prędkościowa figurę(12)

axis([0,10,0,2])    % dopasować zakresy osi

hołd on

[wy,czas]=lsim(system,wymuszenie_predkosciowe, t);

plot(czas,wy,opis (i,:))

xlabel('czas [s]');

ylabel(['T=',num2str(T)])

title('odpowiedź prędkościowa')

%charakterystyki częstotliwościowe Bodego figurę(13)

axis([0.01,100,-90,30])    % dopasować zakresy osi

hołd on

bodę(system,opis(i,:));

title('charakterystyki częstotliwościowe Bodego') end

4. Układ całkujący idealny

Utworzyć i uruchomić skrypt uklad_calkujacy_id.m

%układ całkujący idealny %G(s)=k/Ts

opis= [ 'b—'; ' r :x'; "m. - ' ];

k=l    %przeprowadzić symulacje dla a) k=l; b) k=2; c) k=0.5

T=[0.5,2,5]    %przeprowadzić symulacje dla a) T=0.5,2,5; b) T=l,3,6 c) T=0.25,l,4 for i=l:3

%definicja systemu licznik=[k] mianownik=[T(i),0]

system=tf(licznik,mianownik)

% odpowiedź impulsowa figurę(10)

axis([0,10,-2,2.5])    % dopasować zakresy osi

hołd on

[wy, czas]=impulse(system,max(t)); plot(czas,wy,opis (i,:)) xlabel('czas [s] '); ylabel(['T=',num2str(T)]) title('odpowiedź impulsowa')

-3-

Ćwiczenie 4 - Podstawowe układy dynamiczne