- •Часть 1 Числовые методы
- •Раздел 1
- •Тема 1.1 Точность вычислительного эксперимента на эвм.
- •§ 1.1.1 Приближенные числа
- •§ 1.1.2. Погрешность приближенных чисел.
- •§ 1.1.3. Погрешность математических операций над приближенными числами.
- •§ 1.1.4 Характеристики методов вычислений.
- •§ 1.2.1 Табличное представление функций. Система обозначений. Операторы.
- •§ 1.2.2. Разностные операторы. Таблицы конечных разностей.
- •§ 1.2.3. Разностные уравнения
- •Раздел 2
- •Тема 2.1. Прикладные задачи интерполяции и аппроксимации и их алгоритмизация. Методы и алгоритмы решения линейных алгебраических уравнений
- •§ 2.1.1. Основные понятия теории приближения функций.
- •Раздел 1
- •Тема 2. Исчисление конечных разностей. Разностные уравнения.
- •§ 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.1. Основные понятия и определения
Нелинейные функции, представляющие математические модели технической системы, обычно характеризуют статические режимы системы (зависимость выходных параметров от входных). Например, коэффициент полезного действия топочного устройства зависит от коэффициента избытка воздуха не линейно.
Как правило, выходной параметр зависит от многих факторов - выходных параметров
y = f(x1, x2, ..., xn)
Прежде всего, при моделировании, интересуются наилучшим (наибольшим, наименьшим) значением выходного параметра. Его поиск называют процедурой параметрической оптимизации. Оптимизировать можно самые разнообразные параметры технической системы (наилучший вариант конструкции, наилучшее распределение ресурсов и т.д.).
Величины x1, x2, ..., xn называют проектными или управляемыми параметрами (их в процессе процедуры можно менять). Число параметров характеризует размерность задачи оптимизации.
Выбор оптимального решения или сравнение двух альтернативных решений проводится с учетом некоторых правил предпочтения на основе установленных критериев, скалярной меры качества.
В основе построения правил предпочтения лежит функция, определяемая проектными параметрами. Такую функцию называют целевой, а также функцией качества или критерием оптимальности.
В зависимости от содержательного смысла проектных параметров и выбранного смысла проектных параметров и выбранного способа их сочетания в целевой функции (ЦФ) качество технической системы будет тем выше, чем больше ее значение (максимизация) или чем
меньше ее значение (минимизация).
Выбор ЦФ носит субъективный характер и поэтому система может быть оптимальна только в смысле данного критерия.
Целевую функцию можно записать в виде:
F(x) = f(x1, x2, ..., xn)
В качестве скалярной величины целевой функции может выступать физический параметр (температура), экономический показатель (прибыль, себестоимость) и др. ЦФ может быть представлена в виде таблицы, алгоритма. Целевых функций может быть несколько и тогда приходится выбирать определяющий критерий или производить их свертку.
Если в технической системе есть n проектных параметров, образующих вектор X = (x1, x2, ..., xn). Целевая функция F(x) определена в некоторой области XO, то максимумом или минимумом F(x) называют ее значение F(X*). При этом для любой точки X Se(X), и достаточно малом >0 (за исключением самой точки X*) выполняются неравенства:
F(x)-F(X*)< 0 или F(x)-F(X*)> 0 соответственно.
Здесь Se(X*) - - окрестность точки X* , которая является множеством точек (векторов), которые находятся от этой точки на расстоянии не превышающих число >0.
Se(X*) = {X | X - X | }, где
|X – X* | - норма вектора X - X* , отождествляемая с расстоянием между точками X и X* .
Точку X* называют точкой локального экстремума. Глобальным экстремумом называют точку, в которой ЦФ имеет наибольшее (наименьшее) значение среди всех экстремумов области определения.
Функцию F(X) называют одно-экстремальной (унимодальной), если она имеет один экстремум и многоэкстремальной, если она имеет более одного максимума (минимума). Если изменить знак целевой функции на обратный, то задача оптимизации станет инверсной.
Если поиск экстремума F(X) производится в неограниченной области изменения вектора X, то экстремум называют безусловным (в противном случае - условным).
Чтобы точка X* была точкой локального безусловного экстремума функции F(X), которая является скалярной от векторного аргумента размерности n и дважды дифференцируема, необходимо, чтобы все её частные производные по элементам Xi вектора X в точке X* были равны 0, т.е.
F(X)/ Xi =0 , X = X* при 1 [1:n]
Это, так называемое, условие стационарности может быть записано в виде
F(X* ) = 0, где
- градиент целевой функции в точке
X*
.
Градиентом функции F(X) в какой-либо точке Xn называют
n-мерный вектор, компоненты которого равны частным производным F(X), вычисленным в точке Xk по всем Xi. Характер экстремума F(X) в стационарной точке определяется
исследованием поведения F(X) в окрестности X* . Для этого формируется матрица вторых производных (матрица Гессе). Если она определена положительно, то в точке X функция F(X) имеет безусловный минимум. Если Ge F(X* ) < 0, то это безусловный максимум. В остальных случаях это седловая точка.
Для n=1 и n=2 задача оптимизации имеет геометрическую интерпретацию. Одномерная оптимизация представляется кривой на плоскости, а двухмернаяповерхность в трёхмерном пространстве. Для двухмерной функции часто используют её представление в виде линий равного уровня, т.е. когда F(X) = C. Эти линии получаются пересечением поверхности F(X) плоскостью, параллельной осям координат X1, X2, что позволяет представлять характерные особенности ЦФ.
Если область определения ЦФ ограничена, то экстремум такой функции называют условным. Различают прямые и функциональные ограничения.
Прямые ограничения имеют вид
Xni < Xi < Xbi при i [1:n],
где Xni, Xbi минимально и максимально допустимые значения i-го проектного параметра. Например, многие проектные параметры не могут быть отрицательными (масса, геометрический размер и т.д.).
Функциональные ограничения представляют собой дополнительные взаимосвязи между проектными параметрами и параметрами, входящими в ЦФ. Они бывают двух видов:
1) тип равенств ψ(X) = 0;
2) тип неравенств ψ(X) > 0.
Прямые и функциональные ограничения формируют допустимую область поиска экстремума.
Постановка задачи оптимизации производится по следующему алгоритму:
1. Определяется вектор проектных параметров X.
2. Формируются ограничения на изменение X.
3. Составляется целевая функция.
Математическая формулировка задачи:
Безусловная оптимизация
min F(X)
X X0
Условная оптимизация
min F(X)
X XA
при ψ (X) = 0; ψ (X) > 0
Xni < Xi < Xbi, i [1:n]
Таким образом, можно различать два типа оптимизации: условная, безусловная и два вида: одномерная, многомерная.
