Numerov yöntemi

Numerov'un yöntemi (aynı zamanda Cowell'in yöntemi olarak da adlandırılır), birinci mertebeden terimin görünmediği ikinci mertebeden adi diferansiyel denklemleri çözmek için sayısal bir yöntemdir . Dördüncü dereceden doğrusal çok adımlı bir yöntemdir . Yöntem örtüktür, ancak diferansiyel denklem lineer ise açık hâle getirilebilir.

Numerov'un yöntemi Rus astronom Boris Vasilyeviç Numerov tarafından geliştirildi.

Numerov yöntemi, diferansiyel denklemleri şu formdaki çözmek için kullanılabilir:

${\displaystyle {\frac {d^{2}y}{dx^{2}}}=-g(x)y(x)+s(x).}$

İçindeki üç değer ${\displaystyle y_{n-1},y_{n},y_{n+1}}$, üç eşit uzaklıktan alınan ${\displaystyle x_{n-1},x_{n},x_{n+1}}$ değerleriyle aşağıdaki gibi ilişkilidir:

${\displaystyle y_{n+1}\left(1+{\frac {h^{2}}{12}}g_{n+1}\right)=2y_{n}\left(1-{\frac {5h^{2}}{12}}g_{n}\right)-y_{n-1}\left(1+{\frac {h^{2}}{12}}g_{n-1}\right)+{\frac {h^{2}}{12}}(s_{n+1}+10s_{n}+s_{n-1})+{\mathcal {O}}(h^{6}),}$

burada; ${\displaystyle y_{n}=y(x_{n})}$, ${\displaystyle g_{n}=g(x_{n})}$, ${\displaystyle s_{n}=s(x_{n})}$ ve ${\displaystyle h=x_{n+1}-x_{n}}$ 'dir.

Aşağıdaki formdaki doğrusal olmayan denklemler için

${\displaystyle {\frac {d^{2}y}{dx^{2}}}=f(x,y),}$

yöntem, şunu verir

${\displaystyle y_{n+1}-2y_{n}+y_{n-1}={\frac {h^{2}}{12}}(f_{n+1}+10f_{n}+f_{n-1})+{\mathcal {O}}(h^{6}).}$

Bu, eğer yukarıda verilen açık yönteme ${\displaystyle f}$ 'i lineer ${\displaystyle y}$ ayarlayarak ${\displaystyle f(x,y)=-g(x)y(x)+s(x)}$ indirgenirse, örtük doğrusal çok adımlı bir yöntemdir ve 4. derece doğruluğuna ulaşır Hairer, Nørsett & Wanner 1993 .

Sayısal fizikte yöntem, keyfî potansiyeller için tek boyutlu Schrödinger denkleminin çözümlerini bulmak için kullanılır. Bunun bir örneği, küresel olarak simetrik bir potansiyel için radyal denklemi çözmektir. Bu örnekte, değişkenleri ayırdıktan ve açısal denklemi analitik olarak çözdükten sonra, radyal fonksiyon ${\displaystyle R(r)}$'ın aşağıdaki denklemi ile kalıyoruz:

${\displaystyle {\frac {d}{dr}}\left(r^{2}{\frac {dR}{dr}}\right)-{\frac {2mr^{2}}{\hbar ^{2}}}(V(r)-E)R(r)=l(l+1)R(r).}$

Bu denklem, aşağıdaki yerine koyma ile Numerov yönteminin uygulanması için gerekli forma indirgenebilir:

${\displaystyle u(r)=rR(r)\Rightarrow R(r)={\frac {u(r)}{r}},}$

${\displaystyle {\frac {dR}{dr}}={\frac {1}{r}}{\frac {du}{dr}}-{\frac {u(r)}{r^{2}}}={\frac {1}{r^{2}}}\left(r{\frac {du}{dr}}-u(r)\right)\Rightarrow {\frac {d}{dr}}\left(r^{2}{\frac {dR}{dr}}\right)={\frac {du}{dr}}+r{\frac {d^{2}u}{dr^{2}}}-{\frac {du}{dr}}=r{\frac {d^{2}u}{dr^{2}}}.}$

