vychisl_matem


Чтобы посмотреть этот PDF файл с форматированием и разметкой, скачайте его и откройте на своем компьютере.
1

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ πЕДЕРАЦИИ

федеральное г
осударственное
бюджетное
образовательное учреждение

высшего профессионального образования



Тольяттинский государственный университет












УЧЕБНО
-
МЕТОДИЧЕСКОЕ ПОСОБИЕ

по проведению практических занятий

по дисциплине


ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА



по направлению

010
3
00.68
Математическое обеспечение и администрирование
информационных систем

















Тольятти 01
5



2

Содержание


1

Планы практических занятий

................................
................................
................................
....

3

1.1

Занятие № 1. Теория погрешности.

................................
................................
...................

3

1.2

Занятие № . Решения алгебраических уравнений. Отделение корней.

......................

11

1.3

Занятие № 3. Метод хорд (секущих) и метод деления пополам.

................................
..

13

1.4

Занятие № 4. Метод Ньютона

(касательных).

................................
................................

17

1.5

Занятие № 5. Комбинированный метод.

................................
................................
.........

20

1.6

Занятие № 6. Метод итераций.

................................
................................
.........................

21

1.7

Занятие № 7. Прямые методы решения систем линейных алгебраических уравнений.

24

1.8

Занятие № 8. Метод итераций для систем уравнений.

................................
..................

29

1.9

Занятие № 9. Равномерное приближение функций многочленами.
Интерполяционный многочлен Лагранжа.

................................
................................
................

35

1.10

Занятие № 10. Интерполяционный многочлен Ньютона. Многочлены Чебышева.

44

1.11

Занятие № 11. Интерполяция с кратными узлами и формула Эрмита.

....................

65

1.12

Занятие № 1. Понятие сплайна. Виды спл
айнов.

................................
.....................

69

1.13

Занятие № 13. Численное дифференцирование.

................................
.........................

76

1.14

Занятие № 14. Численное интегрирование. Квадратурные формулы Ньютона
-
Котеса.

79

1.15

Занятие № 15. Квадратурны
е формулы Гаусса.

................................
.........................

88

1.16

Занятие № 16. Численные методы решения задачи Коши для обыкновенных
дифференциальных уравнений. Классифик
ация методов. Метод Пикара.

............................

93

1.17

Занятие № 17. Методы Эйлера.

................................
................................
....................

96

1.18

Занятие № 18. Методы Рунге
-
Кутта.

................................
................................
.........

105

1.19

Занятие № 19. Пошаговый контроль точности. Метод Кутты
-
Мерсона
.

..............

110

1.20

Занятие № 0. Многошаговые методы Адамса.

................................
.......................

114

1.21

Занятие № 1.
Методы прогноза и коррекции.

................................
.........................

119

1.22

Занятие № . Метод Милна четвертого порядка

................................
....................

122

Приложение 1

................................
................................
................................
................................
.

126

Приложение .

................................
................................
................................
................................

128

Индивидуальное домашнее задание № 1

................................
................................
.....................

129

Индивидуальное задание №.

................................
................................
................................
.......

134

Индивидуальное задание №3

................................
................................
................................
........

139

Литература

................................
................................
................................
................................
......

143




3

1

План
ы

практических занятий

1.1

Заняти
е



1
.
Теория погрешности.

Источники и виды погрешности

На некоторых этапах решения задачи на ЭВМ могут возникать
погрешности, искажающие результаты
вычислений. Оценка степени
достоверности получаемых результатов является важнейшим вопросом при
организации вычислительных работ. Это особенно важно при отсутствии
опытных или других данных для проверки адекватности модели.

Погрешность решения задачи обусл
авливается следующими причинами:

1.

Математическое описание задачи является неточным, в частности,
неточно заданы исходные данные.

2.

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

3.

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

1)

неустранимой погрешностью,

2)

погрешностью метода,

3)

вычислительной погрешностью.

Часто неустранимую погрешность разделяют на две части:

а)
неустранимой погрешностью
называют лишь погрешность,
являющуюся следствием неточности задания числовых данных, входящих в
математическое описание задачи;

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

4

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

условиях и быть
совершенно неприемлемой в других; поэтому важно правильно учитывать
область еѐ применимости.

Численный метод также является источником погрешностей. Это связано,
например, с заменой интеграла суммой, усечением рядов при вычислениях
значени
й функций, интерполированием табличных данных и т. п. Как правило,
погрешность численного метода регулируема, т. е. она может быть уменьшена
до любого разумного значения путем изменения некоторого параметра
(например, шага интегрирования, числа членов усеч
енного ряда и т.п.).
Погрешность метода обычно стараются довести до величины, в несколько раз
меньшей погрешности исходных данных. Дальнейшее снижение
погрешности не приведет к повышению точности результатов, а лишь
увеличит стоимость расчето
в из
-
за необоснованного увеличения объема
вычислении.

При вычислениях с помощью ЭВМ неизбежны
погрешности округлений,
связанные с ограниченностью разрядной сетки машины. Обычно после
выполнения операции производится не округление результата, а простое
отбр
асывание лишних разрядов с целью экономии машинного времени. Правда,
в современных машинах предусмотрена свобода выбора программистом
способа округления; соответствующими средствами располагают и некоторые
алгоритмические языки.

Максимальная относительная
погрешность при округлении есть
δ
max
=0.5
α
1
-
k
, где
α

-

основание системы счисления,
к

количество разрядов
мантиссы числа. При простом отбрасывании лишних разрядов эта погрешность
увеличивается вдвое.

В современных машинах с памятью, измеряемой в байтах, при
нята
шестнадцатеричная система счисления, и любое число с плавающей точкой
5

содержит шесть значащих цифр. Следовательно,
α

=16,
к

=
6, максимальная
погрешность округления α
m
ах
= 0.5* 16
-
5

≈ 0.5* 10
-
8
.

Несмотря на то, что при решении больших задач выполняются

миллиарды
операций, это вовсе не означает механического умножения погрешности при
одном округлении на число операций, так как при отдельных действиях
погрешности могут компенсировать друг друга (например, при сложении чисел
разных знаков). Вместе с тем ин
огда погрешности округлений в сочетании с
плохо организованным алгоритмом могут сильно исказить результаты.

Перевод чисел из одной системы счисления в другую также может быть
источником погрешности из
-
за того, что основание одной системы счисления
не являе
тся степенью основания другой (например, 10 и ). Это может привести
к тому, что в новой системе счисления число становится иррациональным.

Например, число 0.1 при переводе в двоичную систему счисления примет
вид 0.10.00011001100... Может оказаться, что с

шагом 0.1 нужно при
вычислениях пройти отрезок [0,1] от х1 до х0; десять шагов не дадут точного
значения
х0.

Абсолютная и относительная погрешности

Если
А
-

точное значение некоторой величины,
a

a

-

известное
приближение к нему, то абсолютной погрешностью приближения
a

числа
А
называют некоторую величину Δ(
a
) удовлетворяющую условию:


Относительной погрешностью
называют некоторую величину
δ(А),
для
которой выполняется условие:


Относительную погрешн
ость часто выражают в процентах.

Информацию о том, что
а
является приближенным значением числа
А
с
абсолютной погрешностью
Δ(а),
принято записывать в виде:


Числа
а
и
Δ(а)
записываются с одинаковым количеством знаков после
запятой.

6

Информацию о том, что
а

является приближенным значением числа
А
с
относительной погрешностью
δ(а),
записывают в виде:


Число δ
α
, заведомо не меньшее относительной погрешности называют
предельной относительной погрешностью, т.е.


Т.к. на практике A≈ α, то приближенно можно прин
ять, что



Погрешность функции

Пусть искомая величина
Y

является функцией параметров
a
1
, а
2
, …,
a
n
;
т.е.
Y

=
Y
(
a
),
и известна область
G

в пространстве переменных
a
1
, а
2
, …,
a
n
,
которой принадлежат параметры. Необходимо получить приближение
y

к
Y

и
оценить его погрешность.

Если
у
-

приближѐнное значение величины
Y
, то
предельной абсолютной
погрешностью А(у)
называют наилучшую при имеющейся информации оценку
погрешности величины
у:


7

Предельной относительной погрешностью δ(у)
называют величину


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


Для оценки предельной относительной погрешности функции имеют
место выражен
ия


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



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


тогда если слагаемые одного знака,
получаем

8


Т.о., если слагаемые одного знака, то предельная относительная
погрешность их суммы не превышает наибольшей из предельных
относительных погрешностей слагаемых.

Если слагаемые разных знаков, то
предельная относительная погрешность
суммы вычисляется по формуле


Пример
. Оценить абсолютную погрешность функции
считая
абсолютную предельную погрешность аргумента известной.


Таким образом, абсолютная погрешность тангенса всегда больше
абсолютной погр
ешности аргумента.

Пример
. Вычислить и определить погрешность результата


Решение: Произведем вычисления


Определим относительные погрешности аргументов

9






Устойчивость, корректность, сходимость

Устойчивость.
Рассмотрим погрешности исходных данных.

Поскольку
это так называемые неустранимые погрешности и вычислитель но может с
ними бороться, то нужно хотя бы иметь представление об их влиянии па
точность окончательных результатов. Конечно, мы вправе надеяться на то, что
погрешность результатов имеет п
орядок погрешности исходных данных.
Всегда ли это так? К сожалению, нет. Некоторые задачи весьма чувствительны
к неточностям в исходных данных. Эта чувствительность характеризуется так
называемой устойчивостью.

Пусть в результате решения задачи по исходном
у значению величины
х
находится значение искомой величины
у.
Если исходная величина имеет
абсолютную погрешность Δх, то решение имеет погрешность Δу. Если решение
непрерывно зависит от входных данных, т.е. всегда ||Δу||→0 при ||Δх||→0 , то
задача называетс
я
устойчивой
по входным данным.

Отсутствие устойчивости обычно означает, что даже незначительные
погрешности в исходных данных приводят к большим погрешностям в решении
или вовсе к неверному результату. О подобных устойчивых задачах также
говорят, что они
чувствительны
к погрешностям исходных данных.

10

Примером такой задачи является отыскание действительных корней
уравнения вида:
tg
(
x
)
-
х  1.
Изменение аргумента
х
на малую величину
Δх
может привести к существенному изменению функции.

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

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

Вместе с тем отметим, что
в настоящее время развиты методы решения
некоторых некорректных задач. Это в основном так называемые
методы
регуляризации.
Они основываются на замене исходной задачи корректно
поставленной задачей. Последняя содержит некоторый параметр, при
стремлении кото
рого к нулю решение этой задачи переходит в решение
исходной задачи.

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

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

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

11

Рассмотрим понятие
сходимости итерационного процесса.
Этот процесс
состоит в том, что для решения некоторой задачи и нахождения искомого
значения определяемого параметра (например, корня нел
инейного уравнения)
строится метод последовательных приближений. В результате многократного
повторения этого процесса (или
итераций)
получаем последовательность
значении
х
1
, х
2
, …, х
п
,...
эта последовательность
сходится
к точному решению
ха, если при нео
граниченном возрастании числа итераций предел этой
последовательности существует и равен
a
:
limx
n

=
а
.
В этом случае имеем
сходящийся численный метод.

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

Таким образом, для получения решения задачи с необходимой точностью,
ее постановка должна быть корректной, а используемый численный метод
должен обладать устойчивостью и сходимостью.

Задание.
Выполнить задание 1 ИДЗ№1.

1.2

Занятие



2
.

Решения алгебраических уравнений. Отделение корней.

Уравнение
f
(
x
)=0
называется алгебраическим уравнением
n
-
й степени,
если
f
(
x
)
является полиномом,

Всякое неалгебраическое уравнение называется трансцендентным. Любое
значение
A
, обращающее функцию
f
(
x
)
в нуль называется корнем уравнения
f
(
x
)=0.

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

корней и трансце
ндентных уравнений. Поэтому, важное значение приобретают
методы приближенного нахождения корней уравнения и оценка степени их
точности.

Сам процесс вычисления корней состоит из двух операций:

1)

отделение корней, т.е. установление возможно тесных промежутков
[
a
,
b
], в которых содержится один и только один корень уравнения;

2)

уточнение приближенных корней, т.е. доведение их до заданной
степени точности.

Для отделения корней часто используют теорему Больцано
-
Коши.

Теорема:
Если непрерывная функция
f
(
x
) принимает
значение разных
знаков на концах отрезка [
a
,
b
] то внутри этого отрезка содержится по крайней
мере один корень уравнения
f
(
x
)=0.
Корень заведомо будет единственным, если
производная
f
'(
x
) существует и сохраняет знак внутри интервала [
a
,
b
].


Отделение корней

уравнения
f
(
x
)=0
можно выполнить графически,
построив график функции
f
(
x
),
по которому можно судить о том, в каких
промежутках находятся точки пересечения его с осью 0х. В некоторых случаях
целесообразно представить исходное уравнение в эквивалентном виде
:
с таким расчетом, чтобы графики функций
y
=
f
1
(
x
)
и
y
=
f
2
(
x
)
строились проще, чем график
y
=
f
(
x
).
Корень уравнения
f
(
x
)=0
представляет
собой абсциссу точки пересечения графиков
y
=
f
1
(
x
)
и
y
=
f
2
(
x
)
.

Пример.
Отделить корни уравнения

Решение. 1
-
способ.
В данном случае
f
(
x
)=
x

3
+
2
х
-
1
,
.
Поскольку
f
'(
x

�0

при всех
x
,
то

функция
f
(
x
)
возрастает в промежутке

Корень
считается отделенным, если указан конечный промежуток
[
a
,
b
],
на котором
он находится.

Методом проб
находим отрезок
[
a
,
b
],

для которого на концах отрезка
функция
f
(
x
)

принимает значения разных знаков. Для этого вычислим
значения функции при некоторых значениях аргумента

13

f
(
-
1)=
-
40,
f
(0)=
-
10,
f
(1)=�20.

Согласно

теореме

Больцано
-
Коши корень находится на отрезке [0,1],
причем он единственный, т.к. производная
f
'
(
x
�)0
сохраняет знак внутри
интервала [0,1].

2
-
способ.
Корень данного уравнения можно отделить и графически.
Придадим уравнению вид
х
3
=
-
2
х1,

т. е.

вид

f
1
(
x
)=
f
2
(
x
),
и построим графики

функций
ух
3

и
y
=
-
2

х

+1
.

Эти графики пересекаются в точке М,
абсцисса которой принадлежит интервалу
(0,1).

Задания.
Выполнить задание 
.1 ИДЗ№1.

1.3

Занятие



3.
Метод хорд (секущих) и метод деления пополам.

Метод деления пополам
(метод дихотомии).

Наиболее надѐжным алгоритмом нахождения корня уравнения
f
(
x
)=0,
особенно когда о поведении функции
f
(
x
) мало что известно, является метод
половинного деления.

Пусть
f
(
x
)
-

функция действительной переменной и пусть известен
интервал [
a
,
b
], на котором функция меняет знак, следовательно, между а и
b

существует точка, в которой функция обращается в нуль. Если разделить
интервал пополам и определить положительна или отрицательна функция в
точке деления, то тем самым найдѐм подынтервал, в кото
ром функция меняет
знак.

В принципе, повторным применением этого приема (деление интервала
пополам) можно сколь угодно близко "подойти к корню". Так как каждый шаг
делит интервал, в котором лежит корень, пополам, то 10 шагов уменьшают
интервал в 
10
104 р
аз. При заданной абсолютной точности
ε
алгоритм метода
деление пополам состоит из следующих шагов:

1.

Вычислить
f
(
a
)
и
f
(
b
).

2.

Положить с 

+
b
) / . Вычислить
f
(
c
).


14

3.

Если
sign
(
f
(
c
))=
sign
(
f
(
a
)), то заменить
а
на с; в противном случае
заменить
b

на с (функция
sign
(
f
(
c
)) означает знак в точке с).

4.

Если
b

-

a


ε
, то перейти к шагу ; в противном случае прекратить
вычисления, поскольку мы достигли требуемой точности. Любой из концов
отрезка или их полусумма может быть использован в качестве корня уравнения
f
(
x
) = 0
.

Алгоритм деления пополам довольно медлителен, но зато
абсолютно застрахован от неудач.

Пример
. Пользуясь методом половинного деления вычислить с
точностью до ε0,1 корень уравнения
x
3
+3
x
2
-
3=0,
заключенный в отрезке [
-
3,
-
2].

Решен
ие.
Описанный выше процесс решения удобно оформить в виде

вычислительного бланка:

n

a
n

b
n

f(a
n
)

f(b
n
)

c
n

f(c
n
)

(b
n
-
a
n
)/2

0

-
3

-
2

-
3

1

-
2,5

0,125

0.5

1

-
3

-
2,5

-
3

0,125

-
2,75

-
1,11

0.25

2

-
2,75

-
2,5

-
1,11

0,125

-
2,625

-
0,42

0.125

3

-
2,625

-
2,5

-
0,42

0,125

-
2,5625

-
0,129

0.0625

Из таблицы следует, что
A
=
-
,565“0,065. Или округлив результат,
получаем
А
=
-
2,60,1.

Метод пропорциональных частей (метод хорд).

Мы будем предполагать выполнение следующих условий:



функция
f
(
x
)
в промежутке [
a
,
b
] непрерывна вместе со своими
производными



значения
f
(
a
) и
f
(
b
)
функции на концах промежутка имеют разные
знаки: т.е.
f
(
a
)*
f
(
b
)0;



обе производные
f
' (
x
) и
f
'' (
x
)
сохраняют каждая определенный знак во
всем промежутке
[а,
b
].

Тогда возможны следующие
чет
ыре
случая, изображенные на рисунке 1.










1.
Рассмотрим случай
f
(
a
)
f
"(
x
)0 ,
которому
соответствуют графики,
изображенные на рис.1б или 1в. Заменим дугу ММ' кривой (рис.1а)
-

хордой
ММ'. Уравнение хорды может быть написано, например, в виде









(1)

Наше правило, по существу, сводится к тому, что вместо точки А
-

пересечения кривой
f
(
x
)
с осью 0
x
, определяется точка
D

-

пересечение с осью
0
x

соответственно хорды
MM
'. Полагая у0 в формуле (1), получим значение
точки
a
1
:







(2)


Через точки
(
a
1
,
f
(
a
1
))
и
(
b
,
f
(
b
))
новую хорду получим точку
a
2
.
Продолжая этот процесс, получим последовательность точек сходящихся к
корню
А.
πормулы для расчета последовательности
a
n

следующие





(3)


2.
В случае
f
(
a
)
f
"(
x
�)0
(см. рис.1а или 1г) приближение к корню А
происходит со стороны точки
b
,
а последовательность точек определяется по
рекуррентной формуле







(4)

Для
оценки погрешности метода хорд
будем предполагать, что
производная функции
f
(
x
)
непрерывна и сохраняет знак на отрезке
[
a
,
b
].
Тогда,
если
А
-
точный корень, а
x
n

-

приближение к корню, то по формуле конечных
приращений
(теорема Лагранжа), имеем









(
5
)

где точка с лежит между
А
и
x
n
.

Из формулы () получаем







(6)

где
m

наименьшее значение модуля производной функции
f
(
x
)
на отрезке [
a
,
b
],
т.е.

Другая формула для оценки погрешности метода хорд







(7)

где
М
и
m

соответственно наибольшее и наименьшее значение модуля
производной функции
f
(
x
).

Пример.
Найти корень уравнения
x
3

-
2
x
2

-
4
x
-
7=0
интервале [3,4], с
точность до 0,01.

Решение.
Обозначим левую часть уравнения через
f
(
x
)
. Подсчитаем

f
(3)
=
-
10<0и<

f
(4)

= 9�
0. Нетрудно показать, что обе производные
f
'(
x
)=3
x
2
-

4
x
-
4,
f
''(
x

) = 6
x
-

4
в промежутке

[3,4
] положительны.

Действительно, для
f
"(
x
)
это очевидно, а т.к.
f
"(
x
�)
0
,
то
f
'(
x
)
монотонно
возрастает. Так как
f
'(3)
=
11

0
,
то в других точках отрезка [3,4] значения


производной будут заведомо не меньше 11. Попутно мы показали, что
наименьшее значение первой производной равно
m
=
f
'(3)
=
11.

Т.к.
f
(3)
f
"(
x
)0,

то для дальнейших расчетов
следует
воспользоваться
формулами (
3
) и (6). Обозначим

,
тогда получим, что


Описанный выше процесс решения удобно оформить в виде
вычислительного бланка.

Процесс прервался после третьего шага, т. к. мы достигли нужную точность:
ε
n
0,004<0,01. Приближенное
значение корня А3,630“0,01.

Задания
.
Выполнить задание

2
. ИДЗ№1.

1.4

Занятие №

4.
Метод Ньютона (касательных).

Цель
-

ознакомить студентов с методом Ньютона решения
алгебраических уравнений.

Метод касательных (Метод Ньютона)

