Метод Нумерова - Википедия - Numerovs method

Метод Нумерова (также называемый методом Коуэлла) - это численный метод решения обыкновенные дифференциальные уравнения второго порядка, в котором член первого порядка не появляется. Это четвертого порядка линейный многоступенчатый метод. Этот метод является неявным, но его можно сделать явным, если дифференциальное уравнение является линейным.

Метод Нумерова разработан российским астрономом. Борис Васильевич Нумеров.

Метод

Метод Нумерова можно использовать для решения дифференциальных уравнений вида

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

В нем три значения ${ displaystyle y_ {n-1}, y_ {n}, y_ {n + 1}}$ взято в трех равноудаленных точках ${ displaystyle x_ {n-1}, x_ {n}, x_ {n + 1}}$ связаны следующим образом:

{ 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}),}

куда ${ Displaystyle у_ {п} = у (х_ {п})}$ , ${ displaystyle g_ {n} = g (x_ {n})}$ , ${ Displaystyle s_ {n} = s (x_ {n})}$ , и ${ displaystyle h = x_ {n + 1} -x_ {n}}$ .

Нелинейные уравнения

Для нелинейных уравнений вида

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

метод дает

{ 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}).}

Это неявный линейный многоступенчатый метод, который сводится к явному методу, данному выше, если ${ displaystyle f}$ линейно по ${ displaystyle y}$ установив ${ Displaystyle е (х, у) = - г (х) у (х) + s (х)}$ . Достигает точности порядка 4 (Хайрер, Норсетт и Ваннер, 1993 г., §III.10).

Заявление

В численной физике метод используется для нахождения решений одномерных Уравнение Шредингера для произвольных потенциалов. Примером может служить решение радиального уравнения для сферически-симметричного потенциала. В этом примере после разделения переменных и аналитического решения углового уравнения остается следующее уравнение радиальной функции ${ Displaystyle R (r)}$ :

{ 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).}

Это уравнение можно привести к виду, необходимому для применения метода Нумерова, с помощью следующей замены:

{ 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}}}.}

И когда мы делаем замену, радиальное уравнение принимает вид

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

или же

{ 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),}

что эквивалентно одномерному уравнению Шредингера, но с модифицированным эффективным потенциалом

{ 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}. }

Это уравнение мы можем решить так же, как если бы мы решили одномерное уравнение Шредингера. Мы можем переписать уравнение немного иначе и, таким образом, более ясно увидеть возможное применение метода Нумерова:

{ displaystyle { frac {d ^ {2} u} {dr ^ {2}}} = - { frac {2m} { hbar ^ {2}}} (E-V _ { text {eff}} (г)) и (г),}

{ displaystyle g (r) = { frac {2m} { hbar ^ {2}}} (E-V _ { text {eff}} (r)),}

{ displaystyle s (r) = 0.}

Вывод

Дано дифференциальное уравнение

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

Чтобы вывести метод Нумерова для решения этого уравнения, начнем с Расширение Тейлора функции, которую мы хотим решить, ${ Displaystyle у (х)}$ , вокруг точки ${ displaystyle x_ {0}}$ :

{ 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}) + { гидроразрыв {(x-x_ {0}) ^ {4}} {4!}} y '' '' (x_ {0}) + { frac {(x-x_ {0}) ^ {5}} { 5!}} Y '' '' '(x_ {0}) + { mathcal {O}} (h ^ {6}).}

Обозначая расстояние от ${ displaystyle x}$ к ${ displaystyle x_ {0}}$ к ${ displaystyle h = x-x_ {0}}$ , мы можем записать это уравнение в виде

{ 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}) .}

Если мы равномерно дискретизируем пространство, мы получим сетку ${ displaystyle x}$ точки, где ${ displaystyle h = x_ {n + 1} -x_ {n}}$ . Применяя приведенные выше уравнения к этому дискретному пространству, мы получаем связь между ${ displaystyle y_ {n}}$ и ${ displaystyle y_ {n + 1}}$ :

{ 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}).}

В вычислительном отношении это равносильно шагу вперед на сумму ${ displaystyle h}$ . Если мы хотим сделать шаг назад, мы заменяем каждый ${ displaystyle h}$ с ${ displaystyle -h}$ и получить выражение для ${ displaystyle y_ {n-1}}$ :

{ 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}).}

Обратите внимание, что только нечетные степени ${ displaystyle h}$ произошла смена знака. Суммируя два уравнения, получаем, что

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

Мы можем решить это уравнение для ${ displaystyle y_ {n + 1}}$ подставив выражение, данное в начале, то есть ${ displaystyle y '' _ {n} = - g_ {n} y_ {n} + s_ {n}}$ . Чтобы получить выражение для ${ displaystyle y '' '' _ {n}}$ фактор, мы просто должны различать ${ displaystyle y '' _ {n} = - g_ {n} y_ {n} + s_ {n}}$ дважды и снова аппроксимируйте его так же, как мы сделали это выше:

{ 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}).}

Если мы теперь подставим это в предыдущее уравнение, мы получим

{ 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}),}

или же

{ 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}).}

Это дает метод Нумерова, если мы проигнорируем член порядка ${ displaystyle h ^ {6}}$ . Отсюда следует, что порядок сходимости (в предположении устойчивости) равен 4.

Метод Нумерова - Википедия - Numerovs method

Содержание

Метод

Нелинейные уравнения

Заявление

Вывод

Рекомендации

внешняя ссылка