wyklad(1), Politechnika Wrocławska - Materiały, teoria systemow


Podstawy techniki systemów

wg dr inż. Andrzeja Żołnierka

Literatura: „Podstawy informacji zarządzania” Budnicki J.

WYKŁAD I

TEORIA I TECHNIKA SYSTEMÓW

1. Tworzenie modeli i modelowanie.

2. Identyfikacja i rozpoznawanie.

3. Analiza i projektowanie.

4. Sterowanie i kierowanie.

Ad1. Tworzenie Modeli i modelowanie.

Badamy system po to, żeby na niego oddziaływać. Musimy stworzyć opis/model takiego urządzenia, aby z niego korzystać.

Modelowanie: Przebadać model empirycznie, jak coś nie gra zmienić program, uprościć, skomplikować model. Im więcej powiązań z otoczeniem, z jakim dany obiekt funkcjonuje, tym bardziej skomplikowany (skomplikowane algorytmy itp.).

0x08 graphic

Tworzenie modeli:

I Przykład(hipotetyczny.)

x- ilość produkcji y- koszt

y = F(x) - model tego systemu (ekonomicznego)

0x08 graphic

II Przykład: y1 = F(x1,x2)

2 wierze destylacyjne, y2 = F(x1,x2)

które mogą robić Podanie takiej

to samo zależności WE/WY

tj. podanie modelu

y1=a11x1 + a12x2

y2=a11x1 + a22x2

Czy jest to model dobry?

Nie wiem, będę musiał sprawdzić: weryfikacje, modulacje.

Ale uznaje, że ten model jest dobry. Jest to model liniowy. Jest to macierz kwadratowa, lub prostokątna w zależności od ilości WE i WY.

Model Statyczny - nie występuje czas, interesuje mnie obserwacja w danym momencie, jak stany znikną podczas badania, to efekt nas interesuje.

Y=Ax

0x08 graphic
Model dynamiczny - występuje czas, stan aktualny zależy od stanów poprzednich.(np. wzrost produkcji i dynamika wzrostu)

a1,a2 > 0

Przykład : Model matematyczny opisujący prąd przemienny.

0x08 graphic

Do opisów modeli dynamicznych zawsze będzie potrzebne równanie różniczkowe

WYKŁAD II

IDENTYFIKACJA, ROZPOZNAWANIE

(b. Specyficzna identyfikacja)

y1 = a11x1 + a12x2

y2 = a21x1 + a22x2

Mamy nieznane parametry:

a11=? a12=? a21=? a22=?

Algorytm identyfikacji - sposób na wyznaczenie tych parametrów.

Eksperyment:

BIERNY

AKTYWNY

W układzie krwionośnym nie można np. podwyższyć ciśnienia u pacjenta, bo umrze,

Można u pacjenta zaobserwować ciśnienie np. dziś i jutro, mam 2 wyniki i je przetwarzam.

Jeżeli sami możemy ustalić, zaplanować wejścia. Wykonuje doświadczenie. Np. Sam ustawiam WE.

Identyfikacje przeprowadzamy w przypadku, kiedy znamy dokładny model systemu.

Przykład:

y1 = a11x1 + a21x2 X1,X2 - ilość surowców przerabianych w 2 agregatach

y2 = a21x1 + a22x2 Y1,Y2 - ilość produktu i kosztu produkcji

Y = AX

A = YX-1 - algorytm identyfikacji

A = Ψ (x1,x2,y1,y2)

_

y = Ø (x, a) zakładam sobie model z dokładnością parametrów

Ilość obserwacji WY zew. obiektu

0x08 graphic
↓ ↓

n _ n

QI =∑ ( yi - yi )2 =(mają nie być różne od siebie)= ∑ ( yi - Ø(a1 xi ))2

i=1 i=1

WY modelu (wylicze sobie)

Min Q(a) → a*

a

Obliczam I pochodną i sprawdzam I lub II warunek wystarczający istnienia min funkcji.

0x08 graphic
_

y = ax

_

y = ax2 - klasa modeli kwadratowa Obliczanie sumy

odl. od prostej

pomiary i-tego a

n n

QI = ∑ (y1 - Ø(a1xi))2 = ∑ (y1 - ax1)2

i=1 i=1

dQ n

― = ∑ 2(y1 - ax1) (-xi) = 0

da i=1

n

∑ (axi2 - xiyi) = 0

i=1

n n

a ∑ xi2 - ∑ xiyi = 0

i=1 i=1

n

xiyi Algorytm identyfikacji wygląda dokładnie tak.

i=1 Model statyczny, liniowy.

