Sayısal Kontrol Sistemleri Bölüm 1 Ayrık Zaman sinyaller ve Sistemler 1.1 Kapalı çevrim sayısal kontrol sistemi 1.2 Sürekli zaman sinyaller ve sistemler : Hatırlatma 1.3 Ayrık zaman sinyaller ve sistemler 1.4 Fark denklemleri, çözümleri ve uygulama örnekleri 1.5 Sürekli zaman sistemlerin ayrıklaştırılması
165
Embed
Sayısal Kontrol Sistemleri Bölüm 1 Ayrık Zaman sinyaller ...
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Sayısal Kontrol Sistemleri
Bölüm 1
Ayrık Zaman sinyaller ve Sistemler
1.1 Kapalı çevrim sayısal kontrol sistemi
1.2 Sürekli zaman sinyaller ve sistemler : Hatırlatma
1.3 Ayrık zaman sinyaller ve sistemler
1.4 Fark denklemleri, çözümleri ve uygulama örnekleri
1.5 Sürekli zaman sistemlerin ayrıklaştırılması
Referans
giriş
r(k)
Sistem
çıkışı
y(t)
Algılayıcı
+
-
YükselticiSayısal
Denetleyici
Hata
e(k)
Denetim
sinyali
v(t)
Sistem
girişi
u(t)
Sistem
ADC
DAC
Sayısal kontrol sistemlerinde genellikle denetlenen
sistem analog (sürekli zaman) bir sistem iken
denetleyici sayısal (ayrık zaman) bir birimdir.
Tutma
Örnekleme
1.1 Kapalı Çevrim Sayısal Kontrol Sistemi
Sistem ? Algılayıcı ? Yükseltici/Dönüştürücü ?
1.2 Sürekli Zaman Sinyaller ve Sistemler: Hatırlatma
y(t)Analog
Sistem
u(t)tetu 2)( tety 210)(
Zamanda sürekliliği olan analog sinyaller ve analog sistemler
y(t)Analog
Sistem
u(t)
t(s)
y(t)
0
-10?
0
1
t(s)
u(t)
Sürekli zaman sistemler:
• Statik / Dinamik sistem ?
• Doğrusal /doğrusal olmayan sistem ?
• Zamanla değişen / zamanla değişmeyen sistem ?
t
x(t)
0
t
x(t)
0
A
1
t
x(t)
0
A
y(t)Analog
Sistem
u(t)
Temel Test Sinyalleri
)t(x4)t(y5dt
)t(dy2
dt
)t(yd2
2
t
0
dt)t(x10)t(y
1)t(x2)t(y
)t(x)wt(Sin)t(y
1.2 Sürekli Zaman Sinyaller ve Sistemler: Hatırlatma
G(s)U(s) Y(s)
• Doğrusal ve zamanla değişmeyen sistemlerin
matematiksel modelleri
)()()(
)()()(
tDutCxty
tButAxtxdt
d
)t(u5dt
)t(du4)t(y3
dt
)t(dy2
dt
)t(yd2
2
3s2s
5s4
)s(U
)s(Y)s(G
2
y(t)Analog
Sistem
u(t) Diferansiyel denklem
Durum denklemiTransfer fonksiyonu
1.2 Sürekli Zaman Sinyaller ve Sistemler: Hatırlatma
G(s)U(s) Y(s)
Örnek: Verilen sistemin birim rampa cevabını bulunuz.
2
1
)(
)()(
ssU
sYsG
G(s)U(s) Y(s)
Örnek: Verilen sistemin birim basamak cevabını bulunuz.
12
1
)(
)()(
2
sssU
sYsG
R(s)Gc(s)
Y(s)+
-Gp(s)
H(s)
U(s)E(s)
Örnek: Verilen kontrol sistemin birim basamak cevabını bulunuz.
1)(
2
1)(
)1(10)(
sH
ssG
s
ssG
p
c
Ayrık birim basamak fonksiyonu kku )1()(
Ayrık birim impuls fonksiyonu d(k)
k)k(r Ayrık birim rampa fonksiyonu r(k)
Ayrık zaman Temel Test Sinyalleri
1.3 Ayrık Zaman Sinyaller ve Sistemler
Analog Sinyallerin ayrıklaştırılması:
Örnekleme peryodu ve seçimi ?
0 0.5 1 1.5 20
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
300
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
20100
t = k.T
tetf 2)(
ka)k(f
kekf 2.0)(
1.3 Ayrık Zaman Sinyaller ve Sistemler
Sönümlü sinüsoidal sinyalin ayrıklaştırılmışı
)5(10)( 2 tSinetf t
)25.0(10)( 1.0 kSinekf k
0 0.5 1 1.5 2-2
-1
0
1
2
3
4
5
6
0 20 40 60-2
-1
0
1
2
3
4
5
6
)5(10)( 2 tSinetf t
Zaman (t) Örnek sayısı (k)
T=?
ttetf 101)(
1-) Verilen sinyali ayrıklaştırarak grafiğini çizen bir
MATLAB programı yazınız. T=?
Örnekler
)10()( 5 tCostetf t
2-) Verilen sinyali ayrıklaştırarak grafiğini çizen bir
MATLAB programı yazınız. T=?
y(k)Sayısal
Sistem
u(k)
Sayısal (Ayrık zaman) Sistemler
Ayrık zaman sistem tanımları
• Statik / Dinamik sistem ya da bellekli /belleksiz sistem
• Doğrusal / doğrusal olmayan sistem
• Nedensel / nedensel olmayan
• Zamanla değişen / zamanla değişmeyen sistem
Ayrık zaman sistem, ayrık zaman giriş sinyalini alarak
ayrık zamanda tanımlı dinamik ya da statik bir giriş-
çıkış ilişkisi verir.
y(k)=F{ y(k),u(k)) }
)()( 2.0 kueky k?
Örnek
y(k-1)
y(k)
y(k+1)
Sayısal (Ayrık zaman) Sistemler
Örnekler: Verilen sistemlerin bellekli, nedensel, doğrusal ve
zamanla değişmeyen olup olmadıklarını belirleyiniz.
k k+1k-1
k=0,1,2…..N
)1(8.0)()1(8.0)( kukukyky
)1(5.0)()( kukuky
)()10()( 2 kukCosky
2)()( kuky
Ayrık zaman Sistemlerin Cevabını Bulma
)k(ue)k(y k2.0
Verilen sistemlerin,
a-) Öteleme elemanını dikkate alarak
zaman bölgesi blok şemasını çiziniz.
b-) birim basamak, birim impuls ve birim rampa cevaplarını
k=0,1,2 örnek için iteratif çözerek bulunuz.
)1k(u5.0)k(u)k(y
)k(u)1k(y8.0)k(y
Öteleme (geciktirme) elemanı
y(k)Sayısal
Sistem
u(k)
Kararlılık: Sınırlı girişe karşı sınırlı çıkış veren sistemler
SGSÇ kararlıdır.
)1(5.0)()( kukuky
)()10()( 2 kukCosky
bkyiseaku )()(
Verilen sistemlerin kararlı olup olmadıklarını gösteriniz.
)k(ue)k(y k2.0 y(k)Sayısal
Sistem
u(k)
1.4 Doğrusal, zamanla değişmeyen ve dinamik ayrık zaman
Sistemler
Fark denklemleri
• Doğrusal ve zamanla değişmeyen sürekli zaman sistemlerdedurumlar, çıkışlar ve giriş fonksiyonları, t zamanın süreklifonksiyonları olduğu halde,
• Ayrık zaman sistemlerde bu büyüklükler, ayrık zaman kdeğişkeninin ayrık anlarında tanımlıdırlar. Burada k = 0, 1, 2,...N
• Örneklenmiş ayrık zaman sistemlerde bu büyüklükler, sürekli tzamanının ancak kT gibi ayrık anlarında tanımlıdırlar. Burada k =0, 1, 2,...
• Sürekli zaman sistemlerin matematiksel modelleri türevleri ihtivaeden diferansiyel denklemlerle tanımlanırken,
• Ayrık zaman sistemlerin modelleri ise giriş ve çıkışlarınötelenmişlerine göre yazılan fark denklemleri ile tanımlanır
• Fark denklemleri geriye farklar ya da ileriye farklar şeklindeyazılabilir.
n. dereceden bir fark denkleminin geriye farklar
biçiminde yazılışı genel olarak :
)mk(ub)1k(ub...)k(u)nk(ya)1k(ya...)k(y m1n1
Burada, k=0 için y(-1), y(-2),……y(-n) başlangıç koşullarıdır.
Örnek: Blok şemasını çiziniz.
)1(3)()2()1(2)( kukukykyky
)k(ub)1k(ub...)mk(u)k(ya)1k(ya...)nk(y 0101
Burada, k=0 için y(n-1), y(n-2),……,y(1), y(0) başlangıç
koşullarıdır.
n. dereceden bir fark denkleminin ileriye farklarbiçiminde yazılışı genel olarak :
)(3)1()()1(2)3( kukukykyky
Örnek: Blok şemasını çiziniz.
Örnekler:
2)0(y:kosuluBaslangic
)1k(u)k(y5.0)1k(y
2)2(y,4)1(y:kosullariBaslangic
)1k(u3)k(u)2k(y)1k(y2)k(y
Fark denklemleri, başlangıç koşulları hesaplanmak koşuluyla
ileri ya da geri yönde ötelenebilir.
k=k+2 ile öteleyiniz.
k=k-1 ile öteleyiniz.
Fark denklemlerinde öteleme
Sabit katsayılı diferansiyel denklemlerin çözümünde olduğu gibi farkdenklemleri de homojen çözüm ve özel çözümün toplamındanmeydana gelir.
n21 r...rr
0arardaya0ra...ra1 n1n
1nn
n1
1
Homojen çözüm yh(k):
Homojen çözüm: Karakteristik denklem
KD’ in kökleri reel ve ayrık ise
KD’ in kökler reel ve katlı ise
KD’in kökleri kompleks ise
knn
k22
k11 rA...rArA)k(yh
jbar,jbar 21
n21 r,...r,r
k1
1nn
k12
k11 rkA...krArA)k(yh
)/(tan,),()()( 122
21 abbakSinAkCosAkyh kk
)k(u)nk(ya...)1k(ya)k(y n1
Fark Denklemlerinin Analitik Çözümü
Fark denklemlerinin özel çözümü
Diferansiyel denklemlerde olduğu gibi kaynak
fonksiyonunun türüne bağlı olan, bir de özel çözüm vardır.
)k(Cosdaya)k(Sin
k
)a(
1.K
n
k
k
özel çözümKaynak
)k(CosC)k(SinC)k(yö
k.C)k(yö
)a.(C)k(yö
1.C)k(yö
21
n
k
k
Özel çözümün aynısı, homojen çözümde de varsa her
defasında özel çözüm k ile çarpılmalıdır.
olsunyykuBurada
kukykyky
k 0)2(,1)1(,)1()(
)()2(1.0)1(3.0)(
Örnek: Verilen fark denklemini çözünüz.
olsunyykukykyky 2.0)1(,0)0(,)()(0625.0)1(5.0)2(
Örnek: Verilen fark denklemini çözünüz.
Fark Denklemlerinin Ardışıl ya da Bilgisayarla Çözümü
Fark denkleminin herhangi bir k=0 için bulunacak çıkışı y(0) olduğuna göre
önceki çıkışları yani y(-1),y(-2),…y(k-n) gibi başlangıç koşuludur. Buna
göre basit bir döngü mantığı ile fark denklemleri ardışıl-iteratif çözülebilir.
olsunyykukykyky 0)2(,1)1(,)()2(1.0)1(3.0)(
Örnek: Verilen fark denkleminin birim basamak cevabını k=0,1,2 örnek
için ardışıl çözünüz. MATLAB programını yazınız.
olsunyy
kukukykyky
4.0)1(,0)0(
,)1(8.0)()1(2.0)(5.0)1(
Örnek: Verilen fark denkleminin birim rampa cevabını k=0,1,2 örnek
için ardışıl çözünüz. MATLAB programını yazınız.
• Sayısal Türev: Bir sinyalin herhangi bir noktadaki türevisinyalin o noktadaki eğimidir. Buna göre Bir f(t)sinyalinin türevi y(t) ile gösterilirse bu sinyalin herhangi bir t(yada k, kT) noktasındaki eğimi, sinyalin zaman ekseni ileyaptığı açı kullanılarak aşağıdaki gibi yazılabilir.
• Fonksiyonun eğimi farklı yaklaşımlarla bulunabilir.
• Bunlardan ileri fark, geri fark ve merkezi fark olarakbilinen yöntemler yaygın olarak kullanılmaktadır.
tan)()( tfdt
dty
k
f(k)
t
f(t)
1.5 Sürekli zaman sistemlerin ayrıklaştırılması
İleri Fark Yöntemi
• İleri fark yönteminde fonksiyonun bir gelecekteki örneklenmiş değeri kullanılarak fonksiyonun o andaki eğimi yaklaşık olarak hesaplanabilir.
•
• İleri fark yöntemine göre sayısal türev,
T
kfkfky
)()1()(
k
f(k)
t
f(t)
k+1
f(k+1)
• Geri fark yönteminde fonksiyonun bir önceki örneklenmiş değeri kullanılarak fonksiyonun o andaki eğimi yaklaşık olarak hesaplanabilir.
•
• Geri fark yöntemine göre sayısal türev,
T
kfkfky
)1()()(
k
f(k)
t
f(t)
k-1
f(k-1)
Geri Fark Yöntemi
Merkezi Fark Yöntemi
(Bilineer Dönüşüm-Tustin Algoritması)
Merkezi fark yönteminde ise fonksiyonun bir önceki ve bir
gelecekteki örneklenmiş değerleri kullanılarak fonksiyonun o
andaki eğimi yaklaşık olarak hesaplanabilir.
Merkezi fark yöntemi ile sayısal türev,
T
kfkfky
2
)1()1()(
k
f(k)
t
f(t)
k-1
f(k-1)
k+1
f(k+1)
Sayısal İntegral
• Bir sinyalin herhangi bir andaki integrali (belirli-sınırlı integrali) bu sinyalin o ana kadar taradığı alanın toplamıdır. Bir f(t) fonksiyonun integrali y(t) ile gösterilirse t anındaki integral,
• Bir fonksiyonun herhangi bir t anındaki taradığı alanı bulmak için çeşitli yöntemler geliştirilmiştir. Burada en çok kullanılan sol kenar, sağ kenar ve yamuk yöntemleri verilmiştir.
t
0t
)0t(ydt)t(f)t(y
(k-1) k
f(k)
f(k-1)
T
t
f(t)
y(k-1) = y(t0) Alan
)1(. kfTAlan
)k(f.TAlan
)1()(.2
kfkfT
Alan
Sol Kenar Kuralı:
Sağ Kenar Kuralı:
Yamuk Kuralı:
Alan)1k(y)k(y
Örnekler
1-) Bir kondansatöre birim rampa gerilim uygulanmıştır.
Kondansatör akımını bulan bir MATLAB programı yazınız.
C=0.1 F. Çeşitli ayrıklaştırma yöntemleri ile tekrarlayınız.
2-) İL(0)=2 A ile yüklü bir bobine v(t)=10Sin(5t) gerilimi
uygulanmıştır. Bobinin akımını bulan bir MATLAB programı
yazınız. L=0.2 H. Çeşitli ayrıklaştırma yöntemleri ile
tekrarlayınız.
Diferansiyel denklemler ayrıklaştırılarak Fark denklemleri
elde edilebilir.
)t(u3dt
)t(du4)t(y2
dt
)t(dy
)1k(uT21
4)k(u
T21
T34)1k(y
T21
1)k(y
Örnek:
Geri fark yöntemi kullanılarak ayrıklaştırılırsa,
Eğer T=0.1 saniye seçilirse aşağıdaki fark
denklemi elde edilir.
)1k(u333.3)k(u583.3)1k(y833.0)k(y
Diferansiyel denklemlerin ayrıklaştırılması
Dolayısıyla dinamik sistemlerin analog modelleri bir integro-
diferansiyel denklem iken dinamik sistemlerin sayısal modelleri
giriş ve çıkışın gecikmişleri ile ifade edilen bir fark denklemi
Örnek: Ayrık zaman üstel fonksiyonun z dönüşümünü alınız.
Çözüm: z-dönüşüm serisinin açılımından,
az
z
z.a1
1)z(Folursa)a()k(f
ez
z
z.e1
1)z(F
1
k
b1b
1|ze| 1bT
Genel olarak üstel fonksiyonlar
k)a()k(f
Örnek: Birim rampa fonksiyonunun z dönüşümü
Çözüm
0k
321k ...z3z2z0z).k(f)z(F
2)1z(
z)z(F
k)k(f
0 1 2 3 4….
f(k)
k
Önceki seriden farklıdır ve yakınsaklığı ve dolayısıyla
z-dönüşümü aşağıdaki gibi bulunur.
Fonksiyon
f(k)Z Dönüşümü
F(z)
)wk(Cos
)wk(Sin
)a(
k
k
)k(u
)k(
k
2
d
1zCoswa2z
)wCosz(z
1wzCos2z
wzSin
az
z
)1z(
)1z(z
)1z(
z
1z
z
1
2
2
3
2
Bazı Ayrık Fonksiyonların z- Dönüşümleri
Analog Fonksiyon
f(t)
z- Dönüşümü
F(z)
)wkT(Cos
)wkT(Sin
e
kT
bkT
1wTzCos2z
)wTCosz(z
1wTzCos2z
wTzSin
ez
z
)1z(
zT
2
2
bT
2
)wt(Cos
)wt(Sin
e
t
bt
Bazı Örneklenmiş Fonksiyonların z- Dönüşümleri
Örneklenmiş Fonksiyon
f(kT)
n
k
kn zkyzYznky
1
)()()(
1-) Toplama-Çıkarma ve sabit çarpan
)z(bY)z(aY)k(by)k(ay 2121
2) Geciktirme:
)z(Yz)nk(y n
3) Öngörme:
)z(Yz)nk(y n
1
0
)()()(n
k
kn zkyzYznky
İki yönlü dönüşüm
Bir yönlü dönüşüm
İki yönlü dönüşüm
Bir yönlü dönüşüm
Bazı Önemli z – Dönüşüm Teoremleri
6) Karmaşık Dönüşüm (üstel ile çarpma)
a
zF)z(F)k(fa)k(f 11
k
4) İlk Değer Teoremi
5) Son Değer Teoremi
)z(YLim)k(yLim)0(yz0k
)z(Y)1z(Lim)k(yLim)(y1zk
7) Rampa ile çarpma
)z(F)dz
dz()z(F)k(fk)k(f n
1n
)ez/(z)z(F 2
k2k5)k(f
0)(f,1)0(f
k4Sine)k(f k2
Örnek: İlk ve son değerlerini bulunuz.
Örnek: Verilen fonksiyonun z dönüşümünü bulunuz.
Örnek:Verilen fonksiyonun z dönüşümünü bulunuz.
422
2
42
4)(
eCoszez
SinzezF
2)2z(
z25)z(F
2
k
)az(
za)z(Fise)a(k)k(f
Genel hali
Uygulamalar
Örnek: Verilen fark denkleminin z-dönüşümünü kullanarak
çözümünü inceleyiniz.
k)5.0()k(uve0)2(y,1)1(y
)k(u)2k(y64.0)1k(y6.1)k(y
Çözüm
2
3
21
)8.0)(5.0()(
)()(64.0)(6.1)(
zz
zzY
çözümBuradan
zUzYzzYzzY
a-) İki yönlü z-dönüşümü başlangıç koşullarını dikkate almadığına göre,
Fark denklemlerinin z-dönüşümünü alma
2
2
2
3
121
)8.0(
64.06.1
)8.0)(5.0()(
5.0)]2()1()([64.0)]1()([6.1)(
z
zz
zz
zzY
çözümBuradan
z
zyyzzYzyzYzzY
b-) Tek yönlü z-dönüşümünde başlangıç koşulları da
dikkate alınmalıdır.
Çözümdeki ilk terimin kaynağın,
ikinci terimin ise başlangıç koşullarının
ortaya çıkardığı çözüm bileşenleri olduğu görülmektedir.
c-) Önceki örnekte verilen fark denklemi k=k+2 yazılarak
ötelenebilir ve elde edilen denkleme tek ya da iki yönlü z-
dönüşümü uygulanabilir. k=k+2 yazılarak
)2k(u)k(y64.0)1k(y6.1)2k(y
y(0) ve y(1) gibi yeni başlangıç koşullarını da belirlemek gerekir.
02.4)1(yburadan5.0)1(y64.0)0(y6.1)1(y
için1k
6.2)0(yburadan1)2(y64.0)1(y6.1)0(y
için2k
Bu denklemin iki yönlü z-dönüşümü ile çözümünün (a) şıkında
bulunan çözüm ile aynı olduğu görülmektedir.
Yeni başlangıç koşulları ile ötelenmiş fark denklemi,
02.4)1(y,6.2)0(y
)2k(u)k(y64.0)1k(y6.1)2k(y
Olduğundan çift yönlü z-dönüşümü uygulanırsa,
)()(64.0)(6.1)( 22 zUzzYzzYzYz
2
3
)8.0)(5.0()(
zz
zzY
d-) Yeni başlangıç koşulları ile ötelenmiş fark denklemi,
02.4)1(y,6.2)0(y
)2k(u)k(y64.0)1k(y6.1)2k(y
Olduğundan tek yönlü z-dönüşümü uygulanırsa,
)1(zu)0(uz)z(Uz
)z(Y64.0)]0(zy)z(zY[6.1)1(zy)0(yz)z(Yz22
22
Elde edilir. Verilen u(k) kullanılarak u(0)=1 ve u(1)=0.5 bulunur.
Bu denklem düzenlenirse, bu denklemin çözümünün (b) şıkkında
elde edilen çözüm ile aynı olduğu görülür.
2
2
2
3
)8.0(
64.06.1
)8.0)(5.0()(
z
zz
zz
zzY
2.2 Ters z – Dönüşümü
.edilireldea)k(fBuradan
z)k(f)z(F
...zazaaz1)z(F
k
0k
k
33221
Ters z-dönüşümünde genellikle güç serisine açma yöntemi ya dakısmi kesirlere ayırma yöntemi kullanılır. Ters dönüşüm bağıntısındanyararlanarak çıkarılan rezidü yöntemi de basit kesirlere ayırmayöntemine benzer bir anlayışla uygulanabilir.
Kuvvet Serisi Yöntemi: Bölme işlemi yapılarak negatif z kuvvetlerinden
oluşan bir seri elde edilir.
Örnek
az
z)z(F
Fonksiyonunun ters z-dönüşümünü
bulmak için bölme işlemi yapılırsa,
Örnek
2z3z
z)z(F
2
12)k(f
z)k(f)z(F
...z15z7z3z)z(F
k
1k
k
4321
Fonksiyonunun ters z-dönüşümünü bulmak için bölme işlemi
yapılırsa,
Kısmi Kesirlere Ayırma Yöntemi: Verilen karmaşık ifade basit
kesirlerine ayrılarak her bir basit kesrin ters z-dönüşümü
bulunur. Ancak z-dönüşüm tablosuna dikkat edilirse basit
terimlerin z-dönüşümlerinin payında her zaman z çarpanı vardır.
Bu nedenle, Z bölgesinde verilen bir F(z) ifadesi, öncelikle
F(z)/z haline getirilerek basit kesirlerine ayrılmalıdır.
)ez).(1z(
z)e1()z(F
aT
aT
Örnek
aTez
z
1z
z)z(F
,...2,1,0k,e1)k(f akT
Payda kökleri reel ve ayrık olduğundan basit kesirleri,
Olduğuna göre
)2).(1(
1)(
zzzE
2z2
1
1z
1
z2
1
)2z)(1z(z
1
z
)z(E
z
z
)2z).(1z(
1)z(E
k22
1)k(u)k(
2
1)k(e d
Payda kökleri reel ve ayrık olduğundan basit kesirleri,
Örnek
22
1
12
1)(
z
z
z
z
z
zzE
Olduğuna göre
)1.()2(
2)(
2
3
zz
zzzX
2)2(1)1.()2(
12)(
)1.()2(
)12.()(
22
2
2
2
z
C
z
B
z
A
zz
z
z
zX
zz
zzzX
3)1.()2(
12.1 |
12
2
zzz
zzA 9
)1.()2(
12.)2( |
22
22
zzz
zzB
1)1.()2(
12.)2(.
!1
1|
22
22
zzz
zz
dz
dC
2)2(
9
1
3)(
2
1
)2(
9
1
3)(22
z
z
z
z
z
zzX
zzzz
zX
)k(u322k2
9)k(x kk
Örnek
Payda kökleri reel ve çakışık olduğundan basit kesirleri,
Karmaşık Kök DurumuKarmaşık kök durumunda da basit kesirlere ayırma yöntem
kullanılabilir. Ancak z-dönüşüm tablosuna dikkat edilirse karmaşık
kökler sinüsoidal fonksiyonlarda ortaya çıkmaktadır. Bu nedenle,
karmaşık kökü olan z-bölgesindeki fonksiyonlar saf sinüsoidal ya
da sönümlü sinüsoidal fonksiyonlara benzetilebilir. Bu amaçla
aşağıdaki z-dönüşümleri tekrar hatırlanabilir.
22
k
22
k
2
2
awazCos2z
)waCosz(z)z(Fise)wk(Cosa)k(f
awazCos2z
wazSin)z(Fise)wk(Sina)k(f
1wzCos2z
)wCosz(z)z(Fise)wk(Cos)k(f
1wzCos2z
wzSin)z(Fise)wk(Sin)k(f
NOT: wp/2 için payda köklerinin saf sanal olduğu görülmektedir.
)k2
sin()k(fise1z
z)z(F)1
2
p
Örnekler
)k2
cos()k(fise1z
z)z(F)2
2
2 p
)k3
(Sin577.0)k3
(Cos)k(f
Buradan1zz
)3/(zCos
1zz
)3/(zCosz
1zz
)w(zCos)w(zCosz)z(F
.olur3/wvedirSinüsoidal
?)k(fise1zz
z)z(F)3
22
2
2
2
2
2
pp
pp
p
?)k(f4z2z
z)z(F)4
2
2
kSinkCoskf kk
3
2)2(
3
1
3
2)2()(
pp
Üstel ve sinüsoidal çarpımıdır ve a=2, w2p/3 olur. Bu
durumda örnek 3 e benzer şekilde düzenlenirse,
0)k(y)1k(y
Örnek: Başlangıç koşulu y(0) olan aşağıdaki farkdenklemini z dönüşümü ile çözünüz.
0)z(Y)0(y)z(Yz
Çözüm:
)0(1
)( yz
zzy
,...2,1,0k),0(y1)k(y k
Ters z-dönüşümü ile
2.3 Fark Denklemlerinin z–Dönüşümü ile Çözümü
Örnek: Başlangıç koşulları verilen aşağıdaki fark denklemini z
dönüşümü ile çözünüz.
1)1(x,0)0(x,0)k(x2)1k(x3)2k(x
0)z(X2)0(zx3)z(zX3)1(zx)0(xz)z(Xz 22
Çözüm:
2z
z
1z
z
2z3z
z)z(X
2
kk )2()1()k(x
Zamanda kaydırarak k=k-2 ötelemesi yapılırsa
0)2(.2)1(.3)( kxkxkx
75.0)2(x,5.0)1(x
Geriye ötelenmiş bu fark denklemi de z-dönüşümü ile çözülürse aynı
sonuç elde edilir.
Fark denklemi elde edilir ve bu denklem k=1 ve k=0 için ya da ilk
verilen denklem k=-1 ve k=-2 için çözülerek yeni başlangıç koşulları
bulunur.
0k0)k(x,1)2k(x)1k(x2)k(x2
Örnek: Başlangıç koşulları verilen aşağıdaki farkdenklemini z
dönüşümü ile çözünüz.
Çözüm:
1z
z)z(Xz)z(Xz2)z(X2 21
502
1
1 2
2
.zz
zz
z
z)z(X
50
50
2
1
50
50
2
1
1 22
2
.zz
z.
.zz
z.z
z
z)z(X
,...,,k)k(Sin)k(Cos)k(x
kk
21042
1
2
1
42
1
2
11
pp
42
1
2
1 1 p Cosw,a
İkinci terim, üstel ile
sinüsoidal çarpımıdır ve
Başlangıç koşulları sıfır alınmak kaydıyla z-bölgesinde çıkışın
girişe oranı transfer fonksiyonunu belirler.
n
n
1
1
m
m
1
10
m
m0
n
n
1
1
m0n1
za...za1
zb...zbb
)z(U
)z(Y)z(G
)z(Uzb...)z(Ub)z(Yza...)z(Yza)z(Y
)mk(ub..)k(ub)nk(ya)1k(ya)k(y
2.4 Transfer Fonksiyonu
G(z)
Y(z)U(z)
?G(z)
)1k(u2)k(u)1k(y3)k(y4)1k(y
Örnek: Fark denklemi verilen sistemin transfer fonksiyonunu
bulunuz.
Örnek: Transfer fonksiyonu verilen sistemin fark denklemini
çıkarınız.
5z4z
3z
)z(U
)z(Y)z(G
2
2.4.1 Transfer Fonksiyonu Verilen Sistemin
verilen bir giriş için Cevabını Bulma
)z(U
)z(Y)z(G )z(U)z(G)z(Y Olduğuna göre
Olur ve sistemin transfer fonksiyonu ile sisteme uygulanan giriş
verilmişse ters z dönüşümü alınarak çıkış bulunur.
)3z)(2z(
2z)z(G
Örnek:TF verilen sistemin birim basamak cevabını bulunuz.
1)3)(2(
5)(
z
z
zz
zzY
U(z) Y(z)G(z)
kkkky )3(5
4)1()2(
5
1)( Sonuç:
2.4.2 Transfer Fonksiyonunun kutupları, sıfırları,
z-düzlemi, birim çember ve kararlılık
Örnekler
Sürekli zaman sistemlerdeki tanımlar ayrık zamanda da geçerli
olduğundan burada sadece örnekler verilecektir.
az
zzG
)(
2)3.0z)(1z(
)2z(z)z(G
25.0z5.0z
z25.0z)z(G
2
2
2.5 Blok Şemalar
Sürekli zaman kontrol sistemlerinin analizinde geçerli olan blok
şema kuralları, sinyalleri ve sistemleri ayrık zamanda
tanımlamak koşuluyla ayrık zamanda da geçerlidir. Ara sinyaller
yazılıp yok edilerek ya da sinyal kaydırarak indirgenebilir.
Örnekler:
G2
R YG1
+-
G3
H2
H1
+
-
Örnek: Transfer fonksiyonu verilen kapalı çevrim kontrol
sisteminin birim impuls cevabını bularak grafiğini çiziniz.
1.0
8.0)(,
5.0)(
zzH
z
zzG
2.6 Sinyal Akış Şemaları
Sürekli zaman kontrol sistemlerinin analizinde geçerli olan
sinyal akış şeması kuralları, sinyalleri ve sistemleri
ayrık zamanda tanımlamak koşuluyla ayrık zamanda da
geçerlidir. Kazanç formülü ile indirgenebilir.
Örnekler:
R 1 Y
X2
o oo oo oo oo oX1 X3 X4
G1 G2 G3 G4
-H1 -H3
-H2
YR 1 Yo oo oo o
G1 G2
-H2-H1
o o1
N
1k
kkT
sR
sYsT
.
)(
)()(
Örnek: Verilen sinyal akış şemasını indirgeyiniz.
R 1 Yo oo oo oo oo o
G1 G2 G3 G4
-H1 -H2
G5
-H3
X
Örnek: Verilen sinyal akış şemasını indirgeyiniz.
2.7 Durum Diyagramları
Sürekli zaman kontrol sistemlerinin analizinde geçerli olan
Durum diyagramı kuralları, sinyalleri ve sistemleri
ayrık zamanda tanımlamak koşuluyla ayrık zamanda da
geçerlidir.
NOT: Başlangıç koşulları da durum diyagramlarında
gösterilebilir ancak burada başlangıç koşulları sıfır alınacaktır.
Örnekler: Durum diyagramlarını a-) zaman bölgesinde b-) z-
dönüşümünü alarak z bölgesinde ve ayrıca sinyal akış
şeması/blok şema gösterimi ile çiziniz.
)k(u5)k(y3)1k(y4)2k(y
)k(u10)k(y)1k(y
2.7.1 Ayrık Zaman Durum Denklemleri ve z-
Dönüşümü ile Çözümü
)()()(
)()()1(
kDukCxky
kBukAxkx
)z(BU)AzI())0(x)AzI(z)z(X
)z(BU)0(zx)z(XAzI
)z(BU)z(AX)0(zx)z(zX
11
z-dönüşümü ile çözümü: Durum değişkenlerinin başlangıç
değeri x(0) vektörü olsun.
u(k) 0
5
)k(x
)k(x
01
23
)1k(x
)1k(x
2
1
2
1
0
2x(0)
k)2()k(u
)2z(
z5)2z(
z5
z2
z2
)z(X
2z
z
0
5
3z1
2z
0
2
)3z(zz
z2z
)z(X
3z1
2z
)AzI(3z1
2z
)2z)(1z(
1)AzI(
22
2
11
Örnek: Verilen durum denklemini z-dönüşümü yardımıyla çözünüz.
Çözüm
2.7.3 Durum Denkleminden
Transfer Fonksiyonuna Dönüşüm
)()()(
)()()1(
kDukCxky
kBukAxkx
Tek giriş ve tek çıkışlı ayrık zaman bir sistemin durum denklemi,
Başlangıç koşulları sıfır alınarak z dönüşümü alınır,
)z(DU)z(CX)z(Y
)z(BU)AzI()z(X
)z(BU)z(AX)z(zX1
DB)AzI(C)z(U
)z(Y)z(G)z(UDB)AzI(C)z(Y 11
ve çıkış/giriş oranı oluşturulursa,
4u(k)(k)x
(k)x 32)k(y
)k(u1
2
)k(x
)k(x
04
23
)1k(x
)1k(x
2
1
2
1
2
1
41
2
)3z(3412z2)z(T4
1
2
3z4
2z
32)z(G
8z3z
5z19z4)z(G
48z3z
37z74
37z74
)3z(34)12z2(2)z(G
2
2
2
Çözüm:
ÖrnekG(z)=?
8z3zBurada3z4
2z1)AzI( 21
Buradan
)(2)(5)1()2(2)3(3 kukykykyky
)3k(y)1k(x
)2k(y)1k(x)k(x
)1k(y)1k(x)k(x
)k(y)k(x
3
23
12
1
)k(u
3
20
0
)k(x
)k(x
)k(x
3
2
3
1
3
5100
010
)1k(x
)1k(x
)1k(x
3
2
1
3
2
1
)k(x
)k(x
)k(x
001)k(y
3
2
1
2.7.4 Fark Denkleminden Durum Denklemine Dönüşüm
Verilen fark denkleminde atanan
değerler yerlerine yazılırsa;
Buna göre Kanonik formda durum denklemi
Üçüncü dereceden bir sistem olduğundan 3 tane durum değişkeni vardır ve
aşağıdaki gibi atamalar yapılabilir.
Aşağıdaki fark denklemi verilmiş olsun. NOT: Çıkış sadece u(k) olursa
)k(u2)k(x5)k(x)k(x2)1k(x3 1233
)k(u4)k(y3)1k(y2)2k(y
)k(x
)k(x 01)k(y),k(u
4
0
)k(x
)k(x
23
10
)1k(x
)1k(x
2
1
2
1
2
1
Örnek:
Buna göre durum denklemi
)2k(y)1k(x
)1k(y)1k(x)k(x
)k(y)k(x
2
12
1
Durum değişkeni atamaları yapılırsa,
Çözüm:
)k(x2)k(x3)k(u4)1k(xBuradan
)k(u4)k(x3)k(x2)1k(x
212
122
2.7.5 Transfer Fonksiyonundan
Durum Denklemlerini Elde Etme
Sürekli zaman kontrol sistemlerinde incelendiği gibi ayrık zaman
kontrol sistemlerinde de verilen bir transfer fonksiyonundan
durum denklemi,
a-) Kanonik
b-) Seri
c-) Paralel
yöntemlerle elde edilebilir. Yöntemler önceden bilindiğinden
burada örnekler verilecektir.
5.02
12)(
23
2
zzz
zzzG
)(
)(.
5.021
2
)(
)(..
5.02
12
)(
)(321
321
3
3
23
2
zE
zE
zzz
zzz
zE
zE
z
z
zzz
zz
zU
zY
)()(2)()( 321 zEzzEzzEzzY
)(5.0)()(2)()( 321 zEzzEzzEzzEzU
)z(Ez5.0)z(Ez)z(Ez2)z(U)z(EBuradan 321
Kanonik Yöntem
Örnek: verilen transfer fonksiyonuna ait işaret akış
diyagramını çiziniz ve durum denklemlerini çıkarınız.
Çözüm
U(z)
1
E(z) z-1
-2
-1
-0,5
z-1
1
2
1
Y(z)
zx1 (z)
zx2 (z)
x1 (z)
z-1
zx3 (z)
x2(z)
x3(z)
)k(x5.0)k(x)k(x2)k(u)1k(x 3211
)k(x)1k(x 12
)()1( 23 kxkx )k(x)k(x2)k(x)k(y 321
)k(x
)k(x
)k(x
121)k(y),k(u
0
0
1
)k(x
)k(x
)k(x
010
001
5.012
)1k(x
)1k(x
)1k(x
3
2
1
3
2
1
3
2
1
Sinyal akış şeması kullanılarak durum değişkenleri ve çıkış
)1z)(8.0z)(5.0z(
2z)z(G
Seri Yöntem
Örnek: verilen transfer fonksiyonuna ait işaret akış diyagramını çiziniz
ve durum denklemlerini çıkarınız.
Çözüm
)1z(
2z
)8.0z(
2
)5.0z(
z)z(G
Paralel Yöntem
Örnek: verilen transfer fonksiyonuna ait işaret akış diyagramını çiziniz
ve durum denklemlerini çıkarınız.
Çözüm
-6
T
2
5
T
3
T
4
-2
2.7.6 Zamandaki Blok Diyagramından Durum Denklemi
)k(x2)k(x)k(x)k(y),k(x4)k(u)1k(x
)k(x10)k(x3)k(u)1k(x),k(x6)k(x2)k(u)1k(x
32133
322211
T T T
-2
-3
-1
1
(k)x(k)xy(k)
(k)x1)(kx
(k)x1)(kx
(k)x-(k)3x-(k)2x-u(k)1)(kx
32
23
12
3211
Örnek
Bölüm 3
Ayrık Zaman (Sayısal ) Kontrol
Sistemlerinin Analizi
3.1 Ayrıklaştırma Yöntemlerinin Laplace Karşılığı ve
Kararlılığı
3.2 Sayısal Kontrol Sistemlerinde Kararlılık Analizi
3.3 Sayısal Kontrol Sistemlerinin Tipi ve Kalıcı Durum
Hataları
3.4 Sayısal PID kontrolörler
3.5 Ayrık Zaman Kontrol Sistemlerinin Analizi
3.1 Ayrıklaştırma Yöntemlerinin Laplace Karşılığı ve
Kararlılığı
T
kxkxtx
)()1()(
)z(XT
1z)s(sX
T
)z(X)z(zX)s(sX
T
zs
1 sTz 1
İleri Fark Yöntemi ya da Sol Kenar Kuralı
olduğuna göre karşılıklı laplace ve z dönüşümleri alınırsa,
Buradan ileri fark yönteminde iki düzlem arasındaki ilişkinin
aşağıdaki gibi olduğu görülür.
Ya da
0)Re( s
01
Re
T
z jyxz
01
Re
T
jyx
01
Re
T
yj
T
x
101
xT
x
Im(s)
Re(s)
Kararlı bölge
Kararlı bölge
Re(z) 1
Im(z)
İleri fark yöntemi ile bir sistemin ayrıklaştırılması neticesinde
kararlılığın nasıl değiştiği aşağıdaki gibi belirlenebilir.
Laplace bölgesinde kararlılık,
Olduğuna göre İleri fark yöntemine göre s yerine ayrık karşılığı yazılırsa z-düzlemindeki kararlı bölge
Burada karmaşık z değişkeni
İle gösterilirse kararlı bölge
Geri Fark Yöntemi ya da Sağ Kenar Kuralı
karşılıklı laplace ve z dünüşümleri alınırsa,
Buradan geri fark yönteminde s ve z düzlemi arasındaki ilişki,
Geri fark yöntemi ile ayrıklaştırılan sistemlerin kararlılığı için,
Bu ifade düzenlenirse,
T
kxkxtx
)1()()(
)z(XT
z1)s(sX
T
)z(Xz)z(X)s(sX
11
Tz
1zs
0Tz
1zRe
0
1.
1Re
jyx
jyx
T
22
2
2
1
2
1
yx
Im(s)
Re(s)
Kararlı bölge
Re(z)
Im(z)
Merkezi Fark Yöntemi ya da Yamuk Kuralı (Bilineer Dönüşüm)
)1()(2
)1()()()( kxkxT
kykydttxty
)s(Xs
1)s(Y )z(Xz)z(X
2
T)z(Yz)z(Y 11
)z(Xz1
z1
2
T)z(Y
1
1
1z
1z
T
2s
Olduğundan karşılıklı laplace ve z dünüşümleri alınırsa,
ve
Olur ve yamuk kuralı ile
ayrıklaştırmada s ve z
düzlemi arasındaki ilişki;
Yamuk kuralı (bilineer dönüşüm) ile ayrıklaştırmanın kararlılığı;
Im(s)
Re(s)
Kararlı bölge
Re(z)
Im(z) 0
1
1.
2Re
z
z
T
01
1.
2Re
jyx
jyx
T
122 yx
Örnek: Verilen sistemleri farklı yöntemlerle ayrıklaştırarak,
a-) s ve z düzlemlerindeki kutup-sıfır yerlerini ve kararlılıklarını
değerlendiriniz.
b-) hem analog hem de ayrıklaştırdığınız sayısal sistemin birim
basamak cevabını bulunuz.
NOT: Uygun örnekleme peryodunu seçiniz.
)1(
2)(
2
1)(
sssG
ssG
3.2 Ayrık Zaman Kontrol Sistemlerinde
Kararlılık Analizi
• Ayrık zaman sistemlerin kararlı olabilmesi için karakteristik
denkleminin köklerinin birim çember içinde olması gerektiği
bilinmektedir.
• Birim çember dışında bulunan herhangi bir kutup sistemi
kararsız yapar.
• Yüksek dereceden sistemlerin köklerini bulma zorluğu
nedeniyle bunun yerine Routh-Hurwitz ya da Jury kriteri ile
ayrık zaman sistemlerin kararlılığı belirlenebilir.
• Routh-Hurwitz kriteri ayrık zaman sistemlerin kararlılık
analizine doğrudan uygulanamaz. Ancak, bilineer dönüşüm
kullanılarak verilen ayrık zaman sistemin sürekli zaman
karşılığı yaklaşık olarak bulunabilir ve bulunan analog
sisteme Routh-Hurwitz kriteri uygulanabilir.
• Burada Jury kriterini açıklamak yeterli olacaktır.
JURY Kriteri ile Kararlılık Analizi
0...)( 11
10
nnnn azazazazP
Jury kriteri ile yüksek dereceden sistemlerin doğrudan z-
bölgesinde kararlılık analizi yapılabilir. Sistemin karakteristik
denklemi P(z) elde edilirse aşağıdaki Jury tablosu oluşturulabilir.
Jury tablosundaki katsayılar, aşağıdaki determinantlarla bulunur.
1,...,2,1,0,10
1
nk
aa
aab
k
knnk 2,...,2,1,0,
10
21
nk
bb
bbc
k
knnk
JURY Kararlılık Kriteri
a0>0 olmak koşuluyla bir n. dereceden bir sistemin kararlı
olabilmesi için aşağıdakilerden sıra ile n+1 adet koşulun
sağlanması gerekir.
NOT: Şartlardan biri sağlanmazsa sistem kararsızdır. P(-1)=0 ya
da P(1)=0 olduğunda diğer şartlar sağlanıyorsa marjinal
kararlıdır (yani z=1 de kutbu vardır), aksi halde kararsızdır.
008.03.007.02.1)( 234 zzzzzP
Örnek: Karakteristik denklemi verilen sistemin kararlılığını
JURY kriteri ile inceleyiniz.
Çözüm:
sistem kararlıdır.
02.01.01.1)( 23 zzzzP
Örnek: Karakteristik denklemi verilen sistemin kararlılığını
JURY kriteri ile inceleyiniz.
Çözüm:
Sistem marjinal kararlıdır.
3679.0z3679.1z
)2642.0z3679.0(K)z(G
2p
Örnek: Kapalı çevrim kontrol sisteminin kararlı olabilmesi için
K’ nın alabileceği değer aralığını JURY kriteri ile belirleyiniz
02642.03679.0)3679.13679.0()( 2 KzKzzP
177.5K3925.21|K2642.03679.0|,a|a|)1 02
0K0K6321.0)1(P
K6321.0K2642.03679.03679.1K3679.01)1(P)2
K383.260K1073.07558.2)1(P)3
03925.2 KSonuç:
Çözüm: Karakteristik denklem,
R(z) E(z) Y(z)+
-
U(z)K Gp(z)
3.3 Ayrık Zaman Kontrol Sistemlerinin Tipi
ve Kalıcı Durum HatalarıSayısal kontrol sistemlerindeki kalıcı durum hataları da
analog kontrol sistemlerindeki İncelemelere benzer şekilde
yapılabilir.
Şekildeki birim geri beslemeli kontrol sisteminde ileri besleme
transfer fonksiyonu G(z), kutup ve sıfırları cinsinden aşağıdaki
gibi yazılabilir.
)pz()z(
)zz(K
)z(G
i
r
j
1
Burada, zj- açık çevrim sistemin sıfırları, pi-açık çevrim sistemin
kutuplarıdır ve r ise sistemin tipini gösterir. Analog sistemlerde
olduğu gibi sayısal kontrol sistemlerinde de kalıcı durum hataları,
sistemin tipine göre değerlendirilebilir.
R(z) +
-
Y(z)
G(z)
Birim rampa giriş için olduğundan21)z/(z)z(R
)z(G)z(limKBuradaK)z(Glim)z()z(
z
)z(G
)z(lime
zv
vz
zss 1
1
10
1
11
11
1
21
)z(GlimKBuradaK)z(Glimz
z
)z(G
)z(lime
zP
pz
zss
1
1
1 1
1
1
1
11
1
Birim basamak giriş için R(z)=z/(z-1) olduğundan
)z(R)z(G
)z()z(E)z(lime
zss
1
11
1
Kalıcı durum hatası (ess), son değer teoreminden,
Bu ifadeden yararlanarak basamak ve rampa girişler için konum
ve hız hata katsayıları aşağıdaki gibi hesaplanabilir.
Buna göre, Tip 0 sistemlerde sonlu basamak hatası
Tip 1 sistemlerde ise sonlu rampa hatası ortaya çıkar.
)z(R)z(G1
1)z(Y)z(R)z(E
Şekilde Hata fonksiyonu,
Örnek: Verilen kapalı çevrim kontrol sisteminde her bir G(z)
için sistemin tipini ve ilgili kalıcı durum hatalarını bulunuz.
9.0z2.0z
1.0z)z(G
5.0z
z)z(G
2
3.4 Sayısal PID Kontrolörler
Analog PID kontrolörler, çeşitli ayrıklaştırma yönemleri kullanılarak
ayrık zamana dönüştürülebilir. Dönüştürme zaman bölgesinden ve
Laplace bölgesinden yapılabilir.
PI Kontrolör: PI kontrolörlerin, genel olarak sistemin tipini
artırarak kalıcı durum hatalarını iyileştirdiği bilinmektedir.
Analog PI kontrolör,
dtteKteKtu ip )()()(
R(z)Gc(z)
Y(z)+
-Gp(z)
H(z)
U(z)E(z)
Örneğin, integral sağ kenar kuralı ile ayrıklaştırılırsa,
)()1()()()()( kTekukukuKkeKku iiiip
Elde edilir ve buna göre programlanabilir.
MATLAB programını yazınız.
Ayrık zaman PI kontrolörün fark denkleminin Z dönüşümü alınırsa,
(Alınız) PI kontrolörün transfer fonksiyonu,
1
)(
1)(
z
azK
z
zTKKzG c
ipc
Bu ifadeden de sayısal PI kontrolörün, sistemin Tipini bir
artırarak kalıcı durum hatasını İyileştireceği ayrıca, sisteme bir
sıfır ekleyeceği ve sistemin kazancını ayarlayacağı
görülmektedir.
Laplace bölgesinde analog kontrolörlerin ayrık zamana
dönüştürülmesi matematiksel olarak,
|)(
)()(zFs
cc sGzG
şeklinde ifade edilebilir. Burada F(z) “ayrıklaştırma
fonksiyonu” olarak adlandırılır. Önceki bölümlerde ileri
fark, geri fark ve merkezi fark (bilineer dönüşüm) yöntemi
için F(z) fonksiyonları belirlenmişti.
1
)()(
1)( |
/)1(
z
azKzG
z
zTKK
s
KKsG c
cip
Tzzs
ipc
Buradan analog PI kontrolör geri fark ile ayrıklaştırılırsa,
PD Kontrolör: Kontrol sisteminin geçici rejim kriterlerini (aşma-
yerleşme-yükselme zamanı) düzeltir. PD kontrolör
)]1()([)()( kekeT
KkeKku d
p
Olduğuna göre geri fark ile ayrıklaştırılırsa
)()()( tedt
dKteKtu dp
Fark denklemi programlanabilir. MATLAB ile programlayınız.
Ayrık zaman PD kontrolörün fark denkleminin Z dönüşümü alınırsa,
(Alınız) YA DA Laplace bölgesinden PD kontrolörün transfer fonksiyonu,
z
azKzGK
Tz
zKKsKsG c
cpd
Tzzs
pdc
)()(
1)( |
/)1(
PID Kontrolör: Kontrol sisteminin geçici+kalıcı rejim kriterlerini
(aşma-yerleşme-yükselme zamanı) + kalıcı durum hatası düzeltir.
PID kontrolör,)()()()( te
dt
dKdtteKteKtu dip
geri fark yöntemi ile ayrıklaştırılırsa,
)()1()(
)],1()([)()()(
kTekuku
kekeT
KkuKkeKku
ii
diip
Fark denklemi programlanabilir. MATLAB ile programlayınız.
Ayrık zaman PID kontrolörün fark denkleminin Z dönüşümü alınırsa,
(Alınız) YA DA Laplace bölgesinden PID kontrolörün transfer fonksiyonu,
)1(
))(()(
1
1)( 21
/)1(
|
zz
azazKzG
Tz
zK
z
TzKKsK
s
KKsG c
cdi
p
Tzzs
di
pc
3.5 Ayrık Zaman Kontrol Sistemlerinin Analizi
Ayrık zaman transfer fonksiyonu verilen ya da elde edilen bir
sisteme sayısal PI, PD ya da PID bağlanarak çeşitli referans
girişler için cevapları analog kontrol sistemlerdekine benzer
şekilde bulunabilir.
R(z)Gc(z)
Y(z)+
-Gp(z)
H(z)
U(z)E(z)
Örnekler:
1-) Şekildeki ayrık zaman kontrol sisteminde transfer
fonksiyonları verilmiştir. Kontrol sisteminin birim basamak
cevaplarını bularak grafiklerini çiziniz.
R(z)Gc(z)
Y(z)+
-Gp(z)
H(z)
U(z)E(z)
9.0
1.0)(,1)(,
1
)7.0(10)(
zzGzH
z
zzG pc
NOT: Ayrık zamanda verilen bir transfer fonksiyonunun
sönüm oranı ve doğal frekansı belirlenebilir mi???
Örnekler:
1-) Şekildeki ayrık zaman kontrol sisteminde transfer
fonksiyonları verilmiştir. Kontrol sisteminin birim basamak
cevaplarını bularak grafiklerini çiziniz.
R(z)Gc(z)
Y(z)+
-Gp(z)
H(z)
U(z)E(z)
5.0
1.0)(,1)(),475.2(4)(
2
zzzGzHzzG pc
Sayısal kontrolörlerin tasarımında iki yaklaşım
izlenebilir.
• Kontrolör, önce analog olarak tasarlanır ve sonra
ayrıklaştırılır.
• Önce sistem ayrıklaştırılır ve sonra kontrolör ayrık
zamanda tasarlanabilir. (bu yöntem sonraki bölümlerde
açıklanacaktır)
Ayrık zaman kontrolör tasarımı
s
s)s(Gc
210
1
334015
14
120140
1
14
21
1410
210
1
12
z
).z(
z
zz
z
z
z
z
s
s)z(G |
z
z
Ts
c
Örnek: Verilen analog sistem için analog bir PI kontrolör
tasarlanmış olsun. Bilineer dönüşüm yöntemiyle ayrık zaman
eşdeğerini bulunuz ve MATLAB programını yazınız. T=0.5 s.
R(s) E(s) Y(s)Gc(s)
+
-
Gp(s)U(s)
1
1
1
515
1
515
1
)334.0(15
)(
)()(
z
z
z
z
z
z
zE
zUzGc
)1(5)(15)1()(
)(515)(1 11
kekekuku
zEzzUz
Öncelikle transfer fonksiyonundan fark denklemi elde edilir.
• İçler dışlar çarpımı yapılarak sayısal kontrolörlerin fark denklemi,
Sayısal kontrolörlerin bilgisayar programlama dili ile
programlanması
Programlanır. MATLAB programını yazınız.
Örnek: Şekilde verilen sayısal kontrol sisteminde sönüm
oranı 0.8 ve sönümsüz doğal frekansı 10 olacak şekilde PI
kontrolörü tasarlayınız.
8010 .,sn/radwn
Çözüm: Kontrolör sürekli zamanda tasarlanarak
ayrıklaştırılabilir. Analog sistem için analog bölgede PI kontrollü
sistemin karakteristik denklemi, standart 2. dereceden referans
sistemin karakteristik denklemine eşdeğer yapılarak kontrolör
kolayca tasarlanabilir. Sürekli zamanda PI kontrollü sistemin
AÇTF,
2s
1
s
)bs(K)s(G)s(G pc
R(s) Y(s)Gp(s)Gc(s)
PI kontrollü sistemin kapalı çevrim TF, standart 2. dereceden
sisteme eşitlenirse,
2
nn
2
2
n
2
pc
pc
wsw2s
w
Kbs)K2(s
)bs(K
)s(G)s(G1
)s(G)s(G)s(T
s
).s()s(GBuradan c
14714
Tasarlanan sürekli zaman kontrolör bilineer dönüşümle T=0.1 s
ile ayrıklaştırılırsa sayısal PI kontrolör aşağıdaki gibi bulunur.
1z
)474.0z(19
s
)14.7s(14)z(G |
1z
1z
T
2s
c
NOT: Bulunan PI kontrolörün programını yazınız.
Örnek: Verilen analog sistem için analog kontrolör
tasarlayarak sayısala dönüştürünüz.
)2s(s
1)s(Gp
R(s) Y(s)Gp(s)Gc(s)
8010 .,sn/radwn
Bölüm 4
Örneklenmiş Verili
Açık Çevrim Kontrol Sistemleri
• 4.1 Örnekleme ve tutma
• 4.2 Laplace Bölgesinden Yıldızlanmış Dönüşüm
• 4.3 s ve z- Düzlemleri Arasındaki İlişki
• 4.4 Açık çevrim ayrık zaman kontrol sistemleri
Referans
giriş
r(k)
Sistem
çıkışı
y(t)
Algılayıcı
+
-
YükselticiSayısal
Denetleyici
Hata
e(k)
Denetim
sinyali
v(t)
Sistem
girişi
u(t)
Sistem
ADC
DAC
Tutma
Örnekleme
4.1 Örnekleme ve Tutma
Bölüm 1 de blok şeması verilen aşağıdaki kapalı çevrim sayısal
kontrol sistemlerinin analiz ve tasarımında Örnekleme ve Tutma
işlemleri önemli bir yer tutar.
Şekildeki gibi sürekli zaman bir sinyalin T örnekleme peryodu ile
örneklendiğini ve sıfırıncı dereceden bir tutma devresi ile sürekli
zamana geri dönüştürüldüğünü yani orijinal analog sinyalin geri
elde edilmeye çalışıldığını kabul edelim. Tutma devreleri yüksek
dereceden olabilir ancak burada sadece sıfırıncı dereceden
tutma ele alınacaktır.
Örnekleme ve Tutma
.....)]T2t(u)Tt(u)[1(x)]Tt(u)t(u)[0(x)t(x
Sıfırıncı dereceden tutulmuş sinyalin
Denklemi aşağıdaki darbe zincirleri ile
tanımlanabilir.
....es
1e
s
1)T(xe
s
1
s
1)0(x)s(X Ts2TsTs
LD alınırsa
....e)T2(xe)T(x)0(xes
1
s
1)s(X Ts2TsTs
0k
kTsTs e)kT(xes
1
s
1)s(X
Burada iki terimden birincisi
tutma, ikincisi örnekleme
işleminin karşılığıdır.
Buna göre ikinci terim ayrılırsa örnekleme işleminin matematiksel
ifadesi,
0
)()(*
k
kTsekTxsX
Olarak elde edilir. Örneklenmiş sinyallerin Z-dönüşüm ifadesi
hatırlanırsa,
0k
kTsTs e)kT(xes
1
s
1)s(X
zesT
0k
kz)kT(x)z(X
olduğu görülür.
Olduğundan yıldızlanmış sinyal ile z-dönüşümü alınmış sinyal
arasındaki ilişkinin
Dolayısıyla örnekleme işlemi, analog bir sinyali ayrıklaştırma
işlemidir ve buradan z-dönüşümüne kolayca geçilebilir.
Diğer bir ifade ile bu bağıntı yardımıyla s ve z düzlemi arasında
da bir ilişki kurulabilir.
)s(Gs
e1
)s(X
)s(Xh
Ts
*
Sonuç olarak örnekleme ve sıfırıncı dereceden tutma devresini
aşağıdaki blok şema ile gösterilebilir.
Diğer taraftan birinci terim ayrılırsa tutma devresinin transfer
fonksiyonu elde edilir.
0k
kTsTs e)kT(xes
1
s
1)s(X
te)t(f
)1s(TTsT
Ts2T2TsTkTs
0k
kT
e1
1
ee1
1)s(F
...eeee1e)kT(f)s(F
e)kT(f
TT1
sT*
ez
z
ez1
1e)s(F)z(F
Örnek: Zaman bölgesi ifadesi verilen sinyalin yıldızlı
dönüşümünü ve z dönüşümünü bulunuz.
Çözüm
NOT: Hatırlanırsa örneklenmiş sinyalin doğrudan z dönüşümünün
bu ifade olduğu belirlenmişti.
Örnek: Şekildeki gibi bir analog sistemini girişine ADC ve DAC
bağlanırsa sistemin blok şeması nasıl çizilebilir ?
4.2 Laplace Bölgesinden Yıldızlı DönüşümYukarıda, zaman bölgesindeki ifadesi bilinen sürekli zaman bir sinyalin
örneklenmişinin Laplace dönüşümü (yıldızlı dönüşümü) alınabilmekte ve
yıldızlı dönüşümden de z dönüşümüne geçilebilmektedir. Laplace
dönüşümü verilen bir sürekli zaman sinyalin, yıldızlı dönüşümü rezidü
teoremi kullanılarak aşağıdaki gibi alınabilir.
)(1
1)(Rezidü )(
sTeEsE
Sonuç olarak buradaki Rezidü katsayısı, basit kesirlere ayırma işlemindeki
katsayılara karşılık gelmektedir.
)4s)(1s(
1)s(E
)4()2( 1
1
2
1
1
1
2
1)(
sTsT eesE
1412 1
1
2
1
1
1
2
1)(
zezezE
TT
Örnek: Laplace dönüşümü verilen sinyalin
Çözüm
?)z(Eve?)s(E
)2s)(1s(s
1)s(E
)2s(T)1s(TTs
2
)s(T
1
)s(T
0
))s((T
e1
1
2
1
e1
1
e1
1
2
1)s(E
e1
1
)1(
1
e1
1
)1(
1
e1
1
)2)(1(
1)s(E
ez
z
2
1
ez
z
1z
z
2
1)z(E
ze1
1
2
1
ze1
1
z1
1
2
1)z(E
T2T
1T21T1
Örnek:Laplace dönüşümü verilen sinyalin yıldızlı dönüşümünü ve
z dönüşümünü bulunuz.
Çözüm
4.3 s ve z- Düzlemleri Arasındaki İlişkiw js 2,1
sTez
s-düzleminde olarak verilen karmaşık bir noktanın
olduğundan,
z-düzlemindeki yeri,
TjTT)j(2,1 e.eez ww
rerz j.2,1Ter Tw
Bu noktalar genlik ve faz bilgileri ile kutupsal olarak gösterilirse,
Olur. Burada
Im(s)
Re(s)
jω
-jω
-α
s düzlemi
s1 = -α + jω
s2 = -α - jω
Im(z)
Re(z)
b
a
z düzlemi
- b
r θ
θ
z-düzleminde jbaz olarak verilen karmaşık bir noktanın
s-düzlemindeki yeri de yani değerleri de bulunabilir.
jbaz 2,1 rz
a
btanvebar 122
rlnT
1er T
22)r(ln
rln
22)(1
InrT
wn
Burada,
olur. Bu eşitliklerden yararlanarak z-düzlemindeki bir karmaşık kutbun,
sönüm oranı ve doğal frekansı da hesaplanabilir.
w ve
a
btan
T
1
TT 1
ww
T1ww
T
rlnw
2n
n
Eşitlikleri
oranlanarak
düzenlenirse
olur
Örnek:
a-) noktasının z-düzlemindeki yerini bulunuz. T=0.1 s
b-) noktasının s-düzlemindeki yerini bulunuz. T=0.1 s5.0j5.0z
2j1s
632.0zz
264.0z368.0)z(G,
1ss
1)s(G
22
Örnek: verilen sistemlerin ζ, wn ve τ değerlerini karşılaştırınız.
Sayısal sistem, analog sistemin T=1 saniye ile örneklenmiş
karşılığıdır. Bu karşılığın bulunuşu ileride açıklanacaktır.
22
2
221
1)(
nn
n
wsws
w
sssG
2
,5.0
,1wn
0632.0zz2
Rad89.0795.0z618.0j5.0z 2,12,1
Çözüm:
Ayrık zaman sistemin karakteristik denklemi,
z düzleminden s düzlemine dönüşümler kullanılarak
25.0
89.0)795.0(
)795.0(
22
In
In
9191.089.0)795.0(1
1 22 InWn
snIn
36.4)795.0(
1
NOT: Örnekleme peryodu daha küçük seçilerek ayrık zaman
karşılığı bulunsaydı ayrık zaman sistemin kriterleri, analog
sistemin kriterlerine oldukça yakın olacaktı.
Şekilde s-düzlemindeki çeşitli noktaların z-düzlemindeki
karşılıkları verilmiştir. Buradan, s-düzlemindeki sanal eksenin
yani kritik kararlılık sınırının, z-düzleminde birim çembere
karşılık geldiği görülmektedir. Dolayısıyla z düzleminde birim
çemberin içi kararlılık bölgesidir. S-düzleminde kutuplar sola
doğru uzaklaştıkça bağıl kararlılık artarken z-düzleminde