- •А.Ю. Крюков, Б.Ф. Потапов
- •Крюков, А.Ю.
- •СОДЕРЖАНИЕ
- •ПРЕДИСЛОВИЕ
- •ВВЕДЕНИЕ
- •1.1. СУЩНОСТЬ И НАЗНАЧЕНИЕ МОДЕЛИРОВАНИЯ
- •1.1.1. Понятие о моделировании и свойства моделей
- •1.1.2. Цели моделирования
- •1.1.3. Виды моделирования и классификация моделей
- •1.2.1. Общие положения и основные определения
- •1.2.2. Структура математической модели и ее построение
- •1.2.3. Иерархия математических моделей
- •1.2.5. Классификация математических моделей
- •1.2.6. Геометрическое представление математических моделей
- •1.3. ЗАКЛЮЧИТЕЛЬНЫЕ ЗАМЕЧАНИЯ
- •2.1. МОДЕЛИ МИКРОУРОВНЯ
- •2.2. МОДЕЛИ МАКРОУРОВНЯ
- •2.2.1. Общая характеристика моделей, их структура и сущность
- •1. Компонентные уравнения
- •2. Топологические уравнения
- •Компонентные и топологические уравнения электрической системы
- •2.2.5. Задания для самостоятельной работы
- •2.3. ЗАКЛЮЧИТЕЛЬНЫЕ ЗАМЕЧАНИЯ
- •3.2.1. Общие принципы моделирования полей
- •3.2.2. Особенности построения моделей
- •3.2.3. Модели стационарных полей
- •3.2.4. Модели нестационарных полей
- •3.2.7. Задание для самостоятельной работы
- •3.3. ЗАКЛЮЧИТЕЛЬНЫЕ ЗАМЕЧАНИЯ
- •4.1.1. Надежность объектов как комплексное свойство
- •4.1.2. Классификация отказов и временные понятия
- •4.3.1. Основные понятия, определения и положения
- •4.3.2. Основные характеристики случайных величин
- •4.4. ПОКАЗАТЕЛИ НАДЕЖНОСТИ
- •4.4.1. Показатели безотказности
- •4.6.1. Общая характеристика и виды моделей
- •4.6.2. Распределения, используемые в теории надежности
- •4.6.3. Композиция законов распределения
- •4.6.4. Задание для самостоятельной работы
- •4.7 ПОТОКИ ОТКАЗОВ И ВОССТАНОВЛЕНИЙ
- •4.8. РАСЧЕТ НАДЕЖНОСТИ ТЕХНИЧЕСКИХ СИСТЕМ
- •4.8.1. Общие положения
- •4.8.2. Модели параметрических отказов
- •5.1. ОСНОВЫ ТЕОРИИ ГРАФОВ И СЕТЕВЫЕ МОДЕЛИ
- •5.1.1. Общие замечания, основные понятия и определения
- •5.1.2. Решение задачи о максимальном потоке
- •5.1.3. Потоки минимальной стоимости
- •5.1.4. Практические примеры сетевых задач
- •5.1.5. Некоторые обобщения по сетевым задачам
- •5.1.6. Альтернативные методы решения сетевых задач
- •5.2.1. Основные положения
- •5.2.2. Структура принятия решения
- •5.2.4. Пример применения классических критериев
- •6.1.1. Постановка задач оптимизации и критерии
- •оптимальности
- •6.1.2. Многокритериальные задачи оптимизации
- •6.1.3. Классификация методов оптимизации
- •6.2. РЕГУЛЯРНЫЕ МЕТОДЫ ОПТИМИЗАЦИИ
- •6.2.1. Методы математического анализа
- •6.2.2. Понятие о вариационном исчислении
- •6.2.3. Принцип максимума Понтрягина
- •6.2.4. Метод динамического программирования
- •6.3. ПРЯМЫЕ МЕТОДЫ ОПТИМИЗАЦИИ
- •6.3.2. Минимаксные стратегии одномерного поиска
- •6.3.4. Многомерные методы безусловной оптимизации
- •6.3.7. Заключение по прямым методам оптимизации
- •7 1. ОБЩИЕ ПОЛОЖЕНИЯ
- •ЗАКЛЮЧЕНИЕ
- •Учебное издание
- •Учебное пособие
6.2. РЕГУЛЯРНЫЕ МЕТОДЫ ОПТИМИЗАЦИИ
6.2.1. Методы математического анализа
Эти методы, по существу, представляют поиск оптимального значения параметров целевой функции классическими спо собами исследования ее первой и второй производных, с кото рыми студенты подробно знакомятся в курсе математики. По этому для закрепления знаний о данных методах и умений их использовать авторы рекомендуют обратиться к учебникам по математическому анализу.
6.2.2. Понятие о вариационном исчислении
Известно, какое большое значение для развития классиче ского анализа функций одной или нескольких переменных име ли задачи на отыскание экстремумов дифференцируемых функ ций. Однако еще с давних времен известны задачи на экстре мум, которые нельзя отнести к задачам конечномерного матема тического анализа. Древнейшей такой задачей считается клас сическая изопериметрическая задача: среди плоских гладких замкнутых кривых заданной длины найти кривую, ограничи вающую наибольшую площадь.
В качестве другой задачи на экстремум, не относящейся ко классическому математическому анализу, следует назвать задачу, которую в 1696 г. поставил известный математик Иоганн Бернулли —задачу о брахистохроне ~ линии быстрейшего ската. В этой задаче требуется найти гладкую линию, соединяющую две заданные точки А и В вертикальной плоскости, не лежащие на вертикальной прямой, по которой под действием силы тяже сти материальная точка скатится из А в В за кратчайшее время (рис. 6.5, а). При этом время, необходимое для перехода из точ-
239
ки А в точку В на основании закона сохранения энергии, опре деляется по формуле
t - p/l + (у ')2dx, |
(6.15) |
а ^ |
|
где у - неизвестная функция, g - ускорение свободного падения.
Рис. 6.5. Иллюстрации к задачам: а — Бернулли; б —Ферма
Также к классическим задачам вариационного исчисле ния относится задача о преломлении света. Согласно принципу Ферма, луч света, выходящий из точки А и попадающий в точ ку В, избирает путь, время перехода по которому является наименьшим. В однородной среде скорость света постоянна, и свет распространяется по прямым. Если же среда неоднород на, то скорость света изменяется от точки к точке, а траектории лучей света уже не будут прямыми. Пусть средой является ат мосфера. Поскольку плотность воздуха зависит от высоты
у над уровнем моря, то правомерно предположить, что и ско рость света V зависит от у и выражается с помощью известной функции V(y). Определим траекторию луча света из точки А в точку В. В вертикальной плоскости, проходящей через точ ки А и В (рис. 6.5, б), выберем прямоугольную систему коор динат так, что ось Ох горизонтальна и расположена на уровне моря. Нам известны координаты А (а, уа) и В (b, уь). Считаем, что луч света распространяется по кривой, являющейся графи-
ком гладкой функции у(*), определенной на отрезке [а, Ь]. При этом время, необходимое для перехода света из точки А в точку В, выражается интегралом:
W i+ ( / ) 2-dx. |
(6.16) |
! V(y) |
|
Задача состоит в определении такой гладкой функции У= У(х)> Удовлетворяющей условиям
У(р )=Уа АЬ)=Ув, |
(6.17) |
что интеграл (6.15) или (6.16) получает наименьшее значение.
Сравнив (6.15) и (6.16), отметим, что задача о брахисто хроне - частный случай задачи о преломлении света. Этот факт, подмеченный впервые И. Бернулли, представляет собой так на зываемую оптико-механическую аналогию.
Оригинальность сформулированных задач - в том, что неизвестными в них являются функции, которые должны сде лать значение интеграла наименьшим.
Таким образом, вариационное исчисление - это раздел ма тематики, посвященный нахождению наибольших и наимень ших значений переменных величин, зависящих от выбора одной
или нескольких функций.
Для дальнейшего изложения введем некоторые определе ния. Функционалом называется переменная величина, зависящая от выбора одной или нескольких функций. В вариационном ис числении важнейшими являются функционалы, заданные с по
мощью интегралов, например
ь |
|
I(W(x)) = JF(x,W,W)dx, |
(6.18) |
dW(x)
где W =
dx
Подынтегральная функция F(x,W,W) называется интегрантом функционала, предполагается непрерывной и имеющей непрерывные частные производные по всем переменным до
второго порядка включительно.
К функционалам сводятся описания многочисленных бал листических и транспортных задач, задач, связанных с распреде лением ресурсов и капиталовложений, с заменой оборудования и т.д. Как представлено выше, многие законы механики и физики могут быть сформулированы в виде требования максимизации или минимизации некоторого функционала. Эти законы носят название вариационных принципов механики или физики.
Методы классического вариационного исчисления при годны для оптимизации функционалов, определенных на классе гладких функций, у которых в рассматриваемой области непре рывна первая производная, или на классе кусочно-гладких функций, у которых первая производная имеет конечное число разрывов первого рода. Основное соотношение вариационного исчисления - знаменитое уравнение Эйлера - выводится из ана
лиза изменений вариаций функционала I(W(x)), играющих роль производных функции Щх) [3].
Предполагается, что оптимальное (минимальное) значение
функционала достигается на кривой W(x) (рис. 6 .6), которая, та ким образом, является оптимальной среди всех близких функ ций fV(x):
/.in = W ( x ) ) = JF (x,W ,ti)dx . |
(6.19) |
a
Выполнив несложные преобразования, которые представ лены, например, в [3] и [24], можно получить выражение
F# |
= 0, при V хе[а, Ь]. |
(6.20) |
Соотношение (6 .20) и пред ставляет собой уравнение Эйлера, выражающее необходимое условие экстремума и в той или иной форме лежащее в основе всех задач вариа ционного исчисления. Общее реше ние уравнения Эйлера содержит две
Рис. 6.6. К определению
постоянные, для определения кото
экстремалей
рых, как правило, задаются значения функционала функционала в начале и конце ис
следуемого интервала W(a), W(b). Для уравнения Эйлера можно показать [3], что в случае мини мума должны выполняться условия FM, > 0, а в случае макси мума, наоборот, F** < 0 .
В рамках классического вариационного исчисления могут решаться и более сложные задачи. Так, например, для функцио нала вида
I(W(x))= \F{x,W,,...,Wn,Wlv..,Wn)dx |
(6.21) |
функции W,(х), реализующие экстремум, должны удовлетворять системе дифференциальных уравнений:
F - ----- F: = 0, i = 1,2, dx w'
Для функционала вида
ь
I(W(x))= \F(x,W ,W ,W ,...,W (n))dx,
а
(6.22)
(6.23)
где функция F считается дифференцируемой и + 2 раз по всем аргументам, функция Щх), реализующая экстремум, должна быть решением уравнения
F |
—— F- +-^— F- +... + (-1)"-^— Fj, = 0. (6.24) |
|
* |
dx * dx2 * |
* dx" |