Ve yerine koyduğumuzda, radyal denklem şöyle olur

${\displaystyle r{\frac {d^{2}u}{dr^{2}}}-{\frac {2mr}{\hbar ^{2}}}(V(r)-E)u(r)={\frac {l(l+1)}{r}}u(r),}$

veya

${\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}u}{dr^{2}}}+\left(V(r)+{\frac {\hbar ^{2}}{2m}}{\frac {l(l+1)}{r^{2}}}\right)u(r)=Eu(r),}$

bu, tek boyutlu Schrödinger denklemine eşdeğerdir, ancak değiştirilmiş etkin potansiyel ile şöyle olur;

${\displaystyle V_{\text{eff}}(r)=V(r)+{\frac {\hbar ^{2}}{2m}}{\frac {l(l+1)}{r^{2}}}=V(r)+{\frac {L^{2}}{2mr^{2}}},\quad L^{2}=l(l+1)\hbar ^{2}.}$

Bu denklemi, tek boyutlu Schrödinger denklemini çözdüğümüz gibi çözmeye devam edebiliriz. Denklemi biraz farklı bir şekilde yeniden yazabiliriz ve böylece Numerov'un yönteminin olası uygulamasını daha net görebiliriz:

${\displaystyle {\frac {d^{2}u}{dr^{2}}}=-{\frac {2m}{\hbar ^{2}}}(E-V_{\text{eff}}(r))u(r),}$

${\displaystyle g(r)={\frac {2m}{\hbar ^{2}}}(E-V_{\text{eff}}(r)),}$

${\displaystyle s(r)=0.}$

Aşağıda verilmiş olan diferansiyel denklemi ele alalım.

${\displaystyle y''(x)=-g(x)y(x)+s(x).}$

Numerov'un bu denklemi çözme yöntemini türetmek için, çözmek istediğimiz ${\displaystyle y(x)}$fonksiyonunun ${\displaystyle x_{0}}$ noktası etrafında Taylor açılımı ile başlıyoruz:

${\displaystyle y(x)=y(x_{0})+(x-x_{0})y'(x_{0})+{\frac {(x-x_{0})^{2}}{2!}}y''(x_{0})+{\frac {(x-x_{0})^{3}}{3!}}y'''(x_{0})+{\frac {(x-x_{0})^{4}}{4!}}y''''(x_{0})+{\frac {(x-x_{0})^{5}}{5!}}y'''''(x_{0})+{\mathcal {O}}(h^{6}).}$

Uzaklığı belirten ${\displaystyle x}$ ve ${\displaystyle x_{0}}$ ile ${\displaystyle h=x-x_{0}}$'yi işaret ederek, yukarıdaki denklemi şu şekilde yazabiliriz:

${\displaystyle y(x_{0}+h)=y(x_{0})+hy'(x_{0})+{\frac {h^{2}}{2!}}y''(x_{0})+{\frac {h^{3}}{3!}}y'''(x_{0})+{\frac {h^{4}}{4!}}y''''(x_{0})+{\frac {h^{5}}{5!}}y'''''(x_{0})+{\mathcal {O}}(h^{6}).}$

Uzayı eşit olarak ayrıklaştırırsak, ${\displaystyle x}$ noktalarının bir ızgarasını elde ederiz. Burada, ${\displaystyle h=x_{n+1}-x_{n}}$ 'dir. Yukarıdaki denklemleri bu ayrık uzaya uygulayarak, ${\displaystyle y_{n}}$ ve ${\displaystyle y_{n+1}}$arasında bir bağıntı elde ederiz:

