Diferansiyel Denklemler: Çözüm Yöntemleri ve Örnekler

Bir denklemin mertebesi, derecesi ve doğrusal olup olmadığı nasıl okunur?

Genel bir diferansiyel denklem $F(x,y,y',y'',\ldots,y^{(n)})=0$ biçiminde yazılabilir. Burada $y$, bilinmeyen fonksiyon; $y'$, $y''$ ve daha yüksek türevler ise bu fonksiyonun değişim oranlarıdır.

Mertebe, denklemde görülen en yüksek mertebeden türevin mertebesidir. Örneğin $y''+3y'-2y=0$ ikinci mertebedendir. Türev derecesi ise denklem türevler bakımından polinom biçimindeyse, en yüksek mertebeden türevin üzerindeki üssü ifade eder. Mesela $(y'')^2+y=0$ ikinci mertebeden ve ikinci derecedendir. Ancak $\sin(y'')+y=0$ için türevin polinom içindeki derecesinden söz etmek uygun değildir; bu nedenle mertebe ile dereceyi aynı kavram gibi kullanmamak gerekir.

Bir denklem, bilinmeyen fonksiyon ve türevleri birinci dereceden olmalı, birbirleriyle veya kendi aralarında çarpılmamalı ve bunların katsayıları yalnızca bağımsız değişkene bağlı olmalıdır. $y'+P(x)y=Q(x)$ doğrusal birinci mertebe denklemdir. Buna karşılık $y^2+y'=0$, $yy'=x$ veya $\sin y+y'=0$ doğrusal değildir. Doğrusallık, çözüm yöntemini doğrudan etkiler: homojen doğrusal denklemlerin çözümleri için süperpozisyon ilkesi geçerlidir. Homojen olmayan doğrusal denklemlerde ise genel çözüm, bir özel çözüm ile karşılık gelen homojen denklemin genel çözümünün toplamıdır.

Adi diferansiyel denklemde bilinmeyen fonksiyon tek bir bağımsız değişkene bağlıdır: $y=y(x)$ veya $y=y(t)$. Kısmi diferansiyel denklemde ise $u(x,t)$ gibi birden fazla bağımsız değişken vardır ve $\partial u/\partial t$ türünde kısmi türevler kullanılır. $\partial u/\partial t=k\partial^2u/\partial x^2$ ısı yayılımını ifade eden bir kısmi diferansiyel denklem örneğidir. Bu ayrım, hangi çözüm tekniğinin seçileceğini belirleyen ilk kontroldür.

Genel çözümden başlangıç koşullu çözüme geçiş

Denklem tek başına çoğu zaman bir çözüm ailesi tanımlar. Örneğin $dy/dx=3x^2$ eşitliğinin integrali $y=x^3+C$ olur. Buradaki $C$, aynı diferansiyel denklemi sağlayan farklı eğrileri birbirinden ayıran keyfi sabittir. Bu ifade genel çözümdür.

Başlangıç koşulu, fonksiyonun belirli bir noktadaki değerini verir. $y(0)=5$ koşulu, genel çözümde $x=0$ yazılarak $5=C$ sonucunu doğurur ve özel çözüm $y=x^3+5$ olur. İkinci mertebeden bir denklemde çoğunlukla iki bağımsız koşul gerekir; örneğin $y(0)$ ve $y'(0)$ değerleri. Ancak bu, her denklem için koşulsuz bir kural değildir: koşulların denklemin yapısıyla uyumlu olması ve çözümün incelenen aralıkta anlamlı kalması gerekir.

Adi diferansiyel denklemlerde başlangıç koşulları genellikle aynı bağımsız değişken değerinde fonksiyon ve türev değerlerini belirtir. Sınır koşulları ise bir aralığın uçlarında veya bir bölgenin sınırında fonksiyon ya da türevleri kısıtlar. Kısmi diferansiyel denklemlerde başlangıç verisi çoğunlukla bir başlangıç yüzeyi veya zamanındaki uzaysal fonksiyon olarak verilir. Başlangıç koşulu zamanla ilerleyen bir sistemi tanımlarken sınır koşulu, bir aralık veya bölge üzerindeki uç davranışını belirler.

Bir çözümün yalnızca cebirsel olarak yazılması yeterli değildir. Çözümün tanımlı olduğu aralık, paydaların sıfır olup olmadığı ve başlangıç koşulunun sağlanıp sağlanmadığı kontrol edilmelidir. Bir fonksiyonu elde ettikten sonra türevini alıp başlangıç denklemine yerleştirmek, işaret ve sabit hatalarını yakalamanın en güvenilir yoludur.

Değişkenlerine ayrılabilir denklemlerde bölme adımının sınırı

Birinci mertebeden denklem $dy/dx=f(x)g(y)$ biçimine getirilebiliyorsa değişkenler ayrılabilir. Amaç, $y$ içeren ifadeleri bir tarafa, $x$ içeren ifadeleri diğer tarafa almaktır:

$\frac{dy}{g(y)}=f(x)\,dx$

Ardından iki taraf ayrı ayrı integre edilir:

$\int\frac{dy}{g(y)}=\int f(x)\,dx+C.$

Bu yöntemde kritik nokta $g(y)$ ile bölme işlemidir. Eğer $g(y)=0$ olabilen sabit bir $y$ değeri varsa, bu değer bölme sırasında kaybolabilir ve ayrıca test edilmelidir. Örneğin $y'=y$ denkleminde $y=0$ doğrudan bir çözümdür; genel ifade $y=Ce^x$ içinde $C=0$ seçilerek görülse de benzer durumlar her denklemde otomatik olarak korunmayabilir.

İntegralden sonra mutlak değer de önemlidir. $dy/y=2x\,dx$ sonucundan $\ln|y|=x^2+C$ elde edilir. Sabit yeniden adlandırılarak $y=Ce^{x^2}$ yazılabilir. Başlangıç değeri pozitif veya negatifse $C$ buna göre belirlenir; çözümün işareti, katsayı sıfır değilse ilgili aralıkta korunur.

Bu yöntem özellikle büyüme-azalma, radyoaktif bozunma ve basit soğuma modellerinde kullanılır. Fakat denklem $y'+P(x)y=Q(x)$ biçiminde olup terimler doğal biçimde ayrılamıyorsa değişkenlerine ayırma yerine integrasyon çarpanı tercih edilir.

Doğrusal birinci mertebe denklemlerde integrasyon çarpanı

Doğrusal birinci mertebe denklem standart biçimde

$\frac{dy}{dx}+P(x)y=Q(x)$

olarak yazılır. Burada $P(x)$ ve $Q(x)$ bilinen fonksiyonlardır. Denklemin her iki tarafı

$\mu(x)=e^{\int P(x)\,dx}$

integrasyon çarpanıyla çarpılır. Çarpanın seçilme nedeni, sol tarafın bir çarpımın türevi hâline gelmesidir:

$\frac{d}{dx}[\mu(x)y]=\mu(x)Q(x).$

Sonra integral alınır ve $y$ yalnız bırakılır:

$\mu(x)y=\int\mu(x)Q(x)\,dx+C,$

$y=\frac{1}{\mu(x)}\left(\int\mu(x)Q(x)\,dx+C\right).$

Örneğin $y'+2y=4$ denkleminde $P(x)=2$ ve $Q(x)=4$ olduğundan $\mu=e^{2x}$ seçilir. Böylece $(e^{2x}y)'=4e^{2x}$ olur. İntegrasyon sonucunda $e^{2x}y=2e^{2x}+C$ ve $y=2+Ce^{-2x}$ bulunur. $y(0)=5$ verilirse $5=2+C$, yani $C=3$ olur.

Bu yöntemde önce denklemi standart biçime getirmek gerekir. Örneğin $2y'+4y=6$ ifadesinde her iki taraf 2'ye bölünmeden $P(x)$ ve $Q(x)$ yanlış okunabilir. Ayrıca integrasyon çarpanı $e^{\int P(x)dx}$ içinde integral sabitinin ayrı bir çarpan oluşturacağı ve bu çarpanın sonuca etkisiz biçimde yutulabileceği unutulmamalıdır.

Çözümlü örnekler: koşul kullanarak sabiti ve çözüm aralığını bulma

Örnek 1 — Doğrudan integral ve başlangıç koşulu

Verilen: $\frac{dy}{dx}=3x^2, \qquad y(0)=5.$

Çözüm adımları:

1. Sağ tarafı $x$'e göre integre edin: $y=\int3x^2dx=x^3+C$.
2. Başlangıç koşulunu kullanın: $y(0)=0^3+C=5$, dolayısıyla $C=5$.
3. Özel çözümü yazın: $y=x^3+5$.
4. Kontrol edin: $y'=3x^2$ ve $y(0)=5$ sağlanır.

Sonuç: $y=x^3+5$. Burada başlangıç koşulu, çözüm ailesindeki tüm eğriler arasından tek bir eğriyi seçmiştir.

Örnek 2 — Değişkenlerine ayrılabilir denklem

Verilen: $\frac{dy}{dx}=2xy, \qquad y(0)=2.$

Çözüm adımları:

1. $y$ sıfır değilmiş gibi değişkenleri ayırın: $dy/y=2x\,dx$.
2. İntegral alın: $\ln|y|=x^2+C$.
3. Üstel biçime geçin: $y=Ce^{x^2}$; sabit, işareti de kapsayacak şekilde yeniden adlandırılmıştır.
4. $y(0)=2$ koşulunu uygulayın: $2=C$.

Sonuç: $y=2e^{x^2}.$

Kontrol için türev alınır: $y'=4xe^{x^2}$. Ayrıca $2xy=2x(2e^{x^2})=4xe^{x^2}$ olduğundan denklem sağlanır. Bu örnekte $y=0$ da denklemin bir çözümüdür; fakat verilen başlangıç koşulu $y(0)=2$ olduğundan aranan özel çözüm sıfır çözümü değildir.

Örnek 3 — Doğrusal denklem

Verilen: $y'+2y=4, \qquad y(0)=5.$

Çözüm adımları:

1. Standart biçim zaten hazırdır; $P(x)=2$, $Q(x)=4$.
2. İntegrasyon çarpanını bulun: $\mu(x)=e^{2x}$.
3. Denklemi çarpanla çarpın: $(e^{2x}y)'=4e^{2x}$.
4. İntegral alın: $e^{2x}y=2e^{2x}+C$.
5. $y=2+Ce^{-2x}$ yazın.
6. $y(0)=5$ ile $C=3$ bulun.

Sonuç: $y=2+3e^{-2x}$. Çözüm $x$ büyüdükçe 2 değerine yaklaşır; bu yorum yalnızca incelenen değişken aralığı ve modelin geçerli olduğu varsayımlar içinde yapılmalıdır.

Analitik çözüm bulunamadığında yaklaşık sonuç ve davranış analizi

Yüksek mertebeli veya doğrusal olmayan denklemlerde kapalı biçimde bir çözüm bulmak zor, hatta seçilen fonksiyon sınıfında mümkün olmayabilir. Böyle durumlarda denklem doğrudan çözümsüz bırakılmaz; Euler ve Runge-Kutta gibi sayısal yöntemlerle belirli noktalardaki yaklaşık değerler hesaplanır. Bu yöntemler, çözümü sürekli bir formülle değil, seçilen adım büyüklüğü üzerinden ardışık tahminlerle üretir.

Birinci mertebeden $y'=f(x,y)$ ve $y(x_0)=y_0$ için basit Euler yaklaşımı

$y_{n+1}=y_n+h f(x_n,y_n)$

şeklindedir. Burada $h$, adım genişliğidir. $h$ küçültüldüğünde çoğu uygun problemde yaklaşım iyileşebilir; ancak hesap yükü artar ve yöntemin kararlılığı ayrıca değerlendirilmelidir. Bu nedenle tek bir yaklaşık sayı, hata büyüklüğü belirtilmeden tam çözüm gibi sunulmamalıdır.

Sistemlerde birden fazla bilinmeyen fonksiyon bulunabilir. İkinci mertebeden bir denklem, yeni değişkenler tanımlanarak birinci mertebeden denklem sistemine dönüştürülebilir. Doğrusal olmayan sistemlerde sabit noktalar, faz düzlemi ve kararlılık incelenerek uzun vadeli davranış hakkında bilgi edinilir. Sabit nokta, sistemin türevlerinin sıfır olduğu durumdur; fakat sabit noktanın bulunması tek başına kararlı olduğu anlamına gelmez.

Diferansiyel denklemler matematiksel modellemenin dinamik kısmını kurar. Modelleme varsayımları, optimizasyon hedefleriyle birleştirilebilir; sayısal analiz ise elde edilen denklemin yaklaşık çözümünü hesaplamak için kullanılır. Bu bağlantılar için matematiksel modelleme, sayısal analiz yöntemleri ve optimizasyon problemleri içeriklerine geçilebilir.

Sık yapılan hatalar ve karıştırılan sınır durumları

1. Mertebeyi derece sanmak: En yüksek türevin mertebesi ile o türevin üssü aynı kavram değildir. $(y'')^2+y=0$ ifadesinde mertebe 2, türev bakımından derece 2'dir; ancak bu iki bilgi tanım olarak farklıdır.

2. Genel çözümü özel çözüm sanmak: $y=x^3+C$ tüm çözüm ailesidir. Başlangıç koşulu uygulanmadan $C$ belirlenmiş olmaz.

3. Başlangıç koşulunu integrasyondan önce yanlış kullanmak: Önce genel çözüm bulunur, sonra koşul uygulanır. Koşul, integrasyon sabitini belirlemek için kullanılır.

4. Değişkenleri tam ayırmamak: $dy/y=2x\,dx$ biçiminde $y$ ve $x$ farklı taraflarda bulunmalıdır. Bir değişkenin terimi yanlış tarafta bırakılırsa integral sonucu denklemi sağlamaz.

5. Bölme sırasında çözüm kaybetmek: $y$ veya $g(y)$ ile bölmeden önce bu ifadeyi sıfır yapan durumlar kontrol edilmelidir. Sıfır fonksiyonu bazı denklemlerde ayrı bir özel çözüm olabilir.

6. Doğrusal denklemde standart biçimi atlamak: $ay'+by=c$ ifadesinde önce $a$'ya bölmek gerekir. Aksi hâlde integrasyon çarpanı yanlış seçilebilir.

7. Her birinci mertebe denklemi ayrılabilir sanmak: $y'+2y=4$ denkleminde $y$ terimini yalnızca bir tarafa taşıyarak $x$ terimlerinden bağımsız bir ayrılma elde edilemez; integrasyon çarpanı yöntemi uygundur.

8. Çözüm aralığını göz ardı etmek: Logaritma, payda veya kök içeren ifadeler gerçek sayılar kümesinde kısıtlı olabilir. Çözüm, başlangıç koşulunun bulunduğu noktayı içeren ve ifadenin tanımlı kaldığı bir aralıkta değerlendirilmelidir.

9. Sayısal sonucu kesin kabul etmek: Euler veya Runge-Kutta çıktısı yaklaşık değerdir. Adım büyüklüğü, yuvarlama ve kararlılık koşulları belirtilmeden analitik çözümle aynı kesinlikte yorumlanmamalıdır.