Riccati Benzeri Diferansiyel Denklem Çözümü

Diferansiyel denklemler dünyasına hoş geldiniz! Karşımızda oldukça ilginç bir başlangıç değer problemi var:

$$ \frac{dy}{dx}=y^2+x,\qquad y(0)=0. $$

Bu denklem, Riccati denklemlerine benzer bir yapıya sahip olmasına rağmen, genel Riccati formu $y' = q_0(x) + q_1(x)y + q_2(x)y^2$ ile tam olarak örtüşmez çünkü $q_1(x)$ terimi sıfır ve $q_0(x)$ terimi sadece $x$'e bağlı. Ancak bu, onu standart birinci dereceden lineer veya ayrılabilir denklemlerden ayırır. Ne yazık ki, bu tür denklemlerin kapalı formda, yani temel fonksiyonlar (polinomlar, üstel, logaritmik, trigonometrik fonksiyonlar) cinsinden analitik bir çözümü yoktur.

Bu tür denklemler genellikle özel fonksiyonlar veya seri çözümleri kullanılarak incelenir. Ancak sorunuz "kesin çözüm" talep ettiğinden ve "sadece sayısal bir yaklaşım değil" dediğinden, bu durum bizi özel fonksiyonlara yönlendiriyor.

Airy Fonksiyonları ve Riccati Denklemleri Arasındaki İlişki

