Файл: Ермолаева Э.Н. Элементы численного анализа учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 07.07.2024
Просмотров: 134
Скачиваний: 0
Вычислить с точностью до 0,01 следующие определенные интегралы:
1 |
dx |
17. . |
1 + х |
Отв.:о 0,69.
1dx
18.J 1 4- XІ
Отв.:о 0,79
19. |
і |
|
|
dx |
|
||
J |
|
1 + |
xs |
|
|||
|
|
|
|||||
о |
|
|
|
|
|
|
|
Отв.: |
0,84. |
|
|
|
|||
|
2 |
|
х |
|
|
л- dx . |
|
20. |
j" |
l g |
|||||
|
і |
|
|
|
|
|
|
Отв.; |
0,28. |
|
|
|
|||
21. |
Jf |
|
|
x |
**-dx. |
||
|
і |
|
|
|
|
|
|
Отв.: |
0,10. |
|
|
|
|||
22. |
J1" |
|
S |
[ |
x |
x |
dx. |
|
|
|
n |
|
|||
|
0 |
|
|
|
|
|
|
Отв.: |
1,61. . |
|
|
||||
23. |
( |
|
- |
^ |
|
L ^ |
( i x _ |
|
J |
|
|
|
x |
|
|
|
0 |
|
1,85. |
|
|
|
|
Отв.: |
|
|
|
|
|||
0 . |
( 2 |
|
cos |
|
x |
, |
|
24. |
J |
|
|
|
x |
dx . |
|
|
і |
|
|
|
|
|
|
Отв.: |
0,09. |
|
|
|
|||
25. |
ГС |
|
|
|
dx . |
||
і* |
|
- |
—• |
|
|||
|
J |
|
1 |
+ |
* |
|
|
|
о |
|
0,67. |
|
|
|
|
Отв.: |
|
|
|
|
|||
|
і |
|
е~*г |
dx . |
|||
26. |
j |
' |
|||||
|
о |
|
|
|
|
|
|
Отв. |
0,75. |
|
|
|
|
Г л а в а 8
ПРИБЛИЖЕННОЕ ИНТЕГРИРОВАНИЕ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
Очень важно уметь применять основные методы прибли женного решения дифференциальных уравнений, так как даже Б случае достаточно простых обыкновенных дифференциаль ных уравнений не всегда удается выразить решение в квадра турах.
В зависимости от формы, в которой представляется прибли женное решение уравнения, методы приближенного интегриро вания дифференциальных уравнений можно разделить на три группы.
Аналитические методы предусматривают получение при ближенного частного решения уравнения либо в виде квадра тур, либо в виде степенных рядов, равномерно сходящихся, при определенных условиях, к точному решению. Среди ана литических методов основными являются:
—метод последовательных приближений;
—разложение решения в ряд по степеням независимой пе ременной или по степеням некоторых параметров уравнения;
—• разложение в ряд по степеням начальных данных. Графические методы предполагают приближенное построе
ние интегральных кривых.
Численные методы дают приближения для значений точ ного решения лишь в отдельных точках и их оформляют в виде таблицы.
Цель настоящей главы состоит в том, чтобы дать достаточ но полное изложение численных методов решения обыкновен ных дифференциальных уравнений первого порядка, и потому для иллюстрации аналитических методов представлено лишь схематическое изложение первых двух из них.
$ 8.1. М Е Т О Д П О С Л Е Д О В А Т Е Л Ь Н Ы Х П Р И Б Л И Ж Е Н И Й
Изложим этот метод применительно к дифференциальному
уравнению первого порядка |
|
У) |
(8-1) |
с начальным условием |
|
У І-ї=і-„ = Уо • |
(8.2) |
Предположим, что в некоторой окрестности течки М0(х0, |
уQ) |
уравнение (8.1) удовлетворяет условиям теоремы существов„а- шія и единственности решения *.
Интегрируя |
правую и левую части уравнения (8.1) |
в преде- |
|
- лах от х0 ДО х, |
получаем интегральное |
уравнение |
|
|
х |
|
|
|
У (*) = У, + J / (*, |
//) dx . |
(8.3) |
Решение уравнения (8.3) удовлетворяет дифференциально му уравнению (8.1) и начальным условиям (8.2).
Если интервал (х0, х) достаточно мал, то приближенно можно полагать в нем функцию у(х) равной у0, и в качестве первого приближения выбрать функцию
|
|
|
|
|
|
|
х |
|
|
|
|
|
|
|
|
|
Уі |
(х) = у<> |
+ |
\ |
f (х, yv) dx . |
|
|
||||||
Подставив |
затем |
в правую часть уравнения |
(8.3) |
функцию |
||||||||||
У І { Х ) , найдем второе приближение: |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
х |
|
|
|
|
|
|
|
|
|
у, |
(х) |
= |
у0 |
+ |
\ |
f |
{х, yt) |
dx . |
|
|
|
|
Все дальнейшие |
приближения |
строятся |
по |
формуле |
|
|||||||||
|
|
Уп (х) |
- |
|
|
X |
/ |
(*. уп-i) dx |
, |
|
(8.4) |
|||
|
|
Уо + |
J |
|
||||||||||
* В о б л а с т и \х—Хо |
\ < а, \у |
у0, |
<Ь |
с у щ е с т в у е т е д и н с т в е н н о е |
р е ш е н и е |
|||||||||
у р а в н е н и я |
(8 . )), |
у д о в л е т в о р я ю щ е е |
н а ч а л ь н ы м у с л о в и я м |
(8.2), |
если |
в этой |
||||||||
о б л а с т и ф у н к ц и я f ( x , у) у д о в л е т в о р я е т у с л о в и ю Л и п ш и ц а |
|
|
||||||||||||
|
|
/ С*. УО" - |
/ (-г, |
у 2 ) |
I |
, Л' | y L |
- |
у 2 ' |
, |
|
|
|||
где Л'' — |
п о с т о я н н а я |
Л и п ш и ц а |
— |
не |
|
з а в и с и т |
от |
х, у\, уі; |
в |
к а ч е с т в е |
||||
Л' м о ж н о |
в ы б р а т ь |
max |
; f'v |
(х, у) | . |
|
|
|
|
|
|
|
причем |
имеет |
место теорема, |
устанавливающая, |
что, если |
|||
функция |
J (х, у) |
в окрестности |
точки Л10 (х0 , уо) удовлетворяет |
||||
условию |
Липшица, то |
последовательность |
{у„ (х)} |
на |
некото-. |
||
ром достаточно |
малом |
отрезке |
[х0; xo + h] |
равномерно |
сходит |
||
ся к точному решению |
у(х): |
|
|
|
|
у (х) — lim у,г (х) .
Если f(x, у) определена и непрерывна в области
R {0 < |
х — х0 |
< |
а; |
|
\ у — {/„ | < |
Ь] , |
и |
|
|
|
|
|
|
max 1 / |
(х, у) |
| < |
М |
|
при (х, у) |
є R, |
то за величину h можно принять |
|
|
, |
|||
|
h — mm |
х, |
М |
|
||
|
|
|
|
|
|
Погрешность метода последовательных приближений опре деляют по формуле
з „ (х) = | у (х) - уп (х) | ; MN* - ^ П у Г '
Пример 8.1. Методом последовательных приближений най ти приближенное решение дифференциального уравнения у'~х—у, удовлетворяющее начальному условию у\х=о = Ь
Р е ш е н и е . Интегрируя уравнение от 0 до х, получаем
х
|
|
|
|
у — 1 -\- f (х — у) dx , |
|
|
|
|
|
||||||
|
|
|
|
I |
|
о |
|
|
|
|
|
|
|
|
|
откуда при ;/о =1 |
|
|
|
|
|
|
|
|
|
|
|||||
у1 |
= |
1 + |
\{у- |
\)dx=]-x+ |
- |
~ |
- |
|
|
( - 1 ) 2 - 2 Т |
|||||
|
|
|
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
х |
|
|
х- \ |
|
|
, |
.. |
, |
., |
|
х?> |
|
у, |
=r. |
1 + |
\ (х |
— 1 + |
х |
dx |
— |
— |
|||||||
-—- |
1 — |
х |
-f |
х2 |
|
||||||||||
|
|
|
|
|
|
у- |
|
|
|
х |
|
|
|
|
|
|
|
|
|
- x |
+ 2 - j T |
+ (- |
|
I ) 3 |
~у |
|
|
|
и аналогично
- Х - + Х - g Г" -24-
= \ - x - r 2 ( ^ r - - £ r ) + ( - 1 ) *
;/4 = 1 — А - ) - А:2 |
х3 |
XІ |
хь |
|
— |
12 |
120 |
|
|
|
|
|
||
2 ! |
З ! 1 |
4 ! / ' v |
' |
5! |
Оценим погрешность четвертого |
приближения: |
|
||
Функция f(x, у) =x—у |
непрерывна на всей плоскости XOY, |
поэтому за а и ft могут быть взяты любые положительные чис ла:
|
R |
{0 < |
х •< а; |
у |
, < |
ft] . |
|
При |
|
(х, |
у) |
є |
/? |
|
|||
/ |
(х, |
у) |
|
•< | х |
— у |
| |
< ! х | 4- |
| у ; |
< |
л |
4- |
ft -= М ; |
|
||||
h — mm |
[а, |
|
• . . |
— min |
а, |
|
—г- |
— |
|
|
-—г- ; |
|
|||||
|
|
\ |
|
М J |
|
х а 4- ft / |
а +- ft |
|
|||||||||
|
|
|
|
Л/ = max |
| / ' у |
(х, |
//) ! |
= |
1 ; |
|
|
|
|
||||
|
|
|
|
г, (х) |
< |
(а + |
ft) |
|
• |
|
|
|
|
|
|||
Для случая а — Ъ = 1, |
|
Л = |
™ |
получаем |
|
|
|
|
|||||||||
|
|
|
|
|
х 5 |
|
^ |
1 |
|
|
|
|
1 |
|
|
|
|
|
|
£ 4 W < |
б 0 |
|
< 6 0 . 32 |
|
|
|
1 9 2 и |
|
|
|
|||||
Методом |
математической |
индукции |
легко |
убедиться, |
что |
||||||||||||
.'/., |
1 • |
х |
• |
2 |
2 ! |
|
3 ! |
' ' • • + t |
|
^ |
Л і |
|
|||||
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
- М - 1 ) « - |
(л |
+ |
1 ) ! |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
г/(х) = |
lim уп |
= |
1 - |
X + |
2 [е - * — (1 |
- |
х)] |
= |
2е~х |
- (1 - |
х). |