- •Численные методы решения задач строительства
- •Часть 1
- •Предисловие
- •В ведение
- •Общие сведения о вычислительном эксперименте и математическом моделировании
- •Численные методы.
- •Погрешности вычислений
- •Понятия точности, устойчивости и сходимости при численном решении
- •Глава 1 Основные понятия матричного исчисления
- •1.1. Матрицы и векторы.
- •1.2. Матрицы специального вида
- •1.3. Действия над матрицами
- •1.4. Нормы матрицы и вектора
- •1.5. Функции ms Excel для операций над матрицами
- •Контрольные вопросы
- •Г лава 2. Численные методы решения нелинейных уравнений
- •2.1. Первый этап. Отделение корней
- •Второй этап. Этап уточнения корня
- •Iteration (итерация)- повторение, ре-зультат повторного применения какой-либо математической операции.
- •Метод половинного деления (бисекции)
- •Метод хорд
- •Метод Ньютона (метод касательных)
- •Модифицированный метод Ньютона
- •Реализация численных методов решения нелинейных уравнений средствами приложения ms Excel
- •Последовательность действий:
- •2.3.1. Решение нелинейных уравнений с использованием надстройки «Подбор параметра»
- •Последовательность действий
- •Контрольные вопросы
- •3.2. Прямые методы решения систем линейных алгебраических уравнений
- •3.2.1. Метод Гаусса
- •3.2.2. Метод прогонки
- •Алгоритм метода прогонки
- •3.3. Итерационные методы решения систем линейных алгебраических уравнений
- •3.3.1. Метод Якоби (простых итераций)
- •3.3.2. Метод Гаусса – Зейделя.
- •3.3.3. Условия сходимости итерационного процесса
- •3.5. Устойчивость решения слау относительно исходных данных (или обусловленность задач и вычислений)
- •3.6. Примеры решения слау с использованием электронных таблиц ms Excel
- •3.6.1. Реализация метода Гаусса
- •Последовательность действий
- •Прямой ход метода Гаусса.
- •3.6.2. Решение слау с помощью надстройки «Поиск решения»
- •Последовательность действий:
- •3.6.3. Реализация метода Якоби средствами приложения ms Excel
- •Последовательность действий
- •3.6.4. Реализация метода прогонки средствами приложения Excel
- •Последовательность действий
- •Контрольные вопросы
- •Г лава 4. Численное интегрирование
- •Алгоритм метода половинного шага.
- •4.1. Квадратурные формулы прямоугольников
- •4.2. Квадратурная формула трапеций
- •4.3. Квадратурная формула Симпсона
- •4.4. Реализация методов численного интегрирования средствами приложения Ms Excel
- •Последовательность действий:
- •Контрольные вопросы
- •Г лава 5. Аппроксимация
- •5.1. Задачи аппроксимации
- •5.2. Интерполирование функций
- •5.2.1. Постановка задачи интерполирования
- •5.2.2. Интерполяционная формула Лагранжа
- •5.3. Среднеквадратичное приближение функций
- •5.3.1. Постановка задачи
- •5.3.2. Метод наименьших квадратов
- •5.3.3. Линейная эмпирическая формула (линейная регрессия)
- •5.3.4. Коэффициент корреляции
- •5.3.5. Квадратичное (параболическое) приближение
- •5.3.6. Эмпирические формулы с двумя параметрами (метод выравнивания)
- •5.4. Решение задач аппроксимации с помощью электронных таблиц Excel
- •5.4.1. Построение уравнений регрессии методом наименьших квадратов с использованием надстройки «Поиск решения»
- •Последовательность действий
- •5.4.2. Построение линейной эмпирической формулы с использованием встроенных функций линейн и тенденция
- •Последовательность действий
- •Последовательность действий
- •Контрольные вопросы
- •Глава 6. Численные методы оптимизации
- •6.1. Общие сведения
- •6.1.1. Математическая модель задачи оптимизации
- •6.1.2. Классификация задач математического программирования
- •6.2. Постановка задачи оптимального проектирования
- •6.3. Задачи линейного программирования
- •6.3.1. Общая постановка задачи
- •6.3.2. Геометрический смысл системы линейных неравенств
- •Случай n проектных параметров.
- •6.3.3. Геометрический метод решения задач линейного программирования
- •Последовательность действий:
- •6.3.4. Симплекс-метод решения задач линейного программирования
- •Примеры задач линейного программирования в сфере проектирования и управления строительным производством
- •6.4.1. Задача об оптимальном плане выпуска продукции
- •6.4.2. Задача об оптимальном раскрое материалов (о минимизации отходов)
- •6.4.3. Задача о планировании смен на предприятии
- •6.4.4. Задача о покрытии местности при строительстве объектов
- •6.4.5. Транспортная задача
- •Задача о назначениях (проблема выбора)
- •6.5. Решение задач оптимизации с помощью ms Excel
- •6.5.1. Решение задачи планирования производства
- •Последовательность действий:
- •6.5.2. Решение транспортной задачи
- •Последовательность действий:
- •Контрольные вопросы
- •Литература
- •Глава 1 16
- •Глава 2. 27
- •Глава 3. 51
- •Глава 4. 86
- •Глава 5. 100
- •Глава 6. 125
Задача о назначениях (проблема выбора)
Постановка задачи. В распоряжении имеется n механизмов (бригад) и n различных видов работ. Производительность каждого механизма на различных работах, вообще говоря, различна. Обозначим через cij (i,j=1,2,..,n) производительность i-го механизма на j-й работе. Матрица C=[cij] называется матрицей эффективностей или производительностей. Требуется так распределить n механизмов на n работах, чтобы каждый механизм выполнял одну и только одну работу и чтобы при заданной производительности каждого механизма на каждой из работ суммарный эффект был бы максимальным. |
Эта задача является частным случаем транспортной задачи
Обозначим через xij переменную, равную единице, если i-й механизм назначен на j-ю работу, и равен нулю, если он на эту работу не назначен. В качестве целевой функции выбираем суммарную производительность всех механизмов
Ограничения записываем из условия, что каждый механизм выполняет одну работу и что каждая работа обеспечивается одним механизмом:
Пример 6.9. Задача о назначениях.
Имеются три бригады В1, В2, В3, каждая из которых может быть использована на каждой из 3-х строительной площадке с производительностью (в условных единицах), заданной в таблице.
Таблица.
-
Строител. площадка
Производительность бригад (у.е.)
В1
В2
В3
1
1
2
3
2
2
4
1
3
3
1
5
Требуется распределить бригады по одной на каждую строительную площадку так, чтобы суммарная производительность всех бригад была максимальной.
Обозначим xij=1, если на i-ю стройплощадку назначена j-я бригада, и xij=0, если она на эту стройплощадку не назначена.
В качестве функции цели возьмем суммарную производительность всех бригад:
Zmax =x11 +2x12 +3х13+2x21 +4x22 +x23+3x31 +x32.+5x33.
Систему ограничений запишите из условий, что каждая бригада выполняет работу на одной определенной стройплощадке и что каждая площадка обеспечена одной бригадой рабочих.
6.5. Решение задач оптимизации с помощью ms Excel
Решение задач оптимизации в Microsoft Excel производится с помощью надстройки «Поиск решения».
При поиске оптимального решения задачи можно изменять условия и варианты поиска решения, используя диалоговое окно «Параметры поиска решения». Пользователь при этом может сам выбрать алгоритм оптимизации.
При этом принципы построения схемы задачи в таблице Excel и ее решения практически совпадают как для линейного, так и для нелинейного программирования.
Рассмотрим процесс решения задач оптимизации на прмерах двух задач: задачи планирования производства, и транспортной задачи, используя возможности таблиц Excel.
6.5.1. Решение задачи планирования производства
Рассмотрим задачу, приведенную в примере 6.2.
Математическая модель этой задачи, как описано выше имеет вид:
при ограничениях
Последовательность действий:
Подготовим таблицу как показано на рис.6.15. Ячейки, содержащие целевую функцию (А10) и изменяемые ячейки (проектные параметры x1, x2) В10:С10 тонируем. Для контроля счета в ячейки В10:С10 введем единицы. Значения проектных параметров х1=1 и х2=1 можно рассматривать как начальное приближение решения задачи.
Систему ограничений (6.12) запишем в виде, как показано на рис.6.15, т.е.:
ai x1 + bi x2 - di ≥ 0 , i=1, 2, 3, 4, 5.
Обратите внимание на то, что знаки неравенств одинаковы для всех строк системы ограничений. В рациональности этого вы убедитесь ниже.
В ячейки В3:D7 введем коэффициенты системы ограничений.
Рис.6.15. Расчетная схема для примера 6.2
В ячейку E3 введем формулу для вычисления левой части 1-го ограничения, т.е.
E3=СУММПРОИЗВ($B$10:$C$10;B3:C3)+D3,
и после ввода скопируем ее вниз до конца таблицы. Будет не лишним проверить результаты счета для заданных значений х1=1 и х2=1.
Формулу для целевой функции запишем в ячейке A10:
A10=1-В10-С10.
Выберем команду меню Данные\Поиск решения и в появившемся окне сделаем соответствующие установки (рис.6.16). Ограничения устанавливаются с помощью кнопки Добавить, которая вызывает окно Добавление ограничения для ввода этих ограничений (см. рис. 6.16).
Рис.6.16. Окно «Поиск решения»
После щелчка на кнопке Выполнить, при условии, что введено все верно, результат решения будет иметь вид, как показано на рис.6.17.
Рис.6.17. Результаты решения примера 6.2
Результаты можно сохранить или отказаться от них (восстановив исходные значения). Можно получить один из видов отчетов (Результаты, Устойчивость, Пределы). Отчет можно оформить на отдельном листе Книги с соответствующим именем.
