
- •Раздел 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» для получения аналитических выражений линейной и квадратичной аппроксимирующих функций, и построение графиков точной и аппроксимирующих функций.
- •Литература
Тема 6.7. Лабораторная работа «Методы оптимизации функций нескольких переменных»
Вопросы, подлежащие изучению
Постановка задачи многомерной оптимизации.
Основные понятия: выпуклое множество, целевая функция, линия уровня, градиент, локальный и глобальный минимум.
Градиентные методы: метод с дроблением шага, метод наискорейшего спуска аналитический, метод наискорейшего спуска численный.
Средства MathCad для вычисления значений частных производных для функции нескольких переменных, а также для вычисления определителя и угловых миноров.
«Расчет средствами MatCad» координат точки экстремума многомерной функции с использованием функций Minimize и Minеrr.
Средства пакета MathCad для построения трехмерных графиков, графиков линий уровней и траектории спуска.
Задание
Выбрать индивидуальное задание по номеру варианта из табл. 6.7-1 - функцию f(x,y).
Проверить условия существования точки минимума для заданной функции f(x.y).
Провести «ручной расчет» по вычислению координат точки минимума функции f(x,y) аналитическим методом.
Выбрать начальную точку (x0, y0) для применения метода наискорейшего спуска.
Выполнить «ручной расчет» 3-х итераций аналитическим методом наискорейшего спуска. Результаты расчета свести в таблицу
-
k
x
y
λ
g1
g2
f(x,y)
1
2
3
4
Вычислить следующие погрешности после трех итераций:
где x*, y* - координаты точки минимума, найденные аналитическим методом, а f*=f(x*,y*) - значение исследуемой функции в точке минимума.
Решить задачу оптимизации «расчетом средствами MathCad» с использованием встроенных функций Minimize и Minerr.
Построить трехмерный график функции f(x,y).
Построить график линий уровня функции f(x,y) и траекторию спуска по результатам 3-х итераций ручного расчета, изобразив схематически линии уровня, проходящие через точки траектории. На графике указать точку минимума, найденную в п.3.
Варианты задания
Таблица 6.7-1
-
№
Функция
№
Функция
1
2 x2 + 3 y2 – 5 x + 6
16
6 x2 + 4 y2 – 5 x + 3 y –13
2
x2 + 2 y2 – 3 y + 7
17
5 x2 + y2 + x
3
3 x2 + y2 – 15
18
x2 + 4 y2 – 2 x + 3 y + 5
4
3 x2 + 5 y2 + x – 2
19
2 x2 + 5 y2 + 2 y + 3
5
2 x2 + 3 y2 + 2 x – 3 y
20
x2 + 3 y2 – x + 2 y + 7
6
5 x2 + 2 y2 + 3 x + 10
21
3 x2 + y2 – y + 3
7
4 x2 + 3 y2 – 3 y – 7
22
6 x2 + 3 y2 + 10
8
5 x2 + 6 y2 + 3 x – 2 y + 3
23
5 x2 + 4 y2 – 4 x – 11
9
3 x2 + y2 + - 3 x + y – 2
24
x2 + 2 y2 – x – y
10
6 x2 + 5 y2 – 10
25
3 x2 + 2 y2 – 5 y + 1
11
5 x2 + 2 y2 – 2 x
26
3 x2 + 4 y2 – 2 x + 3 y – 5
12
x2 + 2 y2 – 3 x + 5 y + 1
27
4 x2 + 5 y2 + 2 x – 4 y + 12
13
x2 + 4 y2 – 2 x
28
6 x2 + 3 y2 – 4 x + 17
14
4 x2 + 3 y2 + y + 3
29
x2 + 5 y2 – x + 2 y + 10
15
3 x2 + y2 + 3
30
3 x2 + y2 – 10