a = ------------- Najlepsze rozwiązanie w miare minimalizacji danego kryterium

n

∑ xi2

i=1

ANALIZA

Analiza (Systemu) ilościowa:

Pytanie: Jak zachowuje się WE, kiedy zmieni się WE?

Dokonuje analizy działania systemu. Jak zachowuje się WY dla różnych WE, tj. zadanie ANALIZY.

W przypadku prostych układów wyliczenie y znając konkretne wartości parametrów.

DEF.

„Jeżeli znamy wartości liczb parametrów a11,a12,a21,a22 to wstawiamy dane x1,x2 do modelu i wyliczamy y1,y2.”

Analiza jakościowa:

Czy system ma pewne własności np. czy jest stabilny(pewne wielkości nie będą dążyć do nieskończoności).

PROJEKTOWANIE

Projektowanie - Podejmowanie decyzji

Czy mogę ustawić y i dobrać x, żeby y były takie jak ja sobie życzę ?

Zaprogramować system - przewidzieć co dać na WE aby WY osiągnęło taki stan jak ja sobie życzę.

DEF.

„Dla danych wartości a oraz zadanych wartości y1=y1* , y2=y2* można wyliczyć decyzje x1,x2

y1*=a11x1+a12x2 Musze rozwiązać ten układ względem tych 2 niewiadomych, I wtedy y

y2*=a21x1+a22x2 osiągnie na pewno stan taki jak chce.

y1*a22 - y2*a21 y2*a11 - y1*a21

0x08 graphic
0x08 graphic
x1 = x2 =

a11a22 - a21a12 a11a22 - a21a12

STEROWANIE/KIEROWANIE

Co by się stało, gdyby y1 i y2 zmieniło się?

Codziennie ustala się np. rozmiar produkcji i koszty) codziennie się zmieniają i wtedy musze przetestować „system”, musze nim kierować.

DEF.

Reakcja na zmieniające się parametry - kierowanie/sterowanie

**Pyt. Kolokwialne:

Wyjaśnij zadanie projektowania a zadanie sterowania na prostym przykładzie.

ETAPY MODELOWANIA MATEMATYCZNEGO

1.Sformułowanie celów modelowania.

2.Wybór kategorii modelu i określenie jego struktury.

3.Identyfikacja.

4.Algorytmizacja obliczeń.

5.Weryfikacja.

SFORMUŁOWANIE CELÓW

Pytanie po co to robię? Co chce zbadać, wykryć? Do czego będzie mi to służyć?

1.Opis i wyjaśnienie działania mechanizmu systemu - model FENOMELOGICZNY

2.Przewidywanie zachowania systemu w przyszłości i przy różnych warunkach, oddziaływanie na system - model PROGNOSTYCZNY.

3.Wybór właściwych oddziaływań WE spełniających określone warunki - model DECYZYJNY.

4.Wybór struktur, lub parametrów spełniającego określone zadania - model NORMATYWNY.

WYBÓR KATEGORII MODELU I OKREŚLENIE JEGO STRUKTURY

Zależy od natury modelu: statyczny, dynamiczny...

Każdy system, obiekt jest powiązany z otoczeniem. Dlatego trzeba uwzględnić wpływ otoczenia. Trzeba pójść na ustępstwa: np. nie bierzemy pod uwagę pogody. Im więcej pominiemy, tym model będzie prostszy, ale nie może być za prosty.

Trzeba osiągnąć kompromis pomiędzy: komplikacją modelu, a zbytnim upraszczaniem.

IDENTYFIKACJA

Model zbyt skomplikowany, jest złym rozwiązaniem. Może być bardzo dokładny, ale przez tą dokładność mam zbyt mało danych pomiarowych, pomiarowych model ten jest niemożliwy do skonstruowania, więc wracam się i upraszczam go.

ALGORYTMIZACJA OBLICZEŃ

W końcu muszę cos obliczyć.

Np. Więcej niewiadomych niż równań (nieskończenie wiele rozwiązań).

Jak przystępuje do obliczeń może wystąpić kilka różnych rozwiązań:

1.Nieskończenie wiele rozwiązań : Źle

2.Sprzeczny(wróć i wykryj błąd) : Źle

3.Jedno rozwiązanie : Dobrze

WERYFIKACJA

Czy da się policzyć?

Czy mam określony rodzaj zgodności?

Czy wyniki doświadczalne zgadzają się z wynikami rzeczywistymi?

Itp....

Jeżeli któryś z warunków się nie zgadza to wracam do poprzednich etapów

DEF.

„Jest to porównanie wyników modelowania z zachowaniem się systemu

rzeczywistego.”

