- •А.Ф. Хрусталев основы математического моделирования
- •Хрусталев Александр Федорович
- •Предисловие
- •Введение
- •§1. Математическое моделирование как метод познания
- •1.1. Понятие математической модели
- •1.2. О гуманитарной значимости математических моделей
- •1.3. О сущности аксиоматического метода
- •1.4. Пример построения и исследования математической
- •1.5. А были ли американские астронавты на Луне?
- •§2. Общая схема применения математики. Множественность и единство моделей
- •2.1. Общая схема применения математики
- •Множественность и единство моделей
- •§3. Требования, предъявляемые к математической модели
- •3.1. Об адекватности математической модели
- •О полноте математической модели
- •3.3. О непротиворечивости модели
- •3.4. О других требованиях
- •§4. Приближенные числа и действия с ними
- •4.1. Десятичная запись приближенных чисел
- •4.2. О приближенных числах и значащих цифрах
- •4.3. О правилах действий с приближенными числами
- •4.4. Об оценке точности решения
- •§5. Система линейных уравнений как математическая модель
- •5.1. О решении переопределенных систем
- •5.2. Об устойчивости решений относительно погрешностей
- •§6. О построении математических моделей, описывающих результаты экспериментов
- •6.1. Математическая модель прямых измерений
- •6.2. О построении математических моделей косвенных измерений
- •§7. Оптимизационные модели
- •7.1. Примеры простейших задач на экстремум
- •7.2. О типах задач математического программирования
- •7.3. Конкретная транспортная задача
- •8.2. Простейшие задачи
- •8.3. Основная теорема о функциональной зависимости размерных величин
- •8.4. О необходимости соблюдения правила размерностей
- •9.1.О геометрическом подобии
- •На основании подобия треугольников можно решать различные задачи, связанные с определением расстояний, измерить которые непосредственно либо сложно, либо невозможно.
- •9.2. О физическом подобии
- •9.3. Измерение силы тяжести на луне с помощью телевизора
- •§10. Конечные уравнения как математические модели
- •10.1. Химическая задача
- •10.2. Физическая задача
- •10.3. Математическая задача №1151 [20]
- •§ 11. Предел и производная как математические модели
- •11.1. Предел как математическая модель
- •Определение длины окружности
- •Определение числа е
- •3. Определение площади поверхности шара
- •4. Определение вероятности события
- •5. Определение скорости свободного падения тела
- •6. Определение суммы ряда
- •11.2. Производная как математическая модель
- •§ 12. Интеграл как математическая модель
- •12.1. Задачи, приводящие к понятию интеграла
- •Задача об определении площади криволинейной трапеции
- •12.2. Две схемы приложения определённого интеграла
- •§13. О математических моделях, описываемых дифференциальными уравнениями
- •13.1. Задача о росте населения
- •13.2. Задача о колебаниях пружинного маятника
- •13.3. Задача о форме зеркала прожектора
- •§14. Вычислительный эксперимент
- •14.1. Вычислительный эксперимент как метод познания
- •14.2. Приложения метода к решению химических задач
- •14.3. Приложения метода к решению математических задач
- •Н. Винер § 15. Вероятность события как математическая модель
- •15.1. О понятии «вероятность события»
- •15.2. Вероятностные модели конкретных задач
- •Случай всегда приходит на помощь
- •§ 16. Математическая модель опыта. Случайная величина как математическая модель
- •16.1. Математическая модель опыта
- •16.2. Случайная величина (св) как математическая модель
- •16.3. Механическая модель св
- •§ 17. Схема бернулли. Цепь а.А. Маркова
- •17.1. Схема Бернулли
- •17.2. Цепь событий как математическая модель
- •17.3. Простейшая цепь а.А. Маркова
- •§ 18. Понятие о методе монте – карло
- •18.1. Общее представление о методе
- •Вероятностная задача:
- •18.2. Применения метода к решению невероятностных задач
- •18.3. О применении метода к проверке статистических гипотез
- •18.4. О сущности метода
- •§ 19. Простейшие теоремы математической химии и их приложения к построению математических моделей
- •19.1. Взаимно простые числа в химии
- •19.2. Критерии алканов
- •19.3. Теоремы гомологии
- •19.4. Теорема о молекулярной формуле соединения
- •§20. О математических моделях географии
- •20.1. О математических моделях Земли
- •20.2. Об измерениях диаметра Земли
- •20.3. О расширении горизонта в зависимости от высоты
- •20.4. О математических моделях холмов и впадин Земли
- •§21. О математических моделях логических связок
- •21.1. Высказывания
- •Отрицание
- •Конъюнкция – математическая модель логической
- •Дизъюнкция – математическая модель логической связки «или»
- •21.5. Импликация – математическая модель логической связки «если…, то»
- •21.6. Об определении модуля
- •21.7. О задании функции разными формулами
- •21.8. Высказывания и контактные схемы
- •21.9. Алгебра правды и лжи
- •§ 22. Понятие об имитационном моделировании
- •§ 23. Применение эвм к анализу математических моделей
- •23.1. Применение к совершенствованию учебных задач и их
- •23.2. Применение к анализу устойчивости решений
- •§ 24. Упрощение моделей
- •24.1. Химические задачи
- •24.2. Математические задачи
- •Другие задачи
- •§ 25. Вариант расчетно - графической работы
- •Заключение
- •Библиографический список
- •Содержание
- •Основи математичного моделювання
- •Основы математического моделирования
§7. Оптимизационные модели
В мире не происходит ничего,
в чем не был бы виден смысл
какого-либо максимума или минимума.
Л. Эйлер
7.1. Примеры простейших задач на экстремум
История человеческого общества, так или иначе, связана с поиском наилучших, оптимальных решений, принимаемых в различных, реальных ситуациях.
При рассмотрении такого типа практических задач, как и всегда математика применяется не непосредственно к реальным объектам, а к их математическим моделям, которые в этих условиях называются оптимизационными, поскольку предназначены для нахождения оптимальных решений.
Задача.
Требуется
изготовить паровой котел цилиндрической
формы заданного объема
.
Желательно,
чтобы полная поверхность котла была
наименьшей,
так как при этом условии на его изготовление
пойдет минимальное
количество металла; кроме того, чем
меньше поверхность
котла, тем меньше он будет охлаждаться
от соприкосновения
с окружающей средой. Найти наилучшие
размеры
котла.
Обозначим
радиус основания цилиндра через
,
а его высоту через
.
Тогда получим систему
(7.1)
Таким
образом, мы пришли к следующей
математической модели: среди всех
решений смешанной системы (7.1) найти
такое, при котором полная поверхность
цилиндра (функция цели)
принимает наименьшее значение. Подчеркнем,
что в математической модели для численных
значений величин обычно сохраняют
обозначения самих величин, как в
рассматриваемой задаче. При этом
предполагается, что все величины выражены
через одни и те же основные единицы: в
нашей задаче – это единицы длины,
например метр (м). Тогда все величины в
созданной модели будут безразмерными.
Из
(7.1) имеем
, при этом
(7.2)
и, следовательно,
.
Приравнивая производную к нулю, из уравнения
найдем
,
поэтому
.
При
найденном
функция
(7.2) принимает наименьшее значение,
поскольку как при очень малых, так и при
очень больших значениях
площадь поверхности
будет большой. А так как
только при одном значении
,
то ему и соответствует наименьшая
площадь полной поверхности цилиндра.
Химическая задача. Смесь водорода и паров йода находится в замкнутом сосуде, помещенном в термостат, причем созданы условия, при которых реакция между этими веществами протекает практически необратимо. Определите начальную, объемную долю водорода (в процентах), при которой скорость реакции будет максимальной.
Сначала определим концентрацию водорода, при которой скорость реакции
(7.3)
максимальна
в начальный момент времени. Для этого
в соответствии с законом действующих
масс: «При постоянной температуре
скорость химической реакции прямо
пропорциональна произведению концентраций
реагирующих веществ» необходимо
проанализировать функцию
,
где
– константа скорости реакции, зависящая
только от температуры (и, значит,
неизменная по условию задачи), а
и
- концентрации
и
- соответственно, выраженные в объемных
процентах.
Поскольку х+у=100 (%), то у=100 – х и, следовательно,
v = kx(100 – х). (7.4)
Найдем
значение х,
при котором v
имеет максимум в начальный момент
времени. Для этого преобразуем функцию
(7.4) следующим образом:
.
Очевидно, v достигает максимума, когда выражение
,
которое
не может быть положительным, обращается
в нуль, поскольку слагаемое 502k,
не зависит от
.
Таким образом, скорость реакции
максимальна при х=50
%,
у=50
% (стехиометрическое соотношение
реагентов
).
Определим
теперь отношение
,
соответствующее максимуму скорости в
начальный момент, при условии, что в
газовой смеси помимо водорода и йода
содержатся компоненты, инертные по
отношению к реакции (7.3). Обозначив
концентрацию последних через
и учитывая, что
(%),
получаем:
;
(7.5)
.
(7.6)
Приведя,
как и выше, выражение (7.6) к виду
,
приходим к выводу, что v
достигает максимума при фиксированном
,
когда
;
но тогда на основании (7.5),
,
тo
есть стехиометрическое соотношение
сохраняется
вне зависимости от того, присутствуют
ли в смеси инертные компоненты или нет.
Но если при
содержание водорода и йода равно 50 %, то
при
оно соответственно меньше.
Наконец,
рассмотрим реакцию (7.3) в любой момент
времени. Поскольку в условиях необратимости
реакции ее продукт
представляет собой инертный компонент,
то и тогда максимум скорости процесса
будет наблюдаться при стехиометрическом
соотношении компонентов. Заметим, что
полученный результат не зависит от
константы скорости реакции
,
и поэтому вывод справедлив для любой
температуры, поддерживаемой термостатом.
Подчеркнем, что знание условий, при которых реакция протекает с максимальной скоростью, позволяет осуществлять оптимальные управления химическими реакциями на производстве, добиваясь наибольшей его производительности.
Теперь
рассмотрим задачу об оптимальной
температуре процесса производства
некоторого вещества химическим заводом
[37, с.112]. Выход продукта
этого вещества определяется температурой:
,
причём
вид
функции
заранее не известен, он зависит от
используемого сырья.
Температуру
можно варьировать в определенных
пределах:
.
Получив очередную партию сырья, нужно
найти температуру
,
при
которой наиболее выгодно вести
производство, т. е.
функция
достигает
своего наибольшего значения.
Однако в данном случае нет никакой формулы для целевой функции . Поэтому, чтобы определить ее значение при некоторой температуре , нужно провести опыт либо в лаборатории, либо прямо в производственных условиях. Совершенно ясно, что возможно лишь конечное число измерений и каждое измерение требует времени, а задерживать производство нельзя. Значит, необходимо получить ответ на поставленный вопрос после небольшого числа измерений, т. е. по значениям функции в нескольких точках. Сравнивая эти значения, мы найдем, конечно, приближенно, температуру, при которой выход продукта будет максимальным.
