3545336547

Laboratorium Modelowania i symulacji 2008 r.

%odpowiedź skokowa figurę(11)

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

hołd on

[wy,czas]=step(system,max(t)); plot(czas,wy,opis(i,:)) xlabel ('czas [s] '); ylabel(['b0=',num2str(bO)]) 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,wymuśzenie_predkosciowe,t); plot(czas,wy,opis (i, :)) xlabel('czas [s]*); ylabel(['b0=',num2str(bO)]) title('odpowiedź prędkościowa')

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

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

hołd on

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

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

3. Układ inercyjny pierwszego rzędu

Utworzyć i uruchomić skrypt układ _inercyjny_p_rz.m

%układ inercyjny pierwszego rzędu %G(s)=k/Ts+l

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

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]

system=tf(licznik, mianownik)

% odpowiedź impulsowa figurę(10)

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

-2-

Ćwiczenie 4 - Podstawowe układy dynamiczne