Kryteria Wewnętrzne:

- Zgodność FORMALNA - czy jest to układ zdań poprawnych i spójnych

- Zgodność ALGORYTMICZNA - Jak są np. 2 rozwiązania, to trzeba w razie potrzeby mieć zaszyte mechanizmy, które pozwolą wybrać jedno z możliwych rozwiązań.

Kryteria Zewnętrzne:

- Zgodność HEURYSTYCZNA - czy zbudowany model ma walory interpretacyjne/naukowe, że wyjaśnia działanie systemu. Nabieranie wiedzy o systemie. Model oddaje, interpretuje mechanizmy, które w nim zachodzą.

- Zgodność PRAGMATYCZNA

- Zgodność REPLIKATYWNA - Jeżeli zbudowany model obiektu ma podstawy obiektu, to ten model będzie działał podobnie jak system na podobnych warunkach.

- Zgodność PREDYKATYWNA - Aby model dobrze działał , jak pojawią się dane, których nie było podczas sprawdzania (np. w modelach ekonomicznych). Gdybanie, Przewidywanie.

- Zgodność STRUKTURALNA - Aby model dokładnie oddawał mechanizmy, które rządzą tym modelem, obiektem.

TERAZ MOŻNA ZAPEWNIĆ SYMULACJE KOMPUTEROWE, KIEDYŚ WSZYSTKO NA PIECHOTE, B.ŻMUDNE OBLICZENIA.

KONIEC INFORMACJI OGÓLNYCH

POJĘCIE SYSTEMU

1.Idea wyodrębniania systemu z otoczenia.

2.Idea budowy systemu z elementów(Podsystemy).

3.Idea funkcji spełnianej przez system.

4.Idea ograniczonej zmienności systemu.

Ad 1.

System jest pewną całością, w której współdziałają wyodrębnione części składowe. Funkcjonowanie systemu zależy od funkcji części składowych, i związków między nimi, powiązania części składowych określają strukturę systemu.

Np. samochód jest systemem.

Ad.2-3.

System jest pewną całością i spełnia wspólny jakiś określony cel.

Ad.4.

Formuły matematyczne jakie przyporządkujemy dla systemu, tylko przez jakis czas będą aktualne.

WYKŁAD III

ETAPY PROJEKTOWANIA SYSTEMU STEROWANIA

1.Analiza systemu O.S.(modelowanie).

2.Identyfikacja O.S.

3.Opracowanie algorytmu sterowania.

4.Opracowanie komputerowego programu sterującego.

5.Opracowanie systemu realizującego program sterujący.

Analiza systemu ma uzmysłowić, że mamy działać celowo. Sprecyzować problem, zrobić tak,

Aby wszystko działało, a następnie stworzyć.

Modelem matematycznym systemu(obiektu, procesu) jest zestaw wzorów matematycznych. Ogólnie zestaw relacji matematyczno-logicznych, określających zależności pomiędzy wyróżnianymi wielkościami, są to wielkości WE/WY systemu, a system taki nazywany względnie odosobniony. W pojęciach WE i WY zawarte jest domniemanie o związku przyczynowo-skutkowym, między tymi wielkościami. wielkościami tym sensie, że obserwując WE można określić WY. Zdarzenia nie są losowe, tylko zaplanowane.

0x08 graphic

System wyodrębniony z otoczenia:

Otoczenie oddziałowywuje, na system, poprzez WE( przykład

takiego systemu- auto: WE- naciskanie

gazu, a WY tj, prędkość)

DYNAMIKA SYSTEMÓW

Własności dynamiczne systemu wynikają z działania następujących ogólnych zasad fizycznych:

- przy ograniczonych wydajnościach źródeł, każda nieskończenie mała zmiana stanu energetycznego, lub materiałowego wymaga pewnego czasu( bezwładność inercja).

- każde skończone przemieszczanie się w przestrzeni zjawiska materialnego wymaga pewnego czasu - opóźnienie.

DYNAMIKA - (nie szybkość, nie działalność niezrównoważonych sił) Stan obiektu/systemu, nie zależy, od tego, co jest na WE w danym momencie, tylko od historii.

//*Przykład I: Ilość biomasy marchewki, uzależnione jest od słońca i deszczu, czy jest duża zależy od historii, a nie od tego, czy w momencie badania świeci słońce, czy nie.//

//*Przykład II: Odpalenie Windowa, stan komputera, aby wykonywała się dana funkcja windowa, np. logowanie, musi nastąpić szereg instrukcji poprzedzających, bez których nie dojdzie do logowania//

