В первую очередь рассмотрим функциональные методы решения простейших дифференциальных уравнений, применение которых дает решение в виде функции, представленной аналитически.
Интегрирование дифференциальных уравнений с помощью степенных рядов.
Рассмотрим этот метод в общем случае для дифференциального уравнения n-го порядка
при начальных условиях:
Предположим, что правая часть уравнения (1) является аналитической функцией в начальной точке, т. е. в некоторой окрестности этой точки разлагается в степенной ряд вида:
где - целые неотрицательные числа и - некоторые постоянные коэффициенты.
Тогда интеграл уравнения (1), отвечающий начальным условиям (2), является аналитическим в точке и, пользуясь рядом Тейлора, можно положить
Первые n+1 коэффициентов ряда (3) определяются непосредственно из начальных условий (2) и дифференциального уравнения (1). Для нахождения следующего (n+2)-го коэффициента продифференцируем уравнение (1) по правилу дифференцирования сложной функции. В результате получим
где для удобства принято.
где значок "0" означает, что значения соответствующих производных берутся в точке (). Повторяя этот прием шаг за шагом, можно найти и дальнейшие производные,….
Пример: Написать несколько членов разложения в степенной ряд решения уравнения
удовлетворяющего начальным условиям; .
Решение: Полагаем
Из уравнения (4) получаем
Отсюда. Дифференцируя последовательно уравнение (5), будем иметь
Из этих равенств вытекает, что
Следовательно,
Написать общий член ряда не представляет больших затруднений.
Метод разложения решения дифференциального уравнения в степенные ряды часто используется как элемент более практичных методов приближенного интегрирования дифференциальных уравнений. В частности, для некоторых численных методов интегрирования дифференциальных уравнений требуется определить значение искомых функций в нескольких точках. Эти значения при соблюдении известных условий гладкости данного уравнения могут быть с любой степенью точности подсчитаны с помощью степенных рядов.
Метод Пикара (последовательных приближений).
Этот метод позволяет получить приближенное решение дифференциального уравнения в виде функции, представленной аналитически. Метод Пикара возник в связи с доказательством теоремы существования и единственности решения уравнения и является, по сути, одним из применений принципа сжимающих отображений.
Пусть в условиях теоремы существования требуется найти решение уравнения с начальным условием. Проинтегрируем обе части уравнения от до:
Очевидно, решение этого интегрального уравнения будет удовлетворять дифференциальному уравнению и начальному условию. Действительно, при получим:
Вместе с тем интегральное уравнение (1) позволяет применить метод последовательных приближений. Положим и получим из (1) первое приближение:
Интеграл в правой части содержит только переменную; после нахождения этого интеграла будет получено аналитическое выражения приближения как функции переменной. Заменим теперь в уравнении (1) у найденным значением и получим второе приближение:
и т.д. В общем случае итерационная формула имеет вид:
(n=1, 2,…). (2)
Циклическое применение этой формулы дает последовательность функций,…,….
Так как функция непрерывна в области, то она ограничена в некоторой области, содержащей точку, т.е. .
Применяя к уравнению (2) в условиях теоремы существования принцип сжимающих отображений, нетрудно показать, что последовательность (3) сходится (имеется в виду сходимость по метрике
в пространстве непрерывных функций, определенных на сегменте, таких, что). Ее предел является решением интегрального уравнения (2), а следовательно, и дифференциального уравнения с начальными условиями. Это означает, что k-й член последовательности (3) является приближением к точному решению уравнения с определенной степенью точности.
Оценка погрешности k-го приближения дается формулой:
где - константа Липшица, - верхняя грань модуля функции из неравенства (4), а величина для определения окрестности вычисляется по формуле
Геометрически последовательные приближения представляют собой кривые (n=1, 2, …), проходящие через общую точку.
Пример: Методом последовательных приближений найти приближенное решение дифференциального уравнения, удовлетворяющее начальному условию.
Запишем для данного случая формулу вида (2):
Аналогично:
Оценим погрешность третьего приближения. Для определения области, заданной как
примем, например, . В прямоугольнике функция определена и непрерывна, причем:
По формуле (6) находим. Используя оценочную формулу (5), получаем: .
Метод Эйлера.
В основе метода ломаных Эйлера лежит идея графического построения решения дифференциального уравнения, однако этот метод дает одновременно и способ нахождения искомой функции в численной (табличной) форме.
Пусть дано уравнение с начальным условием. Выбрав достаточно малый шаг, построим, начиная с точки, систему равностоящих точек (i=0, 1, 2,…). Вместо искомой интегральной кривой на отрезке рассмотрим отрезок касательной к ней в точке (обозначим ее) с уравнением
При из уравнения касательной получаем: , откуда видно, что приращение значения функции на первом шаге имеет вид: .
Аналогично, проводя касательную к некоторой интегральной кривой семейства в точке, получим:
что при дает
т.е. получается из добавлением приращения
Таким образом, получение таблицы значений искомой функции по методу Эйлера заключается в циклическом применении пары формул:
K=0, 1, 2,….
Метод Эйлера обладает малой точностью, к тому же погрешность каждого нового шага, вообще говоря, систематически возрастает. Наиболее приемлемым для практики методом оценки точности является в данном случае способ двойного счета - с шагом h и с шагом h/2. Совпадение десятичных знаков в полученных двумя способами результатах дает естественные основания считать их верными.
Если правая часть уравнения непрерывна, то последовательность ломаных Эйлера при на достаточно малом отрезке равномерно стремиться к искомой интегральной кривой.
Пример: Применяя метод Эйлера, составить на отрезке таблицу значений интеграла дифференциального уравнения
удовлетворяющего начальному условию, выбрав шаг h=0,1.
Результаты вычислений приведены в таблице. Для сравнения в последнем столбце помещены значения точного решения
|
Точное значение |
|||||
|
|
|
|
Из приведенной таблицы видно, что абсолютная погрешность значения составляет. Отсюда относительная погрешность примерно равна 3%.
Модификации метода Эйлера (метод Эйлера - Коши).
Существуют различные уточнения метода Эйлера, повышающие его точность. Модификации метода обычно направлены на то, чтобы более точно определить направление перехода из точки в точку.
Рассмотрим снова дифференциальное уравнение с начальным условием.
Выбрав шаг h, положим
(i=0, 1, 2,…).
Согласно методу Эйлера последовательные значения искомого решения вычисляются по приближенной формуле
где. Более точным является усовершенствованный метод ломаных, при котором сначала вычисляют промежуточные значения
и находят значение направления поля интегральных кривых в средней точке, т.е.
а затем полагают
Другой модификацией метода Эйлера является усовершенствованный метод Эйлера - Коши, при котором сначала определяется "грубое приближение" решения
исходя из которого находится направление поля интегральных кривых
Затем приближенно полагают
Геометрически это означает, что мы определяем направление интегральной кривой в исходной точке (xi, yi) и во вспомогательной точке (xi+1, y*i+1), а в качестве окончательного берем среднее этих направлений.
Пример: Первым и вторым усовершенствованными методами Эйлера проинтегрировать уравнение
на отрезке .
Примем шаг h=0,2 и.
Приближенные значения искомого решения, определенные с помощью усовершенствованного метода ломаных, помещены в таблице 1.
В таблице 2 приведены результаты вычислений интеграла усовершенствованным методом Эйлера - Коши, причем шаг сохранен прежний h=0,2.
Таблица 1
|
|
|
|
Для сравнения приводим точное решение, откуда …
Таблица 2
 |