Пусть корень уравнения
f
(
x
)=0
отделен на отрезке
[
a
,
b
],
причем
f
'(
x
) и
f
''(
x
)
непрерывны и сохраняют свои знаки
на этом отрезке. Тогда, как было
сказано выше,

возможны
четыре
случая, графически изображенные на рис.1 а)
-
г).

Геометрически метод Ньютона эквивалентен замене дуги кривой ММ'

касательной, проведенной к некоторой точке данной кривой.



1.
Рассмотрим случай
f
(
a
)
f
"(
x
)0 ,
которому соответствуют графики
изображенные на рис.1 б) или 1 в). В этом случае, приближение к корню
А
происходит со стороны точки
b

(см. рис.). Касательные проводятся в точках
M
0
,
M
1
,
M
2
,

... .Пересечением касательных с осью 0х является
последовательность точек
b
0
,
b
1
,
b
2
,
которая сходится к точному корню
А.
Уравнение касательной к функции
f
(
x
)
в точке
M
n

с координатами (
b
n

,
f
(
b
n
))
имеет вид:








(8)

Полагая
y
=0 ,
x
=
b
n
+1

получаем





(9)

Заметим, что если бы мы в данном случае провели касательную в точке
(
a
,
f
(
a
)),
то точка пересечения этой касательной с осью 0х лежала бы вне отрезка
[
a
,
b
],
и требуемой сходимости к корню не наблюдалось бы.


2.
В случае
f
(
a
)
f
"(
x
�)0
(см. рис.1 а или 1г) приближение к корню
А
происходит со стороны точки
a
,
а последовательность точек опр
еделяется по
рекуррентной формуле





(
10)

Для
оценки погрешности метода касательных
будем предполагать, что
обе производные функции
f
(
x
)
непрерывны и сохраняют свой знак на отрезке
[
a
,
b
].
Тогда, если
А
-

точный корень, а
x
n

-

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







(11)





(12)

где


Если

,
то формула (1) упрощается





(12*)

Пример.
Вычислить приближенно по методу касательных с точностью до
пяти десятичных знаков после запятой больший отрицательный корень
уравнения
x
3

-
12
x
-
8 = 0.

Решение.
Графически отделяя корни данного уравнения, заключаем, что
уравнение имеет три действительных кор
ня; больший отрицательный корень
принадлежит отрезку [
-
1,0]. Укажем отрезок меньшей длины, на котором
находится корень; это отрезок [
-
0,7;
-
0,65], поскольку
f
(
-
0,7) = 0,057
,
f
(
-
0,65)=
-
0,474625.

Т.к.
f
"(
x
) = 6
x
,

то очевидно для отрицательных
x

имеем

f
"(
x
)0
,

то
выполнено условие
f
(
a
)
f
"(
x
)0
, поэтому в качестве нулевого приближения
берем
b
0
=
-
0,65,
а последующие приближения вычисляем по формуле (9).

Описанный выше процесс решения удобно оформить в виде
вычислительного бланка.


Процесс прервался, после

третьего шага, т.к. мы достигли нужную
точность, которая в данном случае может быть рассчитана по формуле (1*).
Приближенное значение корня А
-
0,6945930,000001.


Задания.

Выполнить
задание

2
.
2

ИДЗ№1.

1.5

Занятие №

5.
Комбинированный метод.

Цель
-

ознакомить
студентов с комбинированным методом решения
алгебраических уравнений.

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

Преимущество комбинированного метода заключается в «двустороннем
сжатии» рассматриваемого отрезка, что наглядно показано на рис.3. Из рис.3
также видно, что истинный корень
А

лежит между приближениями со стороны
точки
a

и со стороны точки
b
,

т.е.
а

n


A


b
n

.


Также как, и отдельно в методе хорд и в методе касательных, для расчета
приближений мы рассмотрим два случая.

1.
Случай




(13)

2.
Случай





(14)

Если задана допустимая погрешность ε
,
то процесс уточнения корня
прекратится, как только
За значение корня
А
следует взять


Задания
.
Выполнить задание
2
.
2

ИДЗ№1.

1.6

Занятие № 6.
Метод
итераций
.

Цель
-

ознакомить студентов с методом итераций решения
алгебраических уравнений.

Если каким
-
либо способом получено приближенное значение
x
0

корня
уравнения
f
(
x
)=0,
то уточнение корня можно осуществить методом
последовательных приближений или
методом итераций. Для этого уравнение

f
(
x
)=0
представляют в виде


x

 φ (
x
)


(15)

что всегда можно сделать, и притом многими способами, например,


x

=
x

+
c


f
(
x
)


(
16)

где с
-

произвольная постоянная.

Пусть число
x
1

есть результат подстановки
x
0

в правую часть уравнения
(15), т.е.
x
1
φ(
x
0
).

Далее,
x
2
φ(
x
1
),
x
3
φ(
x
2
),...,




x
n
+1

φ (
x
n
)





(17)

Процесс последовательного вычисления чисел
x
n

(
n
1,,3,...). по формуле
(17) называется методом последовательных приближений или методом
итераций.

Утверждение.
Итерационный процесс сходится

(здесь А
-

точный корень исходного уравнения
f
(
x
)=0),

если на отрезке
[
a
,
b
]

содержащем
корень
А

и его последовательные приближения
x
n
,

выполнено условие






(
1
8
)

Замечание.
В качестве первого приближения
x
0

можно взять
произвольное значение из интервала, содержащего корень
A
, причем такой
интервал можно сделать достаточно малым.

Для оценки погрешности метода итераций используется следующая
формула






(19)

При
оценка погрешности упрощается





(20)

Пример 1 .
Методом итераций найти меньший положительный корень
уравнения
x
3

-
5
x
+1=0.

Решение.
Графически отделяя корни данного уравнения, заключаем, что
уравнение имеет три действительных корня, лежащих в отрезках [
-
3;
-
2], [0;1],
[;3]. Найдем корень принадлежащий отрезку [0; I]. Укажем отрезок меньшей
длины, на котором находится корень. Поскольку,

f
(
x
)=
x
3

-
5
x
+1,
f
(0)=1�0,
,
то корень принадлежит отрезку [0; 0,5].
Данное уравнение
приведем к виду (15), разрешая его относительно
x
:

т.е.

Так как,

то условие (18)
выполняется,
следовательно процесс итераций будет сходиться, а для оценки
погрешности
можно использовать (0).

Взяв в качестве начального приближения середину отрезка, т.е. приняв
x
0
0,5, вычисление последующих приближений проведем по формуле



Результаты этих вычислений представлены в таблице_1, из которой
видно, что искомый к
орень
A
=
0,0164, вычислен с заданной точностью
ε0,0001, т.к. согласно формуле (0) имеем



Замечание.
При нахождении двух других корней исходного уравнения
методом последовательных приближений уже нельзя пользоваться формулой


т.е. условие (18) не выполняется. В этом случае данное уравнение следует
представить в другом виде, например,

для функции

условие (18) на отрезках [
-
3;
-
], [;3] будет выполняться.

Пример .
Методом итераций найти отрицательный корень урав
нения

x
4

+
x

-
3 = 0

Решение.
Данное уравнение имеет два действительных корня;
отрицательный корень находится на отрезке [
-
1,5;
-
1,4], так как для его концов
выполняется условие


Уравнение запишем в виде (16):

x

=
x

+
c

• (
x

4

+
x



3
)

где
с
-

произвольная постоянная. Выберем значение
с
так, чтобы для
функции
на отрезке [
-
1,5;
-
1,4] выполнялось условие (18).
В качестве такого значения можно взять с0,1, тогда


φ (
x
) =
x

 0,1 • (
x
4

+
x



3 )

,


т.е. выполнено условие (1
8
).

Возьмем
x
0
=
-
1,45.
Вычисление последующих приближений осуществим
по формуле


и представим результаты в таблице_, из которой видно, что
A
=
-
1,45262
является приближенным корнем данного уравнения.

Таблица_


Задания:

Выполнить задание
2
. и 3

ИДЗ№1.

1.7

Занятие № 7
.
Прямые
методы решения систем линейных
алгебраических уравнений.

Существует несколько методов, позволяющих получить решение
системы
n

линейных уравнений с
n

неизвестными приближенно, с заданной
точностью.

Способы решения систем линейных уравнений делятся на две г
руппы:

1)

точные методы
, представляющие собой конечные алгоритмы для
вычисления корней системы (решение систем с помощью
обратной матрицы, правило Крамера, метод Гаусса и др.),

2)

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

(метод итерации, метод Зейделя и др.).

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

Эффект
ивное применение итерационных методов существенно зависит от
удачного выбора начального приближения и быстроты сходимости процесса.

Метод Гаусса

1. Схема единственного деления.
Рассмотрим сначала простейший
вариант метода Гаусса, называемый
схемой
единственного деления.

Прямой ход

состоит из
n

-

1
шагов исключения.

1
-
й шаг.

Целью этого шага является исключение неизвестного
x
1

из
уравнений с номерами
i

= 2, 3,
n
. Предположим, что коэффициент
a
11



0.
Будем называть его
главным элементом
1
-
го

шага.

Найдем величины


называемые
множителями 1
-
го шага.
Вычтем последовательно из второго,

третьего, …,
n
-
го уравнений системы первое уравнение, умноженное
соответственно на
q
21
,
q
31
, …,
q
n
1
.
Это позволит обратить в нуль коэффициенты
при
x
1

во

всех уравнениях, кроме первого. В результате получим
эквивалентную систему



в
которой
вычисляются по формулам


2
-
й шаг.

Целью этого шага является исключение неизвестного
х
2

из
уравнений с номерами
i

=
3, 4,
n
.
Пусть
a
22
(1)


0,
где
a
22
(1)

-

коэффициент,

называемый
главным
(или
ведущим) элементом
2
-
го шага.
Вычислим
множители
2
-
го шага


и вычтем последовательно из третьего, четвертого,
.
..
,
n
-
го уравнения
системы второе уравнение, умноженное соответственно на
q
32
,
q
42
, …,
q
m
2
.
В
результате
получим систему


Здесь коэффициенты
вычисляются по формулам


Аналогично проводятся остальные шаги. Опишем очередной
k
-
й шаг.

к
-
й шаг.

В предположении, что
главный (ведущий) элемент
k
-
го

шага
a
kk
(
k
-
1)

отличен от нуля, вычислим
множители
k
-
го шага


и
вычтем последовательно из (
k

+ 1
)
-
го, …,
n
-
го уравнений полученной на
предыдущем шаге системы
k
-
e

уравнение, умноженное соответственно на

После
(
n

-

1
)
-
го шага исключения получим систему уравнений


матрица

которой является верхней треугольной. На
этом вычисления
прямого хода заканчиваются.


Обратный ход.

Из последнего уравнения системы находим
x
n
.
Подставляя найденное значение
x
n

в предпоследнее уравнение, получим
x
n
-
1
.
Осуществляя обратную подстановку, далее последовательно находим
x
n
-
1
,
x
n
-
2
, .,
x
1
.
Вычисления неизвестных здесь проводятся по формулам


Необходимость выбора главных элементов.

Заметим, что вычисление
множителей, а также обратная подстановка требуют деления на главные
элементы
a
kk
(
k
-
1)
.
Поэтому, если один из главных элементов
оказывается равным
нулю, то схема единственного деления не может быть реализована. Здравый
смысл подсказывает, что и в ситуации, когда все главные элементы отличны от
нуля, но среди них есть близкие к нулю, возможен неконтролируемый рост
погрешности.


Метод Гаусса с выбором главного элемента по всей матрице (схема
полного выбора).

В этой схеме допускается нарушение естественного порядка исключения
неизвестных.

На
1
-
м шаге метода среди элементов
a
iJ

определяют максимальный по
модулю элемент

Первое ура
внение системы и уравнение с номером
i
1

меняют местами. Далее стандартным образом производят исключение
неизвестного
из всех уравнений, кроме первого.

На
k
-
м шаге метода среди коэффициентов

при неизвестных в
уравнениях системы с номерами
i

=
k
,
.,
n

выбирают максимальный по модулю
коэффициент
Затем
k
-
е уравнение и уравнение, содержащее
найденный коэффициент, меняют местами и исключают неизвестное

из
уравнений с номерами
i

=
k

+ 1, .,
n
.

На этапе обратного хода неизвестные вычисляют в следующем поря
дке:



Пример.
Решить систему методом Гаусса с выбором главного элемента:


Решение
системы удобно оформить в виде таблицы


Используя выделенные строки, получим треугольную систему для
определения неизвестных


Откуда последовательно находим


Задания:

выполнить задание 4

ИДЗ№1.



1.8

Занятие № 8
.
Метод итераций для систем уравнений.

Цель
-

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

Рассмотрим систему линейных алгебраических уравнений (СЛАУ):

a
11

x
1

+ a
12

x
2

 …  
1n

x
n

= b
1

a
21

x
1

+ a
22

x
2

 …  
21

x
n

= b
2

……………………………………. .

a
n1

x
1

+ a
n1

x
2

…
+
a
nn

x
n

=
b
n
,

или в векторно
-
матричном виде:


Ax

=
B
, (1)

где


А 

11


12




1
21


22




2


1


2






,
B
=

1
2


,
x=

1
2



Итерационные методы решения СЛАУ используются для решения СЛАУ
большой размерности с разреженными матрицами, а также для уточнения
решения СЛАУ,
полученного с помощью прямого метода. πормулировка и
применение итерационных методов требует определенных знаний и
определенного опыта. Выбор эффективного итерационного метода решения
конкретной задачи существенно зависит от ее характерных свойств и от
арх
итектуры вычислительной машины, на которой будет решаться задача.
Поэтому никаких общих правил выбора наилучшего итерационного метода
решения не существует. Метод простой итерации приведен здесь как
иллюстрация действия механизма вычисления решения на осно
ве
итерационной процедуры.

Суть метода состоит в следующем. От системы уравнений вида (1)
переходят к системе уравнений


x
=
D
х  С


(2)

Например, от системы у
равнений


а
11
х
1
а
12
х
2
а
13
х
3
в
1



а
21
х
1
+

а
22
х
2
а
23
х
3
в
2



(3)



а
31
х
1
а
32
х
2
а
33
х
3
в
3

можно перейти к виду (), выразив из первого уравнения х
1
, из второго
-

х
2
, из
третьего
-

х
3
:


х
1
=
-

а
12

11
х
2

-

а
13

11
х
3
в
1

11


х
2
=
-

а
21

22
х
1

-

а
23

22
х
3
в
2

22




(4)


х
3
=
-

а
31

33
х
1

-

а
32

33
х
2
 в
3

33

Приведение исходной системы уравнений в виду () можно осуществить
различными способами. Например, в СЛАУ (3) из первого уравнения можно
выразить х
2
, из второго
-

х
1
, из третьего
-

х
3
и, переставив уравнения для
сохранения порядка следования переменных в векторе решения х, снова прийти
к виду (). Естественно, что матрица
D

и вектор с будут уже иными. Возможны
и другие способы преобразования исходных уравнений.

После преобразования (1
) к виду () назначается нулевое приближение
решения х

(0):



х
1

(0)

х

(0)
 х
2

(0)




х
3

(0).

Если приблизительно известны значения х
i

вектора решения х, то они
выбираются в качестве нулевого приближения, если нет, то в
качестве вектора
х

(0)
выбирается любой вектор, например х

(0)
О.

Первый шаг итерационного процесса состоит в вычислении приближения
х
(1)
:

х

(1)
=
Dx

(0)
С.

Например, назначив х

(0)
и подставив его в систему уравнений (4),
получим:

х
1

(1)
=
-

а
12

11
х

(0
)
-

а
13

11
х
3

(0)
в
1

11

х
2

(1)
=
-

а
21

22
х
1

(0)
-

а
23

22
х
3

(0)
в
2

22

х
3

(1)
=
-

а
31

33
х
1

(0)
-

а
32

33
х
2

(0)
в
3

33.

Далее вычисляем:


х

(2)
=
Dx

(1)
+
C

х

(к)
=
D
х


-
1)
С

и т.д.

Достаточное условие сходимости метода итерации заключается в
следующем, если норма матрицы
D

(обозначается ║
D
║) меньше 1, то система
уравнений () имеет единственное решение х
*
и итерации сходятся к этому
решению со скоростью геометрической прогрессии. Ины
ми словами, если


D
║<1,


(5)

то

lim
к



х

(
к
)



х



= 0

и выполняется тождество

х
*
=
D
х
*
С.

В качестве нормы матрицы
D

используются нормы ║
D

1

или ║
D


:


D

1

=
max
1



|

d
ij

|
=
1
,


D


=
max
1



|

d
ij

|
=
1
|.

Аналогично вводятся нормы вектора х:

║х║
1
=

|

x
i

|
=
1

║х║

=
max
1


|
x
i
|.

Из условия сходимости (5) ясно, что не всякое преобразование исходной
системы к виду () позволит получить решение
уравнения на основе
итерационного процесса, а только такое, которое обеспечит выполнение
условия ║
D
║<1. Важно иметь в виду, что при выполнении этого условия
итерационный процесс сходится для любого начального приближения х

(0)
и
выбор х

(0)
=
O

диктуется пр
осто соображениями удобства назначения х

(0).

Условие (5) выполняется, если модули диагональных коэффициентов для
каждого уравнения системы больше суммы модулей недиагональных
коэффициентов, т. е.



Если задана допустимая погрешность вычислений
Δ
, то для
оценки
погрешности к
-

го приближения широко используется следующее неравенство:


║х

(к)
-

х*║≤║
D
║ ∕ (1
-

D
║) •║х

(к)
-

х


-
1)
║<
Δ



(6)

Из этого неравенства следует критерий окончания итерационного
процесса



║х

(к)
-

х


-
1)
║ < (1
-
║ D║) •¨ ∕ ║D║ (7)

Каждый раз при вычислении очередного приближения х

(k)
проверяется
выполнение неравенства (7).

Выполнение неравенства (7) означает выполнение неравенства

║х

(к)
-

х
*
║ < ¨

и, следовательно, прекращение итерационного процесса.

Метод Зейделя.

Метод Зейделя

представляет собой некоторую модификацию метода
итераций. Основная его идея заключается в том, что при вычислении (
k

+ 1)
-
го
приближения неизвестной
x
i

учитываются уже вычисленные ранее (
k

+ 1)
-
е
приближения неизвестных
x
1
,
x
2
, …,
x
i

-

1
.

Пусть
для системы (1)
получена эквивалентная система


Выберем произвольно начальные приближения корней
.
Далее, предполагая, что
k
-
ые приближения

корней известны, согласно
Зейделю будем

строить (
k

+ 1)
-
е приближения корней по формулам:


Метод прогонки.

Метод прогонки является модификацией метода Гаусса для частно
го
случая разреженных систем
-

системы уравнений с трѐхдиагональной
матрицей. Такие системы получаются при моделировании некоторых

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

Запишем систему уравнений в
виде:


На главной диагонали матрицы этой системы стоят элементы

над ней
-

элементы
под ней
-

элементы


При этом обычно все
коэффициенты
b
i

не равны нулю.

Метод прогонки состоит из двух этапов
-

прямой прогонки (аналога
прямого хода метода Гаусса) и
обратной прогонки (аналога обратного хода
метода Гаусса). Прямая прогонка состоит в том, что каждое неизвестное х
i

выражается через

X
I
+1


с помощью прогоночных коэффициентов

AA
I
,


ВВ
I
:

(1)

Из первого уравнения системы найдѐм


С другой стороны, по формуле

(1)


Приравнивая
коэффициенты в обоих выражениях для х
1
, получаем


(2)

Из второго уравнения системы выразим
х
2

через
х
3
,
заменяя
x
i

по формуле
(1):


Отсюда найдѐм


или






Аналогично можно вычислить прогоночные коэффициенты


для любого

номера

i
:


(3)

Обратная прогонка состоит в последовательном вычислении неизвестных
x
i
.

Сначала нужно найти
х
n
.

Для этого воспользуемся выражением (1) при
i
=
n
-
1

и последним уравнением системы. Запишем их:


Отсюда


Далее, используя формулы (1) и выражения для прогоночных
коэффициентов (), (3), последовательно вычисляем все неизвестные
x
i

для
i
=
n
-
1,
n
-
2,.. .,1.

При анализе алгоритма метода прогонки надо учитывать возможность
деления на нуль в формулах (3). Можно
показать, что при выполнении условия
преобладания диагональных элементов, т.е. если


причѐм хотя бы для одного значения
i

имеет место строгое неравенство,
деления на нуль не возникает, и система имеет единственное решение.

Приведѐнное условие преобладания

диагональных элементов
обеспечивает также устойчивость метода прогонки относительно погрешностей
округлений. Последнее обстоятельство позволяет использовать метод прогонки
для решения больших систем уравнений. Заметим, что данное условие
устойчивости прог
онки является достаточным, но не необходимым. В ряде
случаев для хорошо обусловленных систем метод прогонки оказывается

устойчивым даже при нарушении условия преобладания диагональных
элементов.

Задания:

выполнить задание 5 ИДЗ№1.


1.9

Занятие №
9
.
Равномерное

