DOĞRUSAL ZAMANLA DEĞİŞMEZ SİSTEMLERDE
DİFERANSİYEL DENKLEMLER
0 aN
Giriş – çıkış ilişkisi,
) ( )
( )
( )
( 1 ( 1) 1 0
)
( t a y t a y t a y t
y
aN N N N
) ( )
( )
( )
( 1 ( 1) 1 0
)
( t b x t b x t b x t
x
bM M M M
M
N b b b
a a
a0, 1,, ; 0, 1,, katsayılarının hepsi zamana göre sabitse, sistem N. mertebeden doğrusal zamanla değişmez (DZD) bir sistemdir.
olmak üzere
) ( ,
), ( ),
(t0 y t0 y( 1) t0
y N
Çözümdeki N adet sabit,
başlangıç şartlarına bağlıdır.
f (t)) ( )
( )
( )
( )
(t b x( ) t b 1x( 1) t b1x t b0x t
f M M M M
Çözüm Adımları
diyelim.
1) f (t) yerine sıfır, y yerine yh yazılarak homojen çözüm bulunur.
yh(t) , henüz belirlenmemiş N adet sabite bağlı olarak yazılır.
2) Denklemdeki f (t) dikkate alınarak, fakat hiçbir homojen çözüm bileşeni içermeyen özel çözüm yö(t) bulunur. Bu çözümdeki bütün sabitler bu aşamada belirlenir.
3) y(t) yh(t) yö(t) N adet sabit, başlangıç şartlarına göre belirlenir.
Eğer f (t) parçalı tanımlıysa, her tanım aralığı için:
yh(t) sabitler için farklı semboller kullanılarak aynı biçimli olarak yazılır.
f (t)’nin yalnızca değişen bileşenlerine karşılık gelen özel çözüm bileşenleri yeniden bulunur.
Verilen şartlar sadece geçerli olduğu zaman aralığında kullanılarak bulunan çözümden, komşu aralığın başlangıç (veya son) değerlerine geçiş yapılır ve bunlar o aralıkta kullanılır.
f (t) geçiş anında etkiyen bir darbe içermiyorsa
sıçrama yapmaz.
geçiş anında y(t), y(t),, y(N1)(t) Homojen çözüm
0 )
( )
( )
( )
( 1 ( 1) 1 0
)
( t a y t a y t a y t
y
aN hN N hN h h
Çözümün,
e
λt biçiminde bileşeni olduğunu varsayalım.Bu çözümü homojen denklemde yerine yazalım.
0 0
1 1
1
t t t
t e e e
e a a a
aNN N N
0 0
1 1
1
a a a
aNN N N
Buradan da
diferansiyel denklemin ya da sistemin karakteristik denklemi bulunur.
N
1, 2,,
0 0
1 1
1
a a a
aNN N N
Denklemin
diferansiyel denklemin ya da sistemin özdeğeri (eigenvalue), köklerinin her birine ise
ya da karakteristik kökü denir.
Her özdeğer için çeşitli katsayılarla birer homojen çözüm bileşeni bulunması mümkündür. Buna göre:
Çakışık özdeğer yoksa (tüm özdeğerler birbirinden farklıysa)
t
t N
e
e
Nh t A A
y ( ) 1 1
çakışık özdeğer varsa, meselâ λk , m-katlı bir özdeğer
(λk = λk+1 = … = λk+m-1) ise bunlara karşılık gelen homojen çözüm kısmı
t t
t k k
k
e e
e
k k m mk A t A t
A 1 1 1
Ak Ak1t Akm1tm1 e
ktİstenirse bu bileşenlerdeki t’lerin hepsinin yerine meselâ (t-2) gibi kaymalı bir ifade de yazılabilir.
Böyle kaymalı ifadeler bazen hesap ve yorum kolaylığı sağlayabilir.
Sabitler buna göre bulunur.
Örnek:
homojen denkleminin çözüm ifadesini bulunuz.
0 )
( 128 )
( 160
) ( 72
) ( 14
)
( (5) (4) (3)
) 6
( t y t y t y t y t
yh h h h h
Çözüm: 6 145 724 1603 1282 0 1 2 0 , ,
3 2
4 5 6 4
t t
t t
t
t
e e e e e
e
A t A A A t A tA t
yh( ) 1 0 2 0 3 2 4 4 5 4 6 2 4
A1 A2t
A3e
2t
A4 A5t A6t2 e
4t
Not:
Eğer λk,k+1 = σ ± jω gibi eşlenik çiftler halinde karmaşık özdeğerler varsa, bunlara karşılık gelen homojen çözüm bileşenleri,
çakışık özdeğer değilseler
t C
t
B
e
t sin e
t cos biçiminde de yazılabilir.Ayrıca bu özdeğerlerin çakışması da varsa yine bunun da t ’nin uygun kuvvetleriyle çarpılmış biçimleri de gelir.
Örnek:
homojen denkleminin çözüm ifadesini bulunuz.
Çözüm:
0 )
( 5 )
( 7
) ( 3 )
)(
3
( t y t y t y t
yh h h h
0 5
7 3 2
3
1 1 , 2,3 1 j2
Çakışık kök olmadığı için:
t A
t A
A t
yh( ) 1
e
t 2e
t cos2 3e
t sin 2Örnek: yh(4)(t) 18yh(t) 81yh(t) 0
homojen denkleminin çözüm ifadesini bulunuz.
Çözüm: 4 182 81 0 1 2 j3 , 3 4 j3 Eşlenik köklerin her biri 2 katlı olduğu için:
t t
A t
t A t
A t
A t
yh( ) 1 cos3 2 sin3 3 cos3 4 sin3
b t b t
t a
a ) cos3 sin3
( 1 2 1 2
Özel çözüm
Diferansiyel denklemin sağ tarafı f (t) bileşenlerine ayrılır.
Her bileşen için ayrı ayrı özel çözüm bileşenleri bulunur
ve hepsinin toplamı özel çözüm olur.
Özel çözüm, homojen çözüm tarafından kapsanan herhangi bir bileşen içermemelidir.
t
e
pr1 1 gibi bir bileşen için özel çözüm bileşeni:
} ,
, { 1
1 N
p ise yö1(t) c1
e
p1tp1 k gibi m-katlı bir özdeğere eşitse
(bu özdeğer çakışık değilse m = 1 demektir)
kt
e
m
ö t c t
y 1( ) 1
f (t) içindeki sabit terimler için p1 = 0 olduğu unutulmamalıdır.
Özel çözüm bileşenlerinin katsayıları, başlangıç şartları kullanılmadan belirlenmelidir.
yerine yö1 , ve f (t) yerine yalnızca ilgili bileşeni
Bunun için, c1 katsayısı, diferansiyel denklemde y
t
e
pr1 1 yazılarak bulunur.
} ,
, { 1
1 N
p ise bu yol,
0 1
1 1
1 1 1
1
1 a p a p a p a
c N r
N N
N
biçiminde karakteristik denklemde λ yerine p1 kullanılan kısa bir
hale gelir. p1’in bir özdeğere eşit olması halinde bu kısa yol geçersizdir.
Tüm bileşenler için özel çözüm bileşenleri bulunduktan sonra
( ) ( )
)
(t y 1 t y 2 t
yö ö ö
çözüm bileşenlerinde de istenirse t yerine meselâ (t-2) gibi kaymalı bir ifade de yazılabilir; bu sadece katsayıların farklı olmasını gerektirir, ki henüz belirlenmemiş katsayıların buna göre belirlenmesi sorun teşkil etmez. Böyle kaymalı ifadeler bazen hesap ve yorum kolaylığı
sağlayabilir.
Örnek: y(t) 5y(t) 6y(t) x(t) x(t) 4
e
tu(t) ,0 )
0 (
y y(0) 0
sisteminin çıkışını,
girişi ve başlangıç şartları için bulunuz.
Çözüm: 2 5 6 0 1 2 , 2 3
0
t için f (t) tx( ) 0 olduğundan
t
t
e
e
AA t
y t
y
t 0 ( ) h( ) 1 2 2 3
Denklemin sağ tarafında darbe olmadığı için
2
0 1
) 0 ( )
0
( y A A
y
2
1 3
2 0
) 0 ( )
0
( y A A
y A1 A2 0 Yani t 0 y(t) 0 .
t
t
e
e
BB t
y
t 0 h( ) 1 2 2 3
} 3 ,
2 {
1
t
e
t x t
f ( ) ( ) 4 ve yö(t) c
e
t6 2 )
1 ( 5 )
1 (
4
2
c yö(t) 2
e
tt t
t
e e
e
BB t
y
t 0 ( ) 1 2 2 3 2 2
0 )
0
( B1 B2 y
2 3
2 0
) 0
( B1 B2
y B1 4 , B2 2
Sonuçta tüm zamanların çözümü:
2 2 0 iseise 0
) 0
( 2 3
4
tt t
y t t t
e e
e
ya da kısaca y(t)
4 e
2t 2e
3t 2e
t
u(t)Not:
Eğer t < t0 için f (t) = 0 ise yani f (t) = (…)∙u(t-t0) biçiminde yazılabiliyor ve t = t0’daki
) ( ,
), ( ),
(t0 y t0 y( 1) t0
y N
başlangıç değerlerinin hepsi sıfır ise t < t0 için y(t) = 0 olacağı zaten bellidir. Bu yüzden böyle bir durumda yalnız t ≥ t0 için u(t-t0) = 1 koyarak bulunan çözümü u(t-t0) ile çarparak tüm zamanların çözümü elde edilir.
Örnek: y(t) y(t) x(t)
) 2 (
5 )
2 (
4 )
(t ( 2)u t u t
x
e
t, 0 )
2 (
y y(2) 0 sisteminin çıkışını
için bulunuz.
Çözüm: 2 1 0 1,2 j
2
t yh(t) Acos(t 2) Bsin(t 2)
2
t f (t) x(t) 4
e
t( 2) 5 }, {
1 1 2
0 {1, 2}
için homojen çözüm:
için denklemin sağ tarafı
olup hem hem de
) 2
4
e
t( yö1(t) c1e
(t2) 2 1) 1 (
4
1 2
c
) 2 (
5 0
5
e
t yö2 c2e
0(t2) c2 5 10 5
2 2
c için
için 5 2
)
(t
e
t( 2) yö5 2
) 2 sin(
) 2 cos(
)
(t A t B t
e
t( 2) y0 5
2 )
2
( A y
0 2
) 2
( B
y A 7 , B 2 Tüm zamanlar için ise:
7cos( 2) 2sin( 2) 2 5
( 2))
(t t t ( 2) u t
y
e
tÖrnek: sisteminin çıkışını
için bulunuz.
) ( 3 )
( 2 )
( 3 )
(t y t y t x t
y
, ) 3 (
) ( 2 )
(t u t 2( 3)u t
x
e
t y(0) 0 , y(0) 1Çözüm: 2 3 2 0 1 1, 2 2. 0
) (
0
x t
t y(t) yh(t) A1
e
t A2e
2t
Denklemin sağ tarafında darbe olmadığı için
2
0 1
) 0 ( )
0
( y A A
y
2
1 2
1 )
0 ( )
0
( y A A
y A1 1, A2 1
t
t
e
e
t y
t 0 ( ) 2
t
t
e
e
BB t
y
t 3 h ( ) 1 2 2
0 f (t) 3 2 6 6
e
0t} ,
{
0 1 2 yö1(t) c1 yö1 0 , yö1 0
Bunlar diferansiyel denklemde yalnız ilgili bileşen için yerine konursa:
6 2
0 3
0 c1 c1 3 yö1(t) y(t) B1
e
t B2e
2t 3 30 )
0
( B1 B2 y
2
1 2
1 )
0
( B B
y B1 5 , B2 2
) 3 ( 2 )
3 (
2
) 1
(
3
y t K
e
t Ke
tt h
Sağ taraftaki 6 terimi değişmediği için yö1 = 3 de aynıdır.
) 3 (
3 2
e
t için p1 2 2 yö2(t) c2(t 3)e
2(t3)) 3 ( 2 )
3 (
2 2 ( 3)
)
( 2 2
2
c
e
t c te
tt yö
2(t) 4c2
e
2(t3) yö
4 3
2( 3) 3 2( 3)) ( 2
) ( 3
)
( 2 2 2 2
2
y t y t c c
e
te
tt
yö ö ö
c23 yö2(t) 3(t 3)
e
2(t3)) ( )
( )
( )
(t y t y 1 t y 2 t
y h ö ö
) 3 ( 2 )
3 ( 2 )
3
( 3 3( 3)
)
(t K1
e
t K2e
t t e
ty
K1, K2 sabitlerini bulmak için bu zaman bölgesinde iki şarta ihtiyaç vardır. Diferansiyel denklemin sağ tarafında geçiş (t = 3) anında
etkiyen bir darbe olmadığı için bir önceki zaman aralığı çözümünün son değerlerini bu bölgenin başlangıç şartları olarak kullanabiliriz:
3 )
3 ( )
3
( y B1
e
3 B2e
23 y6
3 4
5 )
3 ( )
3
( y
e
e
y
3 2
5 3 6
e
e
K1 K2 3 32 2
1
K K
3 2 6
2
e
K 3
5 3
1
e
KGenel çözümü yazarsak:
Not:
f (t) trigonometrik veya üstel çarpanlı trigonometrik bir bileşen içeriyorsa bunun üstel sinyaller biçimindeki ifadesinin üs katsayıları
j p1,2
3 )
3 (
3 3 3
0 3
2
0 )
(
) 3 ( 2 )
3 ( 6 2
) 3 3 (
2 2
) 2
3 ( )
5 3
( 5
t t
t t
t y
t t
t t t
t t
e e
e e
e e e
e e
gibi eşlenik çiftlere karşılık gelir.
t r
t
r1
e
t cos 2e
t sin }, ,
{ 1 N
j
c1
e
t cost c2e
t sin tDolayısıyla, f (t) içindeki için özel çözüm
ise
Eğer α ± jβ , m-katlı bir özdeğer çiftine eşitse, yukarıdaki ifadenin tm ile çarpılmışı yazılır. α = 0 için bu, zorlanmış rezonans olduğu anlamına gelir.
Örnek: y(t) 3y(t) x(t) , x(t) 6 4
e
3t cos2t , y(0) 2 Çözüm: 3 0 3 yh(t) Ae
3t6 bileşeni için, 0 yö1 c1 2 ( ) 3
0 6
1
1 y t
c ö
t t
e
cos24 3 için 3 j2 yö2 c2
e
3t cos2t c3e
3t sin 2t )( 3
)
( 2
2 t y t
yö ö
3c2 2c3 3c2 e
3t cos2t
2c2 3c3 3c3 e
3t sin 2tt t
e
cos24 3
, 4 2 3
c c2 2 0 c2 0 , c3 2
Genel çözüm: y(t) A
e
3t 2 2e
3t sin 2t
2 2 0
) 0
( A A
y y(t) 2
e
3t 2e
3t sin 2t Örnek (Zorlanmış rezonans):) ( 5 )
( 4 )
(t y t x t
y x(t) u(t)sin 2t
, 0 )
0 (
y y(0) 0
sisteminin çıkışını, girişi ve başlangıç şartları için hesaplayınız.
Çözüm:
0
2 4
1,2 j2
0
t için yh(t) Acos2t Bsin 2t Sağdaki 5sin 2t için, p1,2 j 2 1,2
t t
c t
t c t
yö( ) 1 cos2 2 sin 2
t t
c t
c t
t c t
c t
yö( ) 1 cos2 2 1 sin 2 2 sin 2 2 2 cos2
t c c t c c t t yö( ) 2 2 2 2 cos2 2 1 2 1 sin 2
tt t c
t c
t y t
yö( ) 4 ö( ) 4 2 cos2 4 1 sin 2 5sin 2
(Son türevde ve özel çözümün yerine yazıldığı diferansiyel denklemde çarpanlı kısmın katsayı bulunmasına faydası olmadığı için … diyerek geçilmiştir.)
0
4c2 c4 1 5 c1 5 4 c2 0 t
t t
B t
A t
y cos2
4 2 5
sin 2
cos )
(
ve ve
Toplam çözüm ise:
A y(0) 0
4 2 5
0 )
0
( B
y 8
, 5
0
A B
) ( 2
4 cos 2 5
8 sin ) 5
(t t t t u t
y
( t ≥ 0 çözümünü u(t) ile çarpmak için gereken şartlar sağlandığı için)
bulunur. Buradaki t çarpanlı bileşen, salınım genliğinin gittikçe artmasına neden olarak pratikte sistemin modelinin geçerli olduğu sınırları zorlayarak bozulmaya neden olur. Köprü ve cam eşyalarda yıkıma neden olduğu söylenen rezonans, normal rezonans değil, bu örnekteki gibi zorlanmış rezonanstır. Salıncak, sarkaç gibi aslında doğrusal olmayıp, küçük salınım sınırlarında yaklaşık doğrusal olan sistemlerde de, salınım genliğinin bir yere kadar kolayca artırılabilmesi için, sallama kuvvetinin, sistemin doğal salınım frekansında olması
gerekir. Meselâ salıncağı bir taraftan itekledikten doğal salınım periyodu kadar sonra salıncak aynı tarafa gelecek, o anda tekrar iteklenince salınım genliği kolayca artacaktır. İki taraftan iteklenen salıncakta da bir taraftan uygulanan kuvvet artı, diğer taraftan
uygulanan kuvvet eksi olduğu için kuvvetin periyodu yine aynıdır.
Ancak sistemi doğrusal varsayılabilecek sınırların dışına çıkartan
salınım genliklerine ulaşınca bu genlik sınırlı hale gelecektir. Salıncağı, doğal salınım frekansından başka bir frekansta kuvvetle sallamak
denenirse, genlik artışının ne kadar zor olduğu görülür.
Diferansiyel Denklemin Sağ Tarafında Darbe Varsa Çözüm
) (
) ( )
( )
( )
( )
( )
( 1 ( 1) 1 0 0
)
( t a y t a y t a y t f t g t r t t
y
aN N N N
Burada aN 0 ve g(t) , t0 anında etkiyen bir darbe içermeyen sonlu bir fonksiyon olsun. r sabit olmasa bile t0 ’daki değeri önemli olduğu için sabit olarak ele alınır.
t < t0 ve t > t0 bölgeleri için ayrı ayrı çözüm yapılır. Tabii bu bölgelerde olduğundan bu darbe dikkate alınmaz.
0 )
(t t0
t0
t y(t0) , y(t0) ,, y(N1)(t0 )
Ancak bölgesine ait son değerler ile
t0
t bölgesine ait y(t0) , y(t0) ,, y(N1)(t0 ) başlangıç değerleri arasındaki geçiş şöyle yapılır:
N N
N
a t r
y t
y( 1)( 0) ( 1)( 0 )
) ( )
(
) ( )
(
0 0
0 ) 2 ( 0
) 2 (
t y t
y
t y
t
y N N
N 2
N = 1 ise yalnızca en üstteki uygulanır.
Yani
aN
t r y t
y( 0 ) ( 0)
Diferansiyel denklemin
t0,t0 aralığında integralini alalım:İspat:
0
0 0
0 0
0 0
0
) ( )
( )
( )
( 1 ( 1) 1 0
) (
t
t t
t t
t
N N
t
t
N
N y t dt a y t dt a y t dt a y t dt
a
0
0 0
0
) (
)
( 0
t
t t
t
dt t
t r
dt t
g
Sol tarafta ilki hariç tüm integraller, sağ taraftaki ilk integral, sonlu büyüklüklerin sonsuz küçük zaman aralığında integralleri olduğu için sıfırdır. Bu zaman aralığında yalnızca sağdaki darbe ve onun doğrudan etkilediği y(N)(t) sonsuz olup bunların integralleri eşitlenerek:
y t
a
y t y t
raN (N1)( ) tt0 N (N1)( 0) (N1)( 0)
0
N N
N
a t r
y t
y( 1)( 0) ( 1)( 0) Varsa (N ≥ 2 ise) daha düşük mertebeli türevlerde ise sıçrama olmaz.
Sağ tarafta darbenin de türevleri olsaydı daha düşük mertebeli türevlerde de sıçrama yaptırabilirdi;
ancak bu konuyu burada ele almayacağız. Sıradaki konu bu konunun uygulaması niteliğindedir.
dolayısıyla
Sistemin Diferansiyel Denkleminden Birim Darbe Tepkisinin Bulunması
Normalde x yerine δ , ve y yerine h yazarak diferansiyel denklemi çözüp birim darbe tepkisini bulabiliriz.
Ancak az önceki konuda anlatılan kuralı, aN 0 olmak üzere, )
( )
( )
( )
( )
( 1 ( 1) 1 0 0
)
( t a y t a y t a y t b x t
y
aN N N N
ile verilen nedensel bir sisteme uygularsak birim darbe tepkisini
bulma yöntemi oldukça basitleşir. Çünkü nedensel DZD sistemlerde 0
0 )
(t t h
olduğu için t < 0 bölgesinde h(t) ’nin bütün türevleri de sıfır olacaktır.
Sonsuz küçük zaman farkını göz ardı ederek başlangıç anını 0+ yerine 0 alırsak yöntem şu şekilde özetlenebilir:
Böylece t > 0 bölgesine ait başlangıç şartları (0+ anı için) kolayca belirlenebilecektir.
y yerine h , sağ tarafa da sıfır yazarak t > 0 bölgesi için 0 )
( )
( )
( )
( 1 ( 1) 1 0
)
( t a h t a h t a h t
h
aN N N N
N N
a h( 1)(0) b0
0 )
0 (
0 )
0
)(
2 (
h h N
diferansiyel denklemi şu başlangıç şartları için çözülür:
Bu kısım N ≥ 2 için uygulanır.
N = 1 ise yalnızca en üstteki, yani
aN
h(0) b0 uygulanır.
Bulunan çözüm u(t) ile çarpılarak birim darbe tepkisinin tüm zamanlar için ifadesi elde edilir.
Denklem ve dolayısıyla çözümü homojendir.
Örnek:
Şekildeki devrenin birim
darbe tepkisini bulunuz. x(t) R L
) (t y
Çözüm:
) ( )
( )
(t y t x t R y
L
Direncin akımı + endüktansın akımı = x(t) 0
) ( )
(
0
h t h t
R
t L
L R R
h L1 )
0 denklemini (
başlangıç şartı için çözmeliyiz.
L R R
L 1 0
Lt
e
RA t
h
( )
L A R
h
(0)
) ( )
( u t
L t R
h Lt
e
R Tüm zamanlar için ifadesi ise:
) (t h
t L
R
Örnek:
) ( 2 )
( 6 )
( 11 )
( 6 )
(t y t y t y t x t
y
ile tanımlı nedensel sistemin
birim darbe tepkisini bulunuz.
Çözüm:
0 )
( 6 )
( 11 )
( 6 )
(
0
h t h t h t h t
t 3 62 11 6 0
3 ,
2 ,
1 2 3
1
t t
t
e e
e
A AA t
h( ) 1 2 2 3 3
3 2
0 1
) 0
( A A A
h
3 2
1 2 3
0 )
0
( A A A
h
3 2
1 4 9
2 1
2 )
0
( A A A
h
1 ,
2 ,
1 2 3
1
A A A
2
( ))
(t 2 3 u t
h
e
t e
t e
tSonuç:
Bazen diferansiyel denklemin yalnız sağ tarafı, Not:
) (
) ( )
( )
( )
( 1 ( 1) 1 0 0
)
( t a y t a y t a y t bx t t
y
aN N N N
sonra h(t) ifadesi u(t) yerine u(t-t0) ile çarpılarak tüm zamanlar için geçerli ifade yazılır.
gibi zamanda ötelenmiş biçimli olabilir. Bu durumda daha önce anlatılan başlangıç şartları 0 yerine t0 anı için yazılır ve katsayılar bulunduktan
Bu çözümü (t-t0) ’a göre düzenlemek, hem katsayıların kolay bulunmasını, hem de çizim ve yorum kolaylığı sağlar.
t > t0 bölgesi için aNh(N)(t) aN1h(N1)(t) a1h(t) a0h(t) 0 diferansiyel denklemi şu başlangıç şartları için çözülür:
N N
a t b
h( 1)( 0)
0 )
(
0 )
(
0 0 ) 2 (
t h
t h N
Bu kısım N ≥ 2 için uygulanır.
Sonra da bu çözüm u(t-t0) ile
çarpılarak tüm zamanlar için geçerli ifade elde edilir.
Örnek: 2y(t) 50y(t) 4x(t 3)
ile verilen nedensel sistemin birim darbe tepkisini bulunuz.
0 )
( 50 )
( 2
3
h t h t
t 22 50 0
2 5
,
1 j
h(t) Acos5(t 3) Bsin5(t 3) A
h(3) 0
B h(3) 4 2 2 5
5 2 ,
0
A B
5( 3)
( 3)5 sin ) 2
(t t u t h