|
|
|
|
|
Метод Эйлера и его модификации являются простейшими представителями конечно-разностных методов (шаговых методов) для приближенного решения обыкновенных дифференциальных уравнений вида (k=1, 2,…, n) при заданных начальных условиях (k=1, 2,…, n).
При применении конечно-разностного метода искомое решение yk=yk(x) (k=1, 2,…, n) последовательно строится на системе точек (узлов) xi=x0+ih (I=0, 1, 2,…), где h - выбранный шаг. Процесс вычислений расчленяется на однообразно повторяющиеся циклы, каждый из которых обеспечивает переход от значения yk(xi) к значению yk(xi+1), начиная с начального y(0)k. Поэтому схема вычислений, вообще говоря, легко программируется и удобна для реализации на ЭВМ.
Метод Рунге-Кутта
Метод Эйлера и метод Эйлера-Коши относятся к семейству методов Рунге-Кутта, имеющих следующий вид. Фиксируем некоторые числа:
последовательно вычисляем:
k2(h)=hf(x+2h,y+21k1(h)),
kq(h)=hf(x+qh,y+q1k1(h)+…+qq-1kq-1(h))
и полагаем:
Рассмотрим вопрос о выборе параметров i, pi, ij. Обозначим
(h)=y(x+h)-z(h).
Будем предполагать, что
(0)="(0)=…=(s)(0)=0
при любых функциях f(x,y), а (s+1)(0)0 для некоторой функции f(x,y). По формуле Тейлора справедливо равенство
где 01. Величина (h) называется погрешностью метода на шаге, а s - порядком погрешности метода.
При q=1 будем иметь:
- (h)=y(x+h)-y(x)-p1hf(x,y),
- (0)=0,
"(0)=(y"(x+h)-p1f(x,y))h=0=f(x,y)(1-p1),
”(h)=y”(x+h).
Ясно, что равенство "(0)=0 выполняется для любых функций f(x,y) лишь при условии, что р1=1. Легко видеть, что при этом значении р1 из формулы (1) получается формула
yk+1=yk+yk, k=0, 1, 2,…,
т.е. в этом случае мы получаем метод Эйлера. Для погрешности этого метода на шаге согласно (2) будем иметь:
Рассмотрим случай q=2, тогда
(h)=y(x+h)-y(x)-p1hf(x,y)-p2hf(x*,y*),
x*=x+2h, y*=y+21hf(x,y).
Согласно исходному дифференциальному уравнению
y=f, y”=fx+fyf, y"""=fxx+2fxyf+fyyf2+fyy”. (3)
(Здесь для краткости через y и f обозначены y(x) и f(x,y) соответственно.)
Вычисляя производные функции (h) и подставляя в выражения для (h), "(h) и ”(h) значение h=0, получим (с учетом соотношений (3)):
"(0)=(1-p1-p2)f,
”(0)=(1-2p22)fx+(1-2p221)fyf.
Хорошо видно, что требование
(0)="(0)=”(0)=0
будет выполняться для всех f(x,y) лишь в том случае, если одновременно будут справедливы следующие три равенства относительно четырех параметров:
- 1-p1-p2=0,
- 1-2p22=0, (4)
- 1-2p221=0.
Произвольно задавая значение одного из параметров, и определяя значения остальных из системы (4), мы будем получать различные методы Рунге-Кутта с порядком погрешности s=2. Например, при р1= из (4) получаем: р2=, 2=1, 21=1. Для этих значений параметров формула (1) принимает вид:
(Здесь yi+1 записано вместо y(x+h), yi - вместо y(x), а через y*i+1 обозначено выражение yi+hf(xi,yi).) Таким образом, в рассматриваемом случае мы приходим к расчетным формулам
y*i+1=yi+hfi, f*i+1=f(xi+1, y*i+1,
соответствующим методу Эйлера-Коши. Из (2) следует, что при этом главная часть погрешности на шаге есть, т.е. пропорциональна третьей степени шага h.
В вычислительной практике наиболее часто используется метод Рунге-Кутта с q=4 и s=4.
Приведем без вывода один из вариантов соответствующих расчетных формул:
k2=hf(x+h/2, y+k1/2),
k3=hf(x+h/2, y+k2/2), (5)
k4=hf(x+h, y+k3),
y=z(h)-y(x)=(k1+2k2+2k3+k4).
Отметим, что в этом случае погрешность на шаге пропорциональна пятой степени шага (h5). Отсюда, в частности, следует, что при достаточно малом h и малых погрешностях вычислений решение уравнения y"=f(x,y), полученное методом Рунге-Кутта по формулам (5), будет близким к точному.
Эффективная оценка погрешности метода Рунге-Кутта затруднительна. Поэтому для определения правильности выбора шага h на практике обычно на каждом этапе из двух шагов применяют двойной пересчет. А именно, исходя из текущего верного значения y(xi), вычисляют величину y(xi+2h) двумя способами: один раз с шагом h, а другой раз с двойным шагом H=2h. Если расхождение полученных значений не превышает допустимой погрешности, то шаг h для данного этапа выбран правильно и полученное с его помощью значение можно принять за y(xi+2h). В противном случае шаг уменьшают в два раза.
Такого рода вычислительную схему легко запрограммировать для работы на ЭВМ.
Геометрический смысл использования метода Рунге-Кутта с расчетными формулами (5) состоит в следующем. Из точки (xi, yi) сдвигаются в направлении, определяемом углом 1, для которого tg1=f(xi,yi). На этом направлении выбирается точка с координатами (xi+, yi+). Затем из точки (xi, yi) сдвигаются в направлении, определяемом углом 2, для которого
tg2=f(xi+, yi+),
и на этом направлении выбирается точка с координатами (xi+, yi+). Наконец, из точки (xi, yi) сдвигаются в направлении, определяемом углом 3, для которого
и на этом направлении выбирается точка с координатами (xi+h, yi+k3). Этим задается еще одно направление, определяемое углом 4, для которого tg4=f(xi+h,yi+k3). Четыре полученные направления усредняются в соответствии с последней из формул (5). На этом окончательном направлении и выбирается очередная точка (xi+1, yi+1)= (xi+h, yi+y).
Для вычисления по формулам (5) удобно пользоваться схемой, приведенной в таблице 1.
Таблица 1.
Метод Рунге-Кутта обладает значительной точностью и, несмотря на свою трудоемкость, широко используется при численном решении дифференциальных уравнений с помощью ЭВМ. Кроме того, важным преимуществом этого метода является возможность применения "переменного шага".
Пример: Методом Рунге-Кутта вычислить на отрезке интеграл дифференциального уравнения y"=x+y, y(0)=1, приняв шаг h=0,1.
Покажем начало процесса. Вычисление у1. Последовательно имеем:
k1=(0+1)0,1=0,1;
k2=0,05+(1+0,05)0,1=0,11;
k3=0,05+(1+0,055)0,1=0,1105;
k4=0,1+(1+0,1105)0,1=0,12105.
Отсюда у0=(0,1+20,11+20,1105+0,12105)=0,1103 и, следовательно,
дифференциальный уравнение пикар эйлер
у1=у0+у0=1+0,1103=1,1103.
Аналогично вычисляются дальнейшие приближения. Результаты вычислений приведены в таблице 2.
Таким образом, у(0,5)=1,7974.
Для сравнения приводим точное решение:
откуда у(0,5)=2е1/2-1,5=1,79744…
Таблица 2
|
|
|
||
|
|
|
||
|
|
|
||
|
|
|
||
|
|
|
||
Практическое занятие №25
«Численное интегрирование с помощью формул прямоугольников, трапеций, Симпсона.
Приближенное решение дифференциальных уравнений с помощью формул Эйлера»
1. Цель: Выработать навыки и умения по применению методов приближённого
интегрирования – формул прямоугольников, трапеций и Симпсона, в решении
приближенными методами дифференциальных уравнений
Пояснения к работе
2.1 Краткие теоретические сведения:
Формула прямоугольников
Известно, что не для всякой непрерывной функции ее первообразная выражается через элементарные функции. Кроме того, на практике сталкиваются с необходимостью вычислять интегралы от функций, заданных табличным или графическим способами, а так же интегралы от функций, первообразные которых выражаются через элементарные очень сложно, что требует большой вычислительной работы и с практической точки зрения нерационально. В этих случаях вычисление определенного интеграла по формуле Ньютона – Лейбница затруднительно, поэтому прибегают к различным методам приближенного интегрирования. Наиболее простым методом приближенного вычисления определенного интеграла является метод прямоугольников, основанный на непосредственном определении интеграла:
где есть интегральная сумма, соответствующая некоторому разбиению отрезка и некоторому набору точек отрезка разбиения.
Вычисление определенного интеграла геометрически сводится к вычислению площади криволинейной трапеции, ограниченной неотрицательной функцией f(x) , осью абсцисс и прямыми и .
Для точности численного интегрирования нужно отрезок разбить на несколько частей и для каждой из них вычислить приближенное значение площади криволинейной трапеции, основанием которой является отрезок , а высотой - число , т.е. значение функции в точке
Выбранное из условия минимума ошибки интегрирования. Тогда за приближенное значение
интеграла на отрезке принимают интегральную сумму:
Практически удобно делить отрезок на равные части 
, а точки совмещать с
левыми 
или правыми 
концами отрезков разбиения. Если точку совместить с левым концом отрезка , то приближенное значение интеграла геометрически равно площади заштрихованной нижней ступенчатой фигуры и может быть представлено формулой левых прямоугольников:
(1)
где - шаг разбиения. Если же в качестве точки . выбрать правый конец отрезка , приближенное значение интеграла графически равно площади верхней ступенчатой фигуры, и вычисляется по формуле правых прямоугольников:
(2)
Погрешность вычисления:
, где - максимум на (3)
Пример 1. Используя формулу прямоугольников при , вычислить с тремя десятичными знаками . Оценить допущенную погрешность.
Решение: разделим отрезок на 10 равных частей точками и найдём значения функции в этих точках:
| 1 | 1.1 | 1.2 | 1.3 | 1.4 | 1.5 | 1.6 | 1.7 | 1.8 | 1.9 | 2 | |
| 1.000 | 0.909 | 0.833 | 0.769 | 0.714 | 0.667 | 0.625 | 0.588 | 0.556 | 0.526 | 0.5 |
Тогда получим и по формуле (1) находим
.
Оценим погрешность. Имеем 
; функция монотонно убывает на отрезке , поэтому 
и 
.
Так как допущенная погрешность влияет уже на второй знак после запятой, то третий знак следует округлить. Значит, 
. Если вычислить этот интеграл по формуле Ньютона –
Лейбница, то получим 
. Таким образом, ответ является приближённым значением . Но ; следовательно, при вычислении допущена погрешность, меньшая .
Формула трапеций
Приближенное значение определенного интеграла можно вычислить и иным способом.
Заменим на отрезке дугу АВ графика подынтегральной функции у = f(x) стягивающей ее хордой (рис.2) и вычислим площадь трапеции АВbа. Примем значение определенного
интеграла численно равным площади этой трапеции.
(4)
Это и есть формула трапеций для приближенного вычисления интеграла. Погрешность вычисления 
для формулы трапеций оценивается так:
, (5)
где точка 
. В случае, если , вычисление по формуле (4) даёт значение интеграла с избытком; если , то интеграл вычисляется с недостатком. Точность вычислений возрастает, если отрезок разделить на несколько частей и применить формулу трапеций к каждому отрезку (рис. 3). Тогда
Рис.2 Рис. 3
Для простоты вычислений удобно делить отрезок на равные части, в этом случае длина каждого из отрезков разбиения есть . Тогда, численное значение интеграла на всем отрезке равно
Эта формула называется общей формулой трапеций. Общую формулу трапеций можно переписать в более удобном виде:
, где шаг (6)
Пример 2. Вычислить интеграл с помощью формулы трапеций при .
Решение: составим таблицу значений подынтегральной функции при и 
:
| 0,2 0,4 0,6 | 0,02 0,16 0,36 | 0,0000 0,0400 0,1593 0,3523 | 0, 1,0 1,2 1,4 1,6 | 0,64 1,0 1,44 1,96 2,56 | 0,5972 0,8415 0,9915 0,9249 0,5487 |
Используя формулу ,
Находим:
Примечание. Если данный интеграл вычислить при , то получим . Следовательно, точность вычислений увеличивается с возрастанием .
Формула Симпсона
Точность приближенного интегрирования заметно возрастает, если подынтегральную функцию на отрезке заменить квадратичной функцией (рис.5), принимающей в узлах х 0 = а, х 1 , х 2 = b значения и . В качестве интерполяционного многочлена используется многочлен Ньютона 2 степени. Тогда
(7)
Соотношение (7) называется формулой Симпсона. Формула Симпсона обладает повышенной точностью и является точной не только для многочленов второй степени, но и третьей. Погрешность формулы Симпсона оценивается следующим образом:
,
где точка 
(8)
Для увеличения точности вычислений отрезок разбивают на п
пар участков 
(рис. 4) и к каждому из них применяют формулу (7). Тогда численное значение определенного интеграла на всем отрезке будет равно
, где (9)
Соотношение (9) называется общей формулой Симпсона
.
Пример 3.
Вычислить по формуле Симпсона при .
По формуле (9) имеем 
. Подставляя в подынтегральную функцию значения , получим
Приближенное решение обыкновенных дифференциальных уравнений. Метод Эйлера
Простейшим обыкновенным дифференциальным уравнением является уравнение первого порядка, разрешенное относительно производной:
Основная задача, связанная с этим уравнением, известна как задана Коши: найти решение уравнения (10) в виде функции у(х), удовлетворяющей начальному условию
Геометрически это означает, что требуется найти интегральную кривую у = у(х), проходящую через
заданную точку M 0 (x 0 , .y 0), при выполнении равенства (11) (см. рис.6). С численной точки зрения задача Коши выглядит следующим образом: требуется построить таблицу значений функции у=у(х), удовлетворяющей уравнению (10) и начальному условию (11) на отрезке [a;b ]с некоторым шагом h. Обычно считается, что х 0 = а, т.е. начальное условие задано в левом конце отрезка.
Простейшим из численных методов решения дифференциальных уравнений является метод Эйлера. В основе метода Эйлера лежит идея графического построения решения дифференциального уравнения, однако этот метод дает одновременно и способ нахождения искомой функции в численной (табличной) форме.
Пусть дано уравнение (10) с начальным условием (11) (т.е. поставлена задача Коши). Решим вначале следующую задачу: найти простейшим способом приближенное значение решения в некоторой точке x 1 = х 0 + h , где h - достаточно малый шаг.
Заметим, что уравнение (10) совместно с начальным условием (11) задают направление касательной к искомой интегральной кривой в точке М 0 (х 0 , у 0). Уравнение касательной имеет вид
Двигаясь вдоль этой касательной (рис. 7), учитывая соотношения (10) и (12), получим приближенное значение решения в точке х 1 :
(13)
Располагая приближенным решением в точке М 1 (х 1 ,y 1), можно повторить описанную ранее процедуру: построить прямую, проходящую через эту точку с угловым коэффициентом f (х 1 , y 1) и по ней найти приближенное значение решения в точке х 2 = х 1 + h . Заметим, что в отличие от ситуации, изображенной на рис. 7, эта прямая не есть касательная к реальной интегральной кривой, поскольку точка M 1 , нам недоступна. Однако представляется интуитивно ясным, что если h достаточно мало, то получаемые приближения будут близки к точным значениям решения.
Продолжая эту идею, построим систему равноотстоящих точек
х i = x 0 + ih (i = 0, 1, 2, ..., n) (14)
Получение таблицы значений искомой функции у(х) по методу Эйлера заключается в циклическом
применении формулы
Геометрическая иллюстрация метода Эйлера приведена на рис. 8. Вместо интегральной кривой в реальности получается совокупность прямых (так называемая ломаная Эйлера ).
Рис.8 Ломаная Эйлера
Методы численного интегрирования дифференциальных уравнений, в которых решение получается от одного узла к другому, называются пошаговыми. Метод Эйлера - простейший представитель семейства пошаговых методов.
Отметим, что оценка погрешности метода при таком элементарном рассмотрении невозможна даже на первом шаге. Кроме того, особенностью любого пошагового метода является то, что, начиная со второго шага исходное значение у, в формуле (13) само является приближенным, т.е. погрешность на каждом следующем шаге систематически возрастает.
Наиболее используемым эмпирическим методом оценки точности как метода Эйлера, так и других пошаговых методов приближенного численного интегрирования обыкновенных дифференциальных уравнений является способ двойного прохождения заданного отрезка - с шагом h и с шагом h/2 . Совпадение соответствующих десятичных знаков в полученных двумя способами результатах дает эмпирические основание считать их верными (хотя полной уверенности в этом быть не может).
Одна из принципиальных трудностей всех пошаговых методов численного решения дифференциальных уравнений состоит в возможности столкнуться с неустойчивостью метода. Оценка погрешности неявно предполагает, что ломаная приближенного решения (см. рис. 8) хотя и не совпадает с интегральной кривой, но качественно на нее похожа. Чаще всего это именно так, но иногда (например, при неудачном выборе шага h) приближенное решение может быть качественно непохожим на точное (например, точное монотонно убывает, а приближенное монотонно возрастает).
Для эмпирического контроля того, не имеет ли места неустойчивость, следует численно интегрировать уравнение с несколькими, значительно отличающимися, значениями шага h , сравнивая качественно поведение решений.
Пример 4. Применяя метод Эйлера, составить таблицу значений решения дифференциального уравнения , с начальным условием на отрезке , приняв h= 0,25. Вычисления проводить с 4-мя знаками после запятой.
Для удобства вычислений составим таблицу.
1-й шаг: по начальным условиям заполним первую строку во 2-м и 3-м столбцах;
2-й шаг: из уравнения вычисляем (i = 0, 1, 2, 3, 4, 5) в столбце (4);
3-й шаг: содержимое столбца (4) умножаем на h
(вычисляем 
) и
записываем результат в столбец (5) этой же строки;
4-й шаг: к содержимом столбца (3) прибавляем содержимое столбца (5) этой же строки
(вычисляеми результат записывает столбец (3)следующей
строки. Определяем х i +1 = x i + h и затем шаги 2-4 повторяем до тех пор, пока не будет пройден
весь отрезок .
| i | x i | y i | ||
| (1) | (2) | (3) | (4) | (5) |
| 1,5000 | 1,5000 | 0,3750 | ||
| 0,25 | 1,8750 | 1,6250 | 0,4062 | |
| 0,50 | 2,2812 | 1,7812 | 0,4453 | |
| 0,75 | 2,7265 | 1,9765 | 0,4951 | |
| 1,00 | 3,2206 | 2,2206 | 0,5552 | |
| 1,25 | 3,7758 | 2,5258 | 0,6314 | |
| 1,50 | 4,4072 |
Пример 5. Решить методом Эйлера дифференциальное уравнение с начальным значением у (0) = 1,3 на отрезке , приняв шаг h = 0,2.
знаками: .
с начальным условием у
(2) = 1, 2 на отрезке , приняв шаг h
= 0,1.
Вариант 2
1.По формуле левых прямоугольников для n=12 вычислить значение интеграла: .
2.По формуле трапеций n=8 вычислить значение интеграла с тремя десятичными
знаками: .
3.По формуле Симпсона для n=6 вычислить значение интеграла с тремя десятичными
знаками: .
4. Применяя методом Эйлера, составить таблицу значений решения дифференциального уравнения
с начальным условием у
(2,6) = 1, 8 на отрезке , приняв шаг h
= 0,2.
Сделайте ручную прикидку графика интегральной кривой на бумаге.
Вариант 3
2.По формуле трапеций для n=10 вычислить значение интеграла с тремя десятичными
знаками:
знаками:
4. Применяя методом Эйлера, составить таблицу значений решения дифференциального уравнения
с начальным условием у
(0,6) = 3,4 на отрезке , приняв шаг h
= 0,2.
Сделайте ручную прикидку графика интегральной кривой на бумаге.
Вариант 4
1.По формуле правых прямоугольников для n=12 вычислить значение интеграла:
2. По формуле трапеций для n=10 вычислить значение интеграла с тремя десятичными знаками:
3. По формуле Симпсона для n=6 вычислить значение интеграла с тремя десятичными
знаками:
4. Применяя методом Эйлера, составить таблицу значений решения дифференциального уравнения
с начальным условием у
(3) = 1,7 на отрезке , приняв шаг h
= 0,2.
Сделайте ручную прикидку графика интегральной кривой на бумаге.
4. Контрольные вопросы:
1. Какие методы приближенного вычисления определенных интегралов вы знаете? Назовите
формулы для вычислений. Какой из них дает наиболее точный результат?
2. На чем основан метод Эйлера приближенно решения дифференциальных уравнений?
5.1 Наименование работы
5.2 Цель работы
5.3 Задание
5.4 Формулы для расчета
5.5 Необходимые расчеты. Анализ результатов расчетов
5.6 Выводы по работе
5.7 Ответы на контрольные вопросы
Литература:
1. Колягин Ю.М. , Луканкин Г.Л., Яковлев Г.Н. Математика в 2-х томах Учебное пособие - М.
Новая волна, 2005, ч.1, с.565-571;
2. Богомолов Н.В. «Практические занятия по математике» - Учебное пособие – М.:Высш. школа,
2003, с. 211-212;
3. Лапчик М.П., Рагулина М.И., Хеннер Е.К. Элементы численных методов: учебник для студ. сред.
проф. образования -М.: Издательский центр «Академия», 2007, с.152-184
Похожая информация.
Рассмотрим задачу Коши (5.2), (5.6) для дифференциального уравнения первого порядка: найти решение уравнения y"=f(x,y), удовлетворяющее условию y(x 0)=y 0 . Пусть y(x)- решение поставленной задачи Коши. Подставив это решение в уравнение (5.2), получим тождество y"(x) ≡ f(x,y(x)). Интегрируя это тождество по x, получаем
или, что тоже самое,
. (5.15)
Таким образом, мы показали, что всякое решение задачи Коши (5.2), (5.6) есть решение интегрального уравнения (5.15). С другой стороны, если y(x)- решение интегрального уравнения (5.15), то дифференцируя (5.15) по x, получаем, что y(x)- решение задачи Коши (5.2), (5.6).
Решение интегрального уравнения (5.15) будем искать с помощью метода последовательных приближений. Положим
y 0 (x)=y 0 ,
. (5.16)
Если оператор
- (5.17)
сжимающий , то последовательные приближения (5.16) сходятся к решению интегрального уравнения (5.15), а, следовательно и дифференциального уравнения y" = f(x,y), удовлетворяющего условию y(x 0) = y 0 . Желающие могут познакомиться с доказательством сжимаемости оператора (5.17) в .
Пример №1 . Найдём с помощью метода последовательных приближений решение уравнения y" = y, удовлетворяющее условию y(0)=1. Подставляя y(0)=1 в (5.16), получаем
y 0 =1,
…,
С другой стороны, решая исходную задачу Коши, имеем y = e x .
Таким образом, нами получено разложение функции e x в ряд Тейлора в нуле (ряд Маклорена).
Перейдём теперь к изложению численного метода Эйлера решения задачи Коши (5.2), (5.6). Разобьём отрезок , на котором мы ищем решение, на части точками x 0 = a
в уравнении (5.2), получаем 
, или, что то же самое,
y i +1 = y i + h·f(x i , y i), (5.17)
Соотношение (5.17) является расчётной формулой метода Эйлера численного решения задачи Коши (5.2), (5.6). Вычислив y i , i = 0,1,..,n получим таблицу значений решения в точках x i , i = 0,1,..,n Для оценки погрешности на одном шаге сетки в методе Эйлера разложим точное решение y(x) по формуле Тейлора в окрестности точки xi до членов второго порядка малости
y(x i +1)=y(x i +h)=y(x i)+y"(x i)h+o(h 2)=y i +hf(x i ,y i)+o(h 2).
Сравнивая с (5.17) видим, что погрешность формулы (5.17) равна o(h 2). К сожалению, метод Эйлера накапливает ошибку от шага к шагу. Поэтому на практике пользуются либо модификациями метода Эйлера, например методом прогноза и коррекции , либо другими методами, в частности методом Рунге-Кутта .
Курсовая работа
по дисциплине: Высшая математика
на тему: Приближенные решения дифференциальных уравнений
Введение
При анализе режимов работы электроэнергетических объектов и разработке новых технологических процессов инженеру часто приходится сталкиваться с дифференциальными уравнениями, т.к. большая часть законов электротехники и теплотехники формулируется в виде дифференциальных уравнений. При этом нередко приходится иметь дело с уравнениями, общее решение которых не выражается в квадратурах. Например, общее решение очень простого уравнения нельзя записать в конечном виде через элементарные функции. Класс задач, для которых можно найти явное решение, весьма узок. В связи с интенсивным применением дифференциальных уравнений в качестве математических моделей широкого круга естественнонаучных задач и с появлением высокопроизводительных ЭВМ важное значение приобрели численные методы их решения. Численные методы - это алгоритмы вычисления приближенных значений искомого решения в точках конечного множества значений аргумента (узлах сетки). Решение при этом получается в виде таблицы. Рассмотрим два таких метода: метод Рунге-Кутта и вытекающий из него метод Эйлера.
Метод Рунге-Кутта. Рассматриваем задачу Коши для дифференциального уравнения первого порядка
решение которой находится на отрезке , h > 0. Считаем, что задача (1) имеет единственное решение у(х), определяемое на этом отрезке. Выбираем на отрезке [х0, x0 + H] сетку значений аргумента хд= х0 + nh , n=0,1,...,N; h= h/n. Разложим решение у (х) в ряд Тейлора в окрестности точки хn, полагая при этом уn=у(хn), у`= y`(xn) и т.д.
Подставим в разложение (2) значение х = xn+1, получая равенство
Стоящие в правой части равенства (3) производные можно найти, последовательно дифференцируя уравнение (1):
с учетом формул (5) равенство (3) можно записать последовательно в виде
Пренебрегая в правых частях формул (6) слагаемыми О(h2),О(h3), малыми при малых h, получаем соответственно формулы
Каждая из формул (7),(8),... позволяет по известному значению y0 решения задачи (1) в начальной точке xо последовательно вычислять приближенные значения этого решения в узлах сетки х1,х0,...,xn; в отличие от точных значений обозначим их
Формулу (8) и тем более, формулы с большим числом членов в практических расчетах не используют, так как если функция f(х,у) правой части имеет несложное выражение, то выражения (4) для ее производных могут оказаться громоздкими. Если функция f(х,у) известна лишь приближенно, то процесс вычислений по этим формулам усложняется еще и из-за необходимости использовать формулы численного, дифференцирования. Расчет приближенных значений уn задачи Коши (1) по формуле (7) называется методом Эйлера, или схемой ломаных. Геометрическая интерпретация этой схемы дана на рис.1, где изображено поле интегральных кривых.
Размещено на http://www.сайт/
При удалении от точки (хо,уо) ломаная Эйлера может заметно отклоняться от графика точного решения. Известна следующая оценка погрешности метода Эйлера. Пусть в D={(х,у): |х-хo |<а; |у-уо| Где С1
=(1+М) (еКН-1) При отсутствии ошибок округлений локальная погрешность метода Эйлера, т.е. погрешность на одном шаге h , возникающая за счет перемещения по касательной к интегральной кривой, проходящей через точку (хn,уn), а не по самой интегральной кривой, есть величина O(h2). Глобальная погрешность или, точнее, максимальная погрешность решения на сетке {х1,x2,...,хN} в целом равна O(h), как следует из неравенства (9). В связи с этим говорят, что метод Эйлера имеет первый порядок точности. С другой стороны в учебнике «Краткий курс математического анализа» А.Ф.Бермант, И.Г.Араманович, М. 1973 г этот метод описывается иначе. Известно, что уравнение задает в некоторой области поле направлений. Решение этого уравнения с некоторыми начальными условиями дает кривую, которая касается поля направлений в любой точке. Если взять последовательность точек х0, х1, х2,
…. и заменить на получившихся отрезках интегральную кривую на отрезки касательных к ней, то получим ломаную линию. При подстановке заданных начальных условий (х0, у0
) в дифференциальное уравнение получаем угловой коэффициент касательной к интегральной кривой в начальной точке Заменив на отрезке интегральную кривую на касательную к ней, получаем значение Производя аналогичную операцию для отрезка , получаем: Если последовательность точек хi выбрать так, чтобы они отстояли друг от друга на одинаковое расстояние h, называемое шагом вычисления, то получаем формулу: Следует отметить, что точность метода Эйлера относительно невысока. Увеличить точность можно, конечно, уменьшив шаг вычислений, однако, это приведет к усложнению расчетов. Поэтому на практике применяется так называемый уточненный метод Эйлера или формула пересчета. Суть метода состоит в том, что в формуле вместо значения берется среднее арифметическое значений f
(x
0,
y
0)
и f
(x
1,
y
1)
. Тогда уточненное значение: Затем находится значение производной в точке. Заменяя f
(x
0,
y
0)
средним арифметическим значений f
(x
0,
y
0)
и
, находят второе уточненное значение у1
. Затем третье: и т.д. пока два последовательных уточненных значения не совпадут в пределах заданной степени точности. Тогда это значение принимается за ординату точки М1 ломаной Эйлера. Аналогичная операция производится для остальных значений у
. Подобное уточнение позволяет существенно повысить точность результата. Одним из методов, позволяющих для решения задачи Коши (1) строить вычислительные схемы более высоких порядков точности, является метод, предложенный Рунге и усовершенствованный Кутта и другими математиками. Схемы метода Рунге-Кутта удобны как, для расчетов на ЭВМ, так и для ручных расчетов. С основной идеей метода ознакомимся на примере построения вычислительных схем второго порядка точности. Воспользуемся теперь второй из формул (6) Покажем, что можно правильно передать члены ряда Тейлора, указанные в формуле (10), избежав дифференцирования, функция f(x,y), С этой целью полагаем где -некоторые постоянные По формуле Тейлора первого порядка находим Подставляя это выражение для значения в равенство (11) получаем Выбираем параметры так, чтобы правые части разложений (10) и (12) совпадали с точностью до слагаемых порядка О(h3). Для этого достаточно положить: Эта система трех уравнений с четырьмя неизвестными имеет бесконечно много решении. Выразим через остальные параметры: и подставим их в формулу (12), пренебрегая при этом слагаемыми О(h3). B результате получаем однопараметрическое семейство двучленных схем Локальная погрешность формулы (13) равна О(h3). Для максимальной погрешности на сетке выполняется оценка где С2 - некоторая постоянная, не зависящая от h. Вычислительная схема расчета приближенных значений решения задачи Коши (1) по формуле (13) называется схемой Рунге-Кутта второго порядка точности. Эти схемы нередко используются в практических вычислениях. При этом полагают, либо, либо. В первом случае получается схема особенно простого вида Для схема (13) имеет вид Методом Рунге-Кутта можно строить схемы различного порядка точности. Например, схема ломаных (7) есть схема Рунге-Кутта первого порядка точности, двучленные схемы (13) имеют второй порядок точности. Наибольшее распространение получили схемы четвертого порядка точности, при построении которых в записи ряда Тейлора (2) удерживаются все члены, включая h4. Приведем без вывода ту из них, которая записана в большинстве стандартных программ для ЭВМ: Для схемы (16) выполняется следующая оценка погрешности: если в прямоугольнике d существуют непрерывные частные производные четвертого порядка функции f
(x,у), то Схемы Рунге-Кутта более высокого порядка точности практически не употребляются, так как расчетные формулы становятся слишком громоздкими. Одним из важных достоинств метода Рунге-Кутта является простота алгоритма вычислений. Для начала вычислений достаточно выбрать сетку {хо, х1, ..., xN} и задать начальное значение у (хо) = уо. Далее вычисления производят последовательно по одним и тем же формулам. Это свойство схем метода Рунге-Кутта очень ценно при расчетах на ЭВМ, программирование расчетных формул метода не представляет труда. Выбор шага. Апостериорная оценка погрешности. Правило Рунге
Правильный выбор шага сетки h
является одним из главных практических вопросов, которые возникают при численном решении дифференциальных уравнений. Как получить требуемую точность? Если шаг выбран слишком большим, то значительной будет локальная погрешность и накопившаяся глобальная погрешность может быть недопустимо большой. Если же шаг слишком мал, то расчет потребует неоправданно много времени работы вычислителя или ЭВМ. При этом наблюдается следующий отрицательный эффект: суммарное действие ошибок округления, малых при выполнении каждой операции, может оказаться столь значительным, что полученный ответ становится бесполезным. Априорные оценки типа (9) мало полезны для получения информации о точности вычислений ввиду их сложности (особенно для схем высокого порядка); при этом они, как правило, во много раз превосходят фактическую ошибку расчета. Основным практическим приемом является апостериорная оценка погрешности. Для ее получения проводят вычисления на двух или более сгущающихся сетках и применяют так называемое правило Рунге, которое заключается в следующем. Пусть обозначает приближенное значение решения задача Коши (1) в точке х = x(h) = xo + nh, вычисленное по некоторой схеме Рунге-Кутта p-го порядка точности, а приближенное значение для y(x), вычисленное по той же схеме и в той же
точке, являющейся узлом более густой сетки с шагом h/2, так что. При некоторых предположениях относительно гладкости функции правой части f(x,y) погрешность схемы Рунге-Кутта р-го порядка точности имеет вид где С зависит от точки х, но не от h. Применяя формулу (17) для оценки погрешности на сетке с шагом h/2 , получаем Сохраняем в равенствах (17) и (18) только главную часть погрешности, пренебрегая слагаемыми и, и вычитаем почленно из первого равенства второе. В итоге приходим к приближенному равенству откуда определяем, что погрешность на сетке с меньшимшагом составляет частности, оценка (19) имеет вид: для схемы ломаных (7) для схемы (13) для схемы (16) Информацию о погрешности вычислений в виде (19) можно использовать для того, чтобы уточнить приближенные значения на сетке с меньшим шагом, внося в них поправки следующим образом. В точке, являющейся общим узлом двух сеток, полагаем в соответствии о формулами (18) и (19): дифференциальный уравнение точность погрешность Значения поправок, в узлах с нечетными номерами m=2n-1 находим, применяя линейную интерполяцию; Можно ожидать, что исправленные значения будут более точными, чем. Пример практического расчета
Используя метод Эйлера, а затем схему Рунге-Кутта второго порядка точности составить таблицу приближенных значений решения задачи Коши на отрезке с шагом h=0,2 и h=0,1. Результат округлить до 10-4. Оценить погрешность на сетке с шагом h=0,1 по методу Рунге. Сравнить полученные результаты с результатами вычислений по схеме четвертого порядка точности на сетке с шагом h=0,1. В таблицах 2 и 3: По методу Эйлера (табл.1) Таблица 1 h
=0,2 Таблица 2 Таблица приближенных значений по схеме Рунге-Кутта второго порядка точности шаг h
=0,1 Таблица 3 Оценим погрешность на сетке с шагом h
=0,1 по методу Рунге. В таблице 4: Таблица 4 Таблица приближенных значений по методу Рунге-Кутта четвертого порядка точности h
=0,1. Таблица 5 По результатам таблиц построим 4 графика yn
=
f
(xn
)
Как видно из графика более точные значения решения задачи Коши (полученные по схеме Рунге-Кутта второго и четвертого порядка точности с шагом h
=0,1 графики 2 и 3) лежат между менее точными значениями (полученные методом Эйлера и по схеме Рунге-Кутта второго порядка точности с шагом h
=0,2 графики 1 и 4) Литература
1. Методическое указание к лабораторной работе «Метод сеток. II. Решение задачи Коши для обыкновенного дифференциального уравнения первого порядка». Архангельск 1985 2. «Краткий курс математического анализа» А.Ф.Бермант, И.Г.Араманович, М. 1973 г. Хемминг Р.В., «Численные методы», «Наука». Основные методы Рунге-Кутта: построение класса расчетных формул. Расчетная формула метода Эйлера. Получение различных методов Рунге-Кутта с погрешностью второго порядка малости при произвольном задавании параметров. Особенности повышения порядка точности. реферат , добавлен 18.04.2015 Общая характеристика и особенности двух методов решения обычных дифференциальных уравнений – Эйлера первого порядка точности и Рунге-Кутта четвёртого порядка точности. Листинг программы для решения обычного дифференциального уравнения в Visual Basic. курсовая работа , добавлен 04.06.2010 Изучение методов Рунге-Кутты четвертого порядка с автоматическим выбором длины шага интегрирования для решения дифференциальных уравнений. Оценка погрешности и сходимость методов, оптимальный выбор шага. Листинг программы для ЭВМ, результаты, иллюстрации. курсовая работа , добавлен 14.09.2010 Численное решение уравнения методом Эйлера и Рунге-Кутта в Excel. Программа на языке Turbo Pascal. Блок-схема алгоритма. Метод Рунге-Кутта для дифференциального уравнения второго порядка. Модель типа "хищник-жертва" с учетом внутривидового взаимодействия. курсовая работа , добавлен 01.03.2012 Решение задачи Коши для дифференциального уравнения. Погрешность приближенных решений. Функция, реализующая явный метод Эйлера. Вычисление погрешности по правилу Рунге. Решение дифференциальных уравнений второго порядка. Условие устойчивости для матрицы. контрольная работа , добавлен 13.06.2012 Составление диагональной системы способом прогонки, нахождение решения задачи Коши для дифференциального уравнения на сетке методом Эйлера и классическим методом Рунге-Кутта. Построение кубического сплайна интерполирующей функции равномерного разбиения. практическая работа , добавлен 06.06.2011 Аналитическое и компьютерное исследования уравнения и модели Ван-дер-Поля. Сущность и особенности применения методов Эйлера и Рунге-Кутта 4 порядка. Сравнение точности метода Эйлера и Рунге-Кутта на одном графике, рисуя фазовые траектории из 1 точки. курсовая работа , добавлен 06.10.2012 Понятие, закономерности формирования и решения дифференциальных уравнений. Теорема о существовании и единственности решения задачи Коши. Существующие подходы и методы решения данной задачи, оценка погрешности полученных значений. Листинг программы. курсовая работа , добавлен 27.01.2014 Задачи Коши и методы их решения. Общие понятия, сходимость явных способов типа Рунге-Кутты, практическая оценка погрешности приближенного решения. Автоматический выбор шага интегрирования, анализ брюсселятора и метод Зонневельда для его расчета. курсовая работа , добавлен 03.11.2011 Дифференциальное уравнение первого порядка, разрешенное относительно производной. Применение рекуррентного соотношения. Техника применения метода Эйлера для численного решения уравнения первого порядка. Численные методы, пригодные для решения задачи Коши. Продолжая изучать практические приложения теории рядов, рассмотрим ещё одну распространённую задачу, название которой вы видите в заголовке. И, чтобы не чувствовать себя газонокосилкой на протяжении урока, давайте сразу же разберёмся в сути задания. Три вопроса и три ответа: Что нужно найти?
Частное решение дифференциального уравнения
. Намёк между строк шепчет, что к данному моменту желательно хотя бы понимать, что такое дифференциальное уравнение
и что такое его решение. КАК по условию требуется это решение?
Приближённо – с помощью ряда
. И третий закономерный вопрос: почему приближённо?
Этот вопрос я уже освещал на уроке Методы Эйлера и Рунге-Кутты
, однако повторение не помешает. Будучи сторонником конкретики, вернусь к простейшему дифференциальному уравнению
. В ходе первой лекции по диффурам мы нашли его общее решение (множество экспонент) и частное решение , соответствующее начальному условию . График функции – это самая обычная линия, которую нетрудно изобразить на чертеже. Но то элементарный случай. На практике встречается великое множество дифференциальных уравнений, неразрешимых аналитически точно (по крайне мере, известными на сегодняшний день способами). Иными словами, как ни крути такое уравнение – проинтегрировать его не удастся. А закавыка состоит в том, что общее решение (семейство линий на плоскости) может существовать
. И тогда на помощь приходят методы вычислительной математики. Встречаем нашу радость! Типовая задача формулируется следующим образом: …
, удовлетворяющее начальному условию , в виде трёх (реже – четырёх-пяти)
отличных от нуля членов ряда Тейлора
. Искомое частное решение раскладывается в данный ряд по известной формуле: Единственное, здесь вместо буквы «эф» используется «игрек» (так уж повелось). Идея и смысл тоже знакомы
: для некоторых диффуров и при некоторых условиях (не будем вдаваться в теорию) построенный степенной ряд будет сходиться
к искомому частному решению . То есть, чем больше членов ряда мы рассмотрим, тем точнее график соответствующего многочлена приблизит график функции . Следует отметить, что вышесказанное применимо и к самым простым случаям. Проведём незамысловатое детское исследование на том же горшке: Пример 1
Найти приближённо частное решение дифференциального уравнения , удовлетворяющее начальному условию в виде четырёх первых отличных от нуля членовряда Тейлора. Решение
: в условиях данной задачи , поэтому общая формула Тейлора трансформируется в частный случай разложения в ряд Маклорена
: Немного забегая вперёд, скажу, что в практических заданиях значительно чаще встречается именно этот, более компактный ряд. Занесите обе рабочие формулы в свой справочник. Разбираемся со значениями . Этапы решения удобно занумеровать: 0) На нулевом шаге записываем значение , которое всегда известно из условия. В тетради итоговые результаты пунктов желательно обводить в кружок, чтобы они были хорошо видны и не затерялись в решении. Мне по техническим причинам сподручнее выделять их жирным шрифтом. Кроме того, отмечаем, что данное значение не равно нулю
! Ведь по условию требуется найти четыре отличных от нуля
членовряда. 1) Вычислим . Для этого в правую часть исходного уравнения вместо «игрека» подставляем известное значение : 2) Вычислим . Сначала находим вторую производную
: Подставляем в правую часть найдённое в предыдущем пункте значение : В распоряжении уже три ненулевых члена разложения, нужен ещё один: Пример 2
Найти приближённо частное решение дифференциального уравнения Решение
начинается стандартной фразой: В данной задаче , следовательно: Теперь последовательно находим значения – до тех пор, пока не будут получены три ненулевых
результата. Если повезёт, то отличны от нуля будут Нарезаем пункты решения: 0) По условию . Вот и первый успех. 1) Вычислим . Сначала разрешим исходное уравнение относительно первой производной, то есть, выразим Получена баранка и это не есть хорошо, поскольку нас интересуют ненулевые
значения. Однако ноль – тоже результат
, который не забываем обвести в кружок или выделить каким-нибудь другим способом. 2) Находим вторую производную и подставляем в правую часть известные значения : Второй «не ноль». 3) Находим – производную от второй производной: Вообще, задание чем-то напоминает Сказку про Репку, когда дедка, бабка и внучка зовут на помощь жучку, кошку и т.д. И в самом деле, каждая следующая производная выражается через своих «предшественников». Подставим в правую часть известные значения : Третье ненулевое значение. Вытащили Репку. Аккуратно и внимательно подставляем «жирные» числа в нашу формулу: Ответ
: искомое приближенное разложение частного решения: В рассмотренном примере попался всего один ноль на втором месте, и это не так уж плохо. В общем случае нулей может встретиться сколько угодно и где угодно. Повторюсь, их очень важно выделять наряду с ненулевыми результатами, чтобы не запутаться в подстановках на завершающем этапе. Вот, пожалуйста – бублик на самом первом месте: Пример 3
Найти приближённо частное решение дифференциального уравнения , соответствующее начальному условию , в виде трёх первых отличных от нуля членовряда Тейлора. Примерный образец оформления задачи в конце урока. Пункты алгоритма можно и не нумеровать (оставляя, например, пустые строки между шагами), но начинающим рекомендую придерживаться строгого шаблона. Рассматриваемая задача требует повышенного внимания – если допустить ошибку на каком-либо шаге, то всё остальное тоже будет неверным! Поэтому ваша ясная голова должна работать как часы. Увы, это не интегралы
или диффуры
, которые надёжно решаются и в утомлённом состоянии, поскольку позволяют выполнить эффективную проверку. На практике заметно чаще встречается разложение в ряд Маклорена
: Пример 4
Решение
: в принципе, можно сразу записать разложение Маклорена
, но оформление задачи академичнее начать с общего случая: Разложение частного решения дифференциального уравнения при начальном условии имеет вид: В данном случае , следовательно: 0) По условию . Ну что поделать…. Будем надеяться, что нулей встретится поменьше. 1) Вычислим . Первая производная уже готова к употреблению. Подставим значения : 2) Найдём вторую производную: И подставим в неё : Резво дело пошло! 3) Находим . Распишу очень подробно: Заметьте, что к производным применимы обычные алгебраические правила: приведение подобных слагаемых на последнем шаге и запись произведения в виде степени: (там же). Подставим в всё, что нажито непосильным трудом : Три ненулевых значения рождены. Подставляем «жирные» числа в формулу Маклорена, получая тем самым приближенное разложение частного решения: Ответ
: Для самостоятельного решения: Пример 5
Представить приближенно частное решение ДУ, соответствующее заданному начальному условию , в виде суммы трех первых отличных от нуля членов степенного ряда. Примерный образец оформления в конце урока. Как видите, задача с частным разложением в ряд Маклорена
оказалась даже труднее общего случая. Сложность рассматриваемого задания, как мы только что убедились, состоит не столько в самом разложении, сколько в трудностях дифференцирования. Более того, порой, приходится находить 5-6 производных (а то и больше), что повышает риск ошибки. И в завершении урока предлагаю пару задач повышенной сложности: Пример 6
Решить дифференциальное уравнение приближённо с помощью разложения частного решения в ряд Маклорена, ограничившись тремя первыми ненулевыми членами ряда Решение
:перед нами диффур второго порядка, но это практически не меняет дела. По условию и нам сразу же предложено воспользоваться рядом Маклорена, чем мы не преминем воспользоваться. Запишем знакомое разложение, прихватив на всякий пожарный побольше слагаемых: Алгоритм работает точно так же: 0) – по условию. 1) – по условию. 2) Разрешим исходное уравнение относительно второй производной: . И подставим : Первое ненулевое значение Щёлкаем производные и выполняем подстановки: Подставим и : Подставим : Второе ненулевое значение. 5) – по ходу дела приводим подобные производные. Подставим : Подставим : Наконец-то. Впрочем, бывает и хуже. Таким образом, приближенное разложение искомого частного решения:
Подобные документы
Как найти частное решение ДУ приближённо с помощью ряда?
, удовлетворяющее начальному условию в виде трёх первых отличных от нуля членовряда Тейлора.
– это идеальный случай с минимальным количеством работы.
. Подставим в правую часть известные значения :
