- •А.Ю. Крюков, Б.Ф. Потапов
- •Крюков, А.Ю.
- •СОДЕРЖАНИЕ
- •ПРЕДИСЛОВИЕ
- •ВВЕДЕНИЕ
- •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. ОБЩИЕ ПОЛОЖЕНИЯ
- •ЗАКЛЮЧЕНИЕ
- •Учебное издание
- •Учебное пособие
для решения которых целесообразно использовать рассмотрен ные выше методы. С этой целью вводят штрафную функцию и образуют новую составную целевую функцию
M(x,) = W(xi) + Y , 4{Gl(xi)\, |
(6.133) |
/=1 |
|
где W(Xj) —рассматриваемая ограниченная функция, М(х,) - со ставная неограниченная функция, a (p[G7{x,)] - штрафная функ ция, учитывающая ограничения, заданные неравенствами.
Как видим, новая целевая функция образуется сложением целевой функции рассматриваемой задачи и штрафной функ ции. Штрафная функция ср(д:,) равна нулю во всех точках про странства проектирования, удовлетворяющих условиям G,{x,) > 0, и стремится к бесконечности в тех точках, в которых эти условия не удовлетворяются. Если все условия удовлетво ряются, то функции М(х,) и W(x,) имеют, очевидно, один и тот же минимум. Если же хотя бы одно условие не удовлетворяется, то целевая функция приобретает бесконечно большие значения, весьма далекие от минимального значения W(xi). Тем самым на каждую конструкцию, не удовлетворяющую поставленным ус ловиям, «налагается штраф».
Существуют также различные другие способы ввода и ис пользования штрафных функций, с которыми читателю предлагает ся ознакомиться самостоятельно по материалам источников [20,26].
6.3.7. Заключение по прямым методам оптимизации
Таким образом, процессы многомерной оптимизации представляют собой целенаправленное движение в пространстве управляемых параметров к точке, в которой достигается экстре мум целевой функции.
Существуют и другие методы поиска, не рассмотренные в данном разделе, в частности, метод сопряженных градиентов и метод переменной метрики. Читателю рекомендуется рассмот реть их самостоятельно по материалам источников [14,20, 24,26].
7 1. ОБЩИЕ ПОЛОЖЕНИЯ
При разработке математических моделей часто возникает задача включения в модель экспериментальных данных, полу ченных при испытаниях технических объектов или их отдельных узлов. Результаты таких испытаний всегда представляют собой наборы величин, характеризующих работу объекта при различ ных сочетаниях режимных параметров. Результаты численных исследований отдельных элементов и узлов также представляют собой наборы величин и с точки зрения дальнейшего использова ния не отличаются от экспериментальных данных. Наиболее эф фективным средством представления результатов экспериментов и численных расчетов в системах математического моделирова ния являются модели идентификации, которые, таким образом, могут быть построены на основе как чисто эмпирических данных, так и результатов теоретических исследований.
При построении модели идентификации обычно предпо лагается, что физическая теория работы объекта отсутствует или по гем или иным причинам не может быть использована. Объ ект идентификации представляет собой так называемый «чер ный ящик» с некоторым числом регулируемых (или, по крайней мере, измеряемых) входов х и одним или несколькими наблю даемыми (измеряемыми) выходами W. Задачей идентификации является построение модели объекга по результатам наблюде ний его реакции на возмущения внешней среды.
Математически задача идентификации формулируется следующим образом. Имеется п экспериментов (г = 1,2, ..., «),
Построение модели идентификации начинается с выбора формы -модели, т.е. вида зависимости W При этом на практике могут встретиться два случая.
1. Форма математической модели известна заранее, и за дача идентификации сводится к определению коэффициентов этой модели. Так, описание ряда затухающих или развивающих ся процессов дается зависимостями экспоненциального типа:
W = a - epx; W = a e‘pi; W = |
(7.3) |
изадачей исследования является определение коэффициентов a
ир. При описании ряда физических процессов известна струк тура математических зависимостей. Так, например, при изуче нии термоэмиссионных преобразователей энергии для опреде ления термоэмиссионного тока используется формула Ричард- сона-Дешмана:
-еср |
|
J = A T 2 e kT |
(7.4) |
где J - плотность термоэмиссионного тока; Т - абсолютная тем пература; е - заряд электрона; к - постоянная Больцмана; ср - работа выхода металла; А - постоянная Ричардсона.
Представленная зависимость включает в себя одну вели чину, определяемую экспериментально, - работу выхода ф. За дачей идентификации в данном случае является подбор значе ний этой величины на основе экспериментальных данных.
Такая же ситуация имеет место при изучении теплообмена, течений жидкости, переходных процессов в системах автомати ческого регулирования. Иногда модель может быть сформулиро вана в виде критериальной зависимости, и задача идентификации тогда заключается в определении степенных показателей крите риев. В этих случаях количество подлежащих определению ко эффициентов модели идентификации точно известно.
где величина р и общее число коэффициентов регрессии т оп ределяются размерностью задачи к и степенью полинома d:
m = c txsi = ^ + |
,р = т - 1, при £ = \ m = l + d n p |
= d. (7.10) |
|
kldl |
|
|
|
Использование конечного отрезка ряда вместо бесконеч |
|||
ного разложения |
и определение |
коэффициентов |
регрессии |
на основе ограниченного набора |
экспериментальных данных |
||
не позволяют точно восстановить зависимости (7.7) и (7.9). По |
|||
этому при идентификации, как правило, ставится задача нахож |
|||
дения приближенной модели. |
|
|
При этом возможны три соотношения между количеством экспериментов п и размерностью массива коэффициентов рег рессии т.
1. Количество коэффициентов регрессии равно количест ву экспериментальных точек (т = п). В этом случае задача сво дится к системе т уравнений с т неизвестными. В данном слу чае можно найти поверхность (кривую), точно проходящую че рез все экспериментальные точки.
В качестве примера на рис. 7.2, а показана такая кривая, построенная в виде полинома 3-й степени (k = 1. т = 4) по четы рем экспериментальным точкам. Для опорных точек полученное решение является точным, т. е. обладает нулевой ошибкой. Од нако точность описания процесса в промежуточных точках и достоверность самого характера модели при этом достаточно сомнительны. В частности, для экспериментальных данных, по казанных на рис. 7.2, а, более точной является совсем другая модель (кривая отклика показана пунктиром). Дело в том, что в случае п = т для определения коэффициентов регрессии ис пользуются все экспериментальные точки и не остается инфор мации - степеней свободы - для проверки точности модели. По этому случай п = т с точки зрения идентификации не является оптимальным.