Nyquist Frequency
Nyquist Frequency
Nyquist Frequency
Milan R. Rapaic
5. maj 2014 13. maj 2014
W ( s)
.
1 + W ( s)
(1)
(2)
Tvrdenje
1. (Koijeva teorema argumenta.) Neka su dati zatvo
rena kriva C u kompleksnoj ravni i kompleksno preslikavanjeF :
C C. Neka je F na samoj krivoj C analiticko, a neka unutar
krive poseduje N nula i P polova. Tada, ukoliko kompleksna
c itavu krivu C, kompleksna slika F (s) c e
promenljiva s obide
obici koordinatni pocetak | N P| puta. Pri tome, ukoliko je
N > P smerovi obilaska su jednaki; u suprotnom oni su razliciti.
W (s)
Podsecanje!
Proces c ija je funkcija povratnog prenosa W ( s ) nakon
zatvaranja povratne sprege
opisuje se funkcijom spregnutog prenosa
Gs ( s ) =
W (s)
.
1 + W (s)
se reavanjem
svojstvene jednacine
F(s) = 1 + W (s) = 0 .
Ukoliko su sva reenja svojstvene
jednacine u levoj poluravni (odnosno ukoliko su im realni delovi
negativni), tada i samo tada je
posmatrani proces u zatvorenoj
sprezi asimptotski stabilan.
Im s
D
C
Re s
Re s
A
(a) Kontura C.
Im s
Im s
C
Re s
Re s
B
A
B
(a) Kontura C.
Im s
Im s
D
Re s
C
Re s
B
B
A
(a) Kontura C.
D=B
Im s
Re s
Re s
A=C
(a) Kontura C.
Im s
Im s
Re s
(a) Kontura C.
Re s
Izvodenje
Nikvistovog kriterijuma Principi
Nikvistov kriterijum se izvodi neposrednom primenom Koijeve
teoreme argumenta na svojstvenu funkciju
Im s
R
Re s
F ( s) = 1 + W ( s)
i na konturu koja obuhvata c itavu desnu poluravan kompleksne
ravni. Takozvana Nikvistova kontura CN je prikazana dijagramom
na slici 7. Ova kontura se protee c itavom imaginarnom osom, a
zatvara se preko polukrunog isecka beskonacnog poluprecnika.
Uobicajeno je Nikvistovu konturu obuhvatati u negativnom matematickom smeru (tj. u smeru kretanja kazaljke na satu). Na taj nacin
kretanje je usaglaeno sa orijentacijom same imaginarne ose.
Jasno, ukoliko je proces nakon zatvaranja povratne sprege apsolutno stabilan, svojstvena jednacina nema reenja (nula) u desnoj
Primetimo ...
Svojstvena funkcija
F(s) = 1 + W (s)
jednovremeno sadri informacije o polovima procesa
i pre i posle zatvaranja povratne sprege. Polovi procesa
u otvorenoj sprezi su jednovremeno i polovi F ( s ), polovi
procesa nakon zatvaranja povratne sprege su nule F ( s ).
Na osnovu Koijeve teoreme,
preslikavanjem Nikvistove
konture CN pomocu svojstvene
funkcije F ( s ) vri se prebrojavanje razlike broja nestabilnih
polova pre i posle zatvaranja
povratne sprege. Svaki obrt
preslikane Nikvistove konture
oko koordinatnog pocetka u
pozitivnom matematickom
smeru oznacava uklanjanje
jednog nestabilnog pola procesa. Ovakvih, nestabilnih
polova ima ukupno P. Otuda
je neophodno izbrojati P obrta
preslikane Nikvistove konture.
Citav
kruni isecak beskonacnog poluprecnika preslikava se u
jednu tacku. Da bismo se u to uverili, primetimo da mora postojati
granicna vrednost
lim W (s) R .
s
y(0) = lim sW ( s )
s
1
= lim W ( s ) .
s
s
Z
0
w ( t) e st dt .
w ( t) e st dt
w ( t) e st dt =
w ( t) e st dt = W ( s ) .
Z
0
w ( t) e st dt
Im W (s)
Im s
R
Re s
(a) Kontura CN .
Re W (s)
Tvrdenje
2. (Nikvistov kriterijum stabilnosti) Proces opisan funk
cijom povratnog prenosa W (s) koja ima P polova u desnoj poluravni kompleksne ravni i nema polova na imaginarnoj osi bice
stabilan nakon zatvaranja jedinicne povratne sprege ako i samo
kriticnu tacku (1, 0) P/2
ako Nikvistov vektor W ( j ) obide
puta, odnosno za ugao P, kako se menja od nule do beskonacnosti.
Primer 1. (Ilustracija primene Nikvistovog kriterijuma) Primenom
Nikvistovog kriterijuma ispitati stabilnost procesa opisanog
funkcijom povratnog prenosa W (s) = ( s+k1)3 .
Nakon smene s sa j, te razdvajanja realnog i imaginarnog dela
dobijenog izraza, funkcija povratngo prenosa se moe zapisati u
obliku
2 3
1 3 2
W ( j ) = k
+
jk
.
(3)
( 2 + 1 )3
( 2 + 1 )3
Prilikom crtanja Nikvistove krive, bitno je odrediti tacku u kojoj
kriva pocinje, tacku u kojoj se kriva zavrava, kao i tacke u kojima
kriva preseca koordinatne ose. Pri tome, narocito su bitni preseci
ta vie, moguce veoma lako odrediti ugao pod kojim kriva uvire
u koordinatni pocetak. Za vrlo veliku vrednost argumenta, svaki
Im W ( j )
k 18
k
=
Re W ( j )
3
=0
k 3 8 3
=
3
3
1
.
( s + 1) 3
k
3
.
= 3 arg( j ) =
2
( j )3
k 3 8 3 .
3
3
takvo da je Im W ( j ) =
1 <
k
8
k<8.
tacaka preseka sa
realnom i imaginarnom osom. Preseke
sa realnom osom odredujemo
kao
tacke u kojima je argument funkcije
povratnog prenosa jednak l, za
neko celobrojno l. Slicno, preseke sa
imaginarnom osom odredujemo
kao
tacke u kojima je argument funkcije
povratnog prenosa 2 + l. Kako je
arg W ( j ) = 3 arg( j + 1)
= 3 arctan ,
l
= tan
+
,
6
3
= tan
pr. sa Re -osom
pr. sa Im -osom
Nenegativna reenja
prvog uslova su
= 0 i = 3. Jedino konacno,
pozitivno
reenje drugog uslova je
3
3 .
Drugim recima, posmatrani je sistem stabilan nakon zatvaranja povratne sprege za svako pozitivno pojacanje k < 8. Nikvistove krive
dobijene za razlicito k prikazane su uporedo na slici 10. Valja primetiti da se sa povecanjem pojacanja Nikvistova kriva radijalno iri. Ova
pojava je ilustrovana dijagramom na slici 10.
Kada se Nikvistov kriterijum primenjuje na funkciju povratnog
prenosa koja poseduje podesivo pojacanje k, kao to je to bio slucaj
u prethodnom primeru, c esto je prikladnije crtati Nikvistovu krivu
samo na osnovu normalizovane funkcije povratnog prenosa dobijene
za k = 1. Analiza stabilnosti za proizvoljno k tada se moe vriti
pomeranjem kriticne tacke u poloaj 1k . Zaista, kako je
1 + kW (s) = 0
1
+ W ( s) = 0 ,
k
analiza stabilnosti se moe ravnopravno vriti posmatranjem broja obrta podesive krive kW ( j ) oko nepokretne kriticne tacke (1, 0) ili posmatranjem
obrta normalizovane, nepokretne krive W ( j ) oko podesive kriticne tacke
( 1k , 0). Slika 11 ilustruje ovaj nacin analize stabilnosti.
Vana napomena
Potpuno je svejedno da li crtamo Nikvistovu krivu funkcije
povratnog prenosa kW ( s ), pa
potom analiziramo broj obrta
oko kriticne tacke (1, 0) ili
crtamo Nikvistovu krivu na
osnovu normalizovane funkcije
povratnog prenosa W ( s ), a potom analiziramo broj obrta oko
podesive kriticne tacke ( 1k , 0).
Sa prakticnog stanovita drugi
nacin je najcece pogodniji.
Im W ( j )
Re W ( j )
k = 0.5
k=2
k=4
k=8
k = 10
Im W ( j )
81
1k
Re W ( j )
Podsecamo c itaoca da je
e j5 = cos(5 ) j sin (5 ) .
Im W ( j )
0.8
Re W ( j )
1k
0.4
1
< 0.8
k
k < 1.25 .
Poredenjem
sa rezultatima prethodnog primera, vidimo da je uvo
denje
kanjenja u znacajnoj meri umanjilo opseg stabilnosti procesa
nakon zatvaranja povratne sprege.
1 3 2
2
( 2 + 1)
(3 3 )
( 2 + 1)
,
.
Nikvistova kriva pocinje u tacki W ( j0) = 1, a zavrava se u koordinatnom pocetku kada tei nuli. Preseke sa koordinatnim osama
nalazimo na uobicajeni nacin.
Im W ( j 3) = 23
3
3 ,
Im W ( j )
3
2
Re W ( j )
stabilisati uvodenjem
odgovarajuce povratne sprege. Medutim,
sem
pitanja Da li je sistem stabilan?, postavlja se i pitanje Koliko je sistem
stabilan?. Drugim recima, u kolikoj meri i na koji nacin moemo
promeniti ponaanje razlicitih komponenata a da sistem u celini i
dalje zadri svojstvo stabilnosti. Ova i slicna pitanja spadaju u oblast
Pretek pojacanja
Pretek pojacanja
Definicija 1. (Pretek pojacanja.) Pretek pojacanja procesa opisanog funkcijom povratnog prenosa W (s) jeste minimalno pojacanje d > 0 takvo da je proces kW (s) stabilan u zatvorenoj sprezi
za svako k < d.
1
.
|W ( j )|
(5)
10
11
Definicija 2. (Pretek kanjenja.) Pretek kanjenja procesa opisanog funkcijom povratnog prenosa W (s) jeste minimalno transportno kanjenje c > 0 takvo da je proces W (s)esc stabilan u
zatvorenoj sprezi za svako < c .
Re W ( j )
0.766
Re W ( j )
1
( s +1 )3
(b) W ( s ) =
1
e 2.56s
( s +1 )3
Vidimo da nakon obrtanja za priblino 67.6 (odnosno 1.18 radijana) u smeru kazaljke na satu Nikvistov dijagram procesa opisanog
funkcijom prenosa W (s) = ( s+11)3 prolazi tacno kroz kriticnu tacku,
odnosno postaje granicno stabilan. Za svaki manji ugao obrta, proces
je i dalje stabilan. Za svaki veci ugao obrta, proces je nestabilan.
Definicija 3. (Pretek faze.) Pretek faze procesa opisanog funkcijom povratnog prenosa W (s) jeste minimalan ugao za koji je
neophhodno obrnuti Nikvistovu krivu da bi ona prola kroz
kriticnu tacku.
Postupak odredivanja
preteka faze je sledeci. Najpre se odredi
vrednost za koju je moduo funkcije povratnog prenosa jednak 1.
Takvo se obeleava sa 1 i ne mora biti jednoznacno u optem
slucaju,
|W ( j1 )| = 1 .
(6)
1
( s + 1) 3
0.766
(a) W ( s ) =
67.6
Im W ( j )
1
e 2.56s .
( s + 1) 3
(7)
p f
.
1
(8)
12