${\displaystyle y_{n+1}=y_{n}+hy'(x_{n})+{\frac {h^{2}}{2!}}y''(x_{n})+{\frac {h^{3}}{3!}}y'''(x_{n})+{\frac {h^{4}}{4!}}y''''(x_{n})+{\frac {h^{5}}{5!}}y'''''(x_{n})+{\mathcal {O}}(h^{6}).}$

Hesaplamalı olarak, bu, ${\displaystyle h}$ miktarında ileri adım atmak anlamına gelir. Geriye doğru bir adım atmak istiyorsak, her ${\displaystyle h}$'ı ${\displaystyle -h}$ ile değiştiririz ${\displaystyle y_{n-1}}$ ifadesini elde ederiz:

${\displaystyle y_{n-1}=y_{n}-hy'(x_{n})+{\frac {h^{2}}{2!}}y''(x_{n})-{\frac {h^{3}}{3!}}y'''(x_{n})+{\frac {h^{4}}{4!}}y''''(x_{n})-{\frac {h^{5}}{5!}}y'''''(x_{n})+{\mathcal {O}}(h^{6}).}$

Sadece ${\displaystyle h}$'ın tek güçlerinin işaretinin değiştiğine dikkat edin. İki denklemi toplayarak şunu elde ederiz:

${\displaystyle y_{n+1}-2y_{n}+y_{n-1}=h^{2}y''_{n}+{\frac {h^{4}}{12}}y''''_{n}+{\mathcal {O}}(h^{6}).}$

Başta verilen ifadede yerine koyarak ${\displaystyle y_{n+1}}$ için bu denklemi çözebiliriz, yani ${\displaystyle y''_{n}=-g_{n}y_{n}+s_{n}}$ . ${\displaystyle y''''_{n}}$ faktörü için bir ifade elde etmek için, sadece iki kez ${\displaystyle y''_{n}=-g_{n}y_{n}+s_{n}}$ türevini almamız ve yukarıda yaptığımız gibi tekrar yaklaştırmamız gerekir:

${\displaystyle y''''_{n}={\frac {d^{2}}{dx^{2}}}(-g_{n}y_{n}+s_{n}),}$

${\displaystyle h^{2}y''''_{n}=-g_{n+1}y_{n+1}+s_{n+1}+2g_{n}y_{n}-2s_{n}-g_{n-1}y_{n-1}+s_{n-1}+{\mathcal {O}}(h^{4}).}$

Şimdi bunu önceki denklemin yerine koyarsak,

${\displaystyle y_{n+1}-2y_{n}+y_{n-1}={h^{2}}(-g_{n}y_{n}+s_{n})+{\frac {h^{2}}{12}}(-g_{n+1}y_{n+1}+s_{n+1}+2g_{n}y_{n}-2s_{n}-g_{n-1}y_{n-1}+s_{n-1})+{\mathcal {O}}(h^{6}),}$

veya

${\displaystyle y_{n+1}\left(1+{\frac {h^{2}}{12}}g_{n+1}\right)-2y_{n}\left(1-{\frac {5h^{2}}{12}}g_{n}\right)+y_{n-1}\left(1+{\frac {h^{2}}{12}}g_{n-1}\right)={\frac {h^{2}}{12}}(s_{n+1}+10s_{n}+s_{n-1})+{\mathcal {O}}(h^{6}).}$

Bu, derece terimini göz ardı edersek Numerov'un yöntemini verir. Yakınsama derecesi (kararlılık varsayılarak) 4'tür.

• Solving ordinary differential equations I: Nonstiff problems, Berlin, New York: Springer-Verlag, 1993, ISBN 978-3-540-56670-0 .

This book includes the following references:

• "A method of extrapolation of perturbations", Monthly Notices of the Royal Astronomical Society, 84, 1924, ss. 592-601, doi:10.1093/mnras/84.8.592 .
• "Note on the numerical integration of d²x/dt² = f(x,t)", Astronomische Nachrichten, 230, 1927, ss. 359-364, doi:10.1002/asna.19272301903 .