- •Н. Б. Борковский Математические методы в экологии
- •Часть I. Численные методы
- •Тема 1 9
- •Тема 2 14
- •Тема 3 30
- •3.0 Постановка задачи 30
- •1. Постановка проблемы
- •2. Выбор или построение математической модели
- •3. Постановка вычислительной задачи
- •4. Предварительный анализ свойств вычислительной задачи
- •5. Выбор или построение численного метода
- •6. Алгоритмизация и программирование
- •7. Отладка программы
- •8. Счет по программе
- •9. Обработка и интерпретация результатов
- •10. Использование результатов и коррекция математической модели
- •1.0 Источники приближенных чисел
- •1.1 Источники и классификация погрешностей
- •1.2 Абсолютная и относительная погрешности вычисления
- •1.3 Погрешности арифметических операций
- •1.4 Обратная задача теории погрешностей
- •Постановка задачи
- •Метод половинного деления (метод «вилки», метод дихотомии)
- •Метод Ньютона (метод касательных)
- •Метод хорд (метод пропорциональных частей)
- •Метод итерации (метод последовательных приближений)
- •Постановка задачи
- •Решение систем линейных уравнений методом итерации
- •Решение систем нелинейных уравнений методом итерации
- •Решение систем линейных уравнений – «подводные камни»
- •4.0 Интерполирование – постановка задачи
- •4.1 Интерполирование полиномами
- •4.2 Аппроксимация функций. Метод наименьших квадратов
- •4.3 Интерполирование сплайнами
- •4.4 Интерполирование – примеры и «подводные камни»
- •4.5 Сглаживание данных (скользящее среднее)
- •5.0 Постановка задачи
- •5.1 Формула трапеций
- •5.2 Формула Симпсона
- •5.3 Формула Ньютона–Котеса
- •5.4 Формула Гаусса
- •5.5 Связь формул интегрирования
- •5.6 Численное интегрирование. Сравнение методов
- •5.7 Численное интегрирование – «подводные камни»
- •6.0 Введение
- •6.1 Обыкновенные ду. Метод Эйлера
- •6.2 Обыкновенные ду. Основные методы решения
- •6.3 Обыкновенные ду. Методы прогноза и коррекции
- •6.5 Обыкновенные ду. Сравнение методов
- •6.6 Оду. Влияние неточности исходной информации
- •6.7 Оду. Представление о конечных разностях
- •6.8 Оду. Краевые задачи – введение
- •8.0 Ду в частных производных. Постановка задачи
- •8.1 Ду в частных производных. Эллиптические уравнения
- •8.2 Ду в частных производных. Гиперболические уравнения
- •8.3 Ду в частных производных. Параболические уравнения
- •9.0 Методы безусловной оптимизации. Введение
- •9.1 Методы безусловной оптимизации. Классификация методов
- •9.2 Методы безусловной оптимизации нулевого порядка
- •9.3 Методы безусловной оптимизации первого порядка
- •9.4 Методы оптимизации второго порядка
- •Тема 10
- •10.0 Метод Монте–Карло. Некоторые задачи
- •10.1 Проблема получения случайных чисел
- •10.2 Общая схема метода
- •10.3 Метод Монте–Карло. Пример
- •10.4 Вычисление кратных интегралов
- •Число испытаний n не зависит от размерности интеграла i0
1.1 Источники и классификация погрешностей
Обычно выделяют следующие источники погрешности решения задачи.
Математическое описание задачи является неточным. В частности, не удалось надлежащим образом описать задачу с помощью набора уравнений, неточно заданы исходные данные описания и др. Это неустранимая погрешность.
Применяемый для решения метод часто не является точным: получение точного решения задачи требует неограниченного или неприемлемо большого числа арифметических операций, поэтому вместо получения точного решения приходится прибегать к приближенному. Это погрешность метода.
При выполнении арифметических операций на компьютере или любым другим образом, как правило, производятся округления (это же относится к вводу чисел в память компьютера и выводу полученных результатов). Это вычислительная погрешность.
1.2 Абсолютная и относительная погрешности вычисления
Е
сли
a – точное значение некоторой
величины, а a* – известное приближение
к нему, то абсолютной погрешностью
приближения a* называют величину
Δ(a*), удовлетворяющую неравенству
О
тносительной
погрешностью называют величину
(a*), удовлетворяющую неравенству
Относительную погрешность часто выражают в процентах.
Рассмотрим способы нахождения погрешностей результата при выполнении арифметических операций.
Пусть X и Y – «точные» значения двух величин. В вычислительном устройстве (после округления) они представлены как x и y.
М
ожно
записать
Здесь a и b – положительные числа.
О
чевидно,
что
Вывод: при сложении (вычитании) двух чисел максимальная погрешность суммы (разности) равна сумме максимальных погрешностей исходных двух чисел.
Рассмотрим теперь случаи произведения и деления приближенных величин.
Перепишем соотношения
в виде
Пренебрегая малыми величинами (a/x)(b/y), получаем для XY
Д
ля
частного находим
Вывод: при умножении (делении) двух чисел максимальная относительная погрешность произведения (частного) равна сумме максимальных относительных погрешностей исходных двух чисел.
Представим сказанное в несколько более формализованном виде.
1.3 Погрешности арифметических операций
Погрешность суммы
П
усть
задана функция
Для абсолютной погрешности получаем
Относительная погрешность
Погрешность разности
Пусть задана функция
Для абсолютной погрешности получаем
Относительная погрешность
Если
и
близки, относительная погрешность
может оказаться много больше
и
.
Погрешность произведения
Пусть задана функция
Для абсолютной погрешности получаем
Относительная погрешность
Погрешность частного
Пусть задана функция
Д
ля
абсолютной погрешности получаем
Относительная погрешность
1.4 Обратная задача теории погрешностей
Рассмотренные выше оценки относились к «Прямой задаче погрешностей» – т. е. нахождению погрешности функции по погрешностям аргументов.
Обратная задача теории погрешностей – нахождение допустимой погрешности аргументов, при которой погрешность значений функции будет не более заданной величины .
В общем случае обратная задача математически некорректна. Возможно, однако, рассмотрение некоторых частных случаев.
В
рассмотрении используем неравенство
Накладывается
условие
.
Рассмотрим три частных случая.
Если функция зависит от одной переменной, т. е. n = 1, то
.При n > 1 если погрешности переменных одинаковы, (x1 ) = … = (xi ) = , то
.При n > 1, если вклады в погрешность результата одинаковы,
,
то
.
Тема 2
Решение алгебраических и трансцендентных уравнений
Необходимость решения уравнений возникает при рассмотрении большого числа разнообразных задач. При этом зачастую приходится прибегать к численным методам, поскольку даже в случае алгебраических уравнений порядка n общие решения найдены лишь для n < 5. Что касается трансцендентных уравнений, то общие решения, как правило, отсутствуют вовсе.
