
- •Часть 1 Числовые методы
- •Раздел 1
- •Тема 1.1 Точность вычислительного эксперимента на эвм.
- •§ 1.1.1 Приближенные числа
- •§ 1.1.2. Погрешность приближенных чисел.
- •§ 1.1.3. Погрешность математических операций над приближенными числами.
- •§ 1.1.4 Характеристики методов вычислений.
- •Раздел 1
- •Тема 2. Исчисление конечных разностей. Разностные уравнения.
- •§ 1.2.1 Табличное представление функций. Система обозначений. Операторы.
- •§ 1.2.2. Разностные операторы. Таблицы конечных разностей.
- •§ 1.2.3. Разностные уравнения
- •Раздел 2
- •Тема 2.1. Прикладные задачи интерполяции и аппроксимации и их алгоритмизация. Методы и алгоритмы решения линейных алгебраических уравнений
- •§ 2.1.1. Основные понятия теории приближения функций.
- •§ 2.1.2. Интерполяция многочленами с неравномерными и кратными узлами.
- •Можно и не вычислять коэффициенты Pn(X), а записать интерполяционную формулу, в виде уравнения прямой, проходящей через две координаты точки на плоскости, координаты которых известны.
- •§ 2.1.3. Аппроксимация функций. Метод наименьших квадратов.
- •Раскроем знак суммы относительно неизвестных a0 и a1
- •§ 2.1.4 Методы и алгоритмы решения системы линейных алгебраических уравнений.
- •Раздел 2. Методы и алгоритмы многочленного и немногочленного приближения
- •Тема 2.2.Прикладные задачи численного дифференцирования и интегрирования.
- •§ 2.2.1 Методы численного дифференцирования
- •§ 2.2.2.Методы численного интегрирования
- •Метод трапеций
- •Метод Симпсона
- •Раздел 2
- •Тема 2.3 Системы линейных дифференциальных уравнений и их решение на эвм
- •§ 2.3.1 Основные понятия и определения
- •§ 2.3.2 Алгоритмы решения задачи Коши
- •§ 2.3.3 Алгоритм решения краевой задачи оду.
- •Раздел 2
- •Тема 2.4 Дифференциальные уравнения в частных производных и их численное решение на эвм.
- •§ 2.4.1. Классификация дифференциальных уравнений в частных производных
- •§ 2.4.2 Метод сеток и алгоритмы его реализации
- •§ 2.4.3 Алгоритмы решения параболического уравнения
- •Раздел 3 Численное определение нулей и экстремумов функций
- •Тема 3.1. Методы и алгоритмы нахождение нулей функций
- •§ 3.1.1. Численное решение рациональных и трансцендентных уравнений
- •Пример 1
- •Метод простых итераций
- •Пример 3
- •§ 3.1.2. Нахождение действительных и комплексных корней многочленов
- •§ 3.2.1. Основные понятия и определения
- •Пример 1
- •§ 3.2.2 Методы решения задач оптимизации
- •§ 3.2.3. Задача поиска одномерного, локального безусловного оптимума.
- •§ 3.2.4 Задачи многомерной оптимизации
§ 3.2.4 Задачи многомерной оптимизации
Как правило, целевая функция сложных технических систем зависит от многих проектных параметров,
поэтому отыскание оптимума значительно усложняется. Существует, как и для задачи одномерной безусловной
оптимизации, несколько методов отыскания многомерного оптимума.
Если уравнение целевой функции дважды дифференцируемо, то задача решается аналитически.
Если можно получить значение функции, и ее первых производных, то можно применить метод численного решен-
ния системы нелинейных уравнений.
Если можно вычислять только f(X), f' ′ (X), f " (X) только численно, то
применяют поисковые методы. Различают регулярные и случайные поисковые методы. В этом будут
рассматриваться только регулярные методы поиска.
Аналитический способ отыскания оптимума сводится к определению частных производных
целевой функции. Если приравнять их к нулю, то образуется система нелинейных уравнений.
Когда она поддается аналитическому решению, то находят значения проектных
параметров, при которых они принимают оптимальные значения.
Однако, аналитическое решение нелинейных систем уравнений редко возможно и поэтому используются
численные методы и методы поиска.
Различают 3 типа регулярных методов поиска оптимума многомерных ЦФ.
1. Методы нулевого порядка - для поиска используются только значения целевой функции F(X).
2. Методы первого порядка - кроме F(X) используются и F'(X).
3. Методы второго порядка - используются значения F(X), F'(X), F"(X).
Функцию нескольких переменных f(x) = f(x1,x2,...,xn) трудно изобразить графически уже для
n=2 . При n=2 функцию отображают на плоскости линиями равного уровня. Линией равного
уровня f(x) называют геометрическое место точек, в которых f(x) принимает одно и то же
значение. Так, например, для функции
f(x)
=
(1)
на множестве
значений переменных
и
можно построить линии равного уровня
следующим образом.
Зададимся значением функции fi = const .
Решим уравнение (1) относительно x2:
x2(x1+1)-
x
=
fi
.
Будем задавать значения x1 , вычислять x2 и на пересечении x1 с x2 фиксировать точку, где f = const .
x2
Для f1
= -1
x1 x2
+1 0 -1
0,5 0,5
1 0
Для f2 = 0
x1
x2
-
1
0 +1
x1 0 0
0,5 0,7
1,0 0,5
0,25 0,05
0,75 0,321
- 0,25 0,083
- 0,625 1,042
-1
Рис. 4.1
Линии равного уровня функции позволяют оценить область оптимума и его характер: седло, овраг и т.д.
Для поиска оптимума многомерной функции необходимо знать её градиент и составить (затем вычислить) матрицу Гессе (H).
Градиент f(x0) определяется как
(2)
Матрица Гессе
, где
(3)
Для функции одной переменной градиент равен первой производной, а матрица вырождается в скаляр, равный второй производной.
Найдём
градиент f(x)
(1) в точке
.
Для этой функции легко найти экстремум:
x1 = -1 ; x2
= 3 .
<
0 , это максимум.
Существует несколько численных методов поиска максимума функции нескольких переменных, основанных на построении улучшающей последовательности решений. Любой алгоритм построения последовательности улучшающихся решений состоит из правила перехода от решений xi , полученного на I -том шаге алгоритма, к решению xi+1 и из правила остановки.
Для любого метода точки связаны соотношением
xi+1
= xi +
, где
-
вектор направления поиска;
-
характеризует величину шага;
-
модуль произведения характеризует
длину шага вдоль выбранного направления,
и если
-
нормирован так, что его модуль равен 1,
то
-
скаляр.
Правило остановки основано на неулучшаемости некоторого полученного в ходе поиска решения.
Чаще всего методы расчёта максимума функции нескольких переменных делят на три типа, в зависимости от того, какие характеристики функции используются для выбора и .
В методе нулевого порядка используются только значения функции. В этом методе используются несколько видов алгоритмов.
Для f (x1,x2) = 0 их можно схематически изобразить так:
x2
x2
x2
x*
M
C
2
x*
x2
x2
x’
1
1 1’
x1
x1
0
0 0’
x0
0’ x0
A x0
x1 x1 x1
Рис. 4.2 а) б) в).
На рисунках не нанесены линии равного уровня, чтобы не усложнять картинку.
В первом случае (рис.4.2а) из точки x0 в точку x* можно попасть, двигаясь поочерёдно по направлениям, параллельным осям x1 и x2 .
Это алгоритм поочерёдного изменения переменных (Гаусса-Зейделя). При движении по координате x1 фиксируют координату x2 и ведут одномерный поиск. Алгоритм обобщается на n переменных.
Во втором случае (рис.4.2б) используется идея параллельных касательных (сопряжённых направлений).
Из некоторой начальной точки по направлению координаты x1 осуществляют одномерный поиск и попадают в точку x0 . Из точки x0 ,как из начальной, методом координатного изменения аргументов находят точку x1. Осуществляют дальнейшее движение по направлению
-
такое значение аргумента функции f(x),
при котором оно максимально на множестве
значений. Т.е. переход из точки 1 в точку
2 осуществляется быстрее.
В третьем случае алгоритм называют симплексным. Для n=2 симплекс представляет собой треугольник.. Симплексом в общем случае (при n переменных) называют выпуклый многоугольник с n+1 вершинами, не лежащих в подпространстве размерности меньшей n. Симплекс называют регулярным (правильным), если расстояние между его вершинами одинаково.
Важным свойством симплекса является тот факт, что в результате отображения любой из вершин относительно противоположной грани, получается новый симплекс. Все вершины этого симплекса, кроме отображённой вершины, совпадают с вершинами предыдущего симплекса.
x2
Если обозначить отображаемую
вершину x0 , две
остальные
x2 иx1 , а отображённую , то координату последней
x1 x2
можно определить из выражения
:
x0
,
где
;
>0
множитель.
x1
При
>1
симплекс растянется;
=2 получит зеркальное отображение;
0< <1 сожмётся;
Известно несколько видов симплексных алгоритмов:
Алгоритм регулярных симплексов (рис.4.2в).
Строится начальный симплекс S1.
Вычисляются значения функции f0 в каждой из его вершин.
Находится худшее (минимальное значение f0) из его вершин.
Производится операция отображения худшей вершины относительно противоположной грани.
Затем процедура циклически повторяется до выполнения правила остановки. Остановить поиск можно в случае, когда отображённая вершина в новом симплексе снова окажется худшей. Для уменьшения погрешности можно далее уменьшить симплекс.
Недлером и Мидом в 1964 году был предложен метод деформируемого симплекса, в котором симплекс меняет свою форму от цикла к циклу.
Для функции с сильно взаимосвязанными переменными, когда линии уровня функции имеют характер граблей. Хуком и Дживсом в 1961 году был предложен метод конфигураций.
Каждый цикл алгоритма состоит из 2-х операций:
оценки характера функции f0 вблизи текущей точки поиска;
движения в направлении роста f0.
При n =2
Здесь верхний индекс – номер точки, а нижний – координата.
Из
исходной точки можно задавшись
координатами
и
числом k = h
построить исходный симплекс, определить
координаты точек 1 и 2.
Точка
1 (
+
h;
-h);
Точка
2 (
-h).
При этом симплекс получится нерегулярным, т.к. сторона противоположная x0 будет равна 2h, а две другие
Для
построения регулярного симплекса надо
выбрать длину стороны h,
шаг по x1, равный
h/2, а по x2
.
Тогда сторона противоположная x0
будет h, а две другие равны
ей, т.к.
.
В этом случае
Точка 1 ( + h ; -h)
Точка 2
(
;
)
Например, для функции
f
(x1 ; x2)
= x1x2+x2-
в точке
=1,
=
-1 можно построить вершину симплекса, в
дальнейшем взяв приращение координат
k=0,5.
x2
1
-1 0,5 1 1,5
x1
-0,5 x1 x2
-1
x0
x2
a = 0,5
h
=
=0,87
=1, = -1
=0,5,
=
-0,5
=1,5, = -0,5
1
наихудшая
=
-1,
=0,5 = -1,5
Выберем сторону симплекса: а = 0.25; a / 2 = 0,125, тогда шаг по x1= =0,43
=1 = -1
=0,875 = -0,57
=1,125 = -0,57
Затем снова повторяются вычисления функции в вершинах симплекса, пока не выполнится правило остановки.
Методы деформируемого многогранника и градиента, подробно изложены в методических указаниях по выполнению курсовой работы.
Тестовые задания для самопроверки усвоения темы.
1. А) Xni < Xi < Xbi при i [1:n],
Б) ψ(X) = 0;
В) ψ(X) > 0. Где выражение для прямого ограничения ЦФ?
2. Классифицировать
Метод поиска используемые значения
а) второго порядка; 1) целевой функции F(X).
б) нулевого порядка; 2) F(X) и F'(X).
в) первого порядка; 3) значения F(X), F'(X), F"(X).
3. Неравенство F(x)-F(X*)< 0 или F(x)-F(X*)> 0 характеризует оптимум в точке X*?
4.
,
это выражение для поиска
оптимума симплексным методом?
5.
-
это операция растяжения при симплексном
методе решения задачи оптимизации?