3545336550

Laboratorium Modelowania i symulacji 2008 r.

figurę(11)

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

hołd 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.5])    % 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,-200,100])    % dopasować zakresy osi

hołd on

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

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

6. Układ różniczkujący rzeczywisty

Utworzyć i uruchomić skrypt uklad_rozniczkujacy_rz.m

%układ różniczkujący rzeczywisty

%G(s)=ks/Ts+l

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

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

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 0] mianownik=[T(i),1]

system=tf(licznik,mianownik)

% odpowiedź impulsowa figurę (10)

axis ([0,10,-5,2])    % 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')

-5-

Ćwiczenie 4 - Podstawowe układy dynamiczne