![](/user_photo/2706_HbeT2.jpg)
- •2. Классификация моделей
- •Классификация математических моделей
- •3. Численные методы решения
- •3.1. Решение нелинейных уравнений
- •3.1.1. Отделение корней или нахождение начального приближения к корню.
- •3.1.2. Метод деления отрезка пополам
- •3.1.3. Метод хорд
- •3.1.4. Метод касательных
- •3.2. Решение систем линейных уравнений
- •3.2.1. Метод Крамера
- •3.2.2. Решение слау с использованием MathCad.
- •Решение систем уравнений в среде MathCad
- •3.3. Численное интегрирование
- •3.3.1. Формула прямоугольников
- •3.3.2. Формула трапеций
- •3.3.3. Вычисление интеграла методом Монте-Карло
- •3.4. Численное дифференцирование
- •Задачи Коши для обыкновенного дифференциального уравнения.
- •4. Методы обработки экспериментальных данных
- •4.1. Приближение функций
- •Интерполяция обобщенными многочленами.
- •Полиноминальная интерполяция.
- •4.2. Элементы корреляционного и регрессионного анализа
- •4.3. Приближение функций с помощью метода наименьших квадратов
- •5. Оптимизация
- •5.1. Постановка задач принятия оптимальных решений
- •5.2. Линейное программирование
- •5.2.1. Постановка задачи и применение в лесном деле
- •Математическая статистика
3.4. Численное дифференцирование
Численное дифференцирование – функция трудно (невозможно) продифференцировать аналитически (Ex – функция задана таблицей). Формулы численного дифференцирования используются при решении дифференциальных уравнений; поиск решений нелинейных уравнений; поиск точек экстремума функций и т.д.
1.1 Формулы численного дифференцирования
1.1.1 Вычисление 1ой производной
Пусть f(x) дифференцируема в окрестности точки x. Из определения производной следует
Естественно использовать для вычисления 2 простейшие формулы:
и
(1.1)
Соответственно
выбору фиксированных значений
,
.
Здесь h>0 (шаг).
Для оценки
погрешностей
(*)
Воспользуемся формулами Тэйлора:
(**)
- некоторые точки в [x,x+h]
и [x-h,x]
соответственно подставляя (**) в (*)
(1.2)
То есть формулы (1.1) имеют первый порядок точности по h.
Геометрия интерпретируется f’(x)=tg
,
=
tg
+;
=tg
-
Естественно
предположить, что лучшим приближением
будет f’(x)
Подставляя в выражения для погрешности
.
Соответствующие разложения Тейлора:
,
получим
;
Таким образом центральная разностная производная аппроксимирует f’(x) со вторым порядком точности относительно h.
1.2 Вычисление 2ой производной.
Простая и
применяемая формула
Подставляя в
выражение для погрешности
Соответствующие разложения по формуле Тейлора:
Получим
Формула имеет 2ой порядок точности.
Задачи Коши для обыкновенного дифференциального уравнения.
Решением обыкновенного ДУ первого порядка y’(t)=f(t,y(t)) (*) называется дифференциальная функция y(t), которая, при подстановке в (*), обращает его в тождество.
Чтобы выделить из семейства решений ДУ(*) одно конкретное, задают начальные условие y(t0)=y0 (**).
Задачу нахождения при t>t0 решения y(t) ДУ (*), удовлетворяющего (**), называют задачей Коши.
Простейшие
дискретный аналог ДУ (*) представляет
собой уравнение
Покажем, что метод Эйлера имеет 1ый порядок аппроксимации. Известно, что
, где tn<
<tn+1
Учитывая равенство y’(tn)=f(tn,y(tn)) для погрешности аппроксимации получаем
Поэтому
,
где
,
т.е. метод имеет 1ый порядок
аппроксимации.
Использование формулы Тейлора.
y’(t) – известна = f(t,y(t))
y”(t) = f’t + f’yy’ – дифференцирование сложной функции
y”’(t) = f(2)tt + f(2)tyy’ + (f(2)yt + f(2)yyy’)f + f’y(f’t + f’yy’)
По мере роста порядка (р) усложняются выражения для производных. Недостаток метода Эйлера - значит-я погрешность – на практике редко используется. Желательно поправить расчетную формулу.
Пусть y(t) – решение ДУ y’(t)=f(t,y(t)), удовлетворяет условию y(tn)=yn
Пусть
(1.3)
- угловой
коэффициент секущей, проходящей через
точки (tn,
y(tn))
и (tn+1,
y(tn+1))
графика функции y(t).
Ясно, что «метод» yn+1
= yn +
hKn
имеет нулевую локальную погрешность.
Следовательно нужно научиться вычислять
значение Kn.
Интегрируя и используя формулу Ньютона
– Лейбница
приходим к равенству
(1.4)
Из (1.3) и (1.4)
следует Kn
=
Примечание. Для приближенного вычисления интеграла формулы прямоугольников:
приводит к методу Эйлера.
Но больший порядок точности имеет формула трапеций:
Итого приходим к правилу трапеций:
Если подставим в правую часть значение yn+1 «предсказанное» методом Эйлера, получим в результате метод Эйлера-Коши:
Этот метод относится к методам прогноза и коррекции.