//*Przykład III: charakter narastania temperatury w czasie(PIEC). Im dajemy wyższe napięcie, tym piec się szybciej nagrzewa, ale i tak jakiś czas to trwa. Czas nagrzania nie może być = 0

Zjawiska dynamiczne:

- charakter bezwładny( uczenie się matematyki: 1 impuls(„ucz się”), bezwładność, opóźnienie, potem 2 impuls(„dostajesz pałke”), zaczyna trybić).

- opóźnienia: Jazda autem, zauważam przeszkodę, nie staje od razu, tylko zajmie to trochę czasu: czas reakcji, wytracenie energii. Cały proces jest opóźniony.

DEF.

Miejsce, gdzie inst.. WE i miejsce WY są rozłożone w czasie.

0x08 graphic
0x08 graphic

Pokonanie przestrzeni czasu zawsze wymaga jakiegoś czasu.

Matematycznie - różne stany przejściowe, lub stany nieustalone

PIEC

0x08 graphic

0x08 graphic

Są to tzw. Przebiegi przejściowe

Tutaj już mniejsza bezwładność (to już nie pasuje do pieca)

Układ będzie dążył do jakiegoś układu ustabilizowanego, lub nie

Magistrala - układ dynamiczny(generator) o stabilnej częstotliwości

DEF.

Precyzyjna definicja matematyczna związana z dynamiką:

Rys wyjaśnienie:

0x08 graphic

Historia WE ma wpływ do tego momentu, ale bez tego punktu,

Czyli C ma wpływ bez E.

Wartość WE w B nie ma wpływu na wart WY w A (ukł. Dynamiczny), ale na wart. WY ma wpływ cała Historia - C, lub

Tylko punk w Historii - P.

To że czas występuje w danym systemie, nie oznacza że układ jest układem dynamicznym, natomiast kiedy czas nie występuje to na pewno układ NIE jest dynamiczny.

Stan systemu dynamicznego :

Jest to najmniejsza liczba danych(wektor stanu), których znajomość danej chwili, przy znajomości wielkości WE, począwszy od tej chwili pozwala jednoznacznie określić stan i wielkości WY systemu w przyszłości.

Zdeterminowanie systemu wynika z postulatów: przyczynowości, zupełności.

0x08 graphic

Ilość zmiennych stałych, które określają położenie danego

Obiektu jest różna.

1 i 2 wpływa na 3 i 4

Więc systemy są zdeterminowane, znając stan WE można

Obliczyć WY.

Postulat przyczynowości:

0x08 graphic

Not: Dopóki WE stanu jest takie same, to sterowanie tez takie

same => przyczyna zmiany stanu.

Postulat zupełności:

0x08 graphic

Not: Zmian jest jednoznaczna w odwrotną stronę. Kiedy WE stanu jest takie same, to sterowanie różne, natomiast jak sterowanie takie same, to WE stanu też

Przykład:

Zbiornik wody - stopień zapełnienia:

Dx(t)

------ = u(t) - y(t)

dt

↑ ↑ ↑

WE WY

prędkość zmiany co wpływa do to co wypłynęło, już

stanu(pochodna) wykonania wypływa

(tutaj woda,

deszcz)

y(y) = y(x(t))

wyznaczam x(t):

t t

x(t) = ∫ u(t)dt - ∫ y(t)dt

Całka oznaczona - suma wielkości, sumowanie przebiegu, który zmienia się w czasie

Całka - sumowanie, czyli uwzględnienie tego , co było po drodze.

WYKŁAD IV

METODY OPISU CIĄGŁYCH LINIOWYCH JEDNOMIAROWYCH OBIEKTÓW STEROWANIA

  1. Równania różniczkowe opisujące jednowymiarowe obiekty sterowania ( w szczególności WE-WY).

  2. Odpowiedź jednostkowa (skokowa).

  3. Odpowiedź impulsowa.

  4. Transmitacja operatorowa.

  5. Transmitacja widmowa.

  6. Metoda przestrzeni fazowej i płaszczyzny fazowej.

  7. Metoda zmiennych stanu.

  8. Charakterystyki częstotliwościowe.

2.Przedstawia jak zachowuje się obiekt kiedy na WE pojawia się sygnał skoku (wymuszenie).