Verilen denklem, Airy fonksiyonları ile yakından ilişkilidir. Riccati denklemleri, ikinci dereceden lineer homojen denklemlere dönüştürülebilir. Eğer elimizde $y' = q_0(x) + q_1(x)y + q_2(x)y^2$ şeklinde bir Riccati denklemi varsa, $y = -\frac{u'}{q_2(x)u}$ dönüşümü ile bu denklem $u'' - \left(\frac{q_2'(x)}{q_2(x)} + q_1(x)\right)u' + q_0(x)q_2(x)u = 0$ şeklinde ikinci dereceden lineer bir denkleme dönüşür.

Bizim denklemimiz için $q_0(x) = x$, $q_1(x) = 0$, $q_2(x) = 1$. Bu durumda dönüşüm $y = -\frac{u'}{u}$ olur.

Denklemde yerine koyarsak:

$\left(-\frac{u'}{u}\right)' = \left(-\frac{u'}{u}\right)^2 + x$

$-\frac{u''u - (u')^2}{u^2} = \frac{(u')^2}{u^2} + x$

$-\frac{u''}{u} + \frac{(u')^2}{u^2} = \frac{(u')^2}{u^2} + x$

$-\frac{u''}{u} = x$

$u'' + xu = 0$

Bu, Airy Denklemi olarak bilinen standart bir ikinci dereceden lineer diferansiyel denklemdir. Airy denkleminin çözümleri Airy fonksiyonları, yani $\text{Ai}(x)$ ve $\text{Bi}(x)$ ile ifade edilir.

Dolayısıyla, $u(x) = C_1 \text{Ai}(-x) + C_2 \text{Bi}(-x)$ şeklindedir. (Genellikle Airy fonksiyonları $\text{Ai}(x)$ ve $\text{Bi}(x)$ olarak tanımlanır, ancak $u'' - xu = 0$ denkleminin çözümleri $\text{Ai}(x)$ ve $\text{Bi}(x)$ iken, $u'' + xu = 0$ denkleminin çözümleri $\text{Ai}(-x)$ ve $\text{Bi}(-x)$'tir.)

Şimdi $y = -\frac{u'}{u}$ dönüşümünü kullanarak $y(x)$'i bulabiliriz.

$u'(x) = -C_1 \text{Ai}'(-x) - C_2 \text{Bi}'(-x)$

Böylece, kesin çözüm:

$$ y(x) = -\frac{-C_1 \text{Ai}'(-x) - C_2 \text{Bi}'(-x)}{C_1 \text{Ai}(-x) + C_2 \text{Bi}(-x)} = \frac{C_1 \text{Ai}'(-x) + C_2 \text{Bi}'(-x)}{C_1 \text{Ai}(-x) + C_2 \text{Bi}(-x)} $$

Başlangıç Koşulunu Uygulama

Başlangıç koşulumuz $y(0)=0$. Bu koşulu çözümümüze uygulayalım:

$$ y(0) = \frac{C_1 \text{Ai}'(0) + C_2 \text{Bi}'(0)}{C_1 \text{Ai}(0) + C_2 \text{Bi}(0)} = 0 $$

Bu ifadenin sıfır olması için pay kısmının sıfır olması gerekir:

$$ C_1 \text{Ai}'(0) + C_2 \text{Bi}'(0) = 0 $$

Airy fonksiyonlarının $x=0$'daki değerleri ve türevleri şunlardır:

• $\text{Ai}(0) = \frac{1}{3^{2/3} \Gamma(2/3)}$
• $\text{Ai}'(0) = -\frac{1}{3^{1/3} \Gamma(1/3)}$
• $\text{Bi}(0) = \frac{1}{3^{1/6} \Gamma(2/3)}$
• $\text{Bi}'(0) = \frac{3^{1/6}}{2 \Gamma(1/3)}$

Bu değerleri yerine koyalım:

$$ C_1 \left(-\frac{1}{3^{1/3} \Gamma(1/3)}\right) + C_2 \left(\frac{3^{1/6}}{2 \Gamma(1/3)}\right) = 0 $$

$\Gamma(1/3)$ sıfır olmadığından, her iki tarafı $\Gamma(1/3)$ ile çarpabiliriz:

$$ -C_1 \frac{1}{3^{1/3}} + C_2 \frac{3^{1/6}}{2} = 0 $$

$$ C_1 \frac{1}{3^{1/3}} = C_2 \frac{3^{1/6}}{2} $$

$$ C_1 = C_2 \frac{3^{1/6}}{2} \cdot 3^{1/3} = C_2 \frac{3^{1/6} \cdot 3^{2/6}}{2} = C_2 \frac{3^{3/6}}{2} = C_2 \frac{3^{1/2}}{2} = C_2 \frac{\sqrt{3}}{2} $$

Yani, $C_1 = \frac{\sqrt{3}}{2} C_2$.

Bu ilişkiyi $y(x)$ ifadesine geri koyarsak:

$$ y(x) = \frac{\frac{\sqrt{3}}{2} C_2 \text{Ai}'(-x) + C_2 \text{Bi}'(-x)}{\frac{\sqrt{3}}{2} C_2 \text{Ai}(-x) + C_2 \text{Bi}(-x)} $$

$C_2$ terimi hem pay hem de paydada sadeleşir (eğer $C_2 \neq 0$ ise):

$$ y(x) = \frac{\frac{\sqrt{3}}{2} \text{Ai}'(-x) + \text{Bi}'(-x)}{\frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x)} $$

Bu, verilen başlangıç değer probleminin Airy fonksiyonları cinsinden kesin çözümüdür.

Çözümün Var Olduğu Maksimal Reel Aralık

Çözüm $y(x) = -\frac{u'(x)}{u(x)}$ şeklinde tanımlandığından, çözümün varlığı $u(x)$'in sıfır olmadığı noktalara bağlıdır. Yani, $C_1 \text{Ai}(-x) + C_2 \text{Bi}(-x) \neq 0$ olmalıdır.

Yukarıdaki $C_1 = \frac{\sqrt{3}}{2} C_2$ ilişkisini kullanarak paydayı inceleyelim:

Payda $= C_2 \left( \frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x) \right)$

Airy fonksiyonları $\text{Ai}(x)$ ve $\text{Bi}(x)$ osilatör davranışına sahiptir. Özellikle $\text{Bi}(x)$, $x \to -\infty$ iken sonsuz sayıda sıfır noktasına sahiptir. $\text{Ai}(x)$ ise $x \to -\infty$ iken sonsuz sayıda sıfır noktasına sahiptir, $x \to \infty$ iken ise üstel olarak sıfıra yaklaşır ve pozitif $x$ değerleri için sıfır noktası yoktur.

Bizim durumumuzda, $\text{Ai}(-x)$ ve $\text{Bi}(-x)$ fonksiyonları söz konusu.

• $\text{Ai}(-x)$ fonksiyonu, $x > 0$ için sonsuz sayıda sıfır noktasına sahiptir.
• $\text{Bi}(-x)$ fonksiyonu, $x > 0$ için sonsuz sayıda sıfır noktasına sahiptir.

Payda ifadesi $F(x) = \frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x)$ şeklindedir. Bu, iki lineer bağımsız Airy fonksiyonunun lineer kombinasyonudur. Lineer bağımsız çözümlerin lineer kombinasyonları genellikle sonsuz sayıda sıfır noktasına sahip olabilir.

Ancak, başlangıç koşulumuz $y(0)=0$ olduğundan, $x=0$ noktasında payda sıfır değildir.

Payda($x=0$) $= C_1 \text{Ai}(0) + C_2 \text{Bi}(0) = C_2 \left( \frac{\sqrt{3}}{2} \text{Ai}(0) + \text{Bi}(0) \right)$.

$\text{Ai}(0) = \frac{1}{3^{2/3} \Gamma(2/3)}$ ve $\text{Bi}(0) = \frac{1}{3^{1/6} \Gamma(2/3)}$.

Payda($x=0$) $= C_2 \left( \frac{\sqrt{3}}{2} \frac{1}{3^{2/3} \Gamma(2/3)} + \frac{1}{3^{1/6} \Gamma(2/3)} \right)$

$= \frac{C_2}{\Gamma(2/3)} \left( \frac{\sqrt{3}}{2 \cdot 3^{2/3}} + \frac{1}{3^{1/6}} \right)$

$= \frac{C_2}{\Gamma(2/3)} \left( \frac{3^{1/2}}{2 \cdot 3^{2/3}} + \frac{1}{3^{1/6}} \right)$

$= \frac{C_2}{\Gamma(2/3)} \left( \frac{1}{2 \cdot 3^{1/6}} + \frac{1}{3^{1/6}} \right)$

$= \frac{C_2}{\Gamma(2/3)} \left( \frac{3}{2 \cdot 3^{1/6}} \right) \neq 0$ (eğer $C_2 \neq 0$ ise).

Dolayısıyla, $x=0$ noktasında payda sıfır değildir ve çözüm bu noktada tanımlıdır.

Genel olarak, $u''+xu=0$ denkleminin çözümleri olan $\text{Ai}(-x)$ ve $\text{Bi}(-x)$ fonksiyonları, $x \in (-\infty, \infty)$ aralığında tanımlıdır ve sürekli türevlenebilirdir. Ancak bunların lineer kombinasyonunun sıfır olduğu noktalar, çözümümüzün sonsuz hale geldiği veya tanımsız olduğu noktalardır.

Özellikle, $\text{Bi}(-x)$ fonksiyonu $x>0$ için sonsuz sayıda sıfıra sahipken, $\text{Ai}(-x)$ fonksiyonu $x>0$ için sonsuz sayıda sıfıra sahiptir. Bu sıfırlar, genellikle birbirleriyle çakışmazlar. Lineer kombinasyonları da genellikle sonsuz sayıda sıfıra sahip olacaktır.

Bu durumda, $y(x)$ çözümünün var olduğu maksimal reel aralık, paydanın sıfır olmadığı bir aralıktır ve $x=0$'ı içermelidir. Airy fonksiyonlarının özellikleri incelendiğinde, $\frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x)$ ifadesinin sıfır olduğu noktalar mevcuttur. Bu sıfır noktalarına $x_k$ diyelim. Çözüm, bu $x_k$ noktaları arasında tanımlıdır.

Ancak, başlangıç koşulu $y(0)=0$ ile belirlenen özel çözüm için, paydanın ilk sıfır noktası (yani $x=0$'a en yakın olanı) çözümün var...var olduğu aralığı belirler.

Çözümün Maksimal Reel Aralığına Devam

$F(x) = \frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x)$ fonksiyonunun sıfırlarını incelememiz gerekiyor. Bu fonksiyon, $u''+xu=0$ denkleminin bir çözümüdür. Bu denklemin çözümleri, $x \to \infty$ (yani $-x \to -\infty$) için salınımlı davranış sergiler ve sonsuz sayıda sıfıra sahiptir. Ancak $x \to -\infty$ (yani $-x \to \infty$) için, $\text{Ai}(-x)$ üstel olarak sıfıra yaklaşırken, $\text{Bi}(-x)$ üstel olarak büyür.

Daha spesifik olarak, $x>0$ için (yani $-x<0$ bölgesi), $\text{Ai}(-x)$ ve $\text{Bi}(-x)$ her ikisi de salınımlı davranış gösterir ve sonsuz sayıda sıfıra sahiptir. Bu sıfır noktaları, çözümümüzün paydasını sıfır yapar ve dolayısıyla çözümün tanımsız olmasına neden olur. $x=0$ noktasında payda sıfır olmadığından, çözüm $x=0$ noktasını içeren bir aralıkta tanımlıdır.

Bu aralık, $x=0$'ı içeren ve $F(x)$'in sıfır olmadığı en büyük aralıktır.

$F(x) = \frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x)$ fonksiyonunun ilk sıfırını bulmak için sayısal yöntemlere başvurmak gerekir. Airy fonksiyonlarının bilinen özelliklerinden ve grafiklerinden, $x>0$ bölgesinde bu tür bir lineer kombinasyonun sıfırları olduğu bilinmektedir.

Örneğin, $\text{Ai}(z)$ ve $\text{Bi}(z)$ fonksiyonlarının sıfırları $z_{Ai,k}$ ve $z_{Bi,k}$ ile gösterilir.

$\text{Ai}(z)$'nin ilk üç sıfırı yaklaşık olarak $-2.338, -4.088, -5.521$ iken, $\text{Bi}(z)$'nin ilk üç sıfırı yaklaşık olarak $-1.173, -3.271, -4.830$'dur.

Bizim durumumuzda argüman $-x$ olduğundan, $\text{Ai}(-x)$'in sıfırları $x \approx 2.338, 4.088, 5.521, \dots$ iken, $\text{Bi}(-x)$'in sıfırları $x \approx 1.173, 3.271, 4.830, \dots$'dir.

$F(x)$ fonksiyonunun sıfırlarını tam olarak belirlemek karmaşık bir iştir ve genellikle sayısal analiz gerektirir. Ancak, $x>0$ için $F(x)$'in sıfırları olduğu kesindir. Bu sıfırlardan $x=0$'a en yakın olanı, çözümün sağ tarafındaki üst sınırını belirleyecektir.

$x<0$ için (yani $-x>0$ bölgesi), $\text{Ai}(-x)$ üstel olarak küçülürken, $\text{Bi}(-x)$ üstel olarak büyür.

$\text{Ai}(z)$ için $z>0$ olduğunda sıfır yoktur.

$\text{Bi}(z)$ için $z>0$ olduğunda sıfır yoktur.

Dolayısıyla, $-x>0$ olduğunda, yani $x<0$ olduğunda, $\text{Ai}(-x) > 0$ ve $\text{Bi}(-x) > 0$ olacaktır.

Bu durumda, $F(x) = \frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x)$ ifadesi $x<0$ için her zaman pozitif olacaktır (çünkü $\text{Ai}(z)$ ve $\text{Bi}(z)$ pozitif $z$ için her zaman pozitiftir).

Bu demektir ki, $x<0$ bölgesinde payda hiçbir zaman sıfır olmaz.

Sonuç olarak, çözümün tanımsız olduğu noktalar sadece $x>0$ bölgesinde bulunur. Bu noktalardan $x=0$'a en yakın olanına $x_0 > 0$ diyelim.

$\frac{\sqrt{3}}{2} \text{Ai}(-x_0) + \text{Bi}(-x_0) = 0$ denklemini sağlayan en küçük pozitif $x_0$ değeri, çözümün varlık aralığının sağ sınırını belirleyecektir. Bu $x_0$ değeri sayısal olarak yaklaşık $2.338$ civarında olacaktır, çünkü $\text{Ai}(-x)$'in ilk sıfırı $x \approx 2.338$'dir ve $\text{Bi}(-x)$ bu noktada sıfır değildir.

Bu nedenle, çözümün var olduğu maksimal reel aralık $x \in (-\infty, x_0)$ şeklindedir, burada $x_0$ yukarıda tanımlanan $F(x)$'in en küçük pozitif sıfırıdır.

Maalesef, $x_0$'ı kapalı formda ifade etmek mümkün değildir; sayısal yöntemlerle bulunması gerekir. Ancak, $x<0$ için paydanın sıfır olmadığını ve $x>0$ için sıfırları olduğunu matematiksel olarak gösterebildik.

Özetle:

1. Kesin Çözüm:

$$ y(x) = \frac{\frac{\sqrt{3}}{2} \text{Ai}'(-x) + \text{Bi}'(-x)}{\frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x)} $$

Burada $\text{Ai}(z)$ ve $\text{Bi}(z)$ standart Airy fonksiyonlarıdır.

1. Maksimal Reel Aralık:

Çözüm, $x \in (-\infty, x_0)$ aralığında var olur. Burada $x_0$, $\frac{\sqrt{3}}{2} \text{Ai}(-x) + \text{Bi}(-x) = 0$ denklemini sağlayan en küçük pozitif $x$ değeridir. Bu değerin kapalı formu mevcut değildir ve sayısal olarak belirlenmesi gerekir. $x<0$ için payda sıfır olmadığından, aralık sol taraftan sınırsızdır.

Bu problem, standart analitik yöntemlerle çözülemeyen ancak özel fonksiyonlar aracılığıyla kesin çözümü ifade edilebilen diferansiyel denklemlerin güzel bir örneğidir.