WYKŁAD 02 SterowCyfrowe

background image

WYKŁAD 2

STEROWANIE CYFROWE

JANUSZ KWAŚNIEWSKI AGH Katedra Automatyzacji Procesów

background image

2

Ciągłe a dyskretne

Równanie różniczkowe

Równanie różnicowe

Transformata prosta
i odwrotna Laplace’a

Transmitancja ciągła

Transmitancja dyskretna

Dyskretne równania stanu

Ciągłe równania stanu

DYSKRETYZACJA

DYSKRETYZACJA

DYSKRETYZACJA

Transformata prosta

i odwrotna z

Na tę konwersję ma wpływ:

• wielkość kwantu
• czas próbkowania
• metoda dyskretyzacji

Dotyczy tylko

układów liniowych

Dotyczy tylko

układów liniowych

background image

3

Analogowy, próbkowany, dyskretny

Powyższy przykład dotyczy 3 bitowego przetwornika

czyli o 2

3

=8 poziomach wartości sygnału

Dyskretny sygnał na wyjściu ma największe

różnice po próbkowaniu i kwantowaniu

z

z

2

z

3

z

-1

z

-2

z

-3

background image

4

Dobór czasu próbkowania

1/6 w praktyce 1/10
dominującej stałej czasowej

Efekt stroboskopowy, (aliasing)

background image

5

Część cyfrowa układu regulacji

O dokładności decyduje:

Liczba stanów (kwantów), n-bitowe przetworniki A/C i C/A,

Okres próbkowania T,

Metoda konwersji

Współczynniki transmitancji G

s

(z) (lub równań różnicowych)

zawierają w sobie

parametry konstrukcyjne

okres próbkowania

STEROWNIK

y

(k)

e

(k)

-

T

T

Zegar (T)

REGULATOR

G

s

(z

)

w T

R

y

0

y(t)

m(t) PROCES

G

p

(s)

T

Interpolator ZOH

(1-e

-sT

) / s

Przetwornik A/C

Przetwornik

C/A

Interpolator ZOH

(1-e

-sT

) / s

Przetworn

ik A/C

LP

Filtr

Czas wykonania algorytmu
regulatora T

R

«T powinien być dużo

krótszy od okresu próbkowania

background image

6

Transmitancja dyskretna a równanie
różnicowe

Ogólna transmitancję dyskretną (impulsową funkcję przejścia)
obiektu ma postać

(pamiętaj, że z z

-1

=1 a z

2

=-1):

Z uwzględnieniem czasu opóźnienia (ang. dead - time, transport

-

delay) o nk okresów próbkowania:

Przyjazna do programowania nosi nazwę rekurencyjnego

równania różnicowego (model ARMA):

G z

y

u

B z

A z

b

b z

b z

a z

a z

i

i

nb

nb

na

na

( )

(

)

(

)

 

 

1

1

0

1

1

1

1

1

na

na

na

nb

nb

nb

a

z

a

z

b

z

b

z

b

z

A

z

B

1

1

0

...

...

)

(

)

(

G z

z

B z

A z

nk

( )

(

)

(

)

1

1

y

i

= - a

1

y

i-1

- a

2

y

i-2

-...- a

na

y

i

- nk

-1

+ b

0

u

i

-

nk

+b

1

u

i - nk -

+...+b

nb

u

i - nk - nb

Rzeczywisty system jednowymiarowy opisywany zarówno modelem analogowym jak
cyfrowym (ciągiem czasowym) można nazwać również filtrem lub przetwornikiem
danych. Wynika to z tego, że mamy tu do czynienia z zasadą jednoznacznego
przekształcenia sygnału wejściowego w wyjściowy.

Odpowiedz się nie zmienia
tylko przesuwa o nk
chwil
czasowych

background image

7

Transmitancja z operatorem
opóźnienia i wyprzedzenia

mnożymy licznik i mianownik przez przeciwny najwyższy operator

)

5

,

0

(

)

2

(

)

5

,

0

1

(

)

2

1

(

1

1

z

z

z

z

Zero = - 2 a biegun = - 0,5,

Czyli obiekt jest stabilny ale
nieminimalnofazowy, statyczny

)

