Файл: Matlab Пособие.docx

Задача, требующая нахождения оптимального значения функции m переменных Q(Х) = Q(x₁, x₂, …, x_m), называется задачей многомерной оптимизации. Так же, как и для случая одномерной оптимизации, задача нахождения максимума функции сводится к задаче нахождения минимума путем замены целевой функции Q на -Q.

В постановке задачи безусловной оптимизации для Q(Х)=Q(x₁,x₂,…,x_m) требуется найти хотя бы одну точку минимума Х* и вычислить Q*=f(Х*). Точка Х* R^m называется точкой глобального минимума функции Q на множестве Х, если для всех Х  R^m выполняется неравенство Q(Х*)  Q(Х). В этом случае значение Q(Х*) называется минимальным значением функции Q на R^m. Точка Х*R^m называется точкой локального минимума функции Q, если существует такая  - окрестность U_этой точки (>0), что для всех Х  Х_= Х U_ выполняется неравенство Q(X*)  Q(X).

Для всякой непрерывно дифференцируемой функции Q достаточным условием того, что функция имеет точку минимума, является положительная определенность матрицы вторых частных производных (матрицы Гессе):

Известно, что для того чтобы матрица была положительно определена, необходимо, чтобы все угловые миноры были положительны. Так, для функции двух переменных Q(x, y) матрица Гессе имеет вид:

а достаточным условием существования минимума является выполнение неравенств:

Аналитический метод поиска минимума применяется только для ограниченного круга задач. В основном это связано с необходимостью решения системы нелинейных уравнений, которая, как правило, решается численными методами. Гораздо проще решать задачу многомерной оптимизации численными методами, например, такими как метод градиентного спуска с дроблением шага и методы наискорейшего спуска с аналитическим или численным выбором шага [1].

Вычисление экстремума функции нескольких переменных

z = f(x₁, x₂, …x_n)

в MatLab осуществляет функция:

[X, Z] = fminsearch(name, x0),

где:

name – имя функции, зависящее от n переменных;
x0 – вектор из n элементов, содержащий координаты точки начального приближения;
X – вектор из n элементов, содержащий координаты точки, в которой достигается минимум;
Z – значение функции в точке с координатами X.

Рассмотрим работу функции fminsearch() на примере определения минимума двумерной функции f(x1, x2) = x₁² + x₂² - x₂- 2x₁+ 2.

Рис. 2.6.1-1. Результат выполнения функции ezsurf()

Построим график (рис. 2.6.1-1) с использованием функции ezsurf(), аргументами которой служат: выражение функции, заключенное в одинарные кавычки, вектор изменения первой производной и вектор изменения второй переменной.

>> ezsurf ('x1.^2+x2.^2-x2-2*x1+2', [-2 2], [-2 2])

При нахождении минимума функции оформим многомерную функцию в виде m-функции (рис. 2.6.1-2).

Рис. 2.6.1-2. Целевая функция F(x)

Определим координаты точки минимума и значение функции в этой точке с использованием функции Matlab minsearch() (рис. 2.6.1-3).

Рис. 2.6.1-3. Использование функции fminsearch() для нахождения

минимума многомерной функции

2.6.2. Лабораторная работа по теме

«Технология решения задач многомерной оптимизации»

1. Вопросы, подлежащие изучению

Построение графиков функции от двух переменных F(x₁,х₂) средствами пакета Matlab.
Нахождение координат точки минимум функции F(x₁,х₂) с использованием встроенных функций пакета Matlab - fminsearch().

2. Общее задание

Изучите материал Темы 2.6 (п. 2.6.1).
Выберите индивидуальное задание из табл. 2.6.2-1.
Постройте график функции F(x₁,х₂).
Найдите координаты точки минимум функции F(x₁,х₂) с использованием встроенной функции fminsearch().
Получите значение функции в точке минимума.
Сохраните текст рабочего окна на внешнем носителе.
Представьте результаты работы преподавателю, ответьте на поставленные вопросы.
Выполните команду clear all.
Оформите отчет по выполненной работе.

Варианты индивидуальных заданий

Таблица 2.6.2-1

№	Функции для вычисления минимума
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30

4. Содержание отчета

В форме комментариев:

Название лабораторной работы
ФИО студента, номер группы
№ варианта
Индивидуальное задание

Протокол вычислений (сессии) в окне Command Window, снабженный необходимыми комментариями.

2.6.3. Контрольные вопросы по теме

Какую функцию называют многомерной?
Что является достаточным условием существования минимума для многомерной функции?
Назначение функции fminsearch(name, x0) и ее параметров.
Можно ли с использованием функции fminsearch() вычислить локальный максимум?
Что служит результатом выполнения функции fminsearch()?
Какие средства Matlab используются для описания целевой функции и как это влияет на параметры функции fminsearch()?

Список литературы

Васильев А.Н. Matlab. Самоучитель. Практический подход. – СПб.: Наука и Техника, 2012. – 448 с.
Дьяконов В.П. Полный самоучитель. – М.: ДМК Пресс, 2012.
Половко А.М., Бутусов П.Н. MATLAB для студентов. – СПб.: БХВ-Петербург, 2005. – 320 с., ил.
Кетков Ю.Л., Кетков А.Ю., Шульц М.М. Matlab 6.x: программирование численных методов. – СПб.: БХВ-Петербург, 2004. – 672 с., ил.