приближение функций многочленами.
Интерполяционный многочлен Лагранжа.

В основе большинства численных методов математического анализа
лежит подмена одной функции
f
(
x
)

(известной, неиз
вестной или частично
известной) другой функцией
φ
(х)

близ
кой к
f
(
x
)

и
обладающей «хорошими»
свойствами, позволяю
щими легко производить над нею те или иные
аналитические или вычислительные операции. Будем называть такую подмену
аппроксимацией
или просто
приближением
функции
f
(
x
)

функцией
φ
(х).

Для
того, чтобы построить какую
-
то разумную теорию таких приближений и
предложить конкретные способы получения аппроксимирующих функций
φ

(х)

по заданным тем

ил
и иным образом аппроксимируемым функциям
f
(
x
),

предва
рительно следует ответить на ряд вопросов.

1)

Что известно о функции
f
(
x
)
?
Задана ли она своим ана
литическим
выражением или таблицей своих значений, какова степень ее гладкости и
доступны ли значения ее производных, как расположены точки в
интересующей части области опреде
ления
f
(
x
)
, где известны ее значения, и
можно ли их зада
вать по своему усмотрению, и т.п.

2)
Какому классу {семейству) функций должна принадле
жать функция


φ
(
х
)?

Какие дополнительные требования предъ
являются к
φ
(
х),

выделяющие
ее
из

заданного класса?

3)

Что понимать под близостью между
f
(
x
)

и

φ(х);

ина
че, какой
принять
критерий согласия
между ними?
Говоря
язы
ком функционального
анализа, по метрике какого пространства должно быть малым расстояние
между
f
(
x
)

и
φ(х
)?

Как видим, задача аппроксимации функции
f
(
x
)

функци
ей
φ(х
)

состоит в
построении для заданной функции
f
(
x
)

такой функции
φ(х
),

что








(1)

причем левая часть приближенного равенства (1) должна быть обусловлена
ответами на
вопросы первой группы, правая часть


второй группы, а ответ на
вопрос 3) должен уточнить значе
ние связывающего
f
(
x
)

и
φ(х
)

символа «

.

Прежде всего, определимся с ответом на второй вопрос.
До
говоримся
использовать в качестве аппроксимирующих функ
ций
φ(х)

только многочлены
или функции, составленные из многочленов
в таком случае будем говорить о
полиномиальной аппроксимации
или
кусочно
-
полиномиальной
аппроксимации

соответственно.

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

привлекательны тем, что они
являются
линейными
ф
ункциями своих параметров (коэффициентов), и их вычисление

св
о
дится

к
выполнению конечного числа простейших арифмети
ческих операций


сложения и умножения.

Будем считать, что аппроксимация функции
f
(
x
)

производится с
помощью многочленов степени
п

N
0
. Тогда в зависи
мости от выбора критерия
согласия и, в частности, от количества
точек согласования
f
(
x
)

с
φ(х
)

(будем
называть их
узлами)
,

т.е. точек, в которых известна информация об
f
(
x
)

и,
возможно, ее производных, можно рассмотреть разные конкретные способы
аппроксима
ции.

Таблица 1



И
нтерполяционный многочлен Лагранжа

Пусть в точках
таких, что
известны значения
функции
у
=
f
(
x
),

т.е. на отрезке
[а,

b
]

задана
табличная (сеточная) функция







(2)

πункция
φ
(х)

называется
интерполирующей
или
интер
поляционной
для
f
(
x
)

на
[а,

b
],

если ее значения
в
заданных точках х
0
, х
1
,


,

х
п
,

называе
мых
узлами интерполяции,
совпадают с заданными значениями

функции
f
(
x
),

т.е. с
у
0
,

у
1
,

,

у
n

соответственно.
Геометри
чески факт
интерполирования означает, что график функции
φ
(
х)

проходит так, что, по
меньшей мере, в
n
+1

заданных точ
ках он пересекает или касается графика
функции
f
(х)

(рис.1).


Рис.1.
Геометрическая интерпретация задачи интерполирования

Легко представить, что таких графиков
φ(
х),

проходящих через заданные
точки, можно изобразить сколько угодно, и они могут отличаться от графика

f
(
x
)

сколь угодно сильно, если не накладывать на
φ(
х)

и
f
(
x
)

определенных
ограничений.

Б
удем считать, что интерполяционная функция
φ
(х)
есть много
член
степени
п.

Тогда
задача интерполяции,
точнее, полино
миальной,
алгебраической или параболической интерполяции (поскольку график любого
многочлена называют параболой
)

формулируется так:

для фун
кции
f
(
x
),

заданной таблицей
(
2),
найти много
член
Р
п
(х)

такой,
что выполняется совокупность
условий интерполяции








(
3)

Найти многочлен
Р
п
(
х)


это значит, учитывая его кано
ническую форму








(
4)

найти его
n
+
1

коэффициент
.

Для этого имеется как раз
n
+
1

условие
(
3). Таким
образом, чтобы многочлен
(
4) был интерполяционным для функции
(
), нужно,
чтобы его ко
эффициенты удовлетворяли сис
теме уравнений


Из курса алгебры известно, что определитель этой линейной сис
темы (так
называемый определитель Вандермонда) отличен от нуля, т.е. решение этой
системы существует и единственно. Од
нако практическое построение
интерполяционного многочлена
таким путем малоэффективно. Поэтому
изберем другой путь.

Будем строить многочлен
п
-
й

степени
L
n
(
x
)

в виде линейной
комбинации

многочленов
п
-
й

же степени
l
i
(
x
)
(
i

=

0,1,...,

п)
.

Для того,
чтобы такой многочлен был интер
поляционным для функции
f
(
x
),
достаточно
зафиксировать в качестве коэффициентов с
i
, этой линейной комбинации
заданные в табл. значения
y
i

=
f
(
x
i
)
,

а от базисных многочленов
l
i
(
x
)

потребовать выполнения условия









(5)

В таком случае для
многочлена

в каждом узле
x
j

(
j
=0
,

1
,


,
п
), в силу (5), справедливо
.


т.е. выполняются условия интерполяции (3).

Чтобы конкретизировать базисные многочлены
l
i
(
x
)
, уч
тем, что они
должны удовлетворять условиям (5). Равенство нулю
i
-
го многочлена во всех
узлах, кроме
i
-
го, означает, что
l
i
(
x
)

можно записать в виде


а коэффициент
А
i


этого представления легко получается из со
держащегося в (5)
требования
l
i
(
x
)
1. Подставляя в выраже
ние
l
i
(
x
)
значение хх
i

и
приравнивая

результат единице, по
лучаем


Таким образом,
базисные многочлены Лагранжа
имеют

вид


а искомый
интерполяционный многочлен Лагранжа
есть







(6)

Пример.
Рассмотрим квадратичную интерполяцию функции
y
=
sinx

на
отрезке [0,
π
/2]
по ее трем значениям:
sin
0 = 0,
sin
π
/4
=

,
sin
π
/2=1.


По формуле (8) строим многочлен Лагранжа второй степени


или после преобразований



Подставим в полученный интерполяционный многочлен кон
трольную
точку
Получим приближенное значение
отличающееся от значения
sin
π
/6 0.5 на вели
чину ≈0.017.

Интерполяционная схема Эйткена

Пусть функция
f
(
x
)
и расположение узлов
x
0
,
x
1
, ...
,
x
n

на отрезке интерпо
-
ляции
[
a
,
b
]
таковы, что имеет место сходимость
процесса интерполяции
.
Пусть
требуется найти не общее выражение
L
n
(
x
),
а лишь его значения при
конкретных
x
, т.е. решается частная задача вычисления отдельных
приближенных значений функции
f
(
x
)
с помощью вычисления соот
-
ветствующих им значений интерполяцио
нного многочлена Лагранжа
L
n
(
x
).
Для
построения эффективного способа решения такой частной задачи интерполя
-
ции учтем следующее:

1)

использование многочлена Лагранжа в виде (6) неудобно из
-
за его гро
-
моздкости, что ведет к большим вычислительным затратам;

2)

заранее не известно, многочлены какой степени нужно использовать
для ин
терполирования данной функции с требуемой точностью. А постепенное
улучшение точности за счет вычислений
L
n
(
x
)
с большим показателем степе
ни
n

невыгодно, т.к.
L
n
-
1
(
x
)
плохо перестраив
ается в
L
n
(
x
);

3)

функция
f
(
x
)
задается таблицей своих приближенных значений.
Процесс сходимости
L
n
(
x
)
к
f
(
x
)
при больших значениях
n

будет нарушен
влиянием на результат исходных ошибок.

Построим вычислительную схему для получения приближенного значе
-
ния сеточной функции
f
(
x
)
в заданной точке
, в основу которой будет по
-
ложена интерполяция Лагранжа на сетке узлов

Организация вычис
-
лений по этой схеме будет иметь итерационный характер. Каждый шаг заклю
-
чается в вычислении некоторого определителя втор
ого порядка.


Пусть даны две точки на кривой
Построим функ
-
цию


Т.е.

совпадает с интерполяционным многочленом Лагранжа первой
сте
пени, построенным по двум данным точкам (сравните с 6).

Построим через определитель функцию

для точек


Она тоже является многочленом Лагранжа первой степени, построенным
по двум точкам

Если на кривой
y
=
f
(
x
)
заданы три точки
(
x
o
,
y
o
), (
x
1
,
y
1
), (
x
2
,
y
2
),
то, исполь
-
зуя введенные линейные функции

образуем новую функцию:


Эта функция есть многочлен второй
степени, решающий задачу пара
-
болической интерполяции по трем точкам
(
x
0
,
y
0
), (
x
1
,
y
1
), (
x
2
,
y
2
).
Но такой мно
-
гочлен единственный, следовательно,

где

-

многочлен
Лагранжа.

Продолжая процесс рассуждения, получим рекуррентное задание после
-
довательности
интерполяционных многочленов Лагранжа, которое составляет
суть интерполяционной схемы Эйткена:







(7)

где
i
1, ,…,
n

и по определению
P
0
(
x
)=
y
0
,
P
1
(
x
)=
y
1
.


Схема Эйткена легко реализуется на ЭВМ. Организация вычисле
ний по
формуле (7) должна быть такова, что если заранее неизвестна степень интер
-
поляционного многочлена, который нужно использовать для вычисления
y
(
x
),
то должно происходить постепенное повышение степени
k

интерполирующих
ее многочленов за счет подключен
ия новых узлов. Счет ведется до тех пор, по
-
ка идет уточнение приближенного значения
y
(
x
).

Об этом можно судить по уменьшению величины
при
увеличении
k

и подходящем фиксировании
i
.

Пример.

Пусть некоторая функция
y
=
y
(
x
)
задана таблицей своих
значений,
округленных до двух знаков после запятой:


Рассмотрим процесс вычисления двух значений этой функции по схеме
Эйткена в точках: а)
б)

Результаты промежуточных вычисле
ний
(в которых один знак запасной) сведем в табл. 3, 4. Числа в столбцах,
помеченных п
осредством
представляют собой значения многочленов
Лагранжа
k
-
ой степени, построенных по узлам от
i
-
го до
(
i
+
k
)
-
го

рекуррентно
по формуле:





(8)

где
k
1, , …,

в соответствии с интерполяционной схемой Эйткена.

Порядок получения этих значений показан проставленными в каждой
клетке номерами.

Таблица 3.

Вычисление по схеме Эйткена значения
y
(0.1)



Таблица 4.

Вычисление по схеме Эйткена значения
y
(0.25)


Процесс вычисления значений многочленов Лагранжа ведется до тех пор,
пока идет уточнение приближенного значения
т.е. уменьшается ве
личина
)

при увеличении
k

и подходящем фиксировании
i
.

Например, для подсчета приближенного значения данной функции в точ
-
ке

расположенной между узлами
x
o
=0.0
и
x
i
=0.2,
целесообразно в каче
-
стве основной последовательности значений многочленов Лагранжа брать
строку таблицы
3,
соответствующую значению
i
=0,
т.е. числа


Вычислив разности между последующими и предыдущими числа
ми этой
строки, а именно:

0.005

0.004

0.010,

видим, что дальнейший счет бессмыслен; разность перестала уменьшаться.
Т.е. по данной информации о функции
y
(
x
)
более точное значение
y
(0.1),
чем
y
(0.1)=1.001,
получаемое с помощью
L
3
(0.1),
найти не удастся.

В случае б) вычисление по схеме Эйткена дает следующий результат:
y
(0.5)≈1.038,
полученный с помощью
L
3
(0.25).

Задания:

выполнить задание 1 ИДЗ№.

1.10

Занятие № 1
0
.
Интерполяционный многочлен Ньютона.

Многочлены
Чебышева.


К
ОНЕЧНЫЕ РАЗНОСТИ

Зададимся целью пр
идать интерполяционной формуле бо
лее простой вид,
подобный виду широко используемой в матема
тическом анализе формулы
Тейлора. Если в интерполяционном многочлене Лагранжа (6) все слагаемые
однотипны и играют одинаковую роль в образовании результата, хотел
ось бы
иметь такое представление интерполяционного многочлена, в котором, как и в
многочлене Тейлора, слагаемые располагались бы в по
рядке убывания их
значимости. Такая структура интерполяцион
ного многочлена позволила бы
более просто перестраивать его ст
епень, добавляя или отбрасывая удаленные от
начала его запи
си члены.

Поставленной цели будем добиваться сначала для несколь
ко суженной
постановки задачи интерполяции. А именно, будем считать, что
интерполируемая функция
у

=
f
(
x
)

задана своими значениями
на
системе
равноотстоящих узлов
т.е. таких, что любой узел х
i
, этой
сетки
можно представить в виде

х
i

 х
0

+
ih
,

где
i
0, 1, …,

п,

a

h

0



некоторая постоянная величина, называемая
шагом
сетки
(таблицы).


Прежде чем строить желаемые интерполяционные форму
лы, рассмотрим
элементы теории конечных разностей.

Вычитая из каждого последующего члена конечной после
довательности
из
n
1 чисел
предыдущий, образу
ем
п
конечных разностей
первого порядка


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

ко
нечных разностей второго
порядка,
или
вторых разностей:


Этот процесс построения разностей может быть продолжен, и весь он,
очевидно, описывается одной рекуррентной
формулой,

выражающей
конечную
разность
к
-
го порядка
через раз
ности


-

1)
-
го порядка:


где
к

=

1
,

, …,
п

и

В

некоторых случаях требуется знать
выражения конеч
ных разностей
непосредственно через значения функции
лежащей в их основе. Для
нескольки
х первых порядков разно
стей их можно получить прямой
подстановкой:


Подметив закономерность в коэффициентах рассмотренных
представлений конечных разностей, записываем общую формулу








(
1
)

которая может быть строго обоснована методом математической индукции и
которая напоминает биномиальное разложение для


-

1)
к
.

Привлекая определение производной, можно обнаружить прямую
связь
между конечными разностями и производны
ми.
А
именно, если учесть, что



то можно сказать, что при малых
h

имеет место приближенное равенство


т.е. первые разности характеризуют первую производную функ
ции
f
(
x
),

по
значениям которой они составлены. Пользуясь этим, имеем для вторых
разностей:


т.е.

и, вообще,








(2)

Таким образом, на конечные разности можно смотреть как на некоторый
аналог производных. Отсюда справедливость многих их свойств, одинаковых
со свойствами производных.

Отметим лишь
простейшие свойства конечных разностей:

1)

конечные разности постоянной равны нулю
;

2)

постоянный множитель у функции можно выносить за

знак конечной разности.

3)

конечная разность от суммы двух функций равна сумме

их конечных разностей в одно
й и той же точке.

Свойства  и 3 характеризуют операцию взятия конечной разности как
линейную операцию.

Учитывая роль, которую играют многочлены в теории интерполирования,
посмотрим, что представляют собой
конечные разности многочлена.

Поскольку многочлен
в своей канонической форме есть ли
нейная
комбинация степенных функций, положим сначала
ух
п
.

Используя
биномиальное разложение
п
-
й

степени дву
члена, получим:



т.е. первая конечная разность степенной функции
ух
п

есть многочлен степени
п
-
1

со старшим членом

Если взять теперь конечную разность от
функции








(3)

то, в силу линейных свойств
Δу,

можно записать


Первое слагаемое в этой сумме, как выяснено, есть многочлен (
n
-
1)
-
й
степени, второе, аналогично,


многочлен степе
ни
п
-
2,

и т.д. Следовательно,
первая конечная разность много
члена (3) в точке
х
с шагом
h

есть тоже
многочлен со стар
шим членом
a
0
nhx
n
-
1
, вторая конечная разность


многочлен
со старшим членом
a
0
n
(
n
-
1)
h
2
x
n
-
2
,…, к
-
я

разность


многочлен со
старшим членом
а
0
п(п
-
1)...(
n
-
к
+ 1)
h
k
x
n
-
k
.

При
к  п
получаем
постоянную

разность
п
-
го порядка


для многочлена (3); конечные разности более высоких поряд
ков, естественно,
равны нулю.

Итак, главный вывод из предыдущих рассуждений:
п
-
е

ко
нечные
разности многочлена
п
-
й

степени постоянны, а (
п
1
)
-
е
и все последующие
равны нулю.

Более важным для понимания сути полиномиального интерполирования
является утверждение, обратное сделанному выше
выводу. А именно, что
если
конечные разности п
-
го порядка некоторой функции у
у(х)

постоянны в
любой точке х при различных фиксированных шагах
h
,
то эта функция
у(х)

есть многочлен степени п.

Для функции
у

=

f
(
x
)
y

заданной таблицей своих значений

в
узлах

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



Прим
ер
.
Составим таблицу конечных разностей для функции
у
=
xln
2
x

по
ее значениям, вычисленным с тремя знаками после запятой в точках

x
j
=
0.4 +
0.2
j
,
где
j
=0, 1,
...,

10
.


Абсолютные величины конечных разностей сначала убывают с
увеличением их порядка, а затем

начинают увеличи
ваться. Это
типичное

поведение конечных разностей при ограничен
ной точности задания значений
сеточной функции.

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


достаточно гладкая, то сначала

с

увеличением
k

в си
лу
упомянутой связи (). Когда эти величины становятся дос
тато
чно малыми,
большую роль начинают играть продукты взаимодействия исходных ошибок
округления (так называемый
шум округлений).


Что происходит с одной отдельно взятой ошибкой величи
ны
ε
у значения
y
i
, можно проследить по таблице. Как видим, с ростом порядка
разностей она
«расползается» по таблице и уве
личивается по абсолютной величине.


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

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

где
ε



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

Вспоминая о том, что многочлен
k
-
й степени имеет
k
-
е

разности
постоянными, а все последующие


нулевыми, прихо
дим к заключению, что
если
k
-
е

разности таблицы конечных разностей некоторой

функции
практически постоянны, то эта функция ведет себя в рассматриваемой
области, как многочлен
k
-
й

степени;
эту степень и следует применять для
интерполиро
вания с наибольшей для данных реалий точностью.


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

0.0005



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

узла

X
0
=

0.4
, нужно применять квадрати
чную
интерполя
цию.

К
онечноразностные интерполяционные формулы

Пусть функция
у
f
(
x
)

задана на сетке равноотстоящих узлов
х
i
х
0
+
ih
,
где
i
 0, 1, …,
п,

и для нее построена таблица конечных разностей.

Будем строить интерполяционный много
член
P
n
(х) в форме







(4)

Его
п
1
коэффициент
а
0
,

а
1
, ….,а
п

будем находить последова
тельно из
п
+1
интерполяционных равенств


А именно, полагая
i
0, т.е. х
х
0

в (4) имеем

а по условию
интерполяции
Р
п

0
)
у
0
,
следова
тельно,
а
о
у
0
.

Далее, при
i

 1 аналогично получаем равенство


в которое подставляем уже найденное значение

Разре
шая это равенство
относительно
а
1

и используя обозначение ко
нечной разности, получаем


Следующий шаг, при
i

=
, дает:



Полной
индукцией можно показать справедливость выражения









(5)

Подставляя найденные коэффициенты (4), получаем многочлен







(6)

который называют
первым интерполяционным
многочленом Ньютона.

Учитывая, что каждое слагаемое многочлена (6), начи
ная со второго,
содержит множитель
X
-
X
0
, естественно предпо
ложить, что этот многочлен
наиболее приспособлен для интер
полирования в окрестности узла х
0.
Будем
называть узел
X
0

базовым
для многочлена (6), и упростим (6) введением новой
переменной
q

равен
ством
или (что то же) равенством х
=
X
0
+
qh
. Так как
h

при любых
i
=0, 1,
...,
n


то в результате подстановки этих разностей в (6) приходим к
первой
интерполяционной формуле Ньютона
в виде







(7)

где обозначение
указывает не только на
п
-
ю

степень многочлена, но
и на базовый узел
x
0

и связь переменных х и
q
.

Первая формула Ньютона (7) обычно применяется
при значениях |q|<1, а
именно,
для интерполирования вперед,
(при

т.е. при
)

и
экстраполирования назад
(при х < х
0
, т.е. при
q

0
).

