Задача коши для разностных уравнений
Данко П. Е., Попов А. Г., Кожевникова Т. Я. Высшая математика в упражнениях и задачах…
Часть II. Глава IV. Обыкновенные дифференциальные уравнения
§ 1. Дифференциальные уравнения первого порядка
1. Основные понятия. Дифференциальным уравнением называется уравнение, связывающее независимые переменные, их функцию и производные (или дифференциалы) этой функции. Если независимая переменная одна, то уравнение называется обыкновенным; если же независимых переменных две или больше, то уравнение называется дифференциальным уравнением в частных производных.
Наивысший порядок производной, входящей в уравнение, называется порядком дифференциального уравнения. Например:
1) х²у’ + 5xy = у² – обыкновенное дифференциальное уравнение первого порядка;
2) – обыкновенное дифференциальное уравнение второго порядка;
3) y’³ + y»y»’ = х – обыкновенное дифференциальное уравнение третьего порядка;
4) F (х, у, у’, у») = 0 – общий вид обыкновенного дифференциального уравнения второго порядка;
5) – уравнение в частных производных первого порядка.
В этом параграфе рассматриваются обыкновенные дифференциальные уравнения первого порядка, т. е. уравнения вида F (х, у, у’) = 0 или (в разрешенном относительно у’ виде) y’ = f(х, у).
Решением дифференциального уравнения называется такая дифференцируемая функция у = φ (x), которая при подстановке в уравнение вместо неизвестной функции обращает его в тождество. Процесс нахождения решения дифференциального уравнения называется интегрированием дифференциального уравнения.
Общим решением дифференциального уравнения первого порядка у’ = f(x, у) в области D называется функция у = φ(x, C), обладающая следующими свойствами: 1) она является решением данного уравнения при любых значениях произвольной постоянной С, принадлежащих некоторому множеству; 2) для любого начального условия у(х0) = у0 такого, что (x0; y0) ∈ 0, существует единственное значение С = С0, при котором решение у = φ(x, C0) удовлетворяет заданному начальному условию.
Всякое решение у = φ(x, C0), получающееся из общего решения у = φ (x, C) при конкретном значении С = С0, называется частным решением.
Задача, в которой требуется найти частное решение уравнения y’ = f(х, у) удовлетворяющее начальному условию у(х0) = y0, называется задачей Коши.
Построенный на плоскости хОу график всякого решения у = φ(х) дифференциального уравнения называется интегральной кривой этого уравнения. Таким образом, общему решению у = φ(х, С) на плоскости хОу соответствует семейство интегральных кривых, зависящее от одного параметра – произвольной постоянной С, а частному решению, удовлетворяющему начальному условию y(x0) = y0, – кривая этого семейства, проходящая через заданную точку М0(x0; у0).
Если функция f(х, у) непрерывна и имеет непрерывную производную в области D, то решение дифференциального уравнения у’= f (х, у) при начальном условии у(х0) = у0 существует и единственно, т. е. через точку (x0; y0) проходит единственная интегральная кривая данного уравнения (теорема Коши).
Особым решением называется такое решение, во всех точках которого условие единственности не выполняется, т. е. в любой окрестности каждой точки (х; у) особого решения существуют по крайней мере две интегральные кривые, проходящие через эту точку.
Особые решения не получаются из общего решения дифференциального управления ни при каких значениях произвольной постоянной С (в том числе и при С = ± ∞).
Особым решением является огибающая семейства интегральных кривых (если она существует), т. е. линия, которая в каждой своей точке касается по меньшей мере одной интегральной кривой.
Например, общее решение уравнения записывается в виде у = sin (х + С). Это семейство интегральных кривых имеет две огибающие: у = 1 и у = -1, которые и будут особыми решениями.
2. Дифференциальные уравнения с разделяющимися переменными. Дифференциальное уравнение вида
относится к типу уравнений с разделяющимися переменными. Если ни одна из функций f1(x), f2(y), φ1(x), φ2(y) не равна тождественно нулю, то в результате деления исходного уравнения на f2 (x) φ1 (y) оно приводится к виду
Почленное интегрирование последнего уравнения приводит к соотношению
которое и определяет (в неявной форме) решение исходного уравнения. (Решение дифференциального уравнения, выраженное в неявной форме, называют интегралом этого уравнения.)
507. Решить уравнение х(у²-4)dx + y dy = 0.
△ Разделив обе части уравнения на у² – 4 ≠ 0, имеем
x² + ln|у² – 4| = ln|C|, или у² – 4 = Сe -λ²
Это общее решение данного дифференциального уравнения.
Пусть теперь у² – 4 = 0, т. е. у = ± 2. Непосредственной подстановкой убеждаемся, что у = ±2 – решение исходного уравнения. Но оно не будет особым решением, так как его можно получить из общего решения при С = 0. ▲
508. Найти частный интеграл уравнения у’ cos х = у / ln у, удовлетворяющий начальному условию y(0) = l.
△ Полагая , перепишем данное уравнение в виде
Проинтегрируем обе части уравнения:
, или
Используя начальное условие у = 1 при х = 0, находим С = 0. Окончательно получаем
▲
509. Найти общий интеграл уравнения у’ = tg x tg y.
△ Полагая и разделяя переменные, приходим к уравнению ctg у dy = tg х dx. Интегрируя, имеем
, или ln|sin у| = -ln|cos x| + ln С.
Отсюда находим sin y = C/cos x, или sin y / cos x = С (общий интеграл). ▲
510. Найти частное решение дифференциального уравнения (l + x²)dy + y dx = 0 при начальном условии у(1) = 1.
△ Преобразуем данное уравнение к виду . Интегрируя, получим
, или ln |y| = – arctg x + С
Это и есть общий интеграл данного уравнения.
Теперь, используя начальное условие, найдем произвольную постоянную С; имеем ln 1 = — arctg 1 + С, т. е. С = π/4. Следовательно,
ln у = – arctg х + π/4,
откуда получаем искомое частное решение y = e π/4 – arctg x . ▲
Данко П. Е., Попов А. Г., Кожевникова Т. Я. Высшая математика в упражнениях и задачах… Ч. II. Стр. 117-119.
Задача Коши онлайн
Данная задача возникает при поиске частного решения дифференциального уравнения. Наш онлайн калькулятор, построенные на основе системы Wolfram Alpha, позволяет найти решение задачи Коши для различных типов дифференциальных уравнений. Чтобы начать работу, необходимо ввести данные своей задачи (дифференциальное уравнение и начальные условия) в калькулятор.
Найти решение задачи Коши для дифференциального уравнения:
при заданных начальных условиях:
При постановке задачи Коши, указываются так называемые начальные условия, позволяющие однозначно выделить искомое частное решение из общего. Эти условия включают в себя значения функции и всех её производных до включительно (где -порядок дифференциального уравнения), заданные в одной и той же точке .
Поясним вышесказанное на конкретном примере. Пусть нам требуется найти частное решение дифференциального уравнения:
удовлетворяющее начальным условиям:
Первым делом, используя различные методы (Бернули, вариации произвольной постоянной Лагранжа), сначала находим общее решение данного дифференциального уравнения:
Теперь, для поиска частного решения, нам необходимо использовать заданные начальные условия. Для этого, находим производную функции полученной ранее:
Далее, поставляем начальные условия в функцию и её производную :
Решая полученную систему уравнений получаем значения произвольных постоянных и :
Подставляем полученные результаты в общее решение дифференциального уравнения, в результате получаем искомое частное решение:
Другие полезные разделы:
Оставить свой комментарий:
Мы в социальных сетях:
Группа ВКонтакте | Бот в Телеграмме
Разностные схемы для решения задачи Коши
Рассмотрим основные принципы конструирования разностных схем.
1. Производная [math]\left.<\left(\frac
2. Дифференциальное уравнение заменяется разностным аналогом на одношаговом (двухточечном) шаблоне [math](x_,x_)[/math] (возможно использование дополнительных точек внутри шаблона) путем согласования с разложением функции [math]y=y(x)[/math] при [math]x=x_[/math] no формуле Тейлора относительно точки [math]x=x_[/math] . По этому принципу строятся одношаговые схемы, в частности семейство схем Рунге-Кутты.
3. Осуществляется интегрирование обеих частей дифференциального уравнения вдоль решения на выбранном шаблоне. При этом правая часть уравнения предварительно заменяется интерполяционным многочленом, степень которого должна соответствовать порядку точности схемы. Данный принцип называется интегрально-интерполяционным и используется для построения одношаговых и многошаговых явных и неявных схем, в том числе схем Адамса.
4. Разностные схемы строятся на основе подобия соотношений теории приближения конструкциям схем. С этой целью используются следствия из интегрально-дифференциальных сплайнов — параметрические соотношения, формулы аппроксимации производных, квадратурные формулы, записанные в общем случае на нерегулярном шаблоне. Данный принцип, который можно назвать принципом аналогий, далее будет применен для конструирования двух- и трехшаговых явных и неявных схем.
Принцип аппроксимаций
Рассмотрим проблему нахождения численного решения задачи Коши (6.9): [math]\frac
y(x_0)=y_0[/math] . Для аппроксимации производной [math]\left.<\left(\frac
и формулу (5.10), записанную на трехточечном нерегулярном шаблоне [math]H_<3,i>= (x_
\delta_= \frac
При [math]\delta_=1[/math] шаблон становится регулярным, а среднее слагаемое с входящим в него [math]y_[/math] — равным нулю. При этом формула (6.17) сводится к традиционной:
Аппроксимация формульных функций, входящих в правую часть (6.9), осуществляется путем перехода к их сеточному представлению.
Методика построения разностных схем с помощью аппроксимаций производной
1. Вводится в общем случае неравномерная сетка [math]\Omega_n= (x_0, x_1,\ldots, x_, x_,\ldots, x_
2. Для получения рассматриваемых здесь двух схем (одношаговой и двух-шаговой) выбирается соответственно двухточечный (одношаговый) и трехточечный (двухшаговый) шаблоны [math]H_<2,i>= (x_,x_)[/math] и [math]H_<3,i>= (x_
3. Выполняется подстановка аппроксимаций (6.16), (6.17) и [math]f(x_,y_)[/math] в дифференциальное уравнение:
После отбрасывания остаточных слагаемых получаются разностные схемы: — явная схема Эйлера первого порядка (явный метод Эйлера):
– обобщенная на нерегулярный шаблон двухшаговая явная схема Эйлера второго порядка (где [math]\Delta \widehat
Обозначим схему (6.20) символами 2Y2B (суть этого обозначения поясняется далее). Анализ остаточного слагаемого [math]\left(\frac
При [math]\delta_=1[/math] (регулярный шаблон) схема (6.20) соответствует методу Эйлера-Коши:
Для начала расчетов по формулам (6.20),(6.21) требуется иметь две «разгонные» точки [math]\widehat
Приведем геометрическую интерпретацию явного метода Эйлера (6.19), которая показана на рис. 6.5.
Пусть известна точка [math](x_,\widehat
Следующая точка [math]\widehat
Схема (6.19) иногда называется схемой ломаных (рис. 6.6), так как очередное значение [math]\widehat
Для анализа устойчивости явного метода Эйлера применим формулу (6.19) при [math]h_=h[/math] к решению задачи (6.12). В результате разностное уравнение будет иметь вид [math]\widehat
Дадим геометрическую интерпретацию метода Эйлера-Коши (6.21) применительно к решению задачи (6.9) (рис. 6.8). Пусть известны две точки [math](x_
Замечание. Все описанные здесь и далее методы легко переносятся на системы обыкновенных дифференциальных уравнений (6.7). Под [math]y[/math] и [math]f(x,y)[/math] следует подразумевать векторы [math]Y=(y_1,\ldots,y_n)^T[/math] и [math]F(x,Y)= \bigl(f_1(x,Y), \ldots, f_n(x,Y)\bigr)^T[/math] . Например, для системы двух дифференциальных уравнений
явный метод Эйлера (6.19) записывается в форме
Пример 6.3. Исследовать устойчивость явного метода Эйлера на примере решения уравнения [math]Ty’+y=1,[/math] [math]y(0)=y_0[/math] , где [math]T[/math] — действительное положительное число. Найти приближенное решение задачи Коши [math]0,\!1y’+y=1,
y(0)=0[/math] на отрезке [math][0;1][/math] явным методом Эйлера.
Точное решение задачи получено в примере 6.1: [math]y(x)=1-e^<-10x>[/math] , поскольку [math]T=0,\!1[/math] . Уравнение [math]Ty’+y=1[/math] перепишем в виде, аналогичном (6.12): [math]y’=-\frac<1>
Для уравнения [math]0,\!1y’+y=1[/math] постоянная времени [math]T=0,\!1[/math] , поэтому [math]h_<\text
Переписав уравнение в виде [math]y’=10-10y[/math] , получим из (6.19) формулу для нахождения приближенного решения:
6.1>>\\\hline x_& \widehat
(h=0,\!5)& y(x_)\\\hline 0,\!00& 0,\!000000& 0,\!000& 0,\!000& 0,\!000000 \\\hline 0,\!05& 0,\!500000& & & 0,\!393469 \\\hline 0,\!01& 0,\!750000& & & 0,\!632121 \\\hline 0,\!15& 0,\!875000& & & 0,\!776870 \\\hline 0,\!20& 0,\!937500& 2,\!000& & 0,\!864665 \\\hline 0,\!25& 0,\!968750& & & 0,\!917915 \\\hline 0,\!30& 0,\!984375& & & 0,\!950213 \\\hline 0,\!35& 0,\!992187& & & 0,\!969803 \\\hline 0,\!40& 0,\!996094& 0,\!000& & 0,\!981684 \\\hline 0,\!45& 0,\!998047& & & 0,\!988891 \\\hline 0,\!50& 0,\!999023& & 5,\!000& 0,\!993262 \\\hline 0,\!55& 0,\!999511& & & 0,\!995913 \\\hline 0,\!60& 0,\!999755& 2,\!000& & 0,\!997521 \\\hline 0,\!65& 0,\!999877& & & 0,\!998497 \\\hline 0,\!70& 0,\!999938& & & 0,\!999088 \\\hline 0,\!75& 0,\!999969& & & 0,\!999446 \\\hline 0,\!80& 0,\!999984& 0,\!000& & 0,\!999665 \\\hline 0,\!85& 0,\!999992& & & 0,\!999797 \\\hline 0,\!90& 0,\!999996& & & 0,\!999877 \\\hline 0,\!95& 0,\!999998& & & 0,\!999925 \\\hline 1,\!00& 0,\!999999& 2,\!000& -15,\!000& 0,\!999955 \\\hline \max\varepsilon_& 0,\!117879& 1,\!135335& 15,\!99995& \\\hline \end
Результаты расчетов, выполненных для различных значений [math]h[/math] , приведены в табл. 6.1. При [math]h=0,\!05[/math] приближенное решение достаточно близко к точному, при [math]h=0,\!2[/math] «пила» не сходится и не расходится, а при [math]h=0,\!5[/math] «пила», аппроксимирующая точное решение, расходится (рис. 6.9,в).
Замечание. Порядок точности метода, как правило, определяется порядком аппроксимации схем. Например, для схемы (6.19) первого порядка, записанной в виде, соответствующем исходному уравнению (6.9), получаем (первое слагаемое аппроксимирует производную):
то есть [math]y’_=f_+O(h_)[/math] . Здесь были использованы формула Тейлора и свойства, приведенные ранее.
Для двухшаговой схемы (6.21) аналогично получается второй порядок аппроксимации. Действительно, воспользуемся частным случаем формулы Тейлора (см. (В. 18)) при
Подставляя полученные соотношения в схему (6.21), записанную в виде [math]\frac
Также легко проверить, что и схема (6.20) на нерегулярном шаблоне имеет второй порядок при безусловной аппроксимации. Если же [math]h_\leqslant h_^2[/math] (условная аппроксимация), то этот порядок равен трем.
Интегрально-интерполяционный принцип
В данном подходе явные или неявные методы соответствующего порядка получаются путем интегрирования уравнения (6.9) вдоль решения [math]y(x)[/math] . Так, на двухточечном шаблоне
Функция [math]f(x,y(x))[/math] в (6.22) заменяется интерполяционным многочленом соответствующей степени [math]p[/math] ( [math]p=0[/math] для метода первого порядка, [math]p=1[/math] для метода второго порядка, [math]p=2[/math] для метода третьего порядка и т.д.). Подчеркнем, что замена [math]f(x,y)[/math] интерполяционным многочленом в (6.22) становится возможной благодаря тому, что [math]f(x,y)[/math] рассматривается не на всей плоскости [math](x,y)[/math] , а только вдоль интегральной кривой [math]y(x)[/math] . Для вычисления интеграла в правой части могут быть использованы известные квадратурные формулы прямоугольников, трапеций, парабол (Симпсона) и др.
Примем [math]x_[/math] в качестве узловой точки многочлена [math]L_<0>(x)[/math] , то есть [math]L_0(x)= f(x_,\widehat
Принимая для интерполяционного многочлена [math]L_1(x)[/math] в качестве узловых точек [math]x_[/math] и [math]x_[/math] и используя квадратурную формулу трапеций (5.46), получаем неявную одношаговую схему второго порядка (метод трапеций):
где [math]f_= f(x_,\widehat
При реализации алгоритма решения задачи Коши неизвестное значение [math]\widehat
и с дифференцированием функции [math]F(\widehat
Как правило, используется метод простых итераций:
Таким образом, уравнение (6.24) решается с помощью формулы (рекуррентного соотношения):
При применении методов Ньютона и простых итераций вначале задается или находится нулевое приближение решения по формуле [math]y_^<\,(0)>= \widehat
Итерации завершаются при выполнении условия окончания [math]\bigl|\widehat
Для сокращения описания и удобства ссылок получаемые здесь схемы обозначаются буквенно-цифровыми символами. Первая цифра в нем обозначает «шаговость» схемы, следующая за ней буква «Я» или буквы «НЯ» — явный или неявный характер, затем цифра указывает порядок точности схемы и за ней буква — возможную модификацию. Таким образом, схема (6.24) обозначается 1НЯ2, т.е. она является одношаговой, неявной и имеет второй порядок точности. Последняя буква указывается только для тех схем, которые имеют модификации.
Приведем геометрическую интерпретацию неявного метода Эйлера (6.23) (рис.6.10). Пусть известна точка [math](x_, \widehat
Исследуем устойчивость неявного метода Эйлера и метода трапеций. Применим формулу (6.23) при [math]h_=h[/math] к решению задачи (6.12). В результате получим разностное уравнение [math]\widehat
Применим формулу (6.24) метода трапеций при [math]h_=h[/math] к решению задачи (6.12). Получим разностное уравнение
Корень характеристического уравнения [math]\lambda_1= \left(1+ \frac
Продолжим рассмотрение интегрально-интерполяционного метода. Использование квадратурной формулы прямоугольников дает модифицированный метод Эйлера второго порядка:
Дадим геометрическую интерпретацию модифицированного метода Эйлера второго порядка (рис. 6.12). Пусть известна точка [math](x_,\widehat
Исследуем устойчивость модифицированного метода Эйлера второго порядка. Применим формулы (6.28) при [math]h_=h[/math] к решению задачи (6.12). В результате получим разностное уравнение:
Корень соответствующего характеристического уравнения [math]\lambda_1= 1+h\mu+ \frac<(h\mu)^2><2>[/math] . Условие устойчивости (6.13) имеет вид [math]\left|1+h\mu+ \frac<(h\mu)^2><2>\right| . Если [math]\mu[/math] — действительное число, то, раскрывая модуль, получаем [math]-1 . Неравенство 0″>[math]2+h\mu+ \frac<(h\mu)^2><2>>0[/math] выполняется для любых [math]h\mu[/math] , а неравенство [math]h\mu+ \frac<(h\mu)^2> <2>, то есть [math]h\mu\cdot (2+h\mu) , справедливо при [math]-2 . Таким образом, интервал устойчивости [math]h\mu\in (-2;0)[/math] модифицированного метода Эйлера совпадает с интервалом устойчивости явного метода Эйлера (6.19).
Для использования регулярного трехточечного шаблона [math](x_
Отметим, что схема парабол имеет четвертый, а не третий порядок точности из-за регулярности шаблона.
В соответствии с классификацией, схема (6.29) относится к классу неявных схем Адамса. Для получения других схем Адамса следует переписать (6.22) в форме
и заменить [math]y'(x)[/math] соответствующим интерполяционным полиномом, построенным на регулярном шаблоне с [math]h=\text
В силу формулы (4.36) второго интерполяционного многочлена Ньютона получаем (где [math]\widehat= \frac
Подставляя полученное выражение в правую часть (6.30) и учитывая, что [math]dx= h\cdot d\widehat[/math] , вычисляем [math]\textstyle<\int\limits_<0>^ <1>y'(\widehat
)\cdot h\, d\widehat
>[/math] , поскольку при [math]x=x_[/math] и [math]x=x_[/math] справедливо [math]\widehat
=0[/math] и [math]\widehat
=1[/math] соответственно. В результате имеем
или с учетом (6.9) общую формулу семейства явных схем Адамса:
где [math]\eta_= h\cdot f(x_, \widehat
Если взять в правой части (6.31) два слагаемых, то получится явный метод Эйлера (6.19), если три слагаемых, то явный метод, описываемый формулой
а если четыре слагаемых, то явный метод вида
В результате имеем схемы Адамса-Бэшфорта
– второго порядка (2Я2В):
– третьего порядка (ЗЯЗА):
– четвертого порядка (4Я4):
Для начала расчетов по формуле (6.32) требуются две «разгонные» точки: [math]\widehat
Если для построения интерполяционного многочлена использовать шаблоны, включающие точку [math]x_[/math] , то получится следующая формула семейства неявных схем Адамса :
где [math]\eta_= h\cdot f(x_, \widehat
Если взять в правой части (6.35) два слагаемых, то получится неявный метод Эйлера (6.23), если три слагаемых, — метод трапеций (6.24):
если четыре слагаемых, — неявный метод, описываемый соотношением (и т.д.)
– первого порядка (6.23);
– второго порядка (6.24);
– третьего порядка (2НЯЗ):
– четвертого порядка (ЗНЯ4):
Для расчетов по формуле (6.36) требуются две «разгонные» точки: [math]\widehat
Замечание. Аналогичные методы могут быть получены из интегрального уравнения вида [math]\textstyle<\int\limits_
k \geqslant 2>[/math] .
Пример 6.4. Найти приближенное решение задачи Коши [math]\begin
Путем прямой подстановки в систему легко убедиться в том, что точное решение задачи имеет вид
Используя соотношения (6.23),(6.24),(6.28),(6.33), выписываем формулы для нахождения приближенного решения указанными методами (при этом применяется векторная форма записи).
Для неявного метода Эйлера из (6.23) имеем
Разрешая эту систему относительно [math]\widehat
Соотношения (6.28) для модифицированного метода Эйлера принимают вид
Для метода Адамса-Бэшфорта третьего порядка из (6.33) находим
Для определения «разгонных» точек [math](x_0,\widehat
Выпишем соотношения для метода трапеций (6.24) и разрешим их относительно неизвестных:
откуда путем разрешения системы относительно [math]\widehat
6.4>>\\\hline x_& \begin
6.5>>\\\hline x_& \begin
Результаты проведенных расчетов даны в табл. 6.4 для [math]\widehat
Пример 6.5. Исследовать устойчивость метода Адамса-Бэшфорта третьего порядка на примере решения задачи Коши
\mu=-1[/math] и [math]\mu=-100[/math] .
Точное решение задачи: [math]y(x)=x^3[/math] . При [math]\mu=-1
\left(h\mu=-\frac<1> <8>\right)[/math] дифференциальное уравнение имеет вид [math]y’=x^3-y+3x^2[/math] , а при [math]\mu=-100
Формула (6.33) при [math]h=1\!\!\not<\phantom<|>>\,8[/math] принимает форму
Положим [math]f_0=y’\! \left(-\frac<1><4>\right)=0,\!1875;
\widehat
6.6>>\\\hline i& \widehat
(h\mu=-100\!\!\not<\phantom<|>>\,8)& \text
Из анализа решения, приведенного в табл. 6.6, видно, что в случае [math]h\mu=-100\!\!\not<\phantom<|>>\,8[/math] процесс не является устойчивым. Дадим обоснование этому факту.
Запишем формулу (6.33) в виде (6.14):
то есть [math]\alpha_0=1,
Составим уравнение (6.15) при каждом значении [math]h\mu[/math] .
При [math]h\mu=-\frac<1><8>[/math] все корни уравнения [math]\xi^3-\frac<73><96>\xi^2-\frac<16><96>\xi+ \frac<5><96>=0[/math] лежат внутри круга единичного радиуса с центром в начале координат, поэтому процесс устойчив.
При [math]h\mu=-\frac<100><8>[/math] один корень уравнения [math]\xi^3+ \frac<2204><96>\xi^2-\frac<1600><96>\xi+ \frac<500><96>=0[/math] , равный [math]\frac<1><3>\cdot \frac<2204><96>[/math] , очевидно не лежит внутри указанного круга. Поэтому процесс неустойчив. Эти выводы подтверждают результаты численного эксперимента. Ошибка вычисления значении [math]\widehat
Принцип согласования с разложением по формуле Тейлора
Рассмотрим задачу Коши (6.9) в общем случае на неравномерной сетке [math]\Omega_
Входящие в правую часть производные можно определить, дифференцируя уравнение (6.9) требуемое число раз:
Если функция [math]f(x,y)[/math] имеет p-е непрерывные производные по обоим аргументам, то в разложении (6.38) можно учесть слагаемые до [math]O(h_^
)[/math] . Чем больше членов разложения используется для вычисления, тем точнее будет приближение и тем выше порядок схемы. Однако рассчитывать производные путем дифференцирования правой части уравнения (6.9) невыгодно. Поэтому Рунге и Кутта независимо друг от друга предложили более удобную и рациональную форму представления численного решения, при использовании которой дифференцировать функцию [math]f(x,y)[/math] не требуется.
В общем случае методы Рунге-Кутты имеют следующую структуру (где [math][/math] ):
где [math]s[/math] — число стадий (этапов), [math]K_[/math] — значения коэффициентов схемы Рунге-Кутты, вычисленные на основе правой части уравнения (6.9), [math]c_
\ell= \overline<2,s>;[/math] [math]m=\overline<1,s-1>;[/math] [math]b_k,
k=\overline<1,s>[/math] . Первый индекс в обозначениях коэффициентов является порядковым номером, а второй соответствует индексу точки [math]x_[/math] — началу отрезка [math][x_,x_][/math] , на котором производится расчет.
В некоторых методах кроме вычисления приближенного решения [math]\widehat
порядок которого, как правило, на единицу больше или меньше обеспечиваемого выражением для [math]\widehat
Коэффициенты схемы подбираются так, чтобы выражение для [math]\widehat
[/math] , где степень [math]p[/math] меняется в зависимости от метода и соответствует порядку схемы.
Наибольшее распространение в вычислительной практике нашел метод Рунге-Кутты четвертого порядка :
Схема (6.41) является четырехчленной, первый коэффициент [math]K_<1,i>[/math] относится к точке [math]x_[/math] , второй и третий — к средней точке [math]x_+ \frac
Приведем геометрическую интерпретацию метода Рунге-Кутты четвертого порядка (рис. 6.13). Пусть известна точка [math](x_,\widehat
1. Схемы Рунге-Кутты легко обобщаются на системы дифференциальных уравнений. Так, применительно к системе двух уравнений
схема Рунге-Кутты четвертого порядка имеет вид [math](i=\overline<0,n-1>)[/math]
2. Соотношения модифицированного метода Эйлера (6.28) можно записать в форме схемы Рунге-Кутты:
Покажем, что эта схема второго порядка. Разложение (6.38) с учетом равенства [math]y»= f_x+ f_\cdot f_
Разложение функции [math]f(x,y)[/math] относительно точки [math](x_,y_)[/math] можно записать следующим образом:
где частные производные [math]f_x,f_y[/math] вычисляются в точке [math](x_,y_)[/math] . Тогда при [math]x=x_+\frac
Следовательно, формулу (6.43) можно переписать в виде
Сравнивая последнее соотношение с (6.44), можно сделать вывод о том, что схема модифицированного метода Эйлера согласуется с разложением по формуле Тейлора вплоть до членов степени [math]h_^2[/math] и поэтому является методом Рунге-Кутты второго порядка.
3. Коэффициенты семейства методов Рунге-Кутты (6.39),(6.40) удобно записывать в виде табл. 6.7. Табл. 6.8 соответствует классическому методу Рунге-Кутты четвертого порядка (6.41), а в табл. 6.9 приведены коэффициенты другого метода четвертого порядка (правило [math]3\!\!\not<\phantom<|>>\,8[/math] ). Табл. 6.10 соответствует методу Рунге–Кутты (6.43) второго порядка, табл. 6.11 и 6.12 — методам третьего порядка. В табл. 6.13 даны коэффициенты метода Фельберга четвертого порядка, где формула для [math]\widetilde
4. Для s-стадийных методов Рунге-Кутты справедливо следующее условие устойчивости, определяемое для задачи (6.12):
Если [math]\mu[/math] — действительное число, можно указать интервалы устойчивости: [math]h\mu\in (-2;0)[/math] при [math]s=1,\!2;
h\mu\in(-2,\!51;0)[/math] при [math]s=3;
h\mu\in (-2,\!78;0)[/math] при [math]s=4[/math] и т.д.
Пример 6.6. Для задачи Коши [math]y’=x+y,
x\in[0;0,\!2][/math] выполнить два шага методом Рунге-Кутты второго порядка (6.43) и один (первый шаг) методом четвертого порядка (6.41). Шаг интефирования принять постоянным и равным 0,1. Сравнить оба численных решения с точным [math]y(x)= 2e^x-x-1[/math] .
В данной задаче [math]x_0=0,
y_0=1[/math] . Применим схему (6.43).
Первый шаг [math](i=0)[/math] . В соответствии с постановкой задачи принимается [math]h_1=0,\!1[/math] . Вычислим
Тогда [math]\widehat
Второй шаг [math](i=1)[/math] . Так как шаг постоянный, то [math]h_2=0,!1[/math] , a [math]x_1=0,\!1[/math] и, следовательно,
Тогда [math]\widehat
Выполним теперь первый шаг методом четвертого порядка по формуле (6.41), записанной при [math]i=0\colon[/math]
При этом вычислим [math]K_<1,0>= x_0+ y_0=0+1=1[/math] ;
Точное решение: [math]y(0,\!1)= 1,\!1103418[/math] .
Сопоставляя точное решения [math]y(x_1)[/math] в точке [math]x_1[/math] с результатами [math]\widehat
В заключение данного раздела отметим преимущества и недостатки одношаговых и многошаговых схем:
1. Одношаговые схемы позволяют достаточно легко строить алгоритмы с переменным шагом и заданной точностью, что необходимо при расчете функций в областях их больших изменений. Однако они имеют недостаток — приближенные значения искомой функции, рассчитанные на предьщущих шагах, в этих схемах не используются.
2. В многошаговых схемах, напротив, учитываются значения искомой функции, рассчитанные на предыдущих шагах интегрирования, что повышает порядок точности этих схем без использования дополнительных точек на отрезке [math][x_,x_][/math] . Однако эти схемы требуют усложнения расчетных алгоритмов при реализации переменного шага.
Замечание. В явных одношаговых методах обеспечить заданную точность на каждом очередном шаге из текущей точки [math](x_, \widehat
1. Задать начальную величину шага [math]h^<(0)>= x_^<(0)>-x_[/math] и рассчитать значение [math]\widehat
2. Вычислить величину [math]\widehat
3. Если [math]\bigl|\widehat
4. Сравнить решения, полученные с шагом [math]h^<(1)>[/math] и [math]h^<(2)>[/math] в точке [math]x_^<(1)>[/math] . Если желаемая точность достигнута, расчеты продолжить с шагом [math]h^<(1)>[/math] , а если нет, то процесс дробления шага продолжить до достижения заданной точности.
Принцип аналогий
Данный метод следует из обобщения интегрально-интерполяционного принципа. Он основан на применении интерполяционных интегрально-дифференциальных сплайн-функций, преобразования которых — параметрические соотношения, аппроксимационные формулы для производных и квадратурные формулы — используются для получения разностных схем решения задачи Коши. Приведем некоторые из этих соотношений, соответствующих задаче аппроксимации некоторой произвольной функции [math]\Phi(x)[/math] многочленами второй степени:
а) параметрическое соотношение, связывающее интерполируемую функцию [math]\Phi(x_)[/math] и ее производную [math]\overline
б) формула аппроксимации первой производной в точке [math]x_[/math] на трехточечном нерегулярном шаблоне [math](x_
s,\,p[/math] — целые числа [math]\left(H_ = h_(s+ p\cdot \delta_), в) квадратурная формула экстраполяционного типа, имеющая четвертый порядок: Эта формула при [math]h_= h=\text г) двухинтервальная (трехточечная) обобщенная на нерегулярный шаблон квадратурная формула парабол четвертого порядка на отрезке [math][x_ д) одноинтервальная (трехточечная) функциональная квадратурная формула четвертого порядка на отрезке [math][x_ которая при [math]h_= h=\text Анализ остаточного слагаемого, получающегося при оценке погрешности квадратурной формулы (6.50) (см. мажоранту в скобках справа от данной формулы), указывает на то, что при сгущающейся вправо сетке (условная аппроксимация интеграла) порядок точности формулы может быть повышен на единицу без увеличения количества точек шаблона. Действительно, при [math]\delta_ легко получить приближенное условие повышения порядка: [math]h_ \leqslant h_^2[/math] . Таким образом, в данном случае порядок повышается за счет значительного уменьшения шага, такого, что [math]h_[/math] не превышает квадрата шага [math]h_[/math] . Этот факт может быть использован при решении задачи Коши на нерегулярном шаблоне; е) две одноинтервальные (двухточечные) квадратурные формулы функционально-дифференциального типа: ж) интегрально-функциональное параметрическое соотношение, осуществляющее связь интегралов и производных: Дискретные явные или неявные схемы для решения задачи (6.9) получаются из приведенных соотношений либо посредством аппроксимации производной, либо путем соответствующей замены одних параметров, относящихся к параметрическому соотношению или к квадратурной формуле, на другие параметры. Этот способ построения схем называется способом аналогий . Так, заменяя в (6.46) [math]\Delta\Phi_ (k=i,i-1)[/math] и разрешая относительно [math]\widehat которая при [math]h_= h= \text Комбинация схем (6.54) и (6.24) дает обобщенную на нерегулярный шаблон схему Адамса—Бэшфорта (2Я2В): которая при [math]h_= h= \text Аппроксимационная формула (6.47) определяет неявную двухшаговую схему второго порядка (2НЯ2): которая при [math]h_= h= \text Замечание. Схема (6.57) является одним из методов Гира , получаемых на основе формул дифференцирования назад, связывающих значение производной в точке [math]x_[/math] со значениями функции в предыдущих точках [math]x_,\ldots, x_ – схема третьего порядка: – схема четвертого порядка: – схема пятого порядка: – схема шестого порядка: Известно, что схема (6.58) является A(α)-устойчивой с [math]\alpha\cong 68^<\circ>[/math] . Все рассматриваемые ниже схемы получаются из квадратурных формул (6.48)-(6.52), в которых вместо разности первообразных [math]F_-F_= I_^[/math] подставляется [math]\Delta \widehat Преобразуя таким образом квадратурную формулу (6.48), получаем явную трехшаговую схему третьего порядка (ЗЯЗ): При [math]h_= h= \text Из обобщенной формулы (6.49) следует первая двухшаговая неявная схема третьего порядка (2НЯЗА), которая при [math]h_= h= \text Одноинтервальная (трехточечная) квадратурная формула (6.50) определяет вторую двухшаговую неявную схему третьего порядка (2НЯЗБ) Эта схема в соответствии с вышеприведенной оценкой точности интеграла при [math]\delta_ \leqslant h_[/math] имеет не третий, а четвертый порядок точности. При [math]h=\text И наконец, последние приводимые здесь схемы следуют из квадратурных формул (6.51), (6.52) (неявные одношаговые схемы) и параметрического соотношения (6.53) (неявная двухшаговая схема): где [math]f’=f_ При [math]h_= h= \text Для применения разностных схем, полученных с помощью четырех различных подходов, можно использовать общую методику. 1. В соответствии с требуемой точностью решения задачи Коши выбрать численный метод (схему) определенного порядка аппроксимации и типа (явный или неявный, одношаговый или многошаговый). 2. Задать сетку [math]\Omega_ i=\overline<0,n-1>[/math] , если сетка неравномерная). Если выбранный метод является ограниченно устойчивым, то для задания параметра сетки оценить величину критического шага. 3. Если выбранный метод многошаговый, найти требуемое число «разгонных» точек. 4. Найти приближенное решение задачи Коши во всех точках сетки [math]\Omega_ 5. Оценить достигнутую точность. В случае невыполнения требований повторить расчеты, меняя значение шага или используя более точный метод. http://mathforyou.net/online/calculus/cauchy/ http://mathhelpplanet.com/static.php?p=raznostnyye-skhemy-dlya-resheniya-zadachi-koshiМетодика решения задачи Коши методами сеток