3545336549

Laboratorium Modelowania i symulacji 2008 r.

%odpowiedź skokowa figurę(XI)

axis([0,10,0,2.5])    % 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ź predkosciowa 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,-90,30])    % dopasować zakresy osi

hołd on

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

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

5. Układ całkujący rzeczywisty

Utworzyć i uruchomić skrypt uklad_calkujacy_rz.m

%układ całkujący rzeczywisty

%G(s)=k/s(Ts+1)

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-1,3,6 c) T-0.25,1,4 for i=l:3

%definicja systemu licznik=[k] mianownik=[T(i),1 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')

%odpowiedź skokowa

-4-

Ćwiczenie 4 - Podstawowe układy dynamiczne