Так как
,

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


номеров узлов нет, то за
базовый для формулы (7) узел
x
0

можно принима
ть узел, ближайший к
заданной фиксированной точке
х,

если за ним имеется достаточное число узлов
для по
строения необходимых для (7) разностей. Поскольку в первой формуле

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

Учет этого обстоятельства приводит к потребности в сим
метричной, в
определенном смысле, для (7) формулы, которая была бы пригодной для
интерполирования в конце таблицы. Для

этого, в отл
ичие от (4), форма
интерполяционного многочлена
Р
п
(х)

берется такой, которая предусматривает
поочередное под
ключение узлов в обратном порядке: сначала последний,
потом предпоследний и т.д., т.е.


Коэффициенты

этого многочлена находятся анало
гично
тому, как они находились для многочлена (4), только здесь подстановка
узловых точек вместо
х

и рассмотрение ин
терполяционных равенств
производится тоже в обратном поряд
ке. Полагая
имеем:


и т.д. В общем случае


Таким образом
,

получаем
второй и
нтерполяционный мно
гочлен
Ньютона






(8)

в котором базовым является узел
х
п

и коэффициенты которого определяются
конечными разностями, расположенными на вос
ходящей от
у
п

диагонали.


Положим в (8)
хх
п
+
qh
,
иначе, введем новую переменную
и
преобразуем к ней входящие в (8) разности:



В результате приходим ко
второй интерполяционной формуле
Ньютона
вида







(9)

Ее также целесообразно использовать при значениях
|
q
|
1

т.е.

в
окрестности узла х
n

для интерполирования назад
(при
q

(
-
1, 0)) н
экстраполирования вперед
(при
q
� 0).

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

Будем считать, что узел
X
0

расположен в середине табли
цы, и

нумерация
остальных узлов производится, начинаясь с х
0
, с использованием как
положительных, так и отрицательных ин
дексов, т.е. считаем
х
i
=
х
0
+
ih
,

где
i
= 0,
1,
2,....

Тогда цен
тральная часть таблицы конечных разностей будет
проиндекси
рована так, как эт
о показано в таблице.

Все подчеркнутые в ней конечные разности (находящиеся с
в одной
строке и на полстроки выше и ниже) называются
центральными разностями.




Интерполяционный многочлен ищем в форме







(10)

предполагающей постепенное подключение узлов х
i
: сначала при
i
=
0,

затем
при
i
1, потом при
i
=
-
1 и т.д., т.е. с двух сто
рон от х
0
. При этом здесь и далее
не будем фиксировать степени многочленов и не будем стремиться выписывать
общие
и,

тем бол
ее, последние члены таких многочленов. Как и в предыдущих
случаях, коэффициенты
находим один за другим
последовательной подстановкой в Р(х)
и

в интерполяци
онные равенства
значений


и т.д. Введя новую переменную
и выразив через нее разности
для всех
i

 0, “ 1, “ , ..., в результате подстановки этих разностей
и выражений коэффициентов в шаб
лон (10), приходим к
первой
интерполяционной формуле Гаусса








(11)

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

Таблица

1


Совершенно аналогично, подключая узлы в другом порядке по
сле

x
0

сначала предшествующий, затем последующий и т.д., т.е.
можно построить
вторую интерполяционную формулу Гаусса







(12)

использующую
верхние

центральные разности (подчеркнутые в табл.1
пунктирной линией).

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

Так, полусумма первого и второго интерполяционных мно
гочленов
Гаусса после преобразований приводит к формуле









(13)

называемой
интерполяционной формулой Стирлинга
К.

Если же взять полусумму второго интерполяционного мно
гочлена Гаусса
и такого же многочлена, но с нижними индекса
ми, увеличенными на единицу
(т.е. с базовой точкой
х
1

вместо
X
0
), то придем к
интерполяционной формуле
Бесселя







(14)

В

последней формуле обращает на себя внимание тот факт, что она
сильно упростится, если в нее подставить значение
q
=1/2
,

соответствующее
значению аргумента

Этот частный случай формулы Бесселя называют
формулой ин
-
терполирования на середину:






(15)

Итак,
если точка х, в которой нужно найти приближенное значение
таблично заданной функции
f
(х),
находится в начале или в конце таблицы,
применяется соответственно первая
(7)
или вторая
(9)
формулы Ньютона с
таким выбором ба
зовой точки, чтобы значение |
q
| было как можно меньше.
Если точка х находится в середине таблицы, то всегда можно за
фиксировать

точку
x
0

в таблице центральных разностей так, чтобы

либо было
по модулю меньше
0.25
и тогда при
менять интерполяционную формулу
Стирлинга
(13),
либо что
бы
q


[0.25, 0.75]
и использовать формулу Бесселя
(14).


Пример
.
Пусть требуется для функции
у

=

f
(
x
),

заданной в предыдущем
примере таблицей нескольких своих значений с тремя знаками после запятой,
найти приближенные значения: а) f(0.5); б)
f
(1.22);
в)
f
(0
-
5); г)
f
(1.94); д)
f
(2.5),
записав предварительно соответст
вующие каждому случаю интерполяционны
е
формулы.

Для решения поставленной задачи учитываем, что значения
f
(
x
) заданы
на сетке
равноотстоящих

узлов, поэтому здесь можно применить
конечноразностную интерполяцию. При этом бу
дем пользоваться уже
составленной таблицей конечных разностей и прове
де
нным ранее ее анализом
на выявление оптимальной степени многочле
на. Для случаев б
-
д фиксируем
вторую степень, для а


третью. В каждом случае, т.е. для конкретного
значения аргумента, выбираем базо
вый узел, подсчитываем значение
вспомогательной переменно
й
q

и, в за
висимости от положения базового узла и
значения
q
, пользуясь представ
ленными
в таблице числами, записываем
требуемую интерполяционную формулу. Подстановка в нее значения
q

приводит к искомому значению

а) При
X
0.5 (начало таблицы) полагаем
х
0
=0.4;

тогда

Соответствующую интерполяционную формулу д
ля
аппроксимации
f
(
x
) при
х
0.4 + 0.2
q

с
q

(
-
1,1)

записываем, глядя на первую
интерполяционную формулу Ньютона (7) и таблицу значений:


Отсюда получаем иском
ое значение


б)

Точка
х
1. находится в средней части таблицы. Поэтому здесь

целесообразно применить формулу Стирлинга или Бесселя. Полагая

и
найдя
останавливаемся на формуле Стирлинга (13), которая
в данном случае имеет вид



и при
q

 0.1 приводит к

искомому значению


в)

Здесь, очевидно, напрашивается применение формулы (15)

интерполирования на середину. Полагая
имеем:


г)

Глядя на положение точки х1.94 в заданной системе узлов
,

видим,
что для вычисления
f
(1.94) также возможно применение

центральных интерполяционных формул. Положив
х
0
=1.8

и вычислив


на основе (
14
) записываем интерполяцион
ную формулу
Бесселя


Из нее получаем


д)

Точка
х
 .5 расположена за последним узлом, поэтому для
экстраполяции
f
(
x
)

здесь однозначно следует применить вторую
интерполяционную формулу Ньютона (9). Считая
х
п

.4 (индекс
п
здесь ис
-

пользуется условно, без придания ему конкретного значения), записываем

формулу экстраполяции


откуда при
находим



О
пределение и свойства
многочленов Чебышева

Многочленом Чебышева
называ
ется функция





(15)

где



πункция
Т
п
(х),

представлен
ная с помощью тригонометрических
функций, на самом деле является многочленом при любом
n
0, 1, ,.... и
справедливо равенство




(16)

πормула (16) определяет при
n
=
1, , 3,... последователь
ность функций

начинающуюся с Т
0
(х)1,
T
1
(
x
) =
x
,
рекуррентно;

при этом нужно
иметь в виду, что здесь
x

[
-
l
, 1], как и в (15).

Подставляя
в (16) заданные начальные члены последова
тельности

n
(х)),

найдем несколько ее последующих членов:


и т. д.

Графики нескольких многочленов Чебышева (с первого по четвертый)
изображены на рис. 1.


Рис. 1
.
Графики многочленов


Анализ рекуррентной
формулы (16) позволяет считать очевидными
следующие факты:

1)

все функции
Т
n
(х),

определенные в
(15),

являются многочленами
при любом натуральном
n
;

2)

степени этих многочленов возрастают с увеличением
n
, причем
старший член многочлена
Т
n
(х)

равен 
n
-
1

х
n
;

3)
многочлены
Т
n
(х)

при четных
n

выражаются через степенные
функции только четных степеней, при нечетных


только нечетных.

Наряду с многочленами Чебышева
Т
n
(х),

часто используют многочлены,
получаемые из
Т
n
(х)

делением на старший коэф
фициент, т.е.





многочлены со старшим коэффициентом 1. Будем называть их
нормированными многочленами Чебышева.

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

1.
Многочлен Чебышева
Т
n
(х)

имеет на отрезке

[
-
1, 1]
ровно
n

раз
-
личных действительных корней; все они задаются форму
лой





(17)

2.
Корни многочленов Чебышева перемежаются с точками их
наибольших и наименьших значе
ний, равны
х соответственно

+1
и

-
1 и

для
. а именно, при j0, 1, …,
n

имеют место экстремумы







(18)

3.
(
теорема Чебышева).
Из всех многочленов степени
n

со старшим
коэффициентом

1
нормиро
ванный многочлен
Чебышева

наименее
уклоняется от нуля на отрезке

[
-
1, 1].

Это свойство означает, что
среди всех многочленов степени
n

вида






(19)

именно нормированный многочлен Чебышева

минимизиру
ет
максимальное расстояние от
графика многочлена при

х

[
-
1,1]
до оси абсцисс.


И
нтерполяция по Чебышевским узлам

Желая, чтобы интерполяционный многочлен Лагранжа
L
n
(
х)
в целом
хорошо приближал функцию
у

f
(
x
)

на отрезке [а,
b
], поставим вопрос: как
расположить на нем
n
+1
узлов
интерполяции х
i
,
(
i
 0, 1, …,
n
), чтобы при этом
минимизировать максимальную на
[а,
b
]

погрешность?

Максимальная погрешность интерполирования достаточно гладкой
функции на отрезке
[
-
1, 1]
многочленом
п
-
й

степени будет минимальной, когда

в качестве узлов интерполяции

берутся корни многочлена
Чебышева
T
n
+1
(
t
).

Будем называть их
чебышевскими узлами интерполяции.

Знание экстремальных значений многочлена Чебышева позволяет
уточнить величину максимального отклоне
ния
L
n
(
x
)

от
f
(
x
)

при таком выборе
узлов, т.е. когда точки
t
i

есть корни
a

А
именно,


Эта оценка называется
наилучшей
равномерной

оценкой погрешности
интерполяции.

Если функция
f
(
x
)

бесконечно дифференцируема на


,
b
]

и
в качестве
узлов интерполяции берутся корни многочленов Чебышева (приведенные к
отрезку
[а,

b
], то


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

Теорема

(Вейерштрасса)
.
Для любой непрерывной на

[а,
b
]

функции
f
(
x
)

най
дется многочлен
Q
n
(
x
)

такой, что


Доказано, что для любой функции
f
(
x
)

существует
единственный

многочлен

такой, который из всех многочленов
Q
n
(
x
)

степени
n

наилучшим образом аппроксимирует на [а,
b
]

функцию
f
(
x
)
,

минимизируя
максимальное расстояние между
f
(
x
)

и

Q
n
(
x
).

Этот многочлен, т.е. многочлен
такой, что


называется
многочленом наилучшего равномерного

приближе
ния
для
f
(
x
)

на
[а,
b
]

или ее
чебышевским приближением.


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

Теорема

(
Чебышева).
Многочлен

является многочленом
наилучшего равномерного приближения для функции
f
(
x
)

тогда и только
тогда, когда на
[а,

b
]

существует не менее
n

+
 точек

х
i

таких, что в них
поочередно принимаются наибольшие положительные и отрицательные
отклонения,
т.

е. поочередно разность

f

(
x
i
)
-
равна
Е
или

-

Е,
где


Эта теорема говорит о том, что максимальная ошибка аппроксимации
функции многочленом наилучшего равно
мерного приближения реализуется в
числе точек, большем, по меньшей мере, на , чем степень
многочлена, причем
знаки ошибки чередуются. Точки х
i
, в которых реализуется максимальное
отклонение многочлена
от
f
(
x
)

на
[а,

b
],

назы
ваются
точками
чебышевского альтернанса.

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

Случай А.
Пусть функция
f
(х) непрерывна на [
a
,
b
], и пусть для нее
требуется построить
многочлен наилучшего рав
номерного прибли
жения нулевой
степени.
Обозначим этот приближающий многочлен через

Он
определяется всего одним параметром:
.

Чтобы найти значение этого
па
раметра для заданной функции
f
(х), воспользуемся тем свойст
вом
непрерывной на замкнутом промежутке функции, согла
сно которому на нем

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

Пусть

Тогда совершенно очевидно, что
полагая
т.е. подменяя функцию
f
(
x
)

функцией
будем
иметь максимальное откло
нение


причем точки отрезка
[а,

b
]
,

в которых оно реализуется


это точки, где
принимаются значения
т

и
М
.

В силу непрерывнос
ти
f
(х), локальные
минимумы и максимумы должны чередоваться; по меньшей мере, два из них
определяют точки чебышевского альтернанса (с
чередованием знаков
разностей f(x)
-
+
2

)
. Поэтому не существует другой постоянной, котор
ая
приближала бы
f
(
x
)

на
[а,
и]

лучше, чем постоянная
+
2
.


Случай Б.
Пусть аппроксимируемая функция
f
(
x
)

диффе
ренцируема и
выпукла (в широком смысле) на отрезке [а,
b
],
a

аппроксимирующая ее
функция


многочлен наилучшего равномерного приближения
первой степени.
Чтобы найти его коэффициенты следует изучить раз
ность
между
f
(
x
)

и
φ(х),

т. е. функцию

Так как функ
ция
f
(
x
)

по предположению выпукла, а сдвиг на линейную
функцию
A
0

+

А
1
Х

не изменяет выпуклости, то и функция
u
(х)

выпукла на
[а,
b
]
.

Следовательно,
существует единственная точка
с

[а,
b
],

в которой
u
(х)

имеет минимум; если
u
(х)

выпукла вверх, то в точке
с

должен быть максимум

u
(х).

В любом случае, точка
с

[а,
b
]

есть точка экстремума
u
(х),

и за счет
возможности варьирования коэффи
циентов функции
φ(х)

(точнее,
коэффициента
А
1
)

можно счи
тать, что точка с является внутренней точкой
отрезка [а,
b
]
.

Потребуем, чтобы точки
а,

с

и
b

в указанной последова
тельности
составляли чебышевский альтернанс, т. е. чтобы в них последовательно
принимались значения
Е,

-

Е,

Е

или
-

Е
,

E
,

-
E
, где



максимальная погрешность
аппроксимации функции
f
(
x
)

функцией
φ
(х).

Добавляя к этим требованиям
еще необходимое условие экстремума дифференцируемой функции
u
(х)

в
точке
с

и воз
вращаясь к исходным функциям, приходим к системе четырех
уравнений относительно четырех неизвестных
A
0
,

А
1
,

Е

и
с

(из которых, в
основном, лишь
первые три неизвестные представля
ют интерес):


Получить решение такой системы в общем виде не представля
ется
возможным, поскольку неизвестная величина
с
входит в нее нелинейным
образом (в каждом конкретном случае подобная система без проблем решается
численно).

Пример
.

Построим многочлены наилучшего равномерного приближения
нулевой и первой степеней для функции
y

=
sinx

на отрезке [0,
π
/6].

Сразу заметим, что данная функция всюду дифференцируема (а значит,
непрерывна) и выпукла вверх на заданном отрезк
е. При этом


Следовательно, согласно рассмотренному выше случаю А, найдя


при
х

[0,
π
/6] можно считать, что
sin
х


0.25
с предельной погрешностью
0.25.


Далее, в соответствии со случаем Б, продифференцировав данную
функцию, составляем систему:


Из нее
последовательно находим:


Таким образом, функцию
y

=
sinx

на отрезке [0,
π
/6] можно подменить
линейной функцией
у 
0.0045 + 0.9549
x
, и наибольшая ошибка при этом не
будет превышать величины
≈0.0045.

Задания:

выполнить задание 
и 3
ИДЗ№.

1.11

Занятие № 1
1
.
Интерполяция с кратными узлами и формула Эрмита.

Рассмотрим задачу полиномиальной интерполяции функ
ции
у
f
(
x
)
в
более общей постановке.

Пусть на промежутке
расположены
m
1 не
совпадающих узлов
и пусть в этих точках извест
ны значения
у
0

=
f
(
x
0
), у
1
=
f
(
x
1
,),...,
у
т

=
f
(
x
m
) данной функции, а также некоторые ее производные (максимальный по
-
рядок производных в разных узлах различен; в каких
-
то узлах производные
могут быть вовсе неизвестны). Такие узлы будем называть
кратными узлами.
Конкретнее, будем
считать, что за
даны:







(1)

тогда
кратность
узла
x
0

считается

равной
k
0
, узла
x
1

-

k
1
, …,
узла
x
m

-

k
m
.

Предполагая, что
суммарная кратность узлов
есть







(2)

ставим задачу построения многочлена
Н
п
(х)
степени
n

(не вы
ше
п
) такого, что





(3)

где



заданные посредством (1) значе
ния функции
f
(
x
)

и ее
производных и по определению счи
тается

Многочлен

будем называть
интерполяционным многочленом Эрмита,
а
совокупность требований (3)


условиями эрмитовой интерполяции.

πормально можно считать, что нахождение такого много
члена состоит в
том, чтобы однозначно определить
п
1 коэффициен
тов
a
0
,
а
1
,..., а
п

его
канонического представления








(4)

из условий (3). В силу предположения () о суммарной кратности узлов
эрмитовой интерполяции, совокупность требо
ваний (3) можно рассматривать
как систему из
п1
уравне
ний относительно
п1
неизвестных


коэффициентов
a
k

мно
гочлена (4):


Выявление общего вида интерполяционных многочленов Эрмита
Н
п
(х)
представляет непростую задачу и требует при
влечения определенных

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

Пусть
L
m
(
x
)


интерполяционный многочлен Лагранжа, построенный по
данным
т1
значения
м
y
i

=
f
(
x

i

= 0,1,...,
т.
Б
удем пользоваться обо
значением

Так как по ус
ловию
т
заведомо не превосходит

п
, то по теореме о делении многочлена с остатком искомый многочлен Эрмита
Н
п
(х)
можно представить в виде








(5)

где
H
n
-
(
m
+1)
(
x
)


некоторый неизвестный пока многочлен степени
п
-
т
-
1.

Для построения многочлена
H
n
-
(
m
+1)
(
x
)
будем привлекать информацию о
производных данной функции, т.е. равенства

в тех узлах х
i

где
первые производные, в соответ
ствии с (1), заданы (информация о самих
значениях функции уже полностью исчерпана: в силу

для
всех х
i

от х
0

до х
m
, согласно (5), будет и

при любых
i

{ 0,1,...,

т}).

Продифференцировав равенство (5), имеем







(6)

Поскольку

в тех узлах х
i
, где по условию эрмито
вой
интерполяции справедливо

можно записать


Отсюда выражаем значения многочлена
H
n
-
(
m
+1)
(
x
)
в этих уз
лах:


Правая часть этого равенства может быть вычислена; обо
значим ее через
. Таким образом, в
ряде узлов х
i

известны зна
чения многочлена
H
n
-
(
m
+1)
(
x
)
=
, по которым этот многочлен однозначно восстанавливается обычной
лагранжевой интерполя
цией, если в условиях (1) не содержится производных
поряд
ка, выше первого (т.е. нет ни одного узла кратности б
ольше 1);
подстановка найденного многочлена
H
n
-
(
m
+1)
(
x
)
в (5) приво
дит к искомому
интерполяционному многочлену Эрмита. Если же в исходной информации (1)
об
f
(
x
)
имеются значения производных более высокого порядка, чем первый, то
для вос
становления многочл
ена
H
n
-
(
m
+1)
(
x
)
ставится задача эрмитовой же
интерполяции, для чего наряду с полученными его значения
ми
, находят
значения его производных путем дифференциро
вания равенства (6) (возможно

неоднократного, в зависимости от максимального порядка заданных
пр
оизводных функции
f
(
x
)).
Эта процедура построения интерполяционных
многочле
нов Эрмита все более низких степеней продолжается до исчер
пывания
всей информации (1) о функции и ее производных.

Рассмотрим реализацию описанного процесса эрмитовой интерполяции
на простом примере, демонстрирующем возмож
ность восстановления
многочлена
n
-
й степени по его значениям и значениям некоторых его
производных при суммарной кратно
сти узлов
п
1.

Пример.
Пусть сведения о некоторой функции

у
f
(
x
)
пред
ставлены
следующей
дискретной информацией:


В соответствии с обозначениями (1) здесь:
т;

к
0
-
1 =1,
k
1

-
1
=
2
,
к
2
-
1=1
�=
п1к
0
к
1
к
2
=
7

�=
n
=6.

Таким образом, по данным сведениям о функции
у
=
f
(
x
)
,
сосредоточенным в трех узлах

кратности,
соответственно, , 3, , следует строить интерполяционный многочлен Эрмита
Н
6
( х).

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

таблицы данных, записываем интерполяционный многочлен Лагранжа
второй степени


Далее по формуле (
5) представляем
Н
6
(х)
через
L
2
(
x
),H
3
(
x
) и
П
3
(х):







(7)

и дифференцируем этот многочлен дважды:



Подстановкой в
значений

х
=
-
1
,
х
=
0
и
х1

и

приравниванием

заданным значениям


получаем значения


Учитывая их, из
условия Я


находим


Итак, для выявления многочлена
H
3
(х)
в (7) снова имеем задачу

эрмитовой интерполяции с данными, содержащимися в следующей

таблице:



Здесь:

В соот
-
ветствии с (5), для этого случая записываем:


(где



многочлен Лагранжа, интерполирующий функцию
Н
3
(х))
-
Остается
найти постоянную
H
0
, для чего воспользуемся условием


Имеем:


Следовательно,


Подставив это в (7), получаем окончательное выражение искомого
многочлена:


1.12

Занятие № 1
2
.
Понятие сплайна.
Виды сплайнов.

Цель
-

ознакомить студентов с понятием сплайн, видами сплайнов.

Помимо «классических» интерполяционных полиномов, в XX веке
получила широкое распространение в практике еще одна разновидность

интерполянтов, называемых
сплайнами.
Исторически о
ни возникли в
авиастроении при конструировании обтекаемых профилей.

Сплайном
S
(
x
) является интерполирующая функция, которая обладает
наименьшей кривизной и имеет квадратично интегрируемую вторую
производную. Так как сплайн является интерполянтом, то для не
го должно
выполняться условие (.1):









(1)

Вид функции сплайна
S
(
x
)
находится из условия минимальности
функционала








(2)

Эта функция физически моделируется гибким и упругим стержнем,
который проходит через все заданные точки интерполяции. Минимум
функционала () соответствует минимуму упругой потенциальной энергии
этого стержня.

Особенно час
то используется
естественный кубический сплайн
-

полином
3
-
й степени, удовлетворяющий требованиям (1) и (). Такой естественный
кубический сплайн является непрерывной функцией и имеет непрерывные
первую и вторую производную на интервале интерполя
ции
(
x
0
,
x
n
).

Как и в предыдущих параграфах данной главы, интерполянт строится на
базе таблицы данных типа


Таблица с
n
+1
узлами определяет
n

подинтервалов между этими узлами.
Кубический полином (естественный сплайн) необходимо построить для
каждого поди
нтервала. Таким образом, интерполяция кубическими сплайнами
относится к классу кусочно
-
многочленной интерполяции. Так как каждый
кубический полином характеризуется четырьмя константами, поэтому для
построения сплайна
S
(
x
)
на всем интервале [
x
0
;
x
n
]

необходимо определить 4
n

параметров.


Требование непрерывности функции, ее первой и второй производных во
внутренних
(
n
-
1)
узлах дает 3(
n
-
1) уравнения для определения неизвестных
параметров. Условие (1) для всех
(
n
+1)
узлов дает еще
(
n
+1)
уравнений.
Недоста
ющие два условия для однозначного определения естественного
сплайна задаются приравниванием нулю вторых производных в точках
x
0

и
x
n







(3)

что означает нулевую кривизну с
плайна на концах интервала.

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

После формулировки общих положений займемся конструированием
явного вида коэффициентов к
убического сплайна. Для этого вначале
рассмотрим отдельный подинтервал
(
x
i
,
x
i
+1
).
Введем обозначения












(4)











(5)









(6)

Сплайн на подинтервале
(
x
i
,
x
i
+
1
)
удобно представить формулой:







(7)

где ζ
i
+
1

и
ζ
i

-

константы, которые необходимо выразить через известные
величины
x
i
,
y
i
.
Очевидно, что функция (7)
-

полином 3
-
й степени относительно
аргумента
x
.

Легко убедиться непосредственной подстановкой, что







(8)

Получим выражения для трех производных функции (7):






(9)








(10)








(11)


πормула (9) позволяет получить выражение для первой производной на
концах любого подинтервала. Декларированная непрерывность первой
производной сплайна
S
(
x
) на всем интервале (
x
0

,
x
n
)
позволяет приравнять
правые и левые производные в точках узлов
x
i
.
В результате получается
соотношение







(12)

В уравнении (1) и последующих этого параграфа использовано новое
обозначение:









(13)

Заметим, что в этом параграфе величины

определенные формулой
(13), называются разделенными разностями, в отличие от конечных разностей,
введенных в раннее.

Соотношения (1) получены для всех
n
-
1 внутренних узлов. Индекс
i

пробегает значения 1, ... ,

n
-
1.

Таким образом, получена система
(
n
-
1) линейных уравнений
относительно неизвестных величин
ζ
i
.

В качестве недостающих условий можно взять уравнения, выражающие
требования (3). Можно поступить иначе
-

использовать единственные
кубические полиномы, кото
рые проходят точно через 4 первые и 4 последние
точки из заданных
x
i
,
y
i

(
i

 о,.
. .
,
n
).
Тогда два дополнительных условия на
концах интервала
(
x
0

,
x
n
)
выразятся следующим образом:







(14)

где






(15)


Полученная система линейных уравнений (1) и (14) для нахождения
n
+1
параметров
ζ
i

(
i

= 0, ...,
n
) имеет единственное решение. Матрица
коэффициентов при неизвестных
ζ
i

может быть записана в следующем виде:







(16)

Правые части уравнений (свободные члены) соответственно равны:



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

представятся в
виде рекуррентных выражений:






(17)

где коэффициенты

выражаются так:







(18)







(19)

При необходимости многократно вычислять сплайн его удобнее
представить в несколько иной форме:






(20)

где

x
i


x

x
i+1
, i =
0, …,
n
-

1.

Преобразование функции (7) к форме (0) дает следующие выражения
для искомых коэффициентов:





(21)

где
i

 0,…,
n

-

1.


Подстановка коэффициентов (1) в выражение (0) дает функцию
сплайна в явном виде для аргументов
x
i

x

x
i
+1
, лежащих в
i
-
м поддиапазоне
интерполяции.

Пример.

Пусть необходимо построить сплайн
-
интерполянт по набору
шести точек.


Для вычисления коэффициентов

сплайна целесообразно по формулам (4)
-

(6), (11) и (13) вычислить промежуточные данные и свести их в таблицу:



Пользуясь полученными величинами можно составить систему линейных
уравнений (1 и (14). Запишем эту систему в матричном виде:






(22)

Решением системы () является следующий набор значений


которые мы представили в виде обыкновенных дробей.

Вычисленные значения
ζ
i

и
h
i

(
i

=
0, 1,.
.
5) подставляются в формулы (4)
-

(7) или (0), (1), по которым можно получить явный вид сплайна
S
(
x
) на
каждом
i
-
м поддиапазоне
x
i


x


x
i
+1
.


Например, на интервале х

[3, 6] (т.е. между узлами
i
1 и
i
) сплайн
-
интерполянт может быть выражен аналитически

следующим кубическим
полиномом:

S
12
(
x
)  1,610  0,06 х  0,3701 х
2

-

0,0374 х
3

,


где коэффициенты вычислены с точностью до 0,0001. График полученного
сплайна представлен на рис. 1.

Аналогичным образом можно построить кубические интерполирующие
полиномы на всех интервалах, ограниченных узлами интерполяции. В
результате получается сплайн
-
интерполянт для диапазона 1 <
x

< 9, график
которого изображен на рис. .


Рис. 1. График сплайн
-
и
нтерполянта на интервале х

[3, 6].

Кружки
-

исходные данные, линия
-

график полинома (3)


Рис. . График сплайн
-
интерполянта на интервале х


[1, 9].

Задания:

выполнить задание 4 ИДЗ№.



1.13

Занятие № 1
3
.
Численное дифференцирование.

Численное дифференцирование применяется в тех случаях, когда:
функция
f
(
x
) задана таблично и, следовательно, методы дифференциального
исчисления неприменимы; аналитическое выражение
f
(
x
) столь сложно, что
вычисления производной представляют значительны тру
дности.

В основе численного дифференцирования лежит следующий прием:
исходная функция
f
(
x
) заменяется на рассматриваемом отрезке [
a
,
b
]
интерполяционным полиномом
P
n
(
x
) и считается, что
f
'(
x
) и
P
'
n
(
x
) примерно
равны, т.е.
f
’(
x
)=
P
'
n
(
x
), .
(
a


x


b

)
.

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

Всегда, когда это возможно, для численного дифференцирования
используется интерполяционный многочлен с равноотстоящими узлами, так
как это значител
ьно упрощает формулы численного дифференцирования. При
равноотстоящих узлах в начале строится интерполяционный полином
Ньютона, а затем он дифференцируется.

Если
f
(
x
) задана таблицей с неравноотстоящими узлами, то вначале
строится интерполяционный полино
м Лагранжа, а затем он дифференцируется.

Будем
изучать

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

в

которой требуется вычислить
производную. Если требуется вычислить производную во всех узлах, то
вначале полином Ньютона строится по первым 3
-
5
узлам и в них вычисляется
производная, потом полином Ньютона строится по следующим 3
-
5 узлам и в
них вычисляется производная и т. д. до тех пор, пока не будет просчитана вся
таблица.


Если по узлам
x
0
,
x
b

x
2
,
x
3
,
x
4
пос
троить первый интерполяционный
полином

Ньютона и его продифференцировать с учетом


то получим приближенное выражение для

в виде:






(1)

Вторая производная
f
"(
x
) в результате дифференцирования формулы (1) с
учетом равенства


имеет следующее приближенное выражение:




(2)

Если требуется вычислить
f
'(
x
) в узлах
x
0
,
x
1
, …,
x
n
, то в формулу (1)
необходимо поочередно поставить
q
=0,
q
=1,
q
=2,
q
=3,
q
4 соответственно.