3.Odpowiedź obiektu na krótkotrwały impuls(np. doprowadzam auto do jakiejś prędkości, przychodzi podmuch wiatru i wytrąca, muszę ulepszyć, zaplanować sterowanie.

4.Transformata Laplace'a.

5-6-7-8.Bardzo złożone.

Dlatego tak dużo opisów, bo różnie są wygodne w różnych analizach.

Matematyczna zmiana stanu - wiemy nie tylko co dzieje się na zewnątrz, ale także wewnątrz.

0x08 graphic

Charakterystyką(Odpowiedzią) skokową:

H(t) jednowymiarowego układu (obiektu) liniowego stacjonarnego, nazywać będziemy odpowiedź tego układu na wymuszenie w postaci jednostkowej funkcji skokowej 1(t), przy zerowych warunkach początkowych.

0x08 graphic

←impuls

Charakterystyką( Odpowiedzią) Impulsową

G(t)liniowym układu( Obiektu) liniowego, stacjonarnego nazywać będziemy odpowiedź tego układu na wymuszenie w postaci funkcji Diraca (t), przy zerowych warunkach początkowych.

0x08 graphic
δ(t) ←impuls o nieskończonej amplitudzie, duży impuls,

o dużym czasie trwania.

Transmisją operatorową:

G(s) liniowym układu (Obiektu) liniowego stacjonarnego nazywać będziemy wartość określoną jako stosunek transformaty wymuszenia U(s) tego układu przy zerowych warunkach początkowych. G(s) =Y(s) / U(s)

0x08 graphic
TYPOWE CZŁONY LINIOWE

- Człon proporcjonalny

- Człon inercyjny I rzędu

- Człon całkujący idealny( na WY jest całka)

- Człon całkujący z inercją( to samo, tylko zaburzone)

- Człon różniczkujący( na WY pochodna sygnału WE)

- Człon różniczkujący z inercją( to samo ale zaburzony)

- Człon oscylacyjny( złożenie członów inercyjnych inercyjnych rzędu)

- Człon oscylacyjny II rzędu

- Człon opóźniający( to co na WE opóźnione na WY)

CZŁON BEZINERCYJNY

( proporcjonalny wzmacniający)

Równanie różniczkowe WE/WY: y= ku

k - współczynnik wzmocnienia określony jako stosunek odpowiedzi y do wymuszenia w stanie ustalonym.

T - stała czasowa = 0

Charakterystyka skokowa: g(t) t

h(t) = k1(t)

Charakterystyka impulsowa: h(t)

h(t) = k δ(t) t

Transmitacja:

Gi(s) = k

Kierunek przepływu prądu

0x08 graphic

**Przykład 1: Wybieram wielkości WE i WY CZWÓRNIK czwórnika, zakładam jeszcze, że (Dzielnik jest nieobciążony.

Napięcia)

I(t) - prąd płynie, wiedzieć to musze, aby znać zależności między Uwe(t) i Uwy(t).

0x08 graphic

Uwe(t) = I(t) * (R1 + R2)

Uwy(t) = I(t) * R2 => I(t) = ....... wyliczam I(t) i podstawiam

R1 + R2

Uwe(t) = Uwy(t) Schemat postępowania: Wypisywanie równań opisuje

R2 opisuje zachowanie się układu.

R2

Uwy(t) = Uwe(t)

R1 + R2

**Przykład 2:

PRĄDNICA TACHOMETRYCZNA PRĄDU STAŁEGO

(Jazda autem prędkościomierz przekłada prędkość obrotu kół na napiecie i pokazuje na prędkościomierzu).

Tu przychodzi strumień

0x08 graphic
0x08 graphic
Wzbudzenia(prądnica)

Nieobciążona, poruszam się w liniowym zakresie magnesowania

E(t) = U(t)

E(t) = CE Ø Ω(t)

U(t) = CE Ø Ω(t)

0x08 graphic

y(t) = k * U(t)

Najprostrzy człon proporcjonalny, bezinercyjny

CZŁON INERCYJNY:

(piec)

(Jak funkcja określona w czasie

I - pochodna mówi o prędkości zmian

II - pochodna o przyspieszeniu)

.

Równanie różniczkowe WE/WY: Ty + y = ku

k - współczynnik wzmocnienia określany jako stosunek odpowiedzi y do wymuszenia u w stanie ustalonym.

0x08 graphic
T - stała czasowa.

Charakterystyka skokowa:

h(t) = k( 1 - ex [(-1)/T]) 1(t)

Odpowiedź impulsowa:

0x08 graphic

g(t) = K/T ex [(-1)/T] 1(t)

Transmitacja:

k

0x08 graphic
Gs =

1 + sT

kierunek przepływu prądu

0x08 graphic
**Przykład 1:

CZWÓRNIK RC

Nieobciążony.

Uwe(t) = RI(t) + 1/C ∫ I(t) dt

Uwe(t) = RI(t) + Uwy(t)

Uwy(t) = 1/C ∫ I(t) dt

d Uwy(t) 1

0x08 graphic
= * I(t) tj. funkcja czasu

dt C

RC d Uwy

Uwe(t) = . * + Uwy

T y + y dt

Oznacza to, że jak włącze napięcie na WY nie pojawi się niezwłocznie, tylko będzie narastać.

Im stała czasowa mniejsza, tym szybciej zostanie osiągnięty stan przez WY.

**Przykład 2:

ELEMENT GRZEWCZY:

P - moc jest funkcją czasu P(t); energia dostarczana do grzejnika

dt - mały czas (delta t) energia ta rozkłada się na ciepło

dodawane i wydalane

P(t)dt = c*m *d Θ(t) + α 0x08 graphic
S Θ dt α - stała

S - powierzchnia grzejnika

Θ - ilość ciepła proporcjonalna do

To co się kumuluje to co się rozprasza temp, zmienia się w czasie

w grzejniku

Mamy bilans cieplny:

P(t)dt = c*m*d Θ(t) + αS*Θ(t) dt

1 c*m d*Θ(t) .

0x08 graphic
0x08 graphic
0x08 graphic
* P(t) = * + Θ(t) k*u(t) = T * y + y

α*S α * S dt

**Przykład 3:

0x08 graphic
GENERATOR

PRĄDU STAŁEGO

Nieobciążony

Iu(t) - prąd wzbudzenia

d*I*w(t)

0x08 graphic
L + R*I*w(t) = Uw(t)

dt

d*E(t) d*I(t)

0x08 graphic
0x08 graphic
= K

dt dt

L d*e(t) R

0x08 graphic
0x08 graphic
0x08 graphic
* + E(t) = Uw(t)

K dt K

Podsumowanie:

Co uzyskujemy?

Układy o różnych naturach podlegają pod te same równania jak odrzuce skale i jednostkę MATEMATYKA jest wszędzie.

Wykład V

CZŁON RÓŻNICZKUJĄCY

.

Równanie różniczkujące WE/WY: y = ku

K - współczynnik wzmocnienia określony jako stosunek odpowiedzi y do pochodnej wymuszenia w stanie ustalonym

T - stała czasowa = 0

Charakterystyka skokowa: h(t) = uδ(t)

dδ(t)

0x08 graphic
Odpowiedź impulsowa: g(t) = k (pochodna z odpowiedzi skokowej)

dt

Transmitacja: g(s) = ks

Przykład:

(nie ma idealnych członów różniczkujących)

0x08 graphic
dU(t)

0x08 graphic
KONDENSATOR i(t) = c

dt

Uc(t) = 1/c ∫ i(t) dt

CZŁON RÓŻNICZKUJĄCY Z INERCJĄ

(rzeczywisty człon różniczkujący)

. .

Równanie różniczkujące WE-WY: Ty + y = ku

K - współczynnik wzmocnienia określany jako odpowiedź y do pochodnej wymuszenia u w stanie ustalonym

T - stała czasowa

k 1

0x08 graphic
0x08 graphic
Charakterystyka skokowa: h(t)= exp ( - ) 1(t)

T T

Odpowiedź impulsowa: k k 1

0x08 graphic
0x08 graphic
0x08 graphic
g(t)= δ(t) - exp( - ) 1(t)

T T2 T

ks

0x08 graphic
Transmitacja: g(s) =

1+sT

0x08 graphic
Przykład:

CZWÓRNIK R,L d I(t)

0x08 graphic
Uwe(t) = R * I(t) + L =

dt

0x08 graphic

= R I(t) + Uwy(t)

dI(t)

0x08 graphic
Uwy(t) = L

dt

d Uwe(t) d I(t) d Uwy(t)

0x08 graphic
0x08 graphic
0x08 graphic
= R +

dt dt dt

Równanie końcowe:

d Uwe R d Uwy . .

0x08 graphic
0x08 graphic
0x08 graphic
= * Uwy(t) + => u = k * y + y

dt L dt

Przykład **:

TRANSFORMATOR :

0x08 graphic

0x08 graphic
d I1

0x08 graphic
U1 = RI1 + L1

0x08 graphic
dt .

0x08 graphic
I

d I1

0x08 graphic
U2 = M

dt

. . ..

U1 = R1 I1 + L1 I

. R L1 .

0x08 graphic
0x08 graphic
U1 = U2 + U2 / M/R

M M

0x08 graphic
0x08 graphic

M . L1 .

0x08 graphic
0x08 graphic
U1 = U2 + U2

R R1

K T

CZŁON CAŁKUJĄCY IDEALNY

(wykonuje operację odwrotną) :

.

Równanie różniczkowe WE - WY: y = ku

0x08 graphic
K - współczynnik wzmocnienia prędkościowego, określany jako stosunek pochodnej odpowiedzi y do wymuszenia u w czasie ustalonym.

Charakterystyka skokowa (stała): h(t) = k t 1(t)

Odpowiedź impulsowa (pochodna): g(t) = k 1(t)

K

0x08 graphic
Transmitacja: g(`s) =

S'

Przykłady:

ZBIORNIK Z CIECZĄ

We - objętość cieczy dopływająca jednostce czasu ( g(t) )

0x08 graphic
Wy - poziom cieczy

0x08 graphic
0x08 graphic
G(t) d V(t)

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
G(t) =

dt

str. dopływający Szybkość narastania

objętości

wys. Słupa wody: h(t) * S = V(t)

d h(t)

0x08 graphic
S = g(t)

dt

t

h(t) = 1/s ∫ g(t) dt

0

2. Kondensator idealny.

3. Układ napędowy pozycyjny.

CZŁON CAŁKUJĄCY Z INRCJĄ

.. .

Równanie różniczkowe WE - WY: Ty + y = ku

K - współczynnik wzmocnienia prędkościowego określany jako stosunek pochodnej odpowiedzi y do wymuszenia u w stanie ustalonym

T - stała czasowa

Charakterystyka skokowa:

0x08 graphic
t

0x08 graphic
0x08 graphic
h(t) = (kt - kT ( 1-exp [- ] ) ) 1(t) układ będzie dążył do całki

idealnej.

T

Odpowiedź impulsowa:

t

0x08 graphic
h(t) = k [ 1-exp [- ] ] 1(t)

T

K

0x08 graphic
Transmitacja: g(s) =

S(1+sT)

Przykłady:

1. OBCOWZBUDNY SILNIK PRĄDU STAŁEGO

Założenia: silnik nieobciążony, uwzględniamy tylko bezwładność

We - napięcie twornika

Wy - kąt położenia wału silnika

0x08 graphic

Ω(t) (WY - położenie kątowe)

α(t)

0x08 graphic

0x08 graphic
Ø - strumień jednego z uzwojeń

Const

0x08 graphic
0x08 graphic

DI(t)

0x08 graphic
U(t) = Rt * I(t) + Lt E(t)

0x08 graphic
dt

pomijam(małe)

dΩ(t)

0x08 graphic
Mn(t) - Mz(t) = I Mn(t) = Cm * Ø * I(t)

dt

Chcemy powiązać 2 wielkości: α(t) i U(t)

d α(t)

0x08 graphic
= Ω(t) No i manipulujemy, aby była zależność α(t) do U(t)

dt

dΩ(t)

0x08 graphic
Cm Ø Ω(t) * I(t) = J

dt

U(t) = Rt * I(t) + CE Ø Ω(t)

J dΩ(t)

0x08 graphic
0x08 graphic
I(t) = *

CM Ø dt

Rt * J d2 α dα

0x08 graphic
0x08 graphic
0x08 graphic
U(t) = * + CEØ

CmØ dt2 dt

CZŁON INERCYJNY II RZĘDU

Równanie różniczkowe: T1T2y'' + T1yt' + y = ku

Lub

TATBy'' + (TA+TB) y' + y = ku

K - współczynnik wzmocnienia określony jako stosunek odpowiedzi y do wymuszenia u w stanie ustalonym.

T1T2TATB - stałe czasowe.

K

0x08 graphic
Transmitacja: G(s) =

T1T2S2 + T1S + 1

Lub

K

0x08 graphic
G(s) =

TATBS2 + (TATB) s + 1

T1 - 4T2 >=0 __

____ 1 √T1

0x08 graphic
0x08 graphic
TB = √T1T2 , S = __

2 √T2

_____ TA + TB

0x08 graphic
T0 = √TATB , ζ = _____ , K = KA - KB

2 √TATB

KA

0x08 graphic
G1(s) =

1 + STA

KB

0x08 graphic
G2(s) =

1 + STB

0x08 graphic
Przykład:

CZWÓRNIK R,C,

0x08 graphic
dI(t) 1 t

0x08 graphic
0x08 graphic
Uwe(t) = R * I(t) + L + ∫ I(t)dt

dt c 0

1 t dUwy

0x08 graphic
0x08 graphic
Uwy(t) = ∫ I(t) dt c = I(t) różniczkuje raz jeszcze

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
C 0 dt

dUwy d2Uwy d2Uwy dI

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
Uwe(t) = RC * + L C + Uwy C = =

0x08 graphic
dt dt dt dt

To czy układ będzie inercyjny II rzędu, czy oscylacyjny zależy od parametrów

CZŁON OPÓŹNIAJĄCY

Równanie różniczkowe: y(t) = ku( t - T0)

K - współczynnik wzmocnienia określony jako stosunek odpowiedzi y do wymuszenia u dla

t >T0

T0 - czas opóźnienia

Charakterystyka skokowa: h(t) = k 1(t - T0)

Odpowiedź impulsowa: g(t) = kδ(t - T0)

Transmitacja: g(s) = ex (-sT0)

CZŁON OSCYLACYJNY

Odpowiada na wymuszenie

0x08 graphic

1. Tłumione

2.Nie Tłumione

3. Narastanie.

TRANSMITACJA OPERATOROWA:

Transformata Laplace'a:

F(s) = ∫ f(s) e-st = L[ f(t) ]

0

Orginały Laplace'owskie ( f(t) będzie orginałem, jak spełni 3 warunki):

1. f(t) = 0 t<0 Obserwacje prowadzimy w przód, nie bierzemy pod uwagę

przeszłości.

2. Funkcja ciągła i monotoniczna przedziałami.

3. f(t) <= m* e αt Można zrobić na funkcjach , które rosną, ale nie szybciej jak funkcja wykładnicza( np. et)2 nie da się jej ograniczyć);

Własności transformaty Laplace'a:

1. Złożenie dwóch operacji

L-1[L(f(t))] - f(t) , L[L-1(F(s))] = F(s)

α(1(t)) = 1/s , α-1(1/s) = 1(t)

1 1

0x08 graphic
L( e-αt * 1(t)) = , L-1 ( ) = e-αt * 1(t)

S+α S+α

Funkcje określam dla czasu T = 0

2. Liniowość:

Obowiązuje zasada superpozycji( odpowiedź na wymuszenie jest odpowiedzią na sumę tych wymuszeń.

L [a1f1(t) + a2f2(t)] = a1L[f1(t)] + a2L[f2(t)]

Przykłady:

( ejωt = cos ωt + j sin ωt)

( e-jωt = cos ω - j sin ωt)

1

0x08 graphic
sin ωt = ( ejωt - e-jωt)

2j

1 1 1 1 2jω ω

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
L { sin ωt * 1(t) } = ( - ) = * =

2j s-jω s+jω 2j s2 + ω2 s2 + ω2

ω

0x08 graphic
L { sin ωt } =

s2 + ω2

4

0x08 graphic
sin 4t =

s2 + 16

1 3

0x08 graphic
0x08 graphic
L-1 ( ) = L-1 ( ⅓ * ) = ⅓* sin 3t

s2 + 9 s2 + 9

3. Holomorficzność (funkcja, która ma pochodne):

dF(s) dF(s)

0x08 graphic
0x08 graphic
= -∫ t f(t) e-st dt czyli L[t f(t)] = -

ds. ds.

oraz

n! n!

0x08 graphic
0x08 graphic
L(fn) = , L(ten-αt) =

Sn+1 (s + α)n+1

L { t*1(t)} = 1/s

-1 1

0x08 graphic
0x08 graphic
- =

s2 s2

s+2 1 3 1 1 1 3 1

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
L-1 ( = L-1 (-2 * + * + * ) = (-2 + et+ e-t)1(t)

s(s2 - 1) s 2 s-1 2 s+1 2 2

3. Transformata pochodnej:

ω s

0x08 graphic
0x08 graphic
L (cos ωt) = L (1/ω(sin ωt)') = 1/ω *(s*L(sin ωt) - sin 0 ) = 1/ω * s * =

s2 + ω2 s22

0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
0x08 graphic
Idea rachunku operatorowego

0x08 graphic
Mamy problem w dziedzinie czasu

Równanie L równanie wyznaczanie L wyznacznie

Dla orginału dla obrazu obrazu oryginału

Przykład: dI(t)

0x08 graphic
0x08 graphic
U *1(t) = R*I(t) + L

dt

0x08 graphic
równanie obrazu bo oper, różniczkowemu

odpow. mnożenie przez s.

*U1/s = R * I(s) + L*s*I(t)

u 1

0x08 graphic
0x08 graphic
I(s) = *

0x08 graphic
S R*SL

Rozw w dziedzinie Laplace'a

Faza końcowa - wyznaczenie oryginału.

1 K/R U/R

0x08 graphic
0x08 graphic
0x08 graphic
= 1/s + = -

s+R/L S s + R/L

0x08 graphic

I(t) = U/R * (1 - e-(R/L)*t) * 1(t)

pisane



Wyszukiwarka