
- •Оглавление
- •Предисловие
- •Лекция 1. Понятие погрешности решения систем алгебраических уравнений
- •1.1. Точные и приближенные числа, источники погрешностей, классификация погрешностей
- •Источники погрешностей
- •1.2. Абсолютные и относительные погрешности
- •1.3. Понятие о системе линейных уравнений
- •1.4. Матричные уравнения
- •Лекция 2. Точные методы решения систем линейных алгебраических уравнений (слау)
- •2.1. Формулы Крамера для решения системы линейных уравнений
- •2.2. Решение системы линейных уравнений методом последовательного исключения неизвестных (метод Гаусса)
- •Пусть дана система
- •Обозначим
- •Лекция 3. Приближенные методы решения слау
- •3.1. Понятие предела для векторов и матриц Пусть дана последовательность векторов
- •Аналогично, если имеется последовательность квадратных матриц
- •3.2. Метод простой итерации Пусть дана система линейных уравнений:
- •В матричном виде:
- •Условия сходимости итерационного процесса
- •3.4. Оценка погрешности приближенного процесса методом итерации
- •Лекция 4. Метод зейделя
- •4.1. Условия сходимости процесса Зейделя
- •Пусть дана линейная система
- •4.2. Оценка погрешности процесса Зейделя
- •4.3. Приведение системы линейных уравнений к виду, удобному для итераций
- •Лекция 5. Методы решения нелинейных уравнений
- •5.1. Итерационные методы решения нелинейных уравнений а). Метод простой итерации.
- •Число итераций при использовании этого метода
- •Если корень нас не устраивает, то мы находим
- •Г). Метод Ньютона.
- •5.2. Метод Рыбакова
- •5.3. Метод наискорейшего спуска
- •Лекция 6. Интерполирование и экстраполирование
- •6.1. Интерполяционный многочлен Лагранжа
- •Перепишем этот многочлен в другой форме
- •Окончательно получаем
- •Лекция 7. Конечные разности
- •Рассмотрим некоторые свойства конечных разностей
- •7.1. Первая интерполяционная формула Ньютона для равноотстоящих узлов интерполяции
- •7.2. Вторая интерполяционная формула Ньютона
- •7.3. Оценка погрешностей интерполяционных формул Ньютона
- •Лекция 8. Линейное интерполирование
- •8.1. Интерполирование по Эйткину
- •8.2. Интерполяция и приближение сплайном
- •Приближение линейными сплайнами
- •Лекция 9. Численное интегрирование функций
- •9.1. Постановка задачи
- •9.2. Формула прямоугольников
- •Лекция 10. Численные методы решения задачи коши для обыкновенных дифференциальных уравнений
- •10.1. Задача Коши. Общие замечания. Постановка задачи
- •10.2. Метод Эйлера
- •10.4. Метод Рунге-Кутта
- •Лекция 11. Разностный метод решения уравнений математической физики (метод сеток)
- •11.1. Метод сеток для уравнения параболического типа
- •Лекция 12. Метод сеток для уравнения гиперболического типа
- •12.2. Метод сеток для уравнений Пуассона и Лапласа
- •Образцы выполнения контрольных работ Контрольная работа № 5
- •Контрольная работа № 6
- •Решение уравнения методом Эйлера
- •Решение уравнения модифицированным методом Эйлера
- •Решение уравнения методом Рунге-Кутта
- •Аналитическое решение заданного уравнения
- •Сравнение точного решения и приближенных решений исходного дифференциального уравнения
- •Контрольная работа №7
- •Контрольная работа №8
- •Решение дифференциального уравнения методом Рунге-Кутта
- •Аналитические решения заданного уравнения
- •Примеры заданий для самостоятельной работы Задание 1
- •Задание 2
- •Задание 3
- •Задание 4
- •Задание 5
- •Варианты задания
- •Задание 6
- •Варианты задания
- •Задание 7
- •Варианты задания
- •Задание 8
- •Варианты задания
- •Задача для самостоятельного решения
- •Вопросы для самостоятельной подготовки
- •Библиографический список
Образцы выполнения контрольных работ Контрольная работа № 5
Вычислить значение интеграла
,
,
а = 0,
b
= 0,
.
Выберем формулу для приближенного вычисления заданного определенного интеграла. Для чего найдем по формулам (9.4), (9.6), (9.10) число n точек разбиения отрезка [0,1] на частичные, которые обеспечат требуемую точность при вычислении по формулам прямоугольников, трапеций и парабол соответственно. А затем остановимся на той из приближенных формул, для которой число n будет наименьшим.
Чтобы
воспользоваться формулами (9.4), (9.6),
(9.10), вычислим и оценим первую, вторую и
четвертую производные подынтегральной
функции
на отрезке [0,1]:
,
,
,
и так как функция f(x) и ее производные убывают на отрезке [0,1], то
,
,
.
Найдем n. Для формулы прямоугольников из (9.4) получаем
.
Для формулы трапеций из (9.6) получаем
,
.
Для формулы парабол из (9.10) получаем
;
.
Таким образом, наименьшего объема вычислений при одинаковой точности потребует формула парабол (9.8) – n = 2m = 8 (n должно быть четным), применяя которую и вычислим приближенно заданный интеграл.
По
числу n
= 2m
= 8
найдем шаг
интегрирования
.
Составим
таблицу (табл.1.) значений подынтегральной
функции в точках xi
= ih,
,
записывая ординаты с четными и нечетными
номерами в разные столбцы. В последней
строке таблицы запишем результаты
суммирования по этим столбцам. Вычисление
будем вести с четырьмя знаками после
запятой, а окончательный ответ округлим
до трех знаков после запятой.
Применяя формулу парабол (9.8), получим
Вычислим заданный интеграл по формуле Ньютона–Лейбница
.
Итак, требуемая точность вычислений достигнута.
Таблица 1
i |
xi=a+ih, |
1+2хi |
Значения уi=1/(1+2хi) |
||
При i=0, i=8 |
При четном i |
При нечетном i |
|||
0 |
0 |
1,0 |
|
|
|
|
0,125 |
1,250 |
|
|
0,80 |
|
0,25 |
1,50 |
|
0,6667 |
|
|
0,375 |
1,750 |
|
|
0,5714 |
|
0,50 |
2,00 |
|
0,50 |
|
|
0,625 |
2,250 |
|
|
0,4444 |
|
0,750 |
2,50 |
|
0,40 |
|
|
0,875 |
2,750 |
|
|
0,3636 |
|
1,0 |
3,0 |
0,3333 |
|
|
Суммы |
1,3333 |
1,5667 |
2,1794 |
Контрольная работа № 6
Пусть дано:
;
у(0)=1;
x
[1,2]
.
Шагом интегрирования h = 0,2 отрезок [0,1] разбивается на пять равных частей точками х0 = 0, х1 = 0,2, х2 = 0,4, х3 = 0,6, х4 = 0,8, х5 = 1,0.