Файл: Ермолаева Э.Н. Элементы численного анализа учеб. пособие.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)

Предположим, что в некоторой окрестности течки М00,

у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 -

х).