Processing math: 100%

dersblog



Ders 22

Mühendislerin Laplace Transformunu sevmesinin sebeplerinden biri fonksiyonlarda kesintili / birdenbire zıplayan geçişleri (jump discontinuities) durumlarında bozulmadan işleyebiliyor olmasıdır.

Kesintili fonksiyonlardan biri birim adım fonksiyonudur. Fonksiyonun kendisi tartışma yaratan bir fonksiyon aslında, sıfır noktasında hangi değere sahip olduğu hala kararlaştırılamadı. Bazıları 0 diyor, bazıları 1 diyor, ben bu derste tanımsız bırakacağım. Yani

u(t) birim adım

u(0) tanımsız

Bazen sıfır yerine başka bir noktada zıplama olmasını isteyebilirim. O zaman fonksiyonu kaydırmak mümkün, mesela a kadar.

ua(t)=u(ta)

Birim Kutu

Kutuyu birim adımlar kullanarak temsil edebilirsek iyi olur, çünkü birim adımların Laplace transformunu biliyoruz.

uab:=ua(t)ub(t)

=u(ta)u(tb)

Mantıklı gözüküyor, b'ye gelene kadar birim adım gibi gidiyoruz, b'de birim adımı iptal edecek şekilde birim adımı çıkartıyoruz.

Peki bu fonksiyonlar niye faydalı? Diğer fonksiyonları çarptıklarında o fonksiyonları faydalı şekillerde değiştirebildikleri için. Diyelim ki şöyle bir f(t) var.

uabf(t)

neye benzerdi? Eğer a,b arasında uab 1 değerine eşitse, çarpım bu aralıkta f(t)'yi olduğu gibi alır. Diğer noktalarda sıfır yapar. Alttaki şekil ortaya çıkar (mavi çizgiler).

f(t)'nin bir parçasını çekip çıkartıyoruz yani. Bu çok ise yarayan bir numara / teknik.

Birim adımın Laplace Transformunu alalım.

L(u(t))=0estu(t)dt

Entegral alt sınırı 0'dan başlıyor. Birim adım 0'dan büyükse hep 1, o zaman üstte u(t)'yi 1 alabiliriz. Yani üstteki aslında 1'in Laplace Transformudur.

=1s   s>0(1)

Kaydırılmış birim adımın Laplace'i

L(u(ta))=a1estdt

Alt sınır a oldu çünkü a'dan önce sıfır, sonrası 1. Bizi tek ilgilendiren kısım fonksiyon 1 olduktan sonra, o zaman alt sınırı değiştiririz.

=ests|a=eass

Diğer bir özellik

L(ectf(t))=F(sc)

İspat: Eğer entegral temsilini açarsak

=0ectf(t)estdt=e(cs)tf(t)dt

Bu son ifadeye başka bir yönden erişebilir miyiz? Mesela F(sc) nedir? s yerine sc koyarsak (temel formülden türetelim)

F(sc)=0f(t)e(sc)tft

=0e(s+c)tdt=0e(cs)tdt

Aynı ifadeye eriştik.

(1) formülünde de ters yönde gidersek

L1(1s)=?

Geriye giderken 1'e mi, birim adıma mı döneceğim? Şimdiye kadar geriye giderken hep 1 seçtim. Artık bu yeterli değil. Problem nerede? Çünkü Laplace Transformu entegral alt sınırı 0'dan başlıyor, yani transform 0 öncesine bakmıyor bile. Bir fonksiyon f(t) 0'dan önce envai türde değere sahip olabilir, Laplace Transformu aynı olacaktır.

Eğer hangi tür problem ile uğraştığımızı biliyorsak bu durum bir sorun çıkarmaz. Mühendisler şöyle kategorize ederler; eğer bir problem şimdiki zamandan, t=0, başlayıp ileri gidiyorsa o bir Laplace Transform problemidir, eğer geçmiş zaman verisi gerekiyorsa o bir Fourier Transform problemidir.

Eğer 1 noktasından ( t=0 ) geriye giderken ortaya çıkan tanımsızlıktan (ambiguity) kurtulmak istiyorsak, ters Laplace'dan elde edilen fonksiyonu birim adım ile çarparız, böylece t=0 öncesi sıfır olmasını garanti etmiş oluruz.

f(t)F(s)

L1(F(s))=u(t)f(t)

Böylece ters Laplace özgün (unique) hale gelir.

Taşıma (translation) ve Laplace

L(f(ta))'yi L(f(t)) temel alarak yazabilir miyiz? Bu ne yazik ki mumkun degil.

