Automaty_2 by Endrju i Falu.doc

Akademia Górniczo-Hutnicza

im. Stanisława Staszica

Wydział Inżynierii Mechanicznej i Robotyki

PODSTAWY AUTOMATYKI – LABORATORIUM

Temat ćwiczenia:

REDUKCJA SCHEMATÓW BLOKOWYCH

ORAZ WYZNACZANIE ICH CHARAKTERYSTYK

Wykonali:

Fal Arkadiusz

                            Latacz Andrzej

                            grupa 7A

Przebieg ćwiczenia:

Przykład 1

Model układu w Simulinku

Zredukowana transmitancja (rozwiązanie analityczne):

Listing m-pliku zawierający funkcje wyznaczającą transmitancje zastępczą:

function rozw2

l1=[1];

m1=[4 1];

l2=[2];

m2=[1 0];

l3=[1];

m3=[3 0];

l4=[4];

m4=[2 1];

l5=[10];

m5=[1];

[l10,m10]=parallel(l3,m3,l5,m5)

[l11,m11]=feedback(l2,m2,l4,m4,-1)

[l12,m12]=series(l1,m1,l11,m11)

[l,m]=series(l12,m12,l10,m10)

g=tf(l,m)

[A,B,C,D]=tf2ss(l,m)

impulse(l,m)

grid

Odpowiedź programu:

Transfer function:

      120 s^2 + 64 s + 2

-------------------------------

24 s^4 + 18 s^3 + 99 s^2 + 24 s

A =

   -0.7500   -4.1250   -1.0000         0

    1.0000         0         0         0

         0    1.0000         0         0

         0         0    1.0000         0

B =

     1

     0

     0

     0

C =

         0    5.0000    2.6667    0.0833

D =

     0

Charakterystyka skokowa

(Simulink)

Charakterystyka impulsowa

(Matlab)

Przykład 2

Wpływ zmiany parametrów na kształt charakterystyk logarytmicznych elementu oscylacyjnego opisanego transmitancją.

Listing m-pliku pozwalający na otrzymanie charakterystyk logarytmicznych

function charakterystyka

k=input ('Podaj wartosc wspolczynnika k= ');

E=input ('Podaj wartosc wspolczynnika E= ');

w=input ('Podaj wartosc wspolczynnika w= ');

w1=input ('Podaj wartosc wspolczynnika w1= ');

E1=input ('Podaj wartosc wspolczynnika E1= ');

L=[k.*w.^2];

M=[1 2.*E.*w w.^2];

L1=[k.*w.^2];

M1=[1 2.*E1.*w w.^2];

L2=[k.*w.^2];

M2=[1 2.*E.*w w.^2];

L3=[k.*w1.^2];

M3=[1 2.*E.*w1 w1.^2];

hold on

bode (L, M)

bode (L1, M1)

legend('Wykres dla E=0.2', 'Wykres dla E=0.8');

grid

hold off

pause

hold on

bode (L2, M2)

bode (L3, M3)

legend('Wykres dla w=1.0', 'Wykres dla w=4.5');

grid

hold off

Zmienne wartości stałe

Zmienne wartości stałe

Wnioski:

·         Zwiększenie wartości współczynnika tłumienia w układzie oscylacyjnym drugiego rzędu spowodowało wzrost modułu dla wartości pulsacji drgań własnych w otoczeniu punktu 1.

·         Położenie punktu przegięcia jest niezależne od wartości współczynnika tłumienia. Jest natomiast zależne od wartości pulsacji drgań własnych

·         Zwiększenie wartości  pulsacji drgań własnych jest widoczne na obu wykresach.

2