6

)(

5

,

0

(

)

5

(

3

5

,

6

)

5

(

3

5

,

6

1

5

)

6

1

)(

5

,

0

1

(

)

5

1

(

2

2

1

2

1

1

1

1

1

z

z

z

z

z

z

z

z

z

z

z

z

z

z

Jedno zero znika, a drugie zero z=5

i bieguny nie zmieniają się

ale obiekt jest cały czas niestabilny i
nieminimalnofazowy, statyczny

)

6

)(

1

(

)

5

(

6

5

5

6

5

1

5

1

)

6

1

)(

1

(

)

5

1

(

2

2

2

1

1

1

1

1

z

z

z

z

z

z

z

z

z

z

z

z

z

z

j.w.
astatyczn
y

background image

8

Faktoryzacja transmitancji
dyskretnej

Stabilność i minimalnofazowość

G z

z

B z

A z

b

z z

z

z p

nk

i

i

nb

nk

i

i

na

( )

(

)

(

)

(

)

(

)

1

1

0

1

1

Re
(s)

Re (z)

Im (z)

Im (s)

obszar
niestabilności

obszar
stabilności

obszar
stabilności

T

z

i

s

i

e

T

+j

z =1

-j

T

T





1

1

1

1

1

1

2

1

2

1

4

,

0

1

7

5

,

0

1

1

3

4

,

0

1

5

,

0

1

2

,

0

1

1

3

2

,

0

1

,

0

1

6

,

0

4

,

2

3

)

(

z

z

z

z

z

z

z

z

z

z

z

G

background image

9

Znaczenie lokalizacji bieguna na
płaszczyźnie z

-0,5

Im (z)

Re (z)

-1

1

0,5

X

X

X

X

X

X

X

X

X

X

X

X

Oscylacyj

ne

Aperiodycz

ne

Bieguny w lewej połówce koła dają oscylacyjne odpowiedzi układu

Tutaj

człon astatyczny

Tutaj

człon astatyczny

background image

10

Twierdzenie inżynierskie

Twierdzenia o wartości końcowej

k

z

G z

z

d

 1

1

1

1

lim (

)

( )

)

(

lim

0

s

G

s

k

d

s

Gdzie: d jest stopniem astatyzmu
(liczbą idealnych członów całkujących)

.

G z

y

u

B z

A z

b

b z

b z

a z

a z

i

i

nb

nb

na

na

( )

(

)

(

)

 

 

1

1

0

1

1

1

1

1

na

na

na

nb

nb

nb

a

z

a

z

b

z

b

z

b

z

A

z

B

1

1

0

...

...

)

(

)

(

background image

11

11

)

6

1

)(

3

1

)(

1

(

)

3

1

)(

2

1

(

)

(

1

1

1

1

1

2

z

z

z

z

z

z

G

Obiekty astatyczne

)

6

1

)(

3

1

)(

1

(

)

3

1

)(

2

1

(

18

15

6

1

6

1

)

(

1

1

1

1

1

3

2

1

2

1

1

z

z

z

z

z

z

z

z

z

z

z

G

3

,

0

)

(

)

1

(

lim

1

1

1

1

z

G

z

k

d

z

background image

12

29

,

0

)

(

lim

1

1

1

1

z

G

k

z

91

,

2

)

(

lim

2

1

2

1

z

G

k

z

)

(

)

1

(

lim

1

1

1

z

G

z

k

d

z

)

5

,

0

)(

1

,

0

(

)

2

,

0

)(

6

,

0

(

)

05

,

0

6

,

0

(

)

12

,

0

4

,

0

(

)

05

,

0

6

,

0

1

(

)

12

,

0

4

,

0

1

(

)

(

2

2

2

1

2

1

1

z

z

z

z

z

z

z

z

z

z

z

z

z

G

)

5

,

0

)(

1

,

0

(

)

2

,

0

)(

5

(

)

05

,

0

6

,

0

(

)

1

8

,

4

(

)

05

,

0

6

,

0

1

(

)

8

,

4

1

(

)

(

2

2

2

1

2

1

2

z

z

z

z

z

z

z

z

z

z

z

z

z

G

background image

13

10

7

1

)

5

)(

2

(

)

1

(

)

(

2

1

s

s

s

s

s

s

s

G

Nieminimalnofazowość z
semioscylacjami dla operatora s

)

5

)(

2

(

)

1

(

)

(

2

s

s

s

s

G

)

5

)(

2

(

)

3

)(

1

(

)

(

3

s

s

s

s

s

G

)

(

lim

0

s

G

s

k

d

s

Jedna semioscylacja
i ujemna odpowiedz

Dwie semioscylacje
i dodatnia odpowiedz

background image

14

Nieminimalnofazowość
z semioscylacjami dla operatora z

)

(

)

1

(

lim

1

1

1

z

G

z

k

d

z

)

5

,

0

)(

1

,

0

(

)

2

(

)

(

4

z

z

z

s

G

)

5

,

0

)(

1

,

0

(

)

5

)(

2

(

)

(

6

z

z

z

z

s

G

)

5

,

0

)(

1

,

0

(

)

2

(

)

(

5

z

z

z

s

G

)

5

,

0

)(

1

,

0

(

)

5

)(

2

(

)

(

7

z

z

z

z

s

G

Jedna semioscylacja
i ujemna odpowiedz

Dwie semioscylacje
i dodatnia odpowiedz

al
e

al
e

background image

15

15

Model obiektu

Model obiektu ze sprzężeniem zwrotnym

68

,

0

)

(

lim

1

1

1

1

z

G

k

z

18

,

2

)

(

lim

2

1

1

1

z

G

k

z

)

2

1

)(

6

1

(

)

2

,

0

1

)(

3

1

(

1

4

12

1

2

,

3

6

,

0

)

(

2

2

1

z

z

z

z

z

z

z

z

z

G

)

2

,

0

)(

85

,

0

(

)

2

,

0

1

)(

3

1

(

2

2

,

7

4

,

11

1

2

,

3

6

,

0

)

(

1

)

(

)

(

2

2

1

1

2

z

z

z

z

z

z

z

z

z

G

z

G

z

G

background image

16

Funkcja wymierna

rational function

i

bezpośrednie dzielenie

direct, long division

Klasa funkcji w postaci ilorazu dwóch

wielomianów (funkcja przejścia) A(z)/B(z)

A(z) = c

0

+ c

1

z

-1

+ c

2

z

-2

+ c

3

z

-3

+...

B(z) = 1

wtedy współczynniki są punktami wykresu

 

  

 

reszta

z

z

z

z

z

z

z

z

z

zer

o

uzupe

ł

z

z

z

z

z

z

z

z

z

z

z

B

z

A

7161

,

1

751

,

3

41

,

4

1

,

3

1

1

,

2

1

,

2

5

.

1

1

,

2

1

,

2

1

1

,

1

1

2

3

4

2

3

6

7

2

3

2

2

2

Jest to suma delt
Kroneckera

(nie Diraca)

background image

17

Najprostsze równanie

x’ = ax + bu sX = aX + bU

a

s

b

U

X

s

G

)

(

1

1

)

(

sT

k

s

G

gdzi
e

a

b

k

a

T

 ,

1

1

x t

e

x t

b
a

e

a t t

a t t

( )

( )

(

)

(

)

(

)

0

0

0

1

na skok jednostkowy
u(t) = 1(t-t

0

)

Niech t - t

0

będzie okresem próbkowania T = t

i+1

- t

i

x

i+1

=

x

i

+

u

i

lub x

i

=

x

i - 1

+

u

i - 1

gdzie:

= e

aT

i

=

a

b

(e

aT

- 1)

x(z)=

z

-1

x(z) +

z

-1

u(z)

1

1

1

)

(

)

(

)

(

z

z

z

u

z

x

z

G

background image

18

Jego rozwiązanie

na wymuszenie skokowe u

i

= 1(t)

dla

=1,

= 0,2 i

= 0,5,

= 0,2

0 T 2T 3T

...

1

t

x

i

0 T 2T 3T

...

1
1

t

x

i


0,6


0,2

0,4

x

i

x

i+1

x

i-1

x

i-2

1

1

1

)

(

)

(

)

(

z

z

z

u

z

x

z

G

x(z)=

z

-1

x(z) +

z

-1

u(z)

background image

19

Zastosowanie transformaty
Laplace’a

6

)

(

6

)

(

5

)

(

2

t

y

dt

t

dy

dt

t

y

d

)

2

(

5

)

3

(

4

1

)

2

)(

3

(

6

12

2

)

6

5

(

6

12

2

)

(

2

2

2

s

s

s

s

s

s

s

s

s

s

s

s

s

s

Y

Po zastosowaniu transformaty Laplace’a:

s

2

Y(s) –sy(0

+

) - dy(0+)/dt -5sY(s) -5 y(0

+

)+6Y(s)=6/s

Podstawiając warunki początkowe otrzymujemy:

Stosując odwrotną transformatę Laplace’a:

Warunki początkowe dy(0

+

)/dt =2, y(0

+

)=2

t

t

e

e

t

y

2

3

5

4

1

)

(

Zaistnienie warunków początkowych wprowadza zera i
często podnosi rząd mianownika i tworzy astatyzm

background image

ALGORYTMY

STERUJĄCE

ALGORYTMY

STERUJĄCE

Podział podstawowych
algorytmów sterowania

20

- Regulatory rozmyte

- Regulatory neuronowe

KLASYCZNE

KLASYCZNE

INTELIGENTNE

INTELIGENTNE

LINIOWE

LINIOWE

NIELINIOWE

NIELINIOWE

REGUŁOWE

REGUŁOWE

EKSPERTOWE

EKSPERTOWE

- PID

- LQ (liniowo-kwadratowe)

- Deadbeat

- IMC(odporne z mod.wew)

- Predykcyjne MPC

- Linearyzujace sprzeż zwrotn

- Czasoptymalne

- Ślizgowe

ADAPTACYJNE

ADAPTACYJNE

MIESZANE

MIESZANE

MODELE
STRUKTURALNE

MODELE
BEHAWIORALNE

background image

21

Klasyczny algorytm sterowania
Przetworniki A/C i C/A

Przetworniki A/C, których model w zależności od typu interpolacji
składa się z próbkowania i modelu przebiegu sygnału pomiędzy
próbkami (ZOH-prostokątów, FOH-trójkątów).

Jego model w zamkniętym układzie regulacji jest uwzględniany
w torze głównym jako dodatkowy człon poprzedzający model procesu.
Przetworniki C/A są bardzo proste i ich model jest pomijalny, gdyż
analogowe wyjście jest sumą ważoną z zestawu impulsów wejściowych
i proces ten jest prawie bezzwłoczny. Prawie wszystkie obiekty są
inercyjne (są filtrami dolnoprzepustowymi), dlatego schodkowy wyjściowy
sygnał z C/A dla dobrze dobranego czasu próbkowania jest wystarczająco
wyfiltrowany.

STEROWNIK

y

n

e

n

-

T

T

Zegar (T)

REGULATOR

G

s

(z)

y

0

y(t)

m(t) PROCES

G

p

(s)

T

Interpolator ZOH

(1-e

-sT

) / s

Przetwornik A/C

Przetwornik

C/A

Interpolator ZOH

(1-e

-sT

) / s

Przetworn

ik A/C

background image

22

Regulator PID z algorytmem pozycyjnym
lub prędkościowym (przyrostowym)

u

T

k

e

T

j

e

T

T

k

e

k

k

u

k

j

d

i

p

0

)

(

)

(

)

(

)

(

gdzie: T - czas próbkowania, czas w którym musi być wykonany algorytm

regulatora,

u – sygnał wyjściowy z regulatora (nastawienie) gdy nie ma uchybu

(odchylenia regulacji) w kolejnych chwilach czasu 0, 1, 2, ..., k

e(k)=y

0

(k)- y(k) - uchyb regulacji w k-tej chwili

e(k)/T=[e(k) - e(k-1)] /T – przybliżenie pochodnej ilorazem

różnicowym

Wady np.:

nieograniczone narastanie sygnału sterującego u(k) w wyniku

nasycenia elementu wykonawczego albo przerwa w sprzężeniu
zwrotnym

)

4

(

)

3

(

2

)

2

(

6

)

1

(

2

)

(

[

6

)

(

)

1

(

)

(

)

(

k

e

k

e

k

e

k

e

k

e

T

T

k

e

T

T

k

e

k

e

k

k

u

d

i

p

Nie występuje sumowanie (odpowiednik całkowania) którego funkcję
przejmuje obiekt całkujący lub całkowanie elementu wykonawczego
np. silnika krokowego.

ALGORYTM
POZYCYJNY oblicza
bezwzględną wartość
sygnału sterującego

ALGORYTM PRĘDKOŚCIOWY

oblicza tylko przyrost sygnału sterującego w następnej próbce

background image

23

Dwie metody projektowania
regulatora cyfrowego

Podstawienia, zastosowanie dla modelu regulatora:

o postaci równania różniczkowego metody dwupunktowej (Eulera) dla różnicy

wstecznej lub czteropunktowej różnicy centralnej do przejścia na model

w postaci równań różnicowych,

o postaci transmitancji G

s

(s) transformaty bilinearnej (Tustina) lub bezpośrednio

transformaty Z do przejścia na model G

s

(z)

i dla obu postaci nie ma potrzeby uwzględniania modelu przetwornika A/C,

Bezpośrednia, polegająca na określeniu dyskretnego modelu procesu

(obiektu) poprzedzonego interpolatorem zerowego rzędu ZOH (ang. Zero

Order Holder) według wzoru:

Gdzie:

Z jest transformatą wyznaczaną z tablic a następnie dla dobranej

struktury regulatora G

s

(z), według klasycznych metod dobieramy

nastawy regulatora dla układu ze sprzężeniem zwrotnym.

Wyrażenie w dużym nawiasie odpowiada za odpowiedz skokową
w każdej próbce. Natomiast wyrażenie (1-z

-1

) we wzorze jest

odpowiedzialne za przesunięcie odpowiedzi o okres.

s

s

G

L

Z

z

z

G

s

G

p

p

ZOH

)

(

)

1

(

)

(

)

(

1

1

background image

24

Metoda podstawienia:
Metoda dwupunktowa (Eulera) dla
różnicy wstecz (regresywnej)

Dwupunktowa różnicy wstecz

Czteropunktowa różnica centralna (czwartego rzędu)

T

x

x

dt

dx

n

n

1

2

2

1

1

2

2

2

T

x

x

x

T

dt

x

d

n

n

n

n

dt

dx

n

dt

dx

]

3

3

[

6

1

3

2

1

n

n

n

n

x

x

x

x

T

dt

dx

n

=x

n

-x

n-1

(różnicy wprost ∆

n

=x

n+1

-x

n

dla T   czasami niestabilne)



T

e

e

T

e

e

T

e

e

T

e

e

T

e

k

D

k

k

k

k

k

k

k

k

k

5

,

1

5

,

0

5

,

0

5

,

1

4

1

)

(

3

2

1

4

3

2

1

k

k

k

k

k

e

e

e

e

e

background image

25

MATLAB
[Ld,Md] = c2dm(L,M,T,'method')

G(s)

T[sek] G(z) dla ‘zoh’ G(z) dla ‘foh’

G(z) dla ‘tustin’

G(z) dla ‘matched’

0,4

1

1

02

.

0

1

3

,

56

70

z

z

1

1

02

.

0

1

14

7

,

27

z

z

1

1

33

.

0

1

14

7

,

32

z

z

1

1

02

.

0

1

2

,

11

9

,

24

z

z

0,5

1

1

006

.

0

1

56

70

z

z

1

1

006

.

0

1

2

,

11

1

,

25

z

z

1

1

43

.

0

1

10

30

z

z

1

1

006

.

0

1

8

22

z

z

b

s

a

s

K

0

10

2

70

s

s

1

1

5

1

5

,

4

1

56

70

z

e

z

1

5

1

5

,

4

1

6

,

5

6

,

19

z

e

z

1

66

.

0

1

3

,

23

z

1

5

1

5

,

4

1

2

,

2

2

,

16

z

e

z

Metoda dyskretyzacji wpływa na postać i współczynniki transmitancji

background image

26

Metoda podstawienia

Dla transmitancji regulatora D(s) określmy równanie różnicowe

b

s

a

s

K

s

E

s

U

s

D

0

)

(

)

(

)

(

)

(

)

(

)

(

)

(

0

s

E

a

s

K

s

U

b

s

ae

dt

de

K

bu

dt

du

0





n

n

n

n

n

n

ae

T

e

e

K

bu

T

u

u

1

0

1





1

1

0

1

1

n

n

n

n

n

n

ae

T

e

e

K

bu

T

u

u

n

n

n

n

u

Tb

e

Ta

K

e

K

u

)

1

(

)

1

(

0

1

0

1

1

1

0

0

)

1

(

)

1

(

n

n

n

n

u

Tb

e

Ta

K

e

K

u

Dla parametrów a=2, b=10 i K

0

=70,

i okresu próbkowani T=0,4 sek

n

n

n

n

u

e

e

u

3

14

70

1

1

n

n

n

u

e

u

4

70

1

1

n

n

n

n

n

n

n

u

e

e

u

e

e

u

9

)

(

70

9

70

70

1

1

1

T=0,5 sek

T=1 sek

Współczynniki w opisie dyskretnym zawierają czas próbkowania

background image

27

Metoda bezpośrednia

Obiekt opisanego transmitancją i regulatora D(z)=k

a

s

a

s

G

)

(

)

(

)

1

(

)

(

)

1

(

)

(

1

1

a

s

s

a

Z

z

s

s

G

Z

z

z

G





z

z

e

z

z

e

z

aT

aT

1

1

)(

1

(

)

1

(

)

1

(

1

1

1

1

gdzie podstawiono

=e

-aT

)

1

(

)

1

(

)

(

)

(

1

)

(

)

(

)

(

)

(

k

z

k

z

G

z

D

z

G

z

D

z

X

z

Y

background image

28

Przekształcenie transmitancji
dyskretnej na równanie różnicowe

bezpośrednia,

kaskadowa,

równoległa

M z

E z

a z

b z

j

j

j

n

j

j

j

n

( )

( )

0

1

1

m

a e

b m

i

j i j

j

n

j

i j

j

n

1

1

2

1

2

1

2

,

0

1

,

0

1

6

,

0

4

,

2

3

)

(

)

(

z

z

z

z

z

E

z

M

m

i

= 3 e

i

+ 2,4 e

i -1

- 0,6 e

i -2

- 0,1 m

i -1

+ 0,2 m

i - 2

przy tej samej liczbie poziomów kwantowania
(długości słowa) daje mniejsze błędy numeryczne

ale m

i

= 3 e

i

+ 2,4 e

2

i -1

- 0,6 e

i -2

- 0,1 kT m

i -1

+ 0,2 m

i – 2

nieliniowe i niestacjonarne

background image

29

Przekształcenie transmitancji
dyskretnej na równanie różnicowe

bezpośrednia,

kaskadowa,

równoległa





1

1

1

1

2

1

2

1

4

,

0

1

5

,

0

1

2

,

0

1

1

3

2

,

0

1

,

0

1

6

,

0

4

,

2

3

)

(

)

(

z

z

z

z

z

z

z

z

z

E

z

M

W

1

= 3

W

2

= 1 + z

-1

W

3

= 1 - 0,2 z

-1

W

z

4

1

1

1 05

,

1

5

4

,

0

1

1

z

W

x

i

(1)

= 3 e

i

x

i

(2)

= x

i

(1)

+ x

i -1

(2)

x

i

(3)

= x

i

(2)

- 0,2 x

i -1

(3)

x

i

(4)

= x

i

(3)

- 0,5 x

i

-1

(4)

m

i

= x

i

(5)

= x

i

(4)

+ 0,4 x

i

-1

(5)

background image

30

Przekształcenie transmitancji
dyskretnej na równanie różnicowe

bezpośrednia,

kaskadowa,

równoległa

1

1

2

1

2

1

4

,

0

1

7

5

,

0

1

1

3

2

,

0

1

,

0

1

6

,

0

4

,

2

3

)

(

)

(

z

z

z

z

z

z

z

E

z

M

W

1

= - 3

W

z

2

1

1

1 05



,

W

z

3

1

7

1 04

,

x

i

(1)

= - 3 e

i

x

i

(2)

= - e

i

- 0,5 x

i

-1

(2)

x

i

(3)

= 7 e

i

+ 0,4 x

i

-1

(3)

m

i

= x

i

(1)

+ x

i

(2)

+ x

i

(3)

background image

31

Struktury modeli dyskretnych

FILTR

ZAKŁÓCENIA

H(z)

OBIEKT
G(z)

y

i

e

i

u

i

zakłócenia
niemierzaln
e

i

nk

i

i

nk

i

i

e

z

D

z

C

u

z

F

z

B

e

z

H

u

z

G

y

z

A

)

(

)

(

)

(

)

(

)

(

)

(

)

(

1

1

1

1

1

1

1

Boxa - Jenkinsa (BJ)

ARMAX (ang. AutoRegresive Moving-Average eXtended)

A(z

-1

) y

i

= B(z

-1

) u

i -nk

+ C(z

-1

) e

i

A(z

-1

) y

i

= B(z

-1

) u

i -nk

+ e

i

ARX dla szumu szerokopasmowego czyli szumu zbliżonego do białego o wariancji

2

Jeżeli sygnał wejściowy u

i

= 0 lub e

i

= 0

wtedy mamy do czynienia z ciągiem czasowym
Najpopularniejszy model to ARMA

MA

część

AR

część

nc

i

nc

i

i

nd

nd

i

i

i

e

c

e

c

e

y

d

y

d

y

d

y

1

1

1

2

2

1

1

y

C z

z

D z

e

i

d

i

(

)

(

)

(

)

1

1

1

1

Model ARIMA d = 1 to ciąg
niestacjonarny o zmieniającej
się wartości średniej

background image

32

Praktyczne zasady (1)

T

g

= 1/f

g

t

f(t)

Czas próbkowania powinien być co najmniej
T

g

/2, w praktyce T

g

/8 (osiem próbek)

.

.

• liczba poziomów kwantowania

• czas próbkowania T

Częstotliwości granicznej sygnału f

g

, którą można oszacować

z amplitudowej charakterystyki częstotliwościowej
lub odpowiedzi skokowej dla obiektu pierwszego rzędu

:

g

(3 dB)

z

=

p

T

1

g

(3 dB) - pulsacja graniczna 3 dB,

z

- pulsacja załamania,

T

u

- czas ustalania się (regulacji),

T

p

- stała czasowa obiektu

(procesu).

T

u

3T

p

u

p

g

T

T

f

2

1

2

1

Rozróżnialne poziomy powinny być
oddalone od siebie o wartość równą,
co najmniej amplitudzie szumów.





SZ

S

s

P

P

f

1

log

2

max

=

C





s

bit

Hartleya i Shannona

background image

33

Praktyczne zasady (2)

Podsumowując można stwierdzić, że czas próbkowania T jest

wielkością krytyczną. Typowo dobiera się jako 1/10 dominującej

stałej czasowej obiektu, jeżeli obiekt ma zauważalny czas

opóźnienia L to T≈0,33 L. Jeżeli sterowanie jest w dalszym ciągu

nie zadowalające (np. nieminimalnofazowy proces) to należy

zmniejszyć czas próbkowania a jeżeli to nie pomoże to należy

zmienić prawo sterownia (np. sterowanie predykcyjne)

Jeżeli projektowanie realizowane jest przy użyciu metody

podstawienia (bezpośredniego projektowania regulatora) to

częstotliwość próbkowania powinna być 20 razy większa od

pasma przenoszenia obiektu natomiast, jeżeli stosujemy metodę

bezpośrednią (bezpośredniego projektowania cyfrowego) to

częstotliwość próbkowania powinna być mniejsza i być około

2 razy większa od pasma przenoszenia. W praktyce jednak,

aby móc pozbyć się losowych zakłóceń stosujemy częstotliwość

próbkowania co najmniej 20 razy większa od pasma przenoszenia

dla układu zamkniętego

background image

34

Dziękuję za uwagę

i zapraszam do wykorzystywania

tej wiedzy w laboratorium i praktyce


Document Outline


Wyszukiwarka

Podobne podstrony:
WYKŁAD 02 SterowCyfrowe
WYKŁAD 02 SterowCyfrowe
WYKŁAD 02 SterowCyfrowe
logika wyklad 02
pdf wykład 02 budowa materii, podstawowe prawa chemiczne 2014
013 HISTORIA SZTUKI WCZESNOCHRZEŚCIJAŃSKIEJ I BIZANTYJSKIEJ, WYKŁAD,# 02 10
Młoda Polska WYKŁAD (02 04 2014)
wykład 02 2013
pmp wykład 02 15
wykład& 02 2013
Wykład 2# 02 12
2006C16 wyklad 02
EIE wykład 3 - 02.04.2011 r, Ekonomia integracji europejskiej
Drogi i ulice wyklad 02
wykład 3" 02 12

więcej podobnych podstron