В случае выбора узлов
x
j
,
x
j
+1
, …,
x
j
+4

вместо узлов
x
0
,
x
1
, ...,
x
4

в
формулах (1), () следует заменить конечные разности
y
0

на конечные
y
j
,
которые берутся из строчки с номером
j

таблицы конечных разностей. Если
требуется найти производные функции
f
(
x
) в основных табличных точка
x
j

то
каждое табличное значение считается

за начальное и тогда
x
j
=
x
0
,
q
0, а
формулы (1), () будут соответственно для всех
f
'(
x
) и
f
'
'

(
x
) иметь вид:







(1')







(2')


Если интерполяционный полином
строится не по пяти, а по четырем или
трем узлам, то в формулах (1), (1'), (), (') отбрасываются соответственно одно
или два последних слагаемых.

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



Метод неопределенных коэффициентов

Аналогичные формулы можно получить и для случая произвольного
расположения узлов. Использование многочлена Лагранжа в этом случ
ае
приводит к вычислению громоздких выражении, поэтому удобнее применять
метод неопределенных коэффициентов. Он заключается в следующем. Искомое
выражение для производной
k
-
го порядка в некоторой точке хх
i

представляется в виде линейной комбинации заданны
х значений функции в
узлах х
0
,
x
1
, ..х
n
:





(3)

Предполагается, что эта формула имеет место для многочленов
y
=1,
y
=
x
-
x
i
, у(х
-
x
i
)
n
. Подставляя последовательно эти выражения в (3), получаем
систему
n
1 линейных
алгебраических уравнений для определения
неизвестных коэффициентов с
0
, с
1
, …, с
n
.

Пример.

Найти выражение для производной ух в случае четырех
равноотстоящих узлов (
n

= 3).

Равенство (
3
) запишется в виде





(
4
)

Используем
следующие многочлены:





(
5
)

Вычислим их производные:




(6)

Подставляем последовательно соотношения (5) и (6) соответственно в
правую и левую части (4) при хх
1
:


Получаем окончательно систему уравнений в виде


Решая
эту систему получаем


Подставляя эти значения в равенство (4), находим выражение для
производной:


Задания:

выполнить задание 1 и  ИДЗ№3.


1.14

Занятие № 1
4
.
Численное интегрирование. Квадратурные формулы
Ньютона
-
Котеса.

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

Пусть требуется найти значение
I

интеграла Римана
для
некоторой заданной на отрезке
[а,

b
]

функ
ции
f
(
x
).

Хорошо известно, что
для функций, допускающих на промежутке
[а,

b
]

конечное число точек
разрыва первого рода, такое значение существует, единственно и может быть
фор
мально

получено по определению:





(1)


где



произвольная упорядоченная система точек от
резка
[а,

b
]

такая,
что


а
ξ
i



произвольная точка элементарного промежутка [х
i
-
1
,_ х
i
].

В математическом анализе обосновывается аналитический способ
нахождения значения
I

с помощью знаменитой формулы Ньютона
-
Лейбница

I

=

F
(
b
)
-
F
(
a
),

(2)

где
F
(
x
)



некоторая первообразная для данной функции
f
(
x
).

К сожалению, применение этого весьма привлекательного под
хода к
вычислению

I

наталкивается на несколько серьезных препятствий. Самое
главное из них



это несуществование пер
вообразной среди элементарных
функций для большинства эле
ментарных функций

f
(
x
)
; например, таким
способом не удается вычислить


Если первообразная

F
(
x
)

для заданной функции
f
(
x
)
все же найдена, то
вычисление двух ее значений

F
(
a
)

и

F
(
b
)

может оказаться более трудоемким,
чем вычисление существенно большего количества значений
f
(
x
).

Поскольку в общем случае значения функций находятся лишь
приближенно, использование точной формулы () приводит к
приближенному

результату, который может быть б
олее эффективно получен с помощью какой
-
либо специ
альной приближенной формулы на основе значений подынте
-
гральной функции
f
(
x
).
Такие специальные приближенные формулы для
вычисления определенных интегралов называют
квадратурными формулами
(механическими к
вадратурами)
или
формулами численного интегрирования.
Первый из этих терминов можно

связать с геометрическим смыслом
определенного интеграла: вычисление

при
f
(х)>
0 равносильно
построению квадрата, равновеликого криволинейной трапеции с основанием
[а,
b
]

и

«крышей»
f
(
x
).


Простые квадратурные формулы можно вывести непосред
ственно из
определения интеграла, т.е. из представления (1). Зафиксировав там некоторое
п≥
1, будем иметь



(
3)

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


i
-
1
,
x
i
], а

высотами


ординаты
f
(
ξ
i
),

рис. 1).


Рис.
1.
Геометрическая интерпретация общей формулы прямоугольников
(.3)


Чтобы из общей формулы (3) получить конструктивное правило
приближенного вычисления интеграла, воспользуемся свободой расположения
точек х,, разбивающих промежуток ин
тегрирования

[а,

b
]

на э
лементарные
отрезки [
x
i
-
1
,
x
i
], и свобо
дой выбора точек
ξ
i

на этих отрезках.

Условимся в дальнейшем

пользоваться равно
мерным разбиением отрезка

[а,

b
]

на

п

частей точками х
i


с шагом

полагая


(4)

При таком разбиении формула (3)
приобретает вид



(5)

Теперь дело за фиксированием точек
ξ
i

на элементарных отрез
ах [
x
i
-
1
,
x
i
].
Рассмотрим три случая.

1)

Положим

ξ
i

=
x
i
-
1
. Тогда из (5) получаем






(6)


2)

Пусть в (5)
ξ
i

=
x
i
. Тогда имеем









(7)

πормулы (6) и (7) называются
квадратурными формулами левых и
правых прямоугольников
соответственно. Совер
шенно очевидно (рис. ),

что

дают двусторонние приближения к значению
I

интеграла от
монотонной функции.


Рис.
.2.

Геометрическое оценивание интеграла от монотонной функции с помощью:

Можно рассчитывать на большую точность получения зна
чения
интеграла, если взять точку
ξ
i

посередине между точками
x
i
-
1

и
x
i
. Отсюда
приходим к следующему случаю.

3) πиксируем


В ре
зу
льтате имеем
квадратурную формулу средних прямоугольни
ков
или, иначе,
формулу средней
точки




(8)

Учитывая априорно большую точность формулы (8) по срав
нению с
формулами (6) и (7) при том же объѐме и
характе
ре вычислений, эту
симметричную формулу будем впредь назы
вать просто
формулой
прямоугольников.

Полученные из определения интеграла квадратурные пра
вила (6)

(8) не
содержат в себе сведений, позволяющих сказать, каким нужно взять число

п

элементарн
ых промежутков

[
x
i
-
1
,
x
i
] или, иначе, каким должен быть шаг

h

разбиения (4) отрезка интегрирования

[а,

b
]
, чтобы гарантированно найти зна
-
чение
I

интеграла с наперед заданной точностью

ε.


Ответ на этот вопрос дает

формула остаточного члена (
гло
бальной
погрешности
)
квадратурной формулы
прямоугольни
ков:





(9)

Как видно из формулы (9), при увеличении числа

п
элементарных
отрезков, на которые разбивается промежуток ин
тегрирования

[а,

b
),

ошибка
численного интегрирования по фор
муле средней точки
(8)

убывает
пропорционально квадрату шага
h
.

Погрешность численного ин
тегрирования
непрерывно дифференцируемой функции по фор
мулам левых и правых
прямоугольников
(6), (.7)

убывает лишь по линейному закону.


С
емейство квадратурных форму
л
Нью
т
она
-
К
отеса

Подстановка в интеграл


вместо функции

f
(
x
)

еѐ
интерполяционного многочлена Лагранжа той или иной степе
ни

п
приводит к
семейству квадратурных формул, называемых формулами Ньютона
-
Котеса.

πункция

f
(
x
)
может быть единственным образом
представлена в виде


где




интерполяционный многочлен Лагранжа
,



остаточный член,

Если система узлов интерполирования
совпадает с точками разбиения
(4) отрезка

[а,

b
]

с шагом
h
,

то замена переменной

x
=
x
0
+
qh

трансформирует
многочлен Лагранжа к
виду







(10)

Для того, чтобы использовать такое выражение
L
n
(
x
) вместо
f
(х) в
нужно изменить границы интегрирования (значению

ха


соответствует значение
q
=
0,а
х
=
b



значе
ние

q

=

п)

и учесть, что

dx
=
hdq
.
Таким образом, получаем


Это равенство, переписанное в виде







(11)

и есть

квадратурная формула Ньютона
-
Котеса,
где


(12)




коэффициенты Котеса.

На самом деле,
формулы (11)
-
(1) определяют
се
мейство

квадратурных
формул. Параметром этого семейства является число
п


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

Для практического применения целесообразно рассмотреть некоторые
частные случаи.

Пусть
n

 1. Это значит, что диапазон интегрирования [
a
,
b
] содержит
только два узла:
x
0

=
a

и
x
1

=
b
.

Вычислим коэффициенты Котеса для этого случая. Подставляя в (1)
n

=
1 и
i

 0, 1 получим:







(13)

Подстановка

коэффициентов (13) в (1) дает квадратурную формулу
Ньютона
-
Котеса первого порядка:




(14)

В данном случае величина шага
h

равна длине всего диапазона
интегрирования [
a
,
b
].


При таком приближении подынтегральная функция
f

(
x
)
заменяется

линейной (полиномом Лагранжа 1
-
й степени). Графически значение
определенного интеграла
I
(
f
,
a
,
b
)
изображается площадью криволинейной
трапеции, ограниченной сверху (или снизу) графиком функции
f

(
x
) (см.
пример на рис. 3).


Рис.3. К выводу формулы
трапеций.


Сплошная линия
-

график интегрируемой функции
f

(
x
),

штрихпунктирная
-

график линейного приближения интегрируемой функции на интервале [
x
0
, х
1
].

Крестиками заполнена область
R
,
площадь которой равна погрешности квадратурной формулы (14).


Приб
лиженное значение интеграла, выраженное квадратурной суммой
(14), равно площади обычной трапеции с той же высотой
h
.
Площадь
криволинейного сектора дает погрешность
R

вычисления определенного
интеграла с помощью квадратурной формулы (14).

Когда данный спос
об вычисления определенного интеграла используют
на практике, то обычно диапазон интегрирования
[
a
,
b
]
разбивают на
N
одинаковых поддиапазонов шириной




h

= (
b

-

a
) /N


(15)

и к каждому из них применяют формулу (14). Тогда искомый интеграл
представится суммой, каждый член которой будет содержать общий множитель
(15). Каждое значение интегрируемой функции
y
k

(
k

=
0, 1, …,
N
)
в эту сумму
будет входить дважды, кроме двух крайних
y
0

и
y
N
. В результате интеграл
I
(
f
,
a
,
b
)
выразится следующей суммой:







(16)

Последняя формула называется
формулой трапеций.