Üstte f(t)'yi a kadar sağa kaydırdık. Kesik çizgili yeni, taşınmış fonksiyonu niye düz çizgili fonksiyon bağlamında temsil edemiyorum? Sebep bariz, çünkü f(t)'nin sıfır öncesi değeri. O değer Laplace Transformu sırasında kullanılmaz, ama taşındıktan sonra yeni fonksiyonda kullanılır. Aradaki bu uyuşmazlığı, eksikliği kapatmak mümkün değildir.

Yapabileceğimiz en iyi şey, kaybettiğimiz bölümü u(ta)f(ta) kullanımıyla Laplace Transformundan silmektir. Altta koyu çizgili olan kısmın transformunu alırız yani.

u(ta)f(ta)easF(s)(2)

(2) geçişi Laplace Tablosunda görülebilir.

Fakat f çoğunlukla ta formunda gelmez, mesela sin(t), cos(t) olur, o zaman şu teknik kullanılabilir,

u(ta)f(t)easL(f(t+a))

Bu tekniğe bazıları üstel kaydırma fonksiyonu (exponential shift formula) diyor. Biz bu terimi kullanmayacağız. Biz t ekseni taşınması (t-axis translation) formülü diyeceğiz.

(2)'nin ispatı

0estu(ta)f(ta)dt

Ama formülde F(s) kullanmak istiyorum, o zaman f(ta)'yi bir şekilde f(t) yapmam lazım. Değişken değişimi yaparım.

t1=ta

=aes(t1+a)u(t1)f(t1)dt1

=esaaest1u(t1)f(t1)dt1

u(t1) t1<0 için sıfır olduğuna göre a,0 arasını entegral alt sınırından atabiliriz,

=esa0est1f(t1)dt1

u(t1)'in kendisini niye attık? Çünkü sıfır sonrası 1 değerinde zaten.

"Ama t değişkeni t1 haline geldi, artık bu Laplace Transformu olmaz" diye düşünenler olabilir. Fakat entegrasyonda kullanılan fonksiyonun illa t ismini taşıması gerekmiyor, çünkü o değişken bir yer tutucu (dummy) sadece. Entegre edilen, sınırlar düzgünse, üstteki de bir Laplace Transformudur.

Yani

u(ta)f(ta)easL(f(t))

t yerine t+a geçirdim.

Eğer u(ta)f(t) olsaydı?

u(ta)f(ta+a)easL(f(t+a)))

Örnek

uab=u(ta)u(tb)

eassebss

Örnek

u(t1)t2

Transformu nedir?

esL((t+1)2)=esL(t2+2t+1)=es(2s3+2s2+1s)

Bayağı kalabalık bir sonuç oldu, fakat transformunu aldığımızda elde edilen yukarıdaki fonksiyon da o kadar basit bir fonksiyon değil zaten; altta yeşil renkle işaretli

Fonksiyon kesintili (discontinuous) bir fonksiyon. Kesintinin nerede oluştuğu eas formülünde a, yani es=e1s için 1.

L1(1+eπss2+1)

Cevabın bir kesikli fonksiyon olacağını şimdiden biliyoruz, çünkü formül içinde üstel bir fonksiyon var. İlk yapacağımız iş üstel fonksiyonu yanlız kalacak şekilde ifadeyi parçalamak.

=1s2+1+eπss2+1

Önce suna bakalım

1s2+11  ?

Eğer üsttekine cevap olarak sin(t) dersek, yanlış olur. Eğer üstel fonksiyon formülün geri kalanında olmasaydı, o zaman sin(t) olabilirdi, ama şimdi

1s2+11u(t)sin(t)

kullanmamız lazım. Peki

eπss2+11  ?

Yine (2)'yi kullanırım,

1u(tπ)sin(tπ)

Nihai cevap nedir? İki parçayı toplarım

1s2+1+eπss2+11u(t)sin(t)+u(tπ)sin(tπ)

Bu cevap teknik olarak doğru, ama biraz masajlayıp onu daha düzgün bir hale getirmek lazım [hatta hoca şakayla karışık bu halde bırakırsanız sınavda not kırarım diyor].

O zaman her parçaya teker teker bakalım. Formülün u(t)sin(t) kısmı sadece t0 için anlamlı, çünkü ondan önce birim step fonksiyonu sıfır. u(tπ)sin(tπ) kısmi ise tπ için anlamlı. Elimizde iki tane senaryo var aslında, o zaman bu fonksiyonu parçalı, senaryolu olarak gösterebiliriz.

f(t)={sin(t)0tπsin(t)+(sin(t))tπ

Terim sin(t), sin(tπ)'dan geldi, çünkü sin'i π kadar yana kaydırırsak (altta kesik çizgili)

elde edilen sin(t)'dir. Nihai formül

f(t)={sin(t)0tπ0tπ


Yukarı