Untitled 37

/ T 1t\v«    \/v' \‘li ł/\

/lin_pn/el_liniil/Integrator /lin_pn/el_linii2/Integrator 1 /lin_pn/el_linii2/Integrator /I in__pn/el_linii3/Integrator 1 /I in__pn/el_linii 3 / Integrator /lin_pn/el_linii4/Integrator1

Na tej podstawie możemy przyporządkować kolejne sygnały poszczególnym wyjściom, jak na rysunku 26.

Uwaga: Należy zwrócić uwagę, że SIMULINK nie przypisuje numeracji zmiennym stanu kolejno, ponadto że to przyporządkowanie się zmienia nawet po drobnej zmianie schematu, i dlatego trzeba to sprawdzać przed każdym seansem symulacji.

Teraz możemy przeprowadzić eksperymenty symulacyjne. Polecenie: (t,x,y|=rk23('linia\10,(],(],ster)

spowoduje rozwiązanie układu równań modelu, czyli wprowadzenie do przestrzeni roboczej MATLAB’a wektora czasu /, macierzy zmiennych stanu x oraz macierzy wielkości wyjściowych w przedziale czasu [0, 10], dla macierzy wielkości wejściowych ster ¹.

Rys. 26. Zdyskrę tyzowcmy do pięciu segmentów schemat blokowy pneumatycznej linii długiej, xl ... xl0 - numeracja zmiennych stanu (lin jpn.m)

Polecenie whos wyprowadzi informację o wymiarach uzyskanych macierzy, w szczególności przydatna będzie informacja o liczbie kroków (wymiar wektora t).

Do inicjacji symulacji możemy też użyć programu uruchomieniowego utworzonego graficznie, np. takiego jak pokazano na rysunku 27.

©-

Clock

Mux

Mux

>ULEL

iin_pn

Constant

Rys. 27. Program uruchomieniowy (lin}tru.nt)

—M pQ

To Workspace

>3

Scope

Z kolei polecenie: plot(t,x)

spowoduje narysowanie przebiegów czasowych wszystkich zmienin* stanu (Rys. 28).Przykładowe polecenie: plot(t,x( 1:59,8)) ²

lub prościej, korzystając z notacji dwukropkowej: plot(t,x(:,8))

spowoduje narysowanie przebiegu czasowego tylko 8--ej zmiennej st;»».• > Obecnie wykorzystamy możliwości grafiki 3-D. Załóżmy, że cłu • ■< narysować czasowe przebiegi strumieni m(i). W tym celu należy utwotf nową macierz, w której znajdą się tylk° ^te zmienne stanu, które są • ‘■ mieniami: z poprzednio otrzymanego komunikatu:

1

Pierwszą kolumną macierzy ster musi być kolumnowy wektor czasu

2

Liczba 59 jest w tym przykład/ir liczbą kroków całkowania, czyli wymiarem wektora t