
- •Раздел 6.
- •Модели и алгоритмы решения задач численными методами с использованием математических пакетов Рекомендации по использованию лабораторного практикума
- •Тема 6.2. Лабораторная работа Методы решения нелинейных уравнений
- •1. Вопросы, подлежащие изучению
- •2. Задание
- •3. Варианты задания
- •4. Пример выполнения задания
- •Задание для решения нелинейных уравнений:
- •Отделение корней
- •«Ручной расчет» трех итераций методом половинного деления.
- •«Ручной расчет» трех итераций методом итераций
- •«Ручной расчет» трех итераций методом Ньютона
- •«Ручной расчет» трех итераций методом хорд
- •Уточним отделенный корень уравнения «расчетом средствами MathCad» с использованием функции root.
- •Тема 6.3. Лабораторная работа Интерполяция функций
- •Вопросы, подлежащие изучению
- •Задание
- •Варианты задания
- •Пример выполнения задания
- •Задание для решения задачи интерполяция функций
- •Результаты линейной интерполяция (узлы интерполяции z1 и z3), и квадратичной сплайн-интерполяция (узлы интерполяции z1, z2 и z3) и графики функций
- •Результаты интерполяция по формуле Лагранжа (узлы интерполяции z1, z2 и z3) и графики функций
- •Оценка погрешностей интерполяции функции f(X) в точках
- •Тема 6.4. Лабораторная работа Численное интегрирование
- •Вопросы, подлежащие изучению
- •Задание
- •Варианты задания
- •Пример выполнения задания
- •Задание для численного интегрирования:
- •Выражения первой и второй производной от заданной функции в символьном виде
- •Значение интеграла, полученное «расчетом средствами MathCad»
- •«Ручной расчет» определенного интеграла по заданному численному методу с шагом и ( и ) и оценка погрешности по правилу Рунге:
- •Оценка погрешности результатов интегрирования с шагом , приняв за точное значение интеграла величину
- •Тема 6.5. Лабораторная работа Методы решения обыкновенных дифференциальных уравнений
- •Вопросы, подлежащие изучению
- •Методы Рунге-Кутты различных порядков. Общие свойства. Погрешности методов.
- •Задание
- •Варианты задания
- •Пример выполнения задания
- •Задание для решения оду:
- •Решение аналитическое (точное) заданного оду, найдем методом разделения переменных.
- •Решение оду «расчетом средствами MathCad»
- •Оценка погрешностей приближенных решений и их погрешности относительно точного (аналитического) решения
- •Графики решений оду, полученных с использованием аналитической формулы ( ) и приближенных методов Эйлера ( ) и Рунге-Кутта 4-го порядка ( ).
- •Тема 6.6. Лабораторная работа «Одномерная оптимизация»
- •Вопросы, подлежащие изучению
- •Задание
- •Варианты задания
- •Пример выполнения контрольного задания
- •Задание для решения задачи одномерной оптимизации:
- •Результаты исследования функции:
- •«Ручной расчет» трех итераций
- •Значение координат точки минимума функции y(X) «расчетом средствами MathCad» с использованием функций root, Minimize и Minеrr .
- •Тема 6.7. Лабораторная работа «Методы оптимизации функций нескольких переменных»
- •Вопросы, подлежащие изучению
- •Задание
- •Варианты задания
- •Пример выполнения задания
- •Тема 6.8. Лабораторная работа «Аппроксимация функций»
- •Вопросы, подлежащие изучению
- •Задание
- •Варианты задания
- •Пример выполнения задания
- •Задание для решения задачи аппроксимации.
- •«Ручной расчет» коэффициентов линейной аппроксимирующей функции по методу наименьших квадратов и получение его аналитического выражения.
- •«Расчет средствами MathCad» для получения аналитических выражений линейной и квадратичной аппроксимирующих функций, и построение графиков точной и аппроксимирующих функций.
- •Литература
«Ручной расчет» трех итераций методом Ньютона
Для применения
метода Ньютона необходимо выбрать
начальное приближение к корню. Для
нашего уравнения 1- 3х +
cos(x)
= 0, известно, что
,
а
.
Выберем за начальное приближение к
корню тот конец отрезка, для которого
знак функции совпадает со знаком второй
производной:
.
Для выполнения
итераций по методу Ньютона воспользуемся
следующей рекуррентной формулой:
Получим
рекуррентную формулу для выполнения
трех итераций
Результаты вычислений представим в виде таблицы, структура которой аналогична табл. 6.2-2.
k |
Xk |
f(xk) |
0 |
1 |
-1.4597 |
1 |
0.6200 |
-4.62•10-2 |
2 |
0.6071 |
-6. 7875 •10-5 |
3 |
0.6071 |
-6.7875 •10-5 |
Оценим погрешность
результата, полученного методом Ньютона
после трех итераций, с использованием
формулы 6.2.3-11 в [2], где
а
«Ручной расчет» трех итераций методом хорд
Для применения метода хорд необходимо произвести выбор неподвижной точки, поскольку от этого зависит вид рекуррентной формулы. Неподвижен тот конец отрезка [a;b] , для которого знак функции f(x) совпадает со знаком ее второй производной. Тогда второй конец отрезка можно принять за начальное приближение к корню, то есть точку х0.
Рекуррентная формула метода хорд имеет вид [2]:
где
- неподвижная точка.
Поскольку для
функции f(x)=1–3x+cosx
<0,
то на отрезке [0;1]
неподвижной точкой является точка
x=b=1, так как
f(1)>0.
Тогда за начальное приближение к
корню положим x0=a=0,
и проведем три итерации по приближению
к корню с использованием следующей
рекуррентной формулы:
Результаты вычислений представим в виде таблицы, структура которой аналогична табл. 6.2-2.
n |
Xn |
f(xn) |
0 |
0 |
2 |
1 |
0.5781 |
0.1032549 |
2 |
0.6059 |
4.080772 •10-3 |
3 |
0.6070 |
1.590771•10-4 |
Погрешность
результата, полученного после трех
итераций, проведенных по методу хорд,
оценим по формуле 6.2-3-15 в [2], где
:
Уточним отделенный корень уравнения «расчетом средствами MathCad» с использованием функции root.
Тема 6.3. Лабораторная работа Интерполяция функций
Вопросы, подлежащие изучению
Постановка задачи интерполяции, основные понятия: узлы интерполяции, интерполирующая и интерполируемая функции.
Условие единственности решения задачи интерполяции.
Интерполяционные формулы Лагранжа и Ньютона, области их применения.
Конечные разности, их назначение и использование.
Правило выбора начальных узлов интерполяции для формул Ньютона.
Практическое правило определение степени интерполяционного многочлена.
Способы задания функции в MathCAD: аналитический, табличный, матричный.
Средство линейной интерполяции – функция linterp.
Средство кубической сплайн-интерполяции – функция interp.
Алгоритм интерполяции по методу Лагранжа.