Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Численные методы. Ч.2. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

Лекция 19

Но в нашем распоряжении имеется только конечное число параметров, и выполнить бесконечное число условий (19.6) практически невозможно. Можно ожидать, что если n-достаточ-

но большое число, то функция

Fn (x) будет близкой к нулю:

Fn (x)= εn (x) , max

 

εn (x)

 

ε и ε > 0

достаточно малое число.

 

 

[a,b]

 

 

 

 

будет удовлетворять диффе-

В этом случае, функция yn (x)

ренциальному уравнению

 

 

 

 

F(x, yn , y'n, y''n ) = εn (x) ,

(19.7)

которое мало отличается от заданного уравнения (19.1). Можно ожидать, что для широкого класса функций F из усло-

вия (19.7) будет следовать достаточная малость на [a,b] разности

y(x) − yn (x) , где y(x) – решение граничной задачи (19.1)–(19.3). В методе Галеркина выбор базисных функций ϕk (x) можно

осуществлять следующим образом.

Так как функция ϕ0 (x) должна удовлетворять условиям

La (ϕ0 )= A ,

Lb (ϕ0 )= B ,

то ϕ0 (x) можно определить как многочлен первой степени:

ϕ0 (x) = α + β x

икоэффициенты α и β выбратькакрешениесистемыуравнений

0 (

 

)

 

1

 

 

α

α + β a

 

+ α

 

β = A,

β0 (α + β a) + β1

β = B.

 

 

 

 

 

 

 

Другие функции ϕk (x) можно,

например, когда α1 = β1 = 0 ,

α0 = β0 =1 вычислять по правилу:

 

 

 

ϕk (x) = (x a)k (x b) ,

k =1,2,...

или

 

 

 

 

 

 

ϕk (x) = (x a)(x b)k ,

k =1,2,...

Пример. Методом Галеркина найти приближенное аналити-

ческое решение уравнения

 

y′′ + y + x = 0 ,

(19.8)

удовлетворяющее краевым условиям:

 

- 101 -

ЧИСЛЕННЫЕ МЕТОДЫ

 

(

 

)

= 0

,

( )

= 0

.

(19.9)

y

 

0

 

 

y 1

 

 

Решение.

В качестве базисных функций ϕk (x) выберем следующие функции:

ϕ0 (x) = 0 , ϕ1 (x) = x(1 − x) ,

ϕ2 (x) = x2 (1 − x) .

Они линейно независимы и удовлетворяют краевым условиям. Приближенное решение задачи ищем в виде

y2 (x) = a1 x (1 − x ) + a2 x2 (1 − x) ;

y2 (x) = a1 x + (a2 a1 )x2 a2 x3

или

y2 (x) = a1 x (a1 a2 )x2 a2 x3 ;

y′′(x) = 2

(

a a

− 6a x

.

2

2

1 )

2

Подставляя в уравнение (19.8)

y2 (x) , y2′(x) , y2′′(x) получим

невязку:

2(a2 a1 ) − 6a2 x + a1 x + (a2 a1 )x2 a2 x3 + x =

= 2(a2 a1 ) + (a1 − 6a2 +1)x + (a2 a1 )x2 a2 x3 = F2 (x,a1,a2 )

Учитывая ортогональность функции F2 (x,a1,a2 ) к функциям ϕ1 (x) , ϕ2 (x) получим систему:

1

F2 (x,a1,a2 )(x x2 )dx = 0,

01

F2 (x,a1,a2 )(x2 x3 )dx = 0.

0

Подставляя в этусистемузначение F2 и вычисляяинтегралы, получим систему линейных алгебраических уравнений для определения a1 и a2 :

 

3

a +

3

a =

1

,

 

 

 

 

 

 

 

 

1

 

20

 

2

12

 

 

10

 

 

 

 

 

 

 

3

 

 

 

13

 

 

 

1

 

 

 

a1

+

a2

=

.

 

20

105

20

 

 

 

 

 

 

 

- 102 -

Лекция 19

Решив эту систему, получим

 

 

a1 =

 

71

 

, a2

=

7

,

 

 

 

 

 

 

 

 

369

 

41

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

таким образом, решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) ≈ y2 (x) =

71

x (1 − x )

+

7

 

x2 (1

x) = x (1

x )

71

+

7

 

x .

369

 

 

 

41

 

 

41

 

 

 

 

369

 

 

Известно точное решение данного дифференциального уравнения:

y(x) = sin(sin(1)x) x .

Для сравнения значений, полученного приближенного и точного решения, приведем таблицу:

xi

0

0.25

0.5

0.75

1

yi

0

0.44

0.069

0.060

0

y(xi )

0

0.44

0.070

0.060

0

Выводы

В лекции освещен вопрос, касающиеся приближенного аналитического решения краевой задачи для дифференциального уравнения 2-го порядка. Детально разобран метод Галеркина, для решения краевой задачи. На примере продемонстрировано применение метода Галеркина для решения краевой задачи для дифференциального уравнения 2-го порядка.

Вопросы для самопроверки

1.Приведите формулу невязки для дифференциального уравнения 2-го порядка.

2.Приведите условия, при которых функции образуют полную замкнутую систему.

3.Приведите пример дифференциального уравнения второго порядка с граничными условиями второго рода.

Рекомендованная литература: 1–10.

- 103 -

ЧИСЛЕННЫЕ МЕТОДЫ

Заключение

В курсе лекций рассмотрены наиболее часто используемые на практике методы, прикладные и научно-технические расчёты: теории приближения функций (интерполяционные формулы Лагранжа и Ньютона, интерполирование сплайнами, наилучшие равномерные приближения, метод наименьших квадратов); численное дифференцирование и интегрирование; решение уравнений и систем; решение задачи Коши; решение краевых задач для ОДУ.

Информация из данного курса лекций способствует пониманию студентами сущности численного решения задач алгебры, математического анализа, дифференциальных уравнений и оценки точности полученных результатов.

Для лучшего освоения материалов приведенных в данном курсе лекций, важную роль играют знания из смежных дисциплин предметной подготовки, в первую очередь «Алгебра» «Математический анализ», «Дифференциальные уравнения», «Вычислительный практикум» и др.

Знания и практические навыки, полученные в ходе применения курса, используются далее при освоении таких дисциплин, как «Дифференциальная геометрия и топология», «Использование алгебраических методов в области защиты информации», «Функциональный анализ», «Комплексный анализ», а также при прохождении производственной практики, выполнения курсовых и выпускных квалификационных работ.

- 104 -

Рекомендованная литература

РЕКОМЕНДОВАННАЯ ЛИТЕРАТУРА

Основная

1.АмосовА.А.идр.Вычислительныеметоды:Учебноепособие.– 4-е изд. стер. – СПб.: Изд-во «Лань», 2014, 672 с.

2.Баxвалов H., Жидков H., Кобельков Г. Численные методы. – Москва – Санкт-Петербyрг: ФИЗМАТЛИТ, 2008.

3.Баxвалов H. С., Лапин А. В., Чижонков Е. В. Численные методы

взадачаx и примераx. – М.: ВШ, 2009.

4.Березин И. С., Жидков Н. П. Методы вычислений. Т.1. – М.: Физматгиз, 1966.

5.Березин И. С., Жидков Н. П. Методы вычислений. Т.2. – М.: Физматгиз, 1962.

6.Вержбицкий В. М. Численные методы. – М.: Высшая школа, 2004.

7.Волков Е. А. Численные методы: Учебное пособие. 5-е изд., стер. – СПб.: Изд-во «Лань», 2008, 256с.

8.Демидович Б. П., Марон И. А. Основы вычислительной математики: Учебное пособие. 8-е изд., стер. – СПб.: Изд-во «Лань», 2011, 672 с.

9.Демидович Б. П., Марон И. А., Шyвалова Э. З. Численные методы анализа. Приближение функций, дифференциальные и интегральные уравнения: Учебное пособие. 5-е изд., стер. – СПб.:

Изд-во «Лань», 2010, 400 с.

10.Зализняк В. Е. Основы научных вычислений. Введение в численные методы для физиков. – М.: Едиториал УРСС, 2002, 296с.

11.Копчёнова Н. В., Марон И. А. Вычислительная математика в задачах и примерах. – М.: Изд-во Лань, 2008, 368 с.

12.СамарскийА.А.,ГyлинА.В.Численныеметоды.–М.:Наука,2008.

13.Срочко В. А. Численные методы. Курс лекций: Учебное посо-

бие. – СПб.: Изд-во «Лань», 2010, 208с.

14.Турчак Л. И., Плотников П. В. Основы численных методов: учебное пособие. – М.: ФИЗМАТЛИТ, 2003, 320 с.

15.Формалёв В. Ф., Ревизников Д. Л. Численные методы. – М.: ФИЗМАТЛИТ, 2006, 400 с.

16.Численные методы: лабораторный практикум / сост. Е. Л. Торопцев, Е. О. Тарасенко, П. К. Корнеев, А. В. Гладков. Ставрополь: Изд-во СКФУ, 2015, 112 с.

17.Шевцов Г. С., Крюкова О. Г., Мызникова Б. И. Численные методы линейной алгебры: Учебное пособие. 2-е изд., испр. и доп. –

СПб.: Изд-во «Лань», 2011, 496с.

-105 -

ЧИСЛЕННЫЕ МЕТОДЫ

Дополнительная

18.Азаров А. И. и др. Сборник задач по методам вычислений. – М.: Hаyка, 1994.

19.Дробышевич В. И., Дымников В. П., Ривин Г. С. Задачи по вычислительной математике. – М.: Hаyка, 1980.

20.Дьяченко В. Ф. Основные понятия вычислительной математи-

ки. – М.: Наука, 1977.

21.Завьялов Ю. С., Квасов Б. И., Мирошниченко В. Л. Методы сплайн-функций. – М.: Наука, 1980.

- 106 -

Учебное издание

Корнеев Петр Кириллович Тарасенко Елена Олеговна Гладков Андрей Владимирович Дерябин Максим Анатольевич

ЧИСЛЕННЫЕ МЕТОДЫ

Часть II

УЧЕБНОЕ ПОСОБИЕ

Издается в авторской редакции Компьютерная верстка: И. В. Бушманова

__________________________________________________________________

 

Подписано к печати 19.11.2018

 

Формат 60x84 1/16

Усл. п. л. 6,28

Уч.-изд. л. 5,81

Бумага офсетная

Заказ 97

Тираж 15 экз.

_________________________________________________________________

Отпечатано в Издательско-полиграфическом комплексе ФГАОУ ВО «Северо-Кавказский федеральный университет» 355009, г. Ставрополь, пр-т Кулакова, 2.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]