- •Оглавление
- •Предисловие
- •Лекция 1. Понятие погрешности решения систем алгебраических уравнений
- •1.1. Точные и приближенные числа, источники погрешностей, классификация погрешностей
- •Источники погрешностей
- •1.2. Абсолютные и относительные погрешности
- •1.3. Понятие о системе линейных уравнений
- •2.1. Формулы Крамера для решения системы линейных уравнений
- •2.2. Решение системы линейных уравнений методом последовательного исключения неизвестных (метод Гаусса)
- •Пусть дана система
- •1.4. Матричные уравнения
- •Лекция 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
- •Варианты задания
- •Задача для самостоятельного решения
- •Вопросы для самостоятельной подготовки
- •Библиографический список
Предисловие
Прогресс вычислительной техники и методов программирования позволяет точнее моделировать физические и технические процессы для оптимизации и эффективности производственной деятельности человека, позволяет прогнозировать протекание исследуемых процессов в других условиях и при других параметрах. Как известно, существует всего 8 физических задач, которые решаются точно, остальные решаются лишь приближенно.
Кроме того, данные экспериментов, необходимые для решения практических задач, тоже даны в виде приближенных чисел. Поэтому каждый инженер, который по роду своей деятельности постоянно сталкивается с вычислительными задачами, возникающими при исследовании физических и технических проблем, должен иметь твердые навыки работы с приближенными числами и стандартными численными методами.
В данном учебном пособии рассматриваются численные методы решения задач наиболее часто встречающихся в инженерной и научно-технической практике, а именно: приближенное решение систем линейных алгебраических уравнений и приближенное решение нелинейных уравнений.
Рассматриваются различные методы интерполяции и экстраполяции для приближенного вычисления значений функций. Дан метод интерполяции сплайнами. Авторы планируют продолжить рассмотрение численных методов дифференциального и интегрального исчисления.
Лекции содержат достаточно полное изложение основных вопросов курса вычислительной математики, соответствующих требованиям к минимуму обязательной программы по подготовке дипломированных специалистов.
Лекция 1. Понятие погрешности решения систем алгебраических уравнений
1.1. Точные и приближенные числа, источники погрешностей, классификация погрешностей
Приближенным числом a называется число, незначительно отличающееся от точного числа A и заменяющее последнее в вычислениях. При решении задач вручную или на компьютере, мы получаем числовой результат, который не является точным, так как при постановке задачи и в ходе вычислений возникают погрешности. Поэтому любая задача, связанная с массовыми действиями над числами, может быть решена с той или иной степенью точности. Поэтому при постановке задачи должна быть указана точность ее решения, т.е. задана погрешность, максимально допустимая в процессе всех вычислений.
Источники погрешностей
Неточное отображение реальных процессов с помощью математики, когда рассматривается не сам процесс, а его математическая идеализация. Поэтому не всегда реальные условия могут быть хорошо описаны математически.
Приближенные значения величин, входящих в условия задачи (погрешности исходных данных).
Замена бесконечных процессов, пределами которых являются искомые величины, конечной последовательностью действий.
Округление исходных данных, промежуточных или окончательных результатов из-за малой разрядности машины.
Кроме указанных выше случаев, погрешности могут появляться в результате действий над приближенными числами. Полная погрешность является результатом сложного взаимодействия всех видов погрешности. При решении конкретных задач те или иные погрешности могут отсутствовать или мало влиять на образование полной погрешности. Однако для полного анализа необходимо учитывать все их виды.
Погрешности можно подразделить на три большие группы:
Исходные (или неустранимые) погрешности, к которым относятся погрешности, возникающие в результате приближенного описания реальных процессов и неточного задания исходных данных. Эти погрешности проходят через все вычисления и являются неустранимыми.
Погрешности округления, которые появляются в результате округления исходных данных, промежуточных и окончательных результатов.
Остаточные погрешности возникают в результате замены бесконечных процессов конечной последовательностью действий.
Оценка погрешности может быть произведена с помощью абсолютной погрешности, с помощью относительной погрешности, с помощью остаточного члена, с помощью статических оценок.
При работе с приближенными величинами, вычислитель должен уметь:
- давать математические характеристики точности приближенных величин;
- зная степень точности исходных данных, оценивать степень точности результатов;
- брать исходные данные с такой степенью точности, чтобы обеспечить заданную точность результата;
- правильно организовать вычислительный процесс, чтобы избавить его от тех выкладок, которые не окажут влияния на точные цифры результатов.
