
- •Раздел 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» для получения аналитических выражений линейной и квадратичной аппроксимирующих функций, и построение графиков точной и аппроксимирующих функций.
- •Литература
4. Пример выполнения задания
Задание для решения нелинейных уравнений:
уравнение:
;
методы решения нелинейных уравнений для «ручного расчета» – половинного деления, итерации, Ньютона и хорд (преподаватель задает студенту только один метод уточнения корня, но здесь в качестве примера приведены « ручные расчеты» при использовании всех четырех методов).
Отделение корней
|
Таким образом, на выбранном отрезке [0;1] , уравнение 1 - 3х + cos(x) = 0 имеет единственный корень.
«Ручной расчет» трех итераций (студент задается только один метод, но в качестве примера ниже приведены «ручные расчеты» с использованием методов половинного деления, итерации, Ньютона и хорд).
«Ручной расчет» трех итераций методом половинного деления.
Проведем необходимые
исследования для применения данного
метода. Известно, что метод половинного
деления сходится, если на выбранном
отрезке отделен один корень. Так как на
отрезке [0;1] функция
меняет знак (
)
и монотонна (f’(x)<0),
то условие сходимости выполняется.
Выберем за начальное
приближение середину отрезка
=0.5 и
проведем расчет трех итераций по
методу половинного деления.
-
>0 следовательно
<0 следовательно
<0 следовательно
Результаты вычислений представим в виде таблицы, структура которой аналогична табл. 6.2-1.
n |
a |
b |
f(a) |
f(b) |
(a+b)/2 |
f( (a+b)/2) |
b-a |
1 |
0 |
1 |
2 |
-1.459 |
0.5 |
0.377 |
0.5 |
2 |
0.5 |
1 |
0.377 |
-1.459 |
0.75 |
-0.518 |
0.25 |
3 |
0.5 |
0.75 |
0.377 |
-0.518 |
0.625 |
-0.064 |
0.125 |
4 |
0.5 |
0.625 |
0.377 |
-0.064 |
0.563 |
0.158 |
0.063 |
После трех итераций приближение к корню x3=0.563.
Оценка погрешности
результатов
.
«Ручной расчет» трех итераций методом итераций
Для применения
метода итераций приведем уравнение
f(x)=0 к виду
.
Тогда рекуррентная формула будет иметь
вид
.
Для сходимости процесса итерации
необходимо, чтобы
при
.
Если
то сходимость не обеспечена.
В случае, если
свободный х выразить
не удается, целесообразно воспользоваться
следующим приемом, позволяющим обеспечить
выполнение условий сходимости. Преобразуем
итерирующую функцию к виду
где
-
параметр, который может быть определен
по правилу: если
то
если
то
где
.
Преобразуем уравнение f(x)=0 к виду x = (cos(x)+1)/3 и проведем исследование итерирующей функции.
Выберем начальное
значение (в методе итераций x0
– произвольное значение из отрезка
[a;b]), например,
x0=0,
и с использованием итерирующей функции
выполним три итерации.
Результаты вычислений представим в виде таблицы, структура которой аналогична табл. 6.2-2.
-
к
Xк
f(xк)
0
0
2
1
0.6667
-0.2141
2
0.5953
4.21 • 10-2
3
0.6093
-7.9496 • 10-3
Оценим погрешность
результата, вычисленного методом
итерации после трех итераций, с помощью
выражения 6.2.3-4 в [2], где
:
.