Теперь пусть
n
. Это значит, что диапазон интегрирования [
a
,
b
]
содержит три узла: х
0
=
a
,
x
1
=(
a
+
b
)
/2,
x
2
=
b
.
В данном случае шаг
интегрирования равен половине диапазона интегрирования
h
=(
b
-
a
)
/ 2

Коэффициенты Котеса в данном случае получаются вычислением
интегралов (11) с подстановкой
n

  для
i

= 0, 1, 2:






(7)

Искомый интеграл представится в виде:



(18)

где
(
b

-

a
) = 2
h
.

πормула (18) была получена заменой подынтегральной функции
f
(
x
)
на
квадратичную, проходящую через заданные точки
(
x
i
,
y
i
)
,
i

=
0, 1, 2.
Полученное
значение искомого интеграла графически изображается площадью
под параболой (см. рис. 4).



Рис. 4. К выводу формулы

Симпсона.

Жирная сплошная линия
-

график интегрируемой функции
f

(
x
),

штриховая
-

график квадратичного приближения
L
2
(
x
) интегрируемой
функции на интервале [
x
0
,
x
2
].
Заштрихованная область
R

выражает погрешность квадратурной формулы (18).

Видно, что на интервале [
x
0
,
x
1
]
квадратичное приближение дает заниженный результат, на интервале
[
x
1
,
x
2
]
-

завышенный.

На практике весь диапазон инте
грирования [
a
,
b
] обычно предварительно
разбивается на четное число
N
= 2
M

равных интервалов, которые попарно
объединяются в
M

поддиапазонов. Каждый поддиапазон содержит три узла:
два крайних и один внутренний. К каждому из
M

поддиапазонов применяется
формула (18), т.е. каждой тройке узлов соответствует своя интерполирующая
парабола. Параболы стыкуются в совпадающих крайних узлах отдельных
поддиапазонов. Номера всех таких узлов четные от
i

  до
i

= 2
M

-

, поэтому
при суммир
овании выражений (18) по поддиапазонам следует учесть, что
слагаемые с этими номерами встретятся дважды.

В результате суммирования с учетом значений коэффициентов Котеса
(17) получается рабочая формула








(19)

гд
е для краткости введены обозначения следующих сумм:





(20)

Заметим, что сумма
S
o

содержит
M

слагаемых, а сумма
S
e

состоит из
M
-
1
слагаемых. Шаг интегрирования равен
h
=(
b

-

a
)/(2
M
)=(
b

-

a
)/
N
.

πормула (19) называется
квадратурной
формулой Симпсона.


При
n

 3 аналогичным путем получается квадратурная формула с тремя
узлами на интервале интегрирования:





(21)

где
h

=
(
b

-

a
)
/ 3.

При решении прикладных задач диапазон интегрирования
[
a
,
b
]
предварительно разбивается на равные поддиапазоны, число которых кратно
трем
(
N

= 3
M
),
и к каждому из поддиапазонов применяется формула (0). Тогда
получится рабочая формула вида:







(22)

Шаг интегрирования равен
h
=(
b
-
a
)
/(3
M
) = (
b
-
a
)
/
N
, а отдельные суммы
выражены следующими формулами:







(23)

Сумма () называется
квадратурной формулой Ньютона

или

формулой «трех восьмых».

В практике квадратурная формула Ньютона
употребляется гораздо реже, чем формулы Симпсона и трапеций.

Задания:

выполнить задание 3 ИДЗ№3.


1.15

Занятие №
15
.
Квадратурные формулы Гаусса.

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

В предыдущем параграфе рассматривались различные квадратурные
формулы для равномерного расположения узлов в интервале интегрирования
[
a
,
b
].
Оказывается, при фиксированном количестве
узлов можно добиться
значительного уменьшения погрешности вычисления определенного интеграла,
если отказаться от их равномерного расположения.


Для нахождения оптимального расположения узлов (которое
обеспечивает уменьшение погрешности численного интегриров
ания) прежде
всего сведем интервал интегрирования общего вида
[
a
,
b
]
к стандартному
интервалу [
-
1, 1] с помощью линейной замены переменной интегрирования:




x

 (а 
b
)
/ 2 +
t

(
b

-

а)
/ 2.


(25)

Исходн
ый интеграл приобретет вид:







(26)

где
t

-

новая переменная интегрирования.

В соответствии с общей квадратурной формулой запишем:








(27)

Нашей целью является выбор таких
п
узлов
t
i

внутри интервала
интегрирования [
-
1, 1] и таких
п
коэффициентов
A
i

(
i

 1, п), чтобы формула
(7) была
точной
для всех
полиномов
максимально большой степени. Мы
имеем возможность варьировать 
n

величин
t
i

и
A
i

(
i
=1,
…,

п),
следовательно, эта максимальная степень полинома равна
N
= 2
n

-

1.

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

Лемма.

Для достижения поставленной цели необходимо и достаточно,
чтобы формула (7) выполнялась точно для следующих 
n

подынтегральных
функций

Доказательство
.

По условиям леммы выполняются следующие равенства для степенных
функций
f
(
t
) =
t
k

(
k

 0, 1, …,
п
-

1):




(28)

Как было принято выше, подынтегральную функцию
f
(
t
)
представим в
виде полинома степени 
п
-

1:





(29)

где
C
k

-

коэффициенты полинома.

Проведем интегрирование полинома (9) по интервалу [
-
1, 1], используя
условие (8) и изменя
я порядок суммирования:




(30)

Последнее равенство обусловлено определением (9).

Мы получили, что равенство (7) является точным для полинома вида
(9), при этом в ходе преобразований мы воспользовались условиями (8).

Лемма

доказана.

Эта лемма обеспечивает справедливость следующего важного
утверждения. Пусть нам удастся найти такие числа
t
i

и
A
i

(
i

= 1,
…,

n
),
обращающие формулу (7) в точную, используя в качестве подынтегральных
степенные функции вида
t
k

(
k

= 0, 1,
...,

2
n

-

1). Тогда эти найденные числа
t
i

и
A
i

обеспечат точность формулы (7) и для любой подынтегральной функции,
которая представляется в виде полинома степени 
n

-

1.

Теперь в уравнениях (8) проинтегрируем аналитически левые части:




(31)

Сравнение (31) с правыми частями уравнений (8) дает нам систему 
n

уравнений для 
n

искомых неизвестных
t
i

,
A
i

(
i

= 1, ..
.,

n
):






(32
)

Система является нелинейной относительно неизвестных
t
i
,
A
i

(
i

 1,…,
n
).
Для ее решения целесообразно воспользоваться свойствами полиномов
Лежандра. Краткие сведения о полиномах Лежандра приведены в приложении
1
.

Выберем подынтегральные функции
f
(
x
)
в
виде:




(33)

где
P
n

(
t
)
-

полином Лежандра степени
n
.

Это значит, что все
n

функций (33) являются полиномами степени ≤ 
n
-
1.
Для таких подынтегральных функций, согласно доказанному выше, формула
(7) является точной и коэффициенты
A
i

(
i

 1,…,
n
)
удовлетворяют системе
(32).

Применим квадратурную формулу (7) для функций (33):





(34)

В
силу свойства ортогональности полиномов Лежандра левые части всех
уравнений (34) равны нулю. Следовательно, и правые части (34) будут равны
нулю




(35)

при произвольных значения
х
A
i
,
если положить
P
n

(
t
i
)
= 0.

Таким образом, оказывается, что в качестве узлов
t
i

следует взять точки
нулей полиномов Лежандра степени
n
. Точки нулей полиномов Лежандра
различных степеней рассчитаны с большой точностью и сведены в
специальные таблицы. Х
арактерным свойством узлов
t
i

и коэффициентов
A
i

в

данном методе является симметрия их значений относительно центра таблицы
t

= 0.

Если числа
t
i

известны, то система (3) становится линейной
относительно величин
A
i
. Для ее решения можно применить любой из
описанных методов. Численные значения узлов
t
i

и коэффициентов
A
i

с
точностью до восьми знаков для различных
n

приведены в таблице приложения
2.

Вычислив значения узлов
t
i

и коэффициентов
A
i

для выбранного
n
,
можно

вернуться к исходной переменной интегрирования
x

и записать
квадратурную формулу Гаусса
в следующем виде:





(36)

где





(37)

Точность
интегрирования методом Гаусса резко возрастает с
увеличением числа узлов в интервале интегрирования [
a
,
b
]. Теоретические
оценки и расчеты показывают, что для интервалов небольшой величины
обычно достаточная точность обеспечивается при
n

 8. Соответствующ
ая
таблица узлов
t
i

и коэффициентов
A
i

приведена в приложении .

Если диапазон интегрирования
[
a
,
b
]
имеет значительную ширину, его
можно разбить на несколько равных поддиапазонов и провести процедуру
вычисления интеграла методом Гаусса для каждого поддиап
азона отдельно.
Интеграл по
j
-
му поддиапазону с границами [
a
j

,
b
j

] с помощью линейного
преобразования (37) представляется в виде суммы, аналогичной (36):




(38)

где

Искомое значение определенного интеграла по полному диапазону
[
a
,
b
]
вычисляется как сумме интегралов (38) по всем поддиапазонам.


Задания:

выполнить задание 3 ИДЗ№3.


1.16

Занятие №
16
.

Численные методы решения задачи Коши для
обыкновенных дифференциальных
уравнений. Классификация
методов.

Метод Пикара.

Будем рассматривать обыкновенное дифференциальное уравнение (ОДУ)
первого порядка



(1)

с начальным условием

y

0
)

=

у
0
,


(2)

где
f
(
x
)



некоторая заданная, в общем случае, нелинейная функция двух
переменных. Будем считать, что для данной задачи (1)
-
(), называемой
начальной задачей
или
задачей Коши,
выполняются требования,
обеспечивающие существование и единственность на отрезке [х
0
,
b
] ее решения
уу(х).


Несмотря на внешнюю простоту уравнения (1), решить его аналитически,
т.е. найти общее решение

уу(х, С)
с тем, чтобы затем выделить из него
интегральную кривую

уу(х),

проходящую через заданную точку


0

0
),
удается лишь для некоторы
х специальных типов таких уравнений. Поэтому, как
и в родственной для (1)
-
() задаче вычисления интегралов, приходится делать
ставку на приближенные способы решения начальных задач для ОДУ, которые
можно разделить на три группы:

1)

приближенно
-
аналитические м
етоды
;

2)

графические или машинно
-
графические методы;

3)

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

К методам первой группы относят такие, которые позволя
ют находить
приближение решения
у(х)
сразу в виде некоторой «хорошей» функции
φ
(х).
Например, широко известен

метод степенных рядов,

в одну из реализаций
которого заложено представление искомой функции

у(х)
отрезком ряда
Тейлора, где тейлоровские коэффициенты, содержащие производные высших

порядков, находят последовательным дифференцирова
нием самого уравнения
(1). Другим представи
тел
ем этой группы методов является метод
последовательных приближений, суть которого приведена чуть ниже.

Название

графические методы

говорит о приближенном представлении
искомого решения
у(х)
на промежутке

[
x
0
,
b
]
в
виде графика, который можно
строить по тем
или иным прави
лам, связанным с графическим толкованием
данной задачи. πи
зическая или, возможно, точнее будет сказать,
электротехниче
ская интерпретация начальных задач для определенных видов
уравнений лежит в основе машинно
-
графических методов при
ближен
ного
решения. Реализуя на физико
-
техническом уровне заданные электрические
процессы, на экране осциллографа на
блюдают поведение решений
дифференциальных уравнений, описывающих эти процессы. Изменение
параметров уравнения приводит к адекватному изменению п
оведения решений,
что по
ложено в основу специализированных
аналоговых вычисли
тельных
машин
(АВМ).

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

приближенных значений

y
i

искомого решения

у(х)
на некото
рой сетке

значений аргумента
х.
Этим способам и будет посвящено
дальнейшее
изложение. Что делать с получае
мыми численными значениями решения,
зависит от прикладной постановки задачи. Если речь идет о нахождении только
значе
ния

у(
b
),
тогда точка

b

включается как конечная в систему рас
четных
точек х
i
, и все приближенн
ые значения

y
i


y
(
x
i
)
,
кроме последнего, участвуют
лишь как промежуточные, т.е. не требуют ни запоминания, ни обработки. Если
же нужно иметь прибли
женное решение

у(х)
в любой точке

х,
то для этого к
получае
мой числовой таблице значений

y
i

можно применить

какой
-
либо из
способов аппроксимации табличных функций, рассмотренных ранее, например,


интерполяцию или сплайн
-
интерполя
цию. Возможны и другие использования
численных дан
ных о решении.

Коснемся одного приближенно
-
аналитического способа решения
начальной

задачи (1)
-
(), в котором искомое ре
шение

уу(х)
в некоторой
правой окрестности точки

х
0

явля
ется пределом последовательности
получаемых определенным образом функций

у
п
(х).

Проинтегрируем левую и правую части уравнения (1) в границах от

х
0

до
х:


Отсюда, с учетом того, что одной из первообразных для

у'(х)
служит

у(х)
,
получаем


или, с использованием начального условия (),




(3)

Таким образом, данное
дифференциальное уравнение (1) с на
чальным
условием () преобразовалось в
интегральное урав
нение
(неизвестная
функция здесь входит под знак интеграла).

Полученное интегральное уравнение (3) имеет вид зада
чи о неподвижной
точке


для оператора

πормально

к этой зада
че
можно применить метод простых итераций





(4)

достаточно обстоятельно рассматривавшийся приме
нительно к системам
линейных и нелинейных алгебраических и трансцендентных

уравнений. Беря в
качестве начальной функции
y
0
(х)
заданную в () по
стоянную

y
0
,
по формуле
(4) при

п0
находим первое при
ближение



Его подстановка в (4) при
п
1 дает второе приближение


и т.д. Таким образом, этот приближенно
-
аналитический метод, называемый
методом последовательных приближений
или

методом Пикара
определяется формулой




(5)

где

n
=0,1, 2,...
и

у
0
(х)
y
0
.

Отметим две характеристики метода последовательных приближений
Пикара, которые можно отнести к негативным
.
Во
-
первых, в силу известных
проблем с эффективным нахождением первообразных, в чистом виде метод
(
5)
редко реализуем. Во
-
вторых, как видно из вышепри
веденного утверждения,
этот ме
тод следует считать локальным, пригодным для приближения решения в
малой правой окрестности начальной точки. Большее значение метод Пикара
имеет для доказательства существования и единственности решения задачи
Коши, нежели дл
я его практи
ческого нахождения.


1.17

Занятие №
17
.
Методы Эйлера.

Цель
-

ознакомить студентов с методами Эйлера решения задачи Коши
для обыкновенных дифференциальных уравнений.

М
етод
Э
йлера


разные подходы к построению

Учитывая ключевую позицию, которую
занимает метод Эйлера в теории
численных методов ОДУ, рассмотрим несколь
ко способов его вывода. При
этом будем считать, что вычисления проводятся с
расчетным шагом

расчетными точками (узлами)
служат точки х
i

 х
0

+
ih

(
i
=0
,1,...,
п)
про
-
межутка [
x
0
,
b
]

и

целью является построение таблицы



приближенных значений
y
i

решения
у  у(х)
задачи (1)
-
() в расчетных точках
х
i
.

Геометрический способ.
Пользуясь тем, что в точке
х
0

из
вестно и
значение решения
у(х
0
)  у
0

(согласно ()), и значе
ние его производной
у'(
x
о
)=
f
(
x
0
,
y
0
)
(согласно (1)), можно записать уравнение касательной к графику
искомой функции
y
 у(х)
в точке (х
0
;
y
0
):




(6)

При достаточно малом шаге
h

ордината




(7)

этой касательной, полученная подстановкой в правую часть (6) значения

по непрерывности должна мало от
учаться от ординаты
у(х
1
) решения
у(х)
задачи (1)
-
(). Следоват
ельно, точка

1
, у
1
)
пересечения касательной (6) с
прямой
x
=
x
1

может быть приближенно принята за новую на
чальную точку.
Через эту точку снова проведем прямую


которая уже приближенно отражает поведение касательной к
уу(х)
в точке
(
х
1
,у(х
1
)).
Подставляя

сюда
хх
2
(х
1
+
h
), иначе, пересекая эту «касательную»
прямой
х  х
2
, получим приближение значения
у(х
2
)
значением


и т.д. В итоге этого процесса, определяемого формулой



i

=
0,1,...,
п


(8)

и называемого
методом Эйлера,
график решения
у  у(х)
дан
ной задачи Коши
(1)
-
() приближенно представляется ло
маной, составленной из отрезков
приближенных касатель
ных (рис. 1), откуда происходит другое название
метода (8)


метод ломаных.



Рис.
1.
Геометрическая интерпретация метода Эйлера


Применение формулы Тейлора.
Описываемый здесь спо
соб вывода
метода Эйлера тесно связан с предыдущим. Линеаризуя решение
в

окрестности
начальной точки по формуле

Тейлора, имеем


Отсюда при
х  х
1

получаем




(9)

Точное равенство (9), переписанное в виде


говорит о том, что здесь мы имеем одновременно как саму фор
мулу Эйлера для
вычисления значения

(сравните с формулой (7)), так и ее

остаточный
член




(10)

где ξ
i



некоторая точка интервала
(
x
0
,
x
1
).

Остаточный член (10) характеризует
локальную
(или, иначе,
шаговую)
ошибку
метода Эйлера, т.е. ошибку,

совершае
мую на одном шаге. Очевидно,
что от шага к шагу, т.е. при мно
гократном применении формулы (8), возможно
наложение ошибок. За
п
шагов, т.е. в точке
b
,
образуется
глобальная ошибка
.
П
орядок глобальной ошибки (относительно шага
h
) на единицу ниже,
чем
порядок локальной ошибки, а по
рядком глобальной ошибки и определяется
порядок соответст
вующего численного процесса решения задачи Коши.
Таким

обра
зом, локальная ошибка метода Эйлера, согласно (10), есть

глобальная


0(
h
),
т.е.
метод Эйлера относит
ся к ме
тодам первого порядка.

Разностный способ.
Рассматривая уравнение (1) в точке
x
0

с учетом ()
имеем равенство


Применяя к его левой части аппроксимацию производной пра
вым
разностным отношением первого порядка точности


получаем


что идентично равенству (9), поставляющему формулу для вычисления
у
1

вида
(7) и локальный остаточный член (10). Ясно, что для получения общей
расчетной формулы (8) мож
но было сразу применить аппроксимацию
по
формуле (6.16) в равенстве





(11)

заменив неизвестное точное значение
y
(
x
i
)
известным прибли
женным
значением
y
i
.

Заметим, что
порядок получающегося таким способом ме
тода
численного интегрирования дифференциальной задачи
(1)
-
(2)
совпадает с
порядком аппроксимации производ
ной в левой части уравнения
(1).

Квадратурный способ.
Как было показано на
чальную задачу для ОДУ
(1)
-
() можно заменить эквива
лентным интегральным уравнением (3). При
хх
1

из него получится равенство




(12)


Применение к интегралу в правой части равенства (1) про
стейшей
(одноточечной) формулы левых прямоуг
ольников дает приближенную формулу



правая часть которой, очевидно, совпадает с выражением (7) для подсчета
значения
у
1
.
В общем случае расчетная форму
ла (8) метода Эйлера получается
численным интегрированием посредством простейшей формулы левых
прямоуго
льников в ра
венстве




(13)

в предположении, что на каждом
i
-
м шаге в роли начальной точки
(
x
0
,
y
0
)
выступает точка
(
x
i
,
y
i
).
Зная точность исполь
зуемой в (13) квадратурной
формулы, легко при
йти к тому же выражению локальной ошибки метода
Эйлера, что и при других способах его построения.

Существуют и другие подходы к выводу метода Эйлера. В частности, он
будет возникать далее как частный случай некото
рых семейств численных
методов решения зад
ачи (1)
-
(2).


Н
есколько простых модификаций метода Эйлера

Разовьем последний из подходов к построению метода Эй
лера. Очевидно,
применение к интегральному равенству (13) других простейших квадратурных
формул будет порождать но
вые методы численного
интегрирования задачи
Коши (1)
-
(2).

Так, если в (13) использовать простейшую квадратурную формулу
правых прямоугольников, придем к методу




(14)

Этот метод называют
неявным
(или
обратным) методом Эйл
е
ра,
поскольку для вычисления неизвестного значения

по известному
значению

требуется ре
шать уравнение, в общем случае нелинейное.

Применение к интегралу в (13) простейшей квадратур
ной формулы
трапеций приводит тоже к неявному методу





(15)

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

более высокий (на единицу) порядок точности метода трапеций (15) по сравн
е
-
нию с явным и с неявным методами Эйлера (8) и (14), т.е.
мето
д трапеций
(15)


это метод второго порядка
.

Некоторый интерес представляет совместное применение явного метода
Эйлера и неявного метода трапеций.

По форме равенство (15) представляет собой ска
лярную задачу о
неподвижной точке относительно неизвестного значе
ния
у
i+1
.
Поэтому, если в
правую часть (15) подставить хо
рошее начальное приближение
,
подсчитываемое по форму
ле (14), то тогда само это равенство можно считать
шагом метода простых итера
ций для уточнения этого значения. Таким образом,
получаем гибридный метод



i
=0,1,...,
n
,

(16)

который называют
методом Хойна
.

Ясно, что можно достичь большей точности, если, исходя из того же
начального
приближения


сделать не одну, а несколько итераций по методу трапеций:





(17)

Такой вариант совместного применения метода Эйлера и метода
трапеций называют
усовершенствованным методом Эйлера
-
Коши с
итерационн
ой обработкой.
Делать много итера
ций по формуле (17) не
рекомендуется (обычно их выполняют не более трех
-
четырех). Совпадение
определенного числа разря
дов в полученных числах
говорит о
точности, с ко
торой решено методом простых итераций уравнение (15)

от
-
носительно
у
i
+1
,
а вовсе не о том, что с такой точностью найдено значение
у(х
i
+1
).


Чтобы получить следующую модификацию метода Эйлера,

проинтегрируем уравнение (1) по отрезку

Имеем



откуда следует равенство




(18)

Применяя к последнему интегралу одноточечную квадратурную формулу
средних прямоугольников и заменяя значения
y
(
x
i
-
1
)
и
y
(
x
i
)
известными
приближенными значениями
у
i
-
1

и
у
i

соответственно, из (18) выводим формулу
для подсчет
а приближенного значения
y
(
x
i
+1
)




(19)

которую будем называть
уточненный методом Эйлера
.

Как известно, квадратурная формула прямоугольников (средней точки)
имеет тот же порядок точности, что и квадра
турная формула трапеций, так что
уточненный метод Эйлера
(19) тоже
является методом второго порядка.
Подтвержде
нием этого факта может служить вывод метода (19) на разно
стной
основе. Применив к равенству (11) формулу симмет
ричной аппроксимации
y
'(
x
i
)
вто
рого порядка точности, получим


откуда после приближенной замены
следует(19).

Обратим внимание на одно принципиальное отличие мето
да (19) от всех
других рассмотренных до этого момента мето
дов: метод (19) является
двухшаговым.
Здесь для вычисления знач
ения
у
i
+1

привлекаются два
предыдущих значения
у
i

и
у
i
-
1
.
Двухшаговость накладывает определенные
ограничения, по крайней мере, на начало численного процесса: значение

не может быть найдено непосредственно этим мето
дом с тем же шагом
h
.
Поэтому
недостающую вторую началь
ную для процесса (19) точку
приходится получать другим пу
тем, например, явным методом Эйлера, а чтобы
не сделать сразу большой ошибки, применяя на старте метод более низкого

порядка точности, рекомендуется осуществлять постепенно
е вхож
дение в
процесс (19). Так, «разгон» можно выполнить по фор
мулам




(20)

а далее уже переключаться на счет по формуле (19).

Пример 1
. Рассмотрим простое линейное уравнение


с
начальным условием
y
(0)=1.

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

Сначала сделаем несколько последовательных приближений по ме
тоду
Пикара. Его итерационная формула (5) для
данной начальной за
дачи имеет вид


Подставляя сюда
у
0

=1
при
п
 0,1,  последовательно получаем:


Эти результаты удобно сравнить с точным решением, если в последнем
разложить в ряд по степеням
х
фигурирующую там функцию
e
-
3
x
. Тогда
получим
представление решения в виде ряда


с которым, как видим, хорошо согласуются приближения
y
1
,
y
2
,
y
3
,
оп
ределяемые методом Пикара.


Теперь проведем подсчет приближенных значений решения
у(х)
данной
задачи в точке
х
0. численным методом Эйлера и его модифи
к
ациями,
принимая
h

 0.1 (т.е. за два шага). Результаты этих вычислении и фактические
ошибки, найденные сравнением с точным значением
y
(0.2) = 0.581881...,
отражены в следующей таблице.

Метод




Эйлера (8)

0.7

0.51

≈0.07

Неявный Эйлера (14)

≈0.7846

≈0.6343


-
0.05

Трапеций (15)

≈0.7478

≈0.5788

≈0.003

Хойна(16)

0.755

≈0.5895


-
0
.008

Уточненный Эйлера (19
)
-
(
20)

0.755

0.587


-
0.005


И
справленный метод
Э
йлера

Пусть найдено приближенное значение
решения
у  у(х)
задачи
(1)
-
() и требуется вычислить
где
Запишем разложение
решения по формуле Тейлора
р
-
го
порядка, принимая за базовую точку
x
i

(т.е.
по степеням х
-
х
i
) и положим в этом разложении
x
=
x
i
+1
.
Имеем




(21)

Если ограничится двумя слагаемыми в право
й части разложения (1), то,
получим обычный ме
тод Эйлера (8). Посмотрим, что дает учитывание третьего
сла
гаемого.

При
р
  из (1) следует равенство




(22)

Значение первой производной в точке
x
i

в силу с
вязи (1), приближенно
известно:




(23)



Дифференцируя (1), по формуле полной производной


находим приближенное значение второй производной:





(24)

Подставляя приближенные выражения
)
в равенство
(), получаем следующую формулу для вычисления

при
i
0, 1, …,
n
:




(25)


Определяемый ею метод будем называть
исправленным мето
дом
Эйлера.

Так как при
i
0 формулы (3) и (4) точны, а
y
0
=
y
(
x
0
),
согласно
начальному условию (), то на первом шаге вычислений по формуле (5) будет
совершаться ошибка, связанная только с усечением ряда Тейлора.
Следовательно, ло
кальная ошибка или, иначе,
шаговая погрешность
метода
(5) составляет величину
o
(
h
3
), а это означает, что
исправ
ленный метод Эйлера
относится к методам второго порядка.


1.18

Занятие №
18
.
Методы Рунге
-
Кутта.

О
семействе методов Рунге
-
Кутта.
М
етоды второго порядка

Недостатком исправленного метода Эйлера (5) и других методов более
высоких порядков, основанных на пошаговом представлении решения
у(х)

задачи (1)
-
() по формуле Тейлора и последовательном дифференцировании
уравнения (1) для получения тейлоровых коэффицие
нтов, является
необходимость вычисления на каждом шаге частных производных функции
f
(
x
,у).

Идея построения явных
методов Рунге
-
Кутты
p
-
го по
рядка заключается
в получении приближений к значениям
f

i
+1
) по формуле вида






(26)

где
φ(х, у,
h
)


некоторая функция, приближающая отрезок ря
да Тейлора до
р
-
го
порядка и не содержащая частных производных функции
f
(
х
,
у).

Так, полагая в (6)
φ(х, у, р)
=
f
(
x
,
у)
приходим к методу Эйлера (8),
т.е.
метод Эйлера можно считать простейшим при
мером методов Рунге
-
Кутты,
соответствующим случаю
р 
1.

Для построения методов Рунге
-
Кутты порядка, выше пер
вого, функцию
φ(х, у,
h
)
берут многопараметрической, и подби
рают ее параметры сравнением
выражения

(6) с многочленом Тейлора для
у(х)
соответствующей желаемому
порядку степени.

Рассмотрим случай
р.
Возьмем функцию
φ
в (6) следующей
структуры:



Ее параметры
с
1
, с
2
, а
и
b

будем подбирать так, чтобы записан
ная,
согласно (6), формула



(27)

определяла метод второго порядка, т.е. чтобы максимальная ло
кальная ошибка
составляла величину
o
(
h
3
).

Разложим функцию двух переменных
f
(
x

+
ah
, у 
bhf
(
x
, у))
по формуле
Тейлора, ограничиваясь линейными
членами:


Ее подстановка в (7) дает




(28)

Сравнение последнего выражения с тейлоровским квадратичным
представлением решения
у(х)
() с точностью до о(
h
3
) равносильно сравнению
его с выражением
у
i
+1

по формуле (5), т е. с исправленным методом Эйлера.
Очевидно, чтобы (8) и (5) совпадали с точностью
O
(
h
3
),
от параметров нужно
потребовать выполнение следующей совокупности условий:






(29)

Полученная система условий содержит три уравнения относительно
четырех параметров метода. Это говорит о наличии одного свободного
параметра. Положим
с
2

 α(≠ 0). Тогда из (9) имеем:


В результате подстановки этих значений параметров в формулу (7)
приходим к
однопараметрическому семейству методов Рунге
-
Кутты
второго порядка.




(30)

Выделим из семейства методов (30) два наиболее про
стых и естеств
енных
частных случая:

при α
=1/2

получаем формулу


в котором узнаѐм
метод Хойна
(16), полученный ранее из других
соображений;

при
α
 1 из (30) выводим новый простой метод




(31)

который назовем
методом средней точки.


М
етоды
Р
унге
-
К
утты произвольного и четвертого порядков

Любой метод из семейства методов Рунге
-
Кутты второго порядка (30)
реализуют по следующей схеме. На каждом ша
ге, т.е. при каждом
i
=0, 1, 2,...,
вычисляют значения функции


а затем находят
шаговую поправку



прибавление которой к результату предыдущего шага дает при
ближенное
значение решения
у(х)
в точке
x
i
+1

=
x
i

+
h
:


Метод такой структуры называют
двухэтапным
по количеству
вычислений значений функции


правой части
уравне
ния (1)


на одном
шаге.

Анализ устройства методов Рунге
-
Кутты второго порядка позволяет
представить, в какой форме следует конструировать явный метод Рунге
-
Кутты
произвольного порядка. По аналогии с предыдущим для семейства методов
Рунге
-
Кутты
p
-
го

поряд
ка используется запись, состоящая из следующей
совокупности формул:




(32)

где
к
 , 3, …,
р
(для
р
-
этапного метода). Многочисленные па
раметры
с
k
,

а
к
,
b
kj
,
фигурирующие в формулах

(3), подбираются так
,
чтобы получаемое методом
(3) значение
у
i
+1
совпадало со значением разложения
y
(
x
i
+1
)
по формуле
Тейлора с погрешностью
O
(
h
p
+1
)
(без учета погрешностей, совершаемых на

предыдущих шагах).

Наиболее употребительным частным случаем
семейства методов

(32)
является следующий
метод Рунге
-
Кутты четвертого порядка
относящийся
к
четырехэтапным
и имеющий вид:





(33)


Не пытаясь воспроизвести выкладки, приводящие от о
бщей записи
семейства (3) при
р4
к конкретному методу (33), дадим геометрическое
толкование последнего.

Обратив внимание на то, что шаговая поправка Δу
i
, есть
средневзвешенная величина поправок
каждого этапа (с
весовыми коэффициентами

1/6, /6, /6, 1/6 соответственно), проанализируем,
как получаются эти поправки этапов. На первом этапе создается приращение
соответствующее шаговой поправке Эйлера,


это

очевидно. На рис  ему отвечает отрезок
ВС
вертикали

х
x
i
+1

(точка
В
получена орто
гональным проектированием

точки
А
на эту вертикаль).


Рис. . Геометрическая иллюстрация одного шага методов Рунге
-
Кутта четвертого
порядка

Так как точка
М
, благодаря свойству средней линии треугольника (см.
ΔАВС),
имеет ординату

определяет значение
f
(
M
),
служащее
(согласно связи
у
=
f
(
x
,
у)
и геомет
рическому смыслу производной) тангенсом
угла
А
в новом тре
угольнике с противолежащим этому углу катетом
Далее, аналогично, подсчитав
на вертикали
x
=
x
i
+1

откладываем следующую промежуточную (этапную) поправку
Вычислив величину
f
(
E
)=

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


длин отрезков

ВС,
BD
,
BE

и
BG
. Точка
Н
будет стартовой для следующего,
i
+1
-
го, шага
метода (33).

Заметим, что если первый этап, как уже упоминалось, соот
ветствует
применению явного метода Эйлера, то
четвертый


неявного, а второй и
третий


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


1.19

Занятие №
19
.
Пошаговый контроль точности. Метод Кутты
-
Мерсона
.

Цель
-

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

Нетрудно понять, что выведение надежных и, в то же вре
мя, простых и
эффективных оценок погрешности, гарантирую
щих

получение таблицы
значений решения
уу
(х)
заданной точности, является делом
малоперспективным, особенно для ме
тодов более
-
менее высоких порядков.
Поэтому главным спосо
бом отслеживания точности при реализации численных
процес
c
ов решения задачи Коши остае
тся применение различных
полуэмпирических правил, основанных на принципе Рунге.

Будем считать, что при использовании метода
р
-
го

порядка абсолютная
шаговая погрешность должна находиться в пределах ε>0. Тогда, согласно
принципу Рунге,
осуществляется счет п
о системе узлов

с шагом
h

и по системе узлов



с шагом
h
/. При четных
j

вторая система
будет совпадать с первой, т.е.
Переход от расчетной точки
x
i

с
приближенным значением решения в ней
у
i

к расчет
ной точке один раз
совершается за один шаг длины
h

и приводит к значению

другой раз


за два шага длины
h
/ («транзитом» через точку
со
значением

и дает значение


Поправка Ричардсона
в таком случае будет составлять величину




(34)

Если величина

меньше заданного
ε
, то можно считать, что ошибка
приближенного равенства
не превосходит
ε
. Если же

то следует уменьшить расчетный шаг
h
.

При условии
стоит
попытаться двигаться дальше с более к
рупным шагом (например, удвоить
h
).

Пример .

(продолжение примера 1). Посмотрим, что дает применение
принципа Рунге к нескольку простым методам численного решения того же
уравнения
y
'х
-
3
y

начальным условием
у(0)
 1. Из точки
х
0 перейдем в
точку
х
=0.2
за

один шаг h 0. четырьмя одношаговыми методами: явным и
неявным методами Эйлера, методом трапеций и методом Хойна


частным
случа
ев метода Рунге
-
Кутты второго порядка. С помощью полученных значе
-
ний

и найденных ранее теми же методами в примере 1 значений
подсчитаем поправки Ричардсона


при р1 для методов Эйлера и
р
для методов трапеций и Хойна. Эти
результаты, а также уточненные прибавлением к значениям
попра
вок
Ричардсона

приближенные значе
ния
решения
y
(0.) и их истинные
погрешности

сведем в следующую таблицу.

Метод







Эйлера (8)

0.4

0.11

0.62


-
0.04

≈0.07

Неявный Эй
лера (14)

0.675


-
0.0407

≈ 0.5936


-
0.01


-
0.05

Трапеций (15)

≈0.569

≈0.003

≈0.580


-
0.0001

≈0.003

Хойна
(16)

0.62


-
0.0102

≈0.5793

≈ 0.006


-
0.008

В эту таблицу последним столбцом помещен последний столбец из
таблицы результатов примера 1, содержащий погрешности значе
ний
Сравнение с ним столбца со значениями поправок Ричардсона показывает, что
эти
поправки хорошо отражают поведение погрешностей методов (хотя и не
дают основания считать их модули оценками погрешностей), а предпоследнего


эффективность уточнения по правилу Рунге
-
Ричардсона.

Грубо обозначенная здесь технология пошагового контроля точн
ости
численного интегрирования дифференциальных уравнений и автоматического
выбора расчетного шага при этом на основе двойного счета в такой
непосредственной форме говорит о ее значительной «дороговизне».
Действительно, предположим, что для решения задачи
(1)
-
() применяется
четырехэтапный

метод Рунге
-
Кутты четвертого порядка (33). Тогда вы
-
полнение одного его шага с контролем точности по правилу Рунге потребует 11
вычислений правой части уравнения (1) (по четыре получения каждого из
значений

и

минус одно общее для

весьма затратно.

Более «дешевый», но, возможно, менее строгий способ су
дить о том,
достаточно ли малым выбран шаг
h

расчетов по ме
тоду Рунге
-
Кутты четвертого
порядка (33),


это вычисление при каждом
i
 0,1, ,... величин


Считае
тся, что если величина Θ
i
, не превосходит нескольких со
тых, то
можно продолжить вычисления с данным шагом или пы
таться при переходе от
i

к
i
1 его увеличить; в противном случае шаг следует уменьшить, например,
вдвое.


Стремление повысить вычислительную эф
фективность при
вело к
появлению различных вычислительных версий методов Рунге
-
Кутты, благо для
этого в семействе методов (3) имеет
ся значительное число свободных
параметров. Основные сооб
ражения, положенные в основу этих версий,
таковы: нужно по
лучить

формулы из семейства методов Рунге
-
Кутты (3),
которые использовали бы одни и те же значения функции


правой части
уравнения (1)


и определяли бы разные кон
кретные методы одного порядка
(или смежных порядков, напри
мер, четвертого и пятого); при этом,
чтобы по
разности результа
тов подсчета приближенных значений решения по
выведенным близким формулам (с одним и тем же шагом
h
(!)) можно было
судить о точности одного из них.

Приведем один из таких методов, который называется
методом Кутты
-
Мерсона
или, иначе,
пятиэтапным мето
дом Рунге
-
Кутты четвертого
порядка,
а также
методов вложенных форм
.

На
i
-
ом шаге решения задачи (1)
-
() последовательно вычисляют:


После этого подсчитывают величину


и проводят сравнения. Если значение
R

окажется больше зада
н
ного
допустимого уровня абсолютных погрешностей
ε,
то шаг уменьшают вдвое
и возвращаются к началу второго эта
па, т.е. заново вычисляют
и
т.д. Если
R

ε
,
то считают
с точностью ε. При переходе к

следующему шагу делается проверка на возможность увеличить расчетный
шаг: если

то далее расчет ведется с шагом
h
:=
2
h
.


1.20

Занятие №
20
.
Многошаговые методы Адамса.

Цель
-

ознакомить студентов с многошаговыми методами Адамса
решения задачи Коши для обыкнов
енных дифференциальных уравнений.

Будем строить численные мето
ды решения начальной задачи




(1)





(2)

Будем считать, что уже найдено несколько приближенных значений
(
j

= 0,1,...,
i
) решения
уу(х)
зада
чи (1)
-
() на равномерной сетке
x
j
=
x
0
+
jh
,
и нужно по
лучить правило для вычисления очередного значения
Для вывода таких правил
используем интегро
-
интерполяционный
подход. А именно, проинтегрировав левую и правую части уравнения (1) по
промежутку
[
x
i

,

x
i
+1
],
в полу
ченном равенстве




(3)

под интеграл вместо функции
f
(
x
,
у(х))

подставим интерполирующий ее
многочлен
Хотя выражение функции
f
(
x
,
y
(
x
))
как функции одной
переменной
х
, вообще говоря, неизвестно, ее дискретные приближенные
значения
обозначаемые в дальнейшем для краткости
f
j
, при
j
=
1, 2, ...,
j

можно считать известными. В таком случае, дополняя эти известные
значения пока что неизвестным значением
можно построить
таблицу конечных разностей (табл. 1), служащую основой для образования

интерполяционных многочленов
к
-
й

степени для интерполирования назад из
точек

Таблица 1


Таблица конечных разностей для построения конечноразностных формул
Адамса


При интерполировании назад из узла х
i
, по второй интерпо
ляционной
формуле Ньютона имеем



(4)

(см.конечные разности, подчеркнутые в табл. 1 сплошной ли
нией), а из узла
x
i
+1

по той же формуле получаем многочлен




(5)

(использующий разности, подчеркнутые пунктиром).

Подстановка многочленов

в равенство (3) приводит к
формулам для вычисления очередного значения

вида




(6)

и






(7)

В результате применения к интегралам в (6) и (7) фор
мулы Ньютона
-
Лейбница получается два семейства методов (с параметром
k

N
0
), называемых
многошаговыми методами

Адамса.
Рассмотрим по отдельности каждое из
этих семейств.

Экст
раполяционные методы Адамса
-
Башфорта.
Чтобы подставить в
(6) многочлен (4), зависящий от переменной
сделаем в интеграле
замену переменной
в

соответствии с которой


Тогда формула (6) может быть переписана в виде



(8)

где




(9)

Таким образом, на основе (8) получается следующая конечно
-
разностная
формула, определяющая
экстраполяционный метод Адамса
-
Башфорта:




(10)


Посмотрим, что представляют собой наиболее

простые
частные случаи
метода Адамса
-
Башфорта,
соответствующие нескольким первым значениям
параметра
к
в формуле (8). Сразу заметим, что при фиксировании
к
 0,1, ,... в
(8) тем самым задается степень интерполяционного многочлена (нуле
вая,
первая, втора
я и т.д.) и, соответственно, число слагаемых, равное 1, , 3, в
правой части (9) (или, что то же, в скобках формулы (10)). Конечные разности в
получающихся при этом конкретных формулах будем раскрывать через
значения функ
ции, приводя формулы к виду, назыв
аемому иногда
ординатным.
Имеем:

при
к
=
0





(11)

при
к
= 1




(12)

при
к 




(13)

при
k
= 3








(14)

πормулы
(11), (12), (13)
и
(14)
определяют ме
тоды Адамса
-
Башфорта
соответственно первого, второго, третьего и четвертого порядков.
Относительно порядка метода (11) сомнен
ий нет: мы узнаѐм
метод Эйлера
(8).

Интерполяционные методы Адамса
-
Моултона.
В интеграле,
фигурирующем в формуле (7), делаем замену


x
=
x
i
+1
+
qh

и подставляем в него выражение

определяемое формулой (5).
Приходим к аналогичному (8) равенству


где




(15)

Отсюда следует конечноразностная формула
интерполяционно
го
метода Адамса
-
Моултона




(16)

Аналогично тому, как это делалось для методов Адамса
-
Башфорта, при
к 0,
1,
,3, т.е. фиксированием одного, двух, трех, четырех членов в
представлении (15) интеграла
получаем следующие частные формулы:

при
к
0




(17)

при
к
= 1




(18)

при
к



(19)

при
к
3








(20)

πормулы (17) и (18) определяют уже известные нам методы, а именно,
неявный метод Эйлера
(14) и
метод
тра
пеций
(15), имеющие первый и
второй порядки точности со
ответственно. Заметим, что оба эти метода

являются одношаговыми, а следующие за ними методы Адамса
-
Моултона (19)
и (0) третьего и четвертого порядков относятся, как легко ви
деть,
соответственно
к двухшаговым и трехшаговым методам. Та
ким образом,
для
интерполяционных методов Адамса
-
Моултона порядок шаговости на единицу
ниже порядка точности метода
(за тривиальным исключением, отвечающим
случаю
к
= 0).

Важное различие в экстраполяционных и интерпо
ляцион
ных методах
Адамса заключается в том, что первые из них явля
ются явными, а вторые


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

1.21

Занятие


21
.
Методы прогноза и коррекции
.

Под названием
методы прогноза и коррекции
(иначе
ме
тоды
предсказания и уточнения, предиктор
-
корректорные методы)
понимается
совместное применение явных и неявных методов одинакового или смежных
порядков. По явной формуле значение решения
у(х)
задачи (1)
-
() в текущей
(расчет
ной) точке
х
i
+1

прогнозируется, т.е. находится его, быть может,
достаточно грубое приближение
, а с помощью неявной форму
лы, в правую
часть которой подставляется спрогнозированное значение, оно уточняется
(корректируется). Пример приближен
ного вычисления
у(х
i
+1
)
по такой явно
-
неявной схеме у нас уже есть: парное использование явного ме
тода Эйлер
а для
предсказания и метода трапеций для уточнения (17).

Остановимся подробнее на методах прогноза и коррекции, базирующихся
на парах явных и неявных методов Адамса одинакового порядка. Обозначим
через
приближенное значение решения
y
(
x
i
+1
),
подсчитываемое

по явной
экстраполяционной формуле Адамса
-
Башфорта, и составим несколько пар из
рас
смотренных в предыдущем параграфе частных формул Адамса
-
Башфорта
(11), (1), (13), (14) и Адамса
-
Моултона (17), (18), (19), (0).

Имеем следующие
предиктор
-
корректорные ме
тоды Адамса:


первого порядка
(он
же явно
-
неявный метод Эйлера)


второго порядка


третьего порядка


четвертого порядка



(21)

Одним из главных достоинств методов прогноза и коррек
ции является
возможность
контролировать шаговую погреш
ность сравнением двух
полученных по явной и неявной форму
лам приближений к
y
(
x
i
+1

). Покажем,
как реализуется эта воз
можность для наиболее употребительного предиктор
-
корректор
ного метода Адамса четвертого порядка (1).

Вспо
мним, что первая из формул (1) была получена из общей формулы
Адамса
-
Башфорта (10), а вторая


из общей формулы Адамса
-
Моултона (16), в
которых последними бра
лись разности третьего порядка (подынтегральная
функция в ра
венстве (3) аппроксимировалась инте
рполяционным много
членом
третьей степени). Считая, что расчетный шаг
h

доста
точно мал и конечные
разности с ростом их порядка убывают, главные части шаговых погрешностей
формул Башфорта и Моултона четвертого порядка, в соответствии с (10) и (16),
характеризуются величинами

для явной и
для
неявной формул. Таким образом,

если наряду с введенным обозначением

обозначить через
прибли
женное значение
у(х
i
+1
),
получаемое по формуле
Адамса
-
Моултона четвертого порядка, то можно записать два прибли
женных
представления
y
(
x
i
+1
)
:





(22)

и




(23)

Отсюда видно, что если четвертые разности функции
f
(
x
,
у(х)) в
используемой части табл. 1 конечных разностей практически постоянны ( а это
можно связать с удачным выбором величины шага
h

при достаточном запасе
знаков в значениях
f
(
x
j
,
y
j
)),

то
во
-
первых, значения

и

дают
двусторонние приближения к точному решению
y
(
x
i
+1
),
а во
-
вторых, через
разность между зна
чениями
и

можно оценить точность каждого из них.

Действительно, приравнивая правые части приближенных равенств () и
(3) и отождествляя

имеем:


откуда


Подставляя последнее в (3), получаем приближенное

равенство




(2
4
)

Использование приближенной формулы (
4
) может быть двояким.
Переписав ее в виде


применяем это для пошагового контроля точности:


если
то полагаем
с точностью ε и переход
им
к следующему шагу
(
i:=
i
+1
),
иначе уменьшаем шаг
h

и снова подсчитываем

Другое назначение формулы (4)


это прямое приме
нение ее правой
части для получения уточненного значения:
полагаем





(25)

Наверное, есть смысл контроль точности делать на каждом шаге, а к
уточнению по формуле (5) прибегать при выводе оконча
тельных результатов.

Замечание

1. При выводе формулы (4) под

мы понима
ем значение,
соответствующее «чистому» методу Адам
са
-
Моултона чет
вертого порядка, т.е.



это точная реализация неявной формулы (0). Вторая же формула
предиктор
-
корректорного метода (1) соот
ветствует лишь одному
приближению к
по методу простых итераций, где в качестве начального
приближения берется
.
Поэтому примене
ния формулы (4) к методу
прогноза и коррекции (1) будут убеди
тельны в том случае, если его вторая
формула итерируется хотя бы один
-
два раза. Однако, чем больше таких
итераций, тем ниже вычислительная эффективность этого метода, в цело
м
весьма высокая по сравнению с мно
гоэтапными методами Рунге
-
Кутты.


1.22

Занятие №
22
.
М
етод

М
илна четвертого порядка

Цель

-

ознакомить

студентов с методо
м

Милна четвертого порядка
решения задачи Коши для обыкновенных дифференциальных уравнений.

Рассмотрим ещ
е один широко известный метод прогноза и коррекции


метод Милна.

Для вывода первой формулы Милна (т.е. формулы предска
зания)
проинтегрируем данное уравнение (1) на промежутке [
x
i
-
3
,
x
i
+1
] и в полученном
интегральном равенстве





(26)

подынтегральную функцию
f
(
x
,у(х))
заменим первым интерпо
ляционным
многочленом Ньютона
Р
3
(х), построенным по четы
рем узлам
с
предполагающимися уже извест
ными приближенными
значениями


Тогда, после замены переменной
на основании (6) имеем:


Отсюда, выразив конечные разности через значения функции, получаем
первую формулу Милна
(предсказания)




(27)

которую, очевидно, следует
отнести к экстраполяционным.

Главный член локальной погрешности формулы (7) на
ходим
интегрированием следующего (первого из неучтенных) слагаемого
интерполяционного многочлена Ньютона. Именно:


Считая четвертые разности примерно одинаковыми, опустим ин
д
екс у
функции
f

в записи
в результате получаем сле
дующее приближенное
представление решения в точке
x
i
+1

:





(28)

Вывод второй формулы Милна более прост. Проинтегрируем уравн
ение (1)
теперь на промежутке [
x
i
-
1
,
x
i
+1
] и в полученном равенстве


применим к интегралу простейшую формулу Симпсона. Имеем




(
29
)

Отбрасывая здесь остаточный член и заменяя значения решения
y
(
x
i
-
1
)
и
y
(
x
i
;) известными приближенными значениями
y
i
-
1

и
y
i
,
а стоящее в правой части
под знаком функции
f

неизвест
ное значение
у(х
i
+1
)
тем значением
которое
получается в результате вычислений по явной первой формуле Милна (7),
приходим ко
второй формуле Милна
(ут
очнения)




(3
0
)

являющейся интерполяционной.

Для вывода приближенной оценки шаговой погрешности воспользуемся
приближенным равенством

где
так же, как и в (8),


условная
запись практически по
стоянных
четвертых разностей. Исходя из точного равенст
-
ва (9), локальную погрешность получаемого с помощью формулы (30)
(возможно, с итерационной обработкой, см. за
мечание) приближенного значения
y
i
+1

можно приближенно охарактеризовать величиной
, т.е.





(3
1
)

Сравнение (
8
) и (3
1
) дает:


следовательно,





(3
2
)

Таким образом, при численном интегрировании
начальной задачи (1)
-
(2)
методом Милна четвертого порядка, опреде
ленным формулами (7) и (30), на
каждом
i
-
м шаге следует вычислять величину


и сравнивать ее модуль с величиной
ε
> 0 допустимой шаговой погрешности. Если
то за
у(х
i
+1
)
принимается получен
ное по второй формуле Милна значение
у
i
+1

(или его уточненное значение
); иначе шаг должен быть
уменьшен.

πигурирующая в приближенном равенстве (3) постоянная 1/9 примерно
вдвое меньше постоянной 19/70≈1/14 в аналогичном равенстве (4) для
предикто
р
-
корректорного метода Адамса четвертого порядка (
2
), что
характеризует метод Милна как несколько более точный при одинаковых
вычисли
тельных затратах.



Приложение 1

ПОЛИНОМЫ ЛЕЖАНДРА

Полиномы Лежандра являются специальными функциями, которые
применяются

при решении многих теоретических и прикладных задач. Полином
Лежандра
n
-
й степени можно определить с помощью производной
n
-
го порядка
следующим образом:





(1)

где
z

-

комплексная переменная.

В данном учебном пособии рассматриваются и используются полиномы
Лежандра для действительного аргумента
x
,
лежащего в интервале
x

[
-
1, 1].

С помощью определения (1) легко получить явные выражения полиномов
Лежандра действительного аргумента низших степеней
:




(2)

Графики перечисленных полиномов приведены на рис.1.

Все полиномы Лежандра
P
n
(
x
)
имеют следующие граничные значения:




(3)

Нетрудно
убедиться, что полиномы Лежандра четной степени являются
четными функциями и наоборот.

Важным для практических применений является свойство ортогональности
полиномов Лежандра:





(4)

где
Q
k
(
x
)
-

любой полином степени
k
, меньшей
n

(
k


n
).

Полиномы Лежандра подчиняются рекуррентному соотношению




(5)


которое, в частности, удобно для последовательного вычисления полиномы
высоких
степеней.



Рис.1. Графики полиномов Лежандра а)
n

=
0, 1, , б)
n

=
3, 4.



Приложение 
.

Параметры квадратурных формул Гаусса



Индивидуальное домашнее задание № 1

1. В приведенных задачах числа m, n, k вычислены с некоторой
погрешностью. Необходимо
вычислить и определить погрешность результата для
Х.

1.

2.

3.

4.


5.

6.


7.

8.

9.

10.

2
.
Ознаком
ьт
е
сь

с методами приближенного вычисления корней уравнений.

Найдите один действительный корень уравнения с точностью 10
-
5
. В ходе
решения ос
уществить следующие шаги:

2
.1. Отделить корень уравнения.

2
.2.
Вычислить с помощью программы значение отдельного корня
методами
: деление отрезка пополам, хорд, касательных, комбинированным
методом, методом итераций
. При использовании метода простых итерац
ий найти
решение при разных начальных приближениях. Результаты вычислений

занести в
табл
ицу
.

Вариант задания выбрать из табл. 1.1.



3
.
Найдите действительный корень уравнения с точностью 10
-
4
, на
интервале [
a
,
b
]. На первом этапе решения методом деления
пополам, уменьшать
интервал, содержащий корень, до тех пор, пока его длина не станет меньше 0,.
Потом, применить один из «более» быстрых методов.

1


11


2


12


3


13


4


14


5


15


6


16


7


17


8


18


9


19


10


20



ОТЧЕТ О

РАБОТЕ

Отчет должен содержать:

1. График исследуемой функции с интервалами отделения корней.


. Таблицы пошаговых расчетов корня уравнения.

3. Обоснованное заключение о преимуществах и недостатках использования
исследованных методов решения

применительно
к заданному уравнению

(для
задания 1)
.

4
.
Используя схему Гаусса
(
схема

единственного деления

и
схема полного
выбора
)

решить систему уравнений




5
.
Решить систему уравнений двумя способами


методом итераций и
методом Зейделя. Продолжать итерации до тех
пор, пока точность
приближенного решения не станет меньше 0,01.





Индивидуальное задание №.

1. πункция
y
=
f
(
x
)
задана таблицей значений:


Указания. Для вариантов 10
-

1 значения аргумента
x

предварительно
перевес
ти из градусов в радианы.

Даны контрол
ьные значения аргумента
x
1
=12;
x
2
=26;
x
3
=42.

a)

Написать подходящие для приближенного вычисления значений
y
1
=
f
(
x
1
),
y
2
=
f
(
x
2
),
y
3
=
f
(
x
3
)
интерполяционные многочлены Лагранжа первой и второй
степени. Получить эти значения.

b)

Составить алгоритм и написать программу

на языке высокого уровня, реали
-
зующую схему Эйткена вычисления с максимально возможной точностью
значения
y
=
f
(
x
)
в произвольной точке
x

промежутка
Пользуясь этим алгоритмом, вычислить приближенные зна
чения

c)

Сделать анализ результатов заданий
1, 2.

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

Для
всех вариантов




3.
Вычислить значения данной функции и ее п
р
о
зводной с помощью
интерполяционного полинома Лагранжа L
n
(x). В качестве узлов

интерполяции
взять:

1) равномерно распределенные точки на отрезке [; ];

) чебышевски
й набор узлов на отрезке [; ].

При табулировании функции вычислять ряд с точностью 10
-
6
.





Замечание.
При вычислении ряда
учесть, что каждый
последующий член ряда
a
n
+1

получается из предыдущего члена
a
n

умножением на
величину q
n
, т.е.
Это позволит избежать переполнения при
вычислении факториалов.

4.
Найти приближенные значения функции
при данных промежуточных

значениях аргумента с помощью кубического сплайна и визуализируйте

результаты сплайн
-
интерполяции.




Отчет должен содержать:




постановку задачи и исходные данные;



описание методов решения;



графики, полученных интерполяционных мног
очленов;



листинг программы.


Индивидуальное задание №3

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

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


узел таблицы)

Таблица 1.

Вариант

Задание 1.

Задание .

Таблица


х

Таблица


х

Используемая
формула

1

2

3.0

2

(
взять 5
последних
значений
)

5,3

Лагранжа

2

3

3.5

3
(взять 4
последних
значения)

6,7

Лагранжа

3

4

2.5

5

2,6

Ньютона

4

5

5.8

4

-
3,2

Ньютона

5

2

3.1

3

2,3

Ньютона

6

3

3.9

2

2,1

Ньютона

7

4

3.3

4
(
взять 5
первых
значений
)

-
0,8

Лагранжа

8

5

6.0

5
(
взять 4
первых
значения
)

3,8

Лагранжа

9

2

3.2

2
(
взять 5
первых
значений
)

2,9

Лагранжа

10

3

5.3

4

1,6

Ньютона

11

4

3.9

3

3,4

Ньютона

12

5

7.
2

5
(
взять 4
первых
значения
)

5

Лагранжа

13

2

4.4

5

6,2

Ньютона

14

3

3.6

3
(
взять 5
последних
значений
)

4,5

Лагранжа

15

4

2.2

4
(
взять 5
последних
значений
)

4

Лагранжа

16

5

6.8

2

3,7

Ньютона

17

2

3.4

3

5,6

Ньютона


18

3

3.7

4
(взять 4
последних
значения)

6,4

Лагранжа

19

4

1.8

5
(взять 5
первых
значений)

7,4

Лагранжа

20

5

7.6

2

4,5

Ньютона



Таблица .

x

f(x)=1/
x
lg
x
+
x
2

1,3

1,7776

2,1

4,5634

2,9

8,5694

3,7

13,8436

4,5

20,3952

5,3

28,2267

6,1

37,3387


Таблица 3.

x

f(x)=ln2,3
x
-
0,8/
x

1,2

0,3486

2,3

1,3180

3,4

1,8214

4,5

2,1592

5,6

2,4128

6,7

2,6156

7,8

2,7845




Таблица 4.

X

f(x)=2,1sin0,37
x

-
3,2

-
1,9449

-
0,8

-
0,6126

1,6

1,1718

4

2,0913

6,4

1,4673

8,8

-
0,2397

11,2

-
1,7698


Таблица 5.

x

f(x)=1,7
-
cos(0,4
-
0,7
x
)

2,6

2,1874

3,8

3,2888

5

3,9061

6,2

3,8209

7,4

3,2452

8,6

2,6949

9,8

2,6535


3
.
Вычислить
значения интеграла, используя квадратурные формулы:



левых прямоугольников,



правых прямоугольников,



центральных
прямоугольников,



трапеции,




Симпсона,



Ньютона,



Гаусса с двумя узлами.

Интеграл вычислить с точностью ε10
-
6
. Точность вычисления интеграла
определяется сравнением результатов при различном числе разбиений отрезка
интегрирования. Именно, точность ε считается

достигнутой, если


здесь
-

значение составной квадратурной формулы при разбиении отрезка на
N частей.

Отчет должен содержать:



постановку задачи и исходные данные,



описание методов решения и расчетные формулы,



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



листинг программы.

Варианты заданий.

1.

2.

3.

4.

5.

6.

7.

8.

9.

10.

11.

12.

13.

14.

15.

16.

17.

18.

19.

20.




Литература

1.

Бахвалов Н. С.

Численные методы: учеб. пособие для вузов/ Н. С.
Бахвалов, Н. П. Жидков, Г. М. Кобельков.
-

4
-
е изд.; Гриф МО.
-

Москва: БИНОМ.
Лаборатория знаний, 006.
-

636 с.
-

(Классический университетский учебник).
-

Библиогр.: с. 64
-
628.
-

Предм. указ.: с. 69
-
63
2.

2.

Вержбицкий В. М.

Численные методы: Линейная алгебра и нелинейные
уравнения: учеб. пособие для вузов / В. М. Вержбицкий.
-

2
-
е изд., испр.; Гриф
МО.
-

Москва: ОНИКС 1 век, 005.
-

431 с.: ил.
-

Библиогр.: с. 419
-
424.
-

Предм.
указ. с. 45
-
429.

3.

Демидови
ч Б. П.

Численные методы анализа: приближение функций,
дифференциальные и интегральные уравнения: учеб. пособие / Б. П. Демидович,
В. А. Марон, Э. З. Шувалова; под ред. Б. П. Демидовича.
-

Изд. 4
-
е, стер.
-

Санкт
-
Петербург: Лань, 008.
-

400 с.
-

(Классиче
ская учебная литература по
математике).
-

Библиогр. в конце гл.

4.

Киреев В. И.

Численные методы в примерах и задачах: у
чеб. пособие для
вузов
/ В. И. Киреев, А. В. Пантелеев.
-

Изд. 
-
е, стер.; Гриф УМО.
-

Москва:
Высш. шк., 006.
-

480 с.: ил.
-

(Прикладная

математика для ВТУЗов).
-

Библиогр.:
с. 477
-
480 .

5.

Пантина И. В. Вычислительная математика [Электронный ресурс]

учебник для вузов
/ И. В. Пантина, А. В. Синчуков.
-

2
-
е изд., перераб. и доп.
-

Москва: Синергия, 01.
-

175 с.

6.

Пирумов У. Г.

Численные метод
ы: учеб. пособие для вузов/ У. Г.
Пирумов.
-

2
-
е изд., испр. и доп.; гриф МО.
-

Москва: Дрофа, 003.
-

1 с.: ил.
-

(Высшее образование).
-

Библиогр.: с. 16.
-

Имен. указ.: с. 17.

7.

Срочко В. А.

Численные методы: курс лекций / В. А. Срочко.
-

Гриф
УМО.
-

Санкт
-
Петербург [и др.]: Лань, 010.
-

0 с.
-

Библиогр.: с. 00.

8.

Турчак Л. И
.

Основы численных методов: учеб. пособие / Л. И. Турчак, П.
В. Плотников.
-

Изд. 
-
е, перераб. и доп.; Гриф МО.
-

Москва: πИЗМАТЛИТ,
2005.
-

300 с.: ил.
-

Библиогр.: с. 90
-
292.
-

Прил.: с. 86
-
289.
-

Предм. указ.: с.
293
-
300.


9.

Шевцов Г. С.

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

Гриф УМО.
-

Москва: πинансы и статистика: ИНπРА
-
М, 008.
-

479 с.
-

(πинансы и
статистика).
-

Библиогр.: с. 473
-
474.
-

Предм. указ.: с. 475
-
479.


Приложенные файлы

  • pdf 18014163
    Размер файла: 6 MB Загрузок: 1

Добавить комментарий