- •Введение
- •Раздел I. МЕТОДЫ ПРИБЛИЖЕННЫХ ВЫЧИСЛЕНИЙ
- •1. Приближенное решение уравнения
- •2. Метод последовательных приближений решения уравнения
- •3. Метод Ньютона (метод касательных) решения уравнения
- •4. Метод секущих (метод хорд) решения уравнения
- •8. Приближенные вычисления значений функций с помощью рядов
- •9. Приближенные вычисления пределов с помощью рядов
- •10. Задачи для самостоятельного решения
- •Раздел II. ПРИБЛИЖЕННОЕ ВЫЧИСЛЕНИЕ ОПРЕДЕЛЕННЫХ ИНТЕГРАЛОВ
- •1. Квадратурные формулы
- •3. Формула трапеций
- •4. Метод парабол (метод Симпсона)
- •5. Вычисление определенных интегралов с помощью рядов
- •6. Задачи для самостоятельного решения
- •Раздел III. ЧИСЛЕННОЕ ИНТЕГРИРОВАНИЕ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
- •1. Метод Пикара последовательных приближений
- •2. Метод Эйлера
- •3. Улучшенный метод Эйлера
- •4. Метод Эйлера–Коши
- •5. Метод Рунге–Кутта
- •6. Метод Адамса
- •8. Метод конечных разностей решения краевых задач для обыкновенных дифференциальных уравнений
- •10. Метод наименьших квадратов
- •11. Задачи для самостоятельного решения
- •Раздел IV. УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ ГИПЕРБОЛИЧЕСКОГО ТИПА
- •1. Виды уравнений математической физики
- •2. Вывод уравнения колебания струны
- •7. Задачи для самостоятельного решения
- •Раздел V. УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ ПАРАБОЛИЧЕСКОГО ТИПА
- •4. Распространение тепла в неограниченном стержне
- •6. Задачи для самостоятельного решения
- •Раздел VI. УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ ЭЛЛИПТИЧЕСКОГО ТИПА
- •1. Уравнения эллиптического типа и краевые задачи для них
- •3. Уравнение Лапласа в цилиндрических координатах
- •5. Примеры решения задач
- •6. Задачи для самостоятельного решения
- •Раздел VII. ЧИСЛЕННОЕ ИНТЕГРИРОВАНИЕ УРАВНЕНИЙ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ
- •1. Метод сеток численного решения дифференциальных уравнений с частными производными
- •3. Метод сеток для уравнения гиперболического типа
- •5. Метод прогонки для уравнения теплопроводности
- •6. Решение уравнения движения грунта
- •7. Задачи для самостоятельного решения
- •Раздел VIII. ВАРИАЦИОННЫЕ МЕТОДЫ
- •1. Некоторые сведения из функционального анализа
- •2. Теоретические основы метода Ритца
- •5. Задачи для самостоятельного решения
- •Раздел IX. КОРРЕЛЯЦИЯ
- •1. Понятие корреляции
- •4. Корреляция
- •5. Задачи для самостоятельного решения
- •Библиографический список
Раздел I. МЕТОДЫ ПРИБЛИЖЕННЫХ ВЫЧИСЛЕНИЙ
1. Приближенное решение уравнения
Различают два класса нелинейных уравнений – уравнения ал- |
||
С |
|
|
гебраические и трансцендентные. Алгебраическими уравнениями на- |
||
зывают уравнения, содержащие только алгебраические функции – |
||
целые, рац ональные, иррациональные, в частности многочлен. |
||
Уравнен я, содержащ е тригонометрические, показательные, лога- |
||
рифмическ е |
|
другие неалгебраические функции, называются |
При |
||
трансцендентными. |
||
Методы решен я нелинейных уравнений делятся на две группы: |
||
точные методы |
|
терационные методы. |
способы |
||
решен |
|
уравнений итерационными методами за счет се- |
рии последовательных при лижений получают решение с необходи- |
||
мой точностью. |
Точные методы решения основываются на поиске |
равнос льных прео разований алгебраических выражений. Напри-
мер, перенос слагаемых из одной части уравнения в другую с проти- |
|
А |
|
воположным знаком, деление |
еих частей уравнения на одинаковое |
число, не равное нулю. Точные |
решений позволяют записать |
корни уравнения в виде некоторого конечного выражения или формулы. Точные решения существуют только для некоторых уравнений, например линейных, квадратных, тригонометрических и некоторых других. Для большинства уравнений приходится использовать графи-
ческие или численные методы приближенного решения с заданной |
|
точностью. В первую очередь это относится к большинству трансцен- |
|
дентных уравнений. Отметим, что доказан факт, что нельзя построить |
|
|
И |
формулу для решения произвольного алгебраического уравнения |
|
выше четвертой степени. |
Д |
Например, уравнение x3 cos x 0 нельзя решить путем равно-
сильных алгебраических преобразований. Но это уравнение можно решать приближенно графическими и численными методами приближенных вычислений.
Решение уравнения проводят численно в два этапа. На первом этапе производится отделение корней – поиск интервалов, на которых содержится только по одному корню. Второй этап решения связан с уточнением корня на выбранном интервале, то есть нахождением значения корня с заданной точностью.
5
Отделение корней уравнения может проводиться графически, то есть путем построения графика функции y f x . Для уравнения ви-
да f x 0 , где f x |
– некоторая непрерывная функция, корень или |
||||||||||||||||
корни этого уравнения являются точкой или точками пересечения |
|||||||||||||||||
графика функции с осью абсцисс. |
|
|
|
|
|
|
|||||||||||
С |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
Пусть имеется уравнение вида |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
f x 0 , |
|
|
|
|
|
|
||
чения |
|
|
|
|
|
|
|||||||||||
где |
f x – алгебра ческая или трансцендентная функция. |
|
|
||||||||||||||
|
|
Реш |
ть уравнен |
е – |
значит найти все его корни, то есть те зна- |
||||||||||||
|
|
|
|
x, |
|
которые |
о ращают |
уравнение в |
тождество. |
||||||||
Если уравнен е достаточно сложно, то задача точного определения |
|||||||||||||||||
|
|
|
|
б |
|
|
|
|
|
|
|||||||
корней является в некоторых случаях неразрешимой. Поэтому ста- |
|||||||||||||||||
вится задача найти такое при лиженное значение корня xпр , |
которое |
||||||||||||||||
отличается от точного значения корня x0 на величину, по модулю не |
|||||||||||||||||
превышающую указанной точности |
(малой положительной вели- |
||||||||||||||||
чины), то есть |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
x0 xпр |
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
Величину также называют допустимой ошибкой, которую |
|||||||||||||||
можно задать по своемуАусмотрению. |
|
|
|||||||||||||||
|
|
Рассмотрим задачу решения уравнения |
s i nx |
1 |
, |
где угол x |
|||||||||||
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
задан в |
градусах. |
Это |
уравнение |
можно |
переписать |
в виде |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
И |
|||||
|
|
|
|
|
1 |
0. Для графическогоДотделения корней достаточно |
|||||||||||
sin |
|
|
x |
|
|||||||||||||
|
|
|
|||||||||||||||
|
180 |
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
построить график функции y sin |
|
x |
|
|
(рис. 1). |
|
|
||||
180 |
|
|
x |
|
6
С |
корней основано на следующих тео- |
||
тделениеАнал т ческое о |
|||
|
|
Рис. 1 |
|
|
б |
|
|
Из р с. 1 в дно, что корень уравнения лежит в промежутке |
|||
x 6; 8 . |
|
|
|
ремах. |
А |
||
Теорема 1. Если непрерывная функция |
f x принимает на кон- |
цах отрезка [ a , b ] значения разных знаков, то есть f a f b 0 , то на этом отрезке содержится по крайней мере один корень уравнения.
Теорема 2. Если непрерывная на отрезке [ a , b ] функция f x
принимает на концах отрезка значения разных знаков, а производ- Д
ная f x сохраняет знак внутри указанного отрезка, то внутри отрез-
ка существует единственный корень уравнения f x 0 .
Для уточнения корней может использоваться один из следую-
щих методов: |
И |
|
– метод последовательных приближений (метод итераций);
– метод Ньютона (метод касательных);
– метод секущих (метод хорд);
– метод половинного деления (метод дихотомии).
2. Метод последовательных приближений решения уравнения
Метод последовательных приближений (метод итераций) – численный метод решения математических задач, используемый для приближённого решения алгебраических уравнений и систем. Суть метода заключается в нахождении по приближённому значению ве-
7
личины следующего приближения, являющегося более точным. Метод позволяет получить решение с заданной точностью в виде предела последовательности итераций. Характер сходимости и сам факт сходимости метода зависит от выбора начального приближения ре-
шения. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
С |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Функциональное |
уравнение |
|
может быть |
записано в виде |
||||||||||||||||
f x x . Функцию |
f |
x |
называют сжимающим |
отображением. |
||||||||||||||||
Последовательность |
ч сел |
x0 , x1,…, xn |
|
называется итерационной, |
||||||||||||||||
если для любого номера n 0 элемент xn |
выражается через элемент |
|||||||||||||||||||
ли |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
xn 1 по рекуррентной формуле |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
xn f xn 1 , |
|
|
|
|
|
|
|
|
|
|
|
|
||||
а в качестве x0 |
взято лю ое число из области задания функции |
|
f x . |
|||||||||||||||||
Получ |
рекуррентное соотношение для нахождения корня |
|||||||||||||||||||
уравнен я f x 0 методом итераций. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
А |
|
|
|
|
|
|
|
|||||||||||
Пример. |
Для |
решения |
уравнения |
sin x |
1 |
|
построить |
рекур- |
||||||||||||
рентную формулуб. |
|
|
|
x |
|
|
|
|
|
|
|
|||||||||
Решение. Уравнение перепишем в форме |
|
x |
1 |
|
|
|
. Ре- |
|||||||||||||
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin |
|
|
x |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
Д |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
180 |
|
||||
куррентная формула решения уравнения методом последовательных |
||||||||||||||||||||
приближений имеет вид |
|
|
|
|
|
|
|
И |
||||||||||||
|
|
|
|
xn |
|
1 |
|
|
. |
|||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
sin |
|
|
xn 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
180 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3. Метод Ньютона (метод касательных) решения уравнения |
|
Пусть функция y f x непрерывна и дифференцируема на от- |
|
резке [ a , b ] и выполняется условие |
f a f b 0 . Предположим, что |
известно начальное приближение x0 |
корня уравнения f x 0 ..Пусть |
найдено приближение x корня уравнения f x 0 . Проведем каса-
n 1
тельную к кривой y f x , проходящую через точку xn 1 ; f xn 1 .
8
Уравнение касательной имеет вид
y f xn 1 f xn 1 x xn 1 .
Найдем координату пересечения касательной с осью Оx:
СибАДИ |
|||||||||||||||||||||
|
|
|
|
|
|
x |
n |
x |
n 1 |
|
f |
xn 1 |
|
. |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
xn 1 |
|
|
|
|
||
Получ ли рекуррентное соотношение для нахождения корня |
|||||||||||||||||||||
уравнен я f x |
0 методом Ньютона. |
|
|
|
|
|
|||||||||||||||
Процесс продолжаем до тех пор, пока не будет достигнута тре- |
|||||||||||||||||||||
буемая точность выч сления. |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Граф ческая |
нтерпретация метода касательных представлена |
||||||||||||||||||||
на рис. 2. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 2 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
||
Пример. |
Для уравнения |
|
f x sin |
|
|
|
x |
|
0 производная |
||||||||||||
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
180 |
|
x |
|
|||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
||||
равна f x |
|
|
cos |
|
|
|
x |
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
180 |
|
|
x2 |
|
|
|
|
|
|
|
|
||||||||
180 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Получаем, что рекуррентная формула решения уравнения методом Ньютона имеет вид
9