Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Прикладная математика. Решение задач с применением табличного процессора Excel. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
2 Мб
Скачать
☆
КАЗАНСКИЙ ГОСУДАРСТВЕННЫЙ
АРХИТЕКТУРНО-СТРОИТЕЛЬНЫЙ УНИВЕРСИТЕТ
Ф.Г. Ахмадиев, Р.Ф. Гиззятов
ПРИКЛАДНАЯ МАТЕМАТИКА.
РЕШЕНИЕ ЗАДАЧ С ПРИМЕНЕНИЕМ
ТАБЛИЧНОГО ПРОЦЕССОРА EXCEL
Ай Пи Ар Медиа
Москва
2022
УДК 51 ББК 22.1
А95
Рецензенты:
Шабалин П.Л. — д-р физ.-мат. наук,
проф. кафедры высшей математики КГАСУ;
Шакирова Л.Р. — д-р пед. наук, проф., зав. кафедрой теории
технологии преподавания математики и информатики КФУ
Ахмадиев, Фаил Габдулбарович.
А95 Прикладная математика. Решение задач с применением табличного процессора Excel :
учебное пособие / Ф.Г. Ахмадиев, Р.Ф. Гиззятов ; Казанский государственный архитек­турно-строительный университет. — Москва : Ай Пи Ар Медиа, 2022. — 135 с. — Текст : электронный.
ISBN 978-5-4497-1392-6
В учебном пособии рассмотрены примеры решения задач прикладной математики, в которых применение излагаемых методов иллюстрируется на конкретном практическом материале и доводится до числового результата. Представлены основные теоретиче­ские положения и математические модели, примеры решения типовых задач с выбором численного метода решения и составлением программы на языке VBA или электронной таблицы в Excel с использованием ее надстройки «Поиск решения».
Предназначено для студентов, обучающихся по укрупненным группам направлений подготовки и специальностей 08.00.00 «Техника и технологии строительства», 23.00.00 «Техника и технологии наземного транспорта», изучающих дисциплины «Прикладная математика», «Математическое моделирование».
Учебное электронное издание
ISBN 978-5-4497-1392-6 © Ахмадиев Ф.Г., Гиззятов Р.Ф., 2016
© Казанский государственный архитектурно-строительный университет, 2016 © Оформление. ООО Компания «Ай Пи Ар Медиа», 2022
СОДЕРЖАНИЕ
Введение…………………………………………………………………….... 5
1. Задачи аппроксимации………………….………………………………... 10
1.1. Методы аппроксимации ……………...……………………………10
1.2. Интерполяция алгебраическими полиномами……………………11
1.3. Интерполяция кубическими сплайнами………………………......18
1.4. Сглаживание экспериментальных данных………………………..24
1.5. Применение сглаживающих сплайнов для обработки экспериментальных данных………………………………………. 30
1.6. Приближение рациональными функциями……………………… 36
1.7. Погрешность многочленной аппроксимации……………………. 38 Задание по теме «Задачи аппроксимации»…………………………... 40
2. Линейное программирование……………………………………………. 42
2.1. Постановка задач линейного программирования (ЛП)……........ 42
2.2. Симплекс-метод решения задачи ЛП…………………………….. 43
2.3. Составление оптимального плана выпуска продукции…………. 45
2.4. Двойственная задача ЛП…………………………………………... 46
2.5. Пример решения задачи ЛП………………………………………. 49
2.6. Анализ задачи ЛП на чувствительность…………………………. 53
2.7. Решение задачи ЛП с помощью табличного процессора Excel… 56
2.8. Транспортная задача (ТЗ)…………………………………………. 60
2.9. Решение ТЗ с помощью табличного процессора Excel………..... 69
2.10. Задача о назначениях………………………………….....……... 72
Задание по теме «Линейное программирование»………...………..... 80
3. Расчет устойчивости откосов земляного полотна…………………...… 82
3.1. Понятие устойчивости откосов насыпи и постановка задачи…. 82
3.2. Алгоритм расчета коэффициента устойчивости………………. 86
3.3. Пример расчета устойчивости откосов насыпи с применением электронных таблиц Excel……………………… 90
3.4. Уточнение значения коэффициента устойчивости…………...…94
3.5. Алгоритм метода поочередного изменения переменных……… 95 Задание по теме «Расчет устойчивости откосов земляного
полотна»………………………………………..............…………. 97
4. Модели массового обслуживания в теории транспортных потоков…. 98
4.1. Одноканальная система с отказами…………………………...….100
4.2. Одноканальная система с ожиданием …………………………...103
4.3. Многоканальная система массового обслуживания с отказами……………………………………………………….… 111
4.4. Многоканальная система массового обслуживания c ожиданием……………………………………………….……… 115
4.5. Расчет систем массового обслуживания методом Монте-Карло…………………………………… ……….118
4.6. Расчет систем массового обслуживания с помощью электронных таблиц Excel…………………………………….…. 119
3
Задание по теме «Модели массового обслуживания в теории
транспортных потоков»………...................................................... 126
Литература…………………………………………….…………………... 128
Приложение 1. Программа для построения сглаживающего сплайна...... 129
Приложение 2. Программа для расчета устойчивости откоса
земляного полотна………………………………………..……….….. 132
Приложение 3. Критические значения t – критерия Стьюдента…….... 134
Приложение 4. Равномерно распределенные случайные числа….……. 135
4
Введение
Математические модели составляют основу методов расчета и опти­мизации исследуемых процессов и явлений, поэтому изучению основных этапов и методов математического моделирования уделяется особое внимание.
Математическое моделирование – приближенное описание свойств объекта на математическом языке, а модель – объект-заменитель, который в определенных условиях может заменить объект-оригинал, воспроизводя интересующие нас свойства и характеристики оригинала. Основными эта­пами моделирования являются: построение модели объекта и ее формали­зация; исследование (решение) математических задач, к которым приводит математическая модель; идентификация и проверка адекватности модели, уточнение, интерпретация результатов моделирования.
На первом этапе осуществляется постановка задачи исследования, определяется объект исследования, выявляются основные особенности и факторы, влияющие на объект, формулируются законы, связывающие ос­новные элементы модели, определяются требования к исходной информа­ции и модели (точность, удобство, простота), выдвигаются гипотезы, предположения, описываются ожидаемые результаты. От содержательного описания объекта переходят к его математической модели, к записи на ма­тематическом языке.
На втором этапе исследуется сама математическая модель, устанав­ливается тип решаемой математической задачи, выбираются методы ее решения. Для оптимального проектирования конструкций, технологиче­ского оборудования, для определения оптимальных параметров модели могут рассматриваться оптимизационные модели: модели теории опти­мального управления, математического программирования. Для решения математической задачи составляется алгоритм расчета. Далее проводятся рабочие расчеты, вычислительный эксперимент, который базируется на использовании существующих комплексов и пакетов прикладных про­грамм, например, Mathcad, Mathematica, MATLAB и др. программ на алго­ритмическом языке. Проводится тестирование программ. При проведении вычислительного эксперимента применяются численные методы.
На третьем этапе осуществляется идентификация и проверка адек­ватности модели, а если необходимо, то и уточнение модели. Под иденти­фикацией модели понимается выбор переменных модели и ее параметров с последующей их оценкой в результате наблюдения или эксперимента, при этом определяются неизвестные параметры модели: адекватность модели – соответствие модели моделируемому объекту, ее точность. В качестве ко­личественных показателей оценки адекватности могут быть использованы различные критерии: относительная погрешность, критерий Фишера и т.д.
5
Если при проверке адекватности установлена недостаточная точность мо­дели, то осуществляется уточнение модели. После проведения идентифи­кации и оценки адекватности модели проводятся итоговые расчеты, вы­числительный эксперимент и интерпретируются результаты моделирова­ния, выдаются рекомендации.
В учебном пособии рассмотрено большое количество примеров, в основном расчетного характера, в которых применение излагаемых мето­дов иллюстрируется на конкретном практическом материале и доводится до числового результата. Все рассмотренные задачи доведены до алгорит­мов, пригодных для практической реализации с применением ЭВМ. Для наглядности проводимых расчетов в учебном пособии реализация вычис­лительных алгоритмов проводится в Excel, с применением электронных таблиц и программ в VBA.
Пособие состоит из четырех глав. Каждая глава содержит соответст­вующий теоретический материал и математическую модель, примеры ре­шения типовых задач с выбором численного метода решения и составле­нием программы на языке VBA или электронной таблицы в Excel с ис­пользованием ее надстройки «Поиск Решения». В конце главы приводятся варианты заданий для студентов с целью научить их пользоваться метода­ми вычислений с применением новых информационных технологий Excel.
Первая глава посвящена вопросам аппроксимации табличных зави­симостей, рассмотрены два основных подхода – интерполяция и сглажива­ние. Большинство численных методов решения инженерных задач так или иначе связаны с аппроксимацией функций. Это и собственно задачи ин­терполяции, сглаживания, приближения функций, и задачи численного решения дифференциальных и интегральных уравнений, в которых ап­проксимация присутствует как промежуточный этап исследования.
Типичной задачей приближения является задача интерполяции, и классический метод ее решения состоит в построении интерполяционного многочлена Лагранжа. Хотя согласно теореме Вейерштрасса всякая непре­рывная функция на отрезке может быть как угодно хорошо приближена многочленами, практические возможности применения многочленов Ла­гранжа ограничены. Прежде всего, чрезмерное увеличение числа узлов ин­терполяции усложняет вычисления и приводит к быстрому накоплению ошибок округления с ростом степени многочлена. Также следует отметить: С.Н. Бернштейном было установлено, что последовательность интерполя­ционных многочленов Лагранжа, построенных для некоторых непрерыв­ных функций по равноотстоящим узлам, с возрастанием количества узлов не стремится к интерполируемой функции. Поэтому на практике, для того чтобы достаточно хорошо приблизить функцию, вместо построения ин­терполяционного многочлена высокой степени используют интерполяцию «кусочными» многочленами, в частности, сплайн-функциями или просто сплайнами. Методы сплайн-функций широко применяются в вычисли-
6
тельной математике и инженерной практике. Этим вопросам в учебном по­собии уделено достаточно большое внимание.
Задача аппроксимации опытных данных, результатов наблюдений связана с задачей сглаживания экспериментальных зависимостей по мето­ду наименьших квадратов. Желательно обработать экспериментальные данные так, чтобы по возможности точно отразить общую тенденцию таб­личной зависимости, но в то же время сгладить случайные ошибки на­блюдений. Этот метод дает возможность так выбрать числовые параметры искомой зависимости, чтобы она в известном смысле наилучшим образом отображала экспериментальные данные. В данной главе приведены алго­ритмы построения интерполяционного и сглаживающего полинома, в ча­стности, сплайнов с использованием электронных таблиц Excel.
Во второй главе рассмотрены задачи линейного программирования и методы их решения с применением электронных таблиц Excel и над­стройки «Поиск решения» для выполнения оптимизации. Рассмотрены та­кие задачи, как составление оптимального плана выпуска продукции, оп­тимального использования ресурсов, транспортной логистики, задачи о на­значениях. Все эти задачи решаются сведением к модели линейного про­граммирования. Описание линейного программирования проводится на основе теории двойственности. Согласно этой теории, при рассмотрении задач линейного программирования можно построить прямую и двойст­венную задачи, они тесно связаны между собой и одна задача вытекает из другой. Применение двойственных задач особенно эффективно при анали­зе моделей на чувствительность. Рассмотрены вопросы устойчивости оп­тимального плана при изменении стоимости продукции и определены ин­тервалы устойчивости двойственных оценок по отношению к изменению количества сырья. Для решения задачи линейного программирования при­водится алгоритм симплекс-метода. Также в пособии изложен алгоритм для построения общей электронной таблицы для решения пары двойствен­ных задач, что позволяет наиболее полно изучить связи между ними.
В третьей главе рассматриваются методы расчета устойчивости от­косов земляного полотна. Определение устойчивости откосов является од­ной из основных задач при проектировании автомобильных дорог. От ка­чества проектирования земляного полотна во многом зависят стоимость строительства автомобильных дорог, надежность в эксплуатации и безо­пасность движения. В практике проектирования автомобильных дорог наибольшее распространение получил метод оценки общей устойчивости откосов шведского ученого Феллениуса – метод круглоцилиндрических поверхностей скольжения, который наиболее просто поддается компью­терной формализации.
Общая устойчивость откоса насыпи оказывается обеспеченной лишь в том случае, если сумма всех сил, сдвигающих массив обрушения, оказы­вается меньше суммы сил или их моментов, его удерживающих, т.е. при расчетном коэффициенте устойчивости, большем или равном единице.
7
Расчетный коэффициент устойчивости определяется как отношение суммы сил удерживающих к сумме сил сдвигающих. Однако учитывая некоторую погрешность методов расчета исходных данных, неучет некоторых факти­ческих условий работы и т.д., с инженерной точки зрения общая устойчи­вость откоса насыпи считается обеспеченной, если расчетный коэффици­ент устойчивости оказывается равным или больше нормативного коэффи­циента устойчивости.
Задача оценки общей устойчивости откосов земляного полотна сво­дится к отысканию центра критической кривой скольжения, характеризуе­мой наименьшим коэффициентом устойчивости. Ни один из известных ме­тодов расчета общей устойчивости откосов не дает сразу точного положе­ния центра наиболее опасной кривой скольжения, который может быть найден лишь методом последовательных приближений. В итоге расчет сводится к безусловной задаче нелинейного программирования. В пособии для определения коэффициента устойчивости приводятся расчетные фор­мулы и алгоритм расчета, на основании которых составлен программный код в VBA в Excel и электронная таблица Excel.
Расчет коэффициента устойчивости с помощью электронной таблицы Excel является наиболее важным для понимания алгоритма и методики расчета. Чтобы найти положение наиболее опасной кривой скольжения, рассматривается несколько возможных ее положений. Для этого намечает­ся семейство кривых, проходящих через подошву откоса и выходящих на поверхность земляного полотна в нескольких местах с определенным ша­гом. Для каждого положения кривой скольжения с помощью электронной таблицы Excel вычисляется коэффициент устойчивости, и среди них выби­рается наименьшее.
В четвертой главе рассматриваются модели различных систем мас­сового обслуживания в теории транспортных потоков, определяются эле­менты систем и их характеристики. При изучении систем массового об­служивания важно знать: распределение длины очереди, распределение времени ожидания в очереди, долю времени, когда система обслуживания свободна. Ответ на эти вопросы зависит непосредственно от характера распределения входящего потока и распределения времени обслуживания, которые могут подчиняться различным законам. В теории массового об­служивания могут применяться: закон равномерной плотности, показа­тельный закон, распределение Пуассона, распределение Эрланга.
Системы массового обслуживания, встречающиеся на практике, яв­ляются весьма разнообразными. Для того чтобы к их решению можно бы­ло подходить с некоторых позиций, вводится специальная классификация.
Первым является деление систем массового обслуживания на систе­мы с отказами и системы с очередью. Системы с очередью в зависимости от длины очереди подразделяются на системы с ограниченной и неогра­ниченной длиной очереди. По числу каналов системы подразделяются на одноканальные и многоканальные. Для некоторых случаев систем массо-
8
вого обслуживания для определения вероятностей состояния системы со­ставляются дифференциальные уравнения и определяются характеристики рассматриваемых систем, эффективность их работы. Рассмотрено множе­ство типовых примеров решения задач.
Для расчета систем массового обслуживания также рассмотрен метод Монте-Карло. Метод Монте-Карло применяется в тех случаях, когда по­строение аналитической модели изучаемого явления является трудным или вовсе неосуществимым. Смысл метода Монте-Карло состоит в том, что ис­следуемый процесс моделируется путем многократных повторений его случайных реализаций. Для реализации этого метода для некоторых случа­ев систем массового обслуживания составлены электронные таблицы Excel. С помощью моделирования путем многократных повторений слу­чайных реализаций (статистических испытаний) определены оценки и до­верительные интервалы числовых характеристик систем массового обслу­живания.
Настоящее учебное пособие будет способствовать развитию у сту­дентов навыков, необходимых для анализа и решения инженерных задач с применением ЭВМ. Предлагаемые в издании задачи являются хорошим средством для практического закрепления знаний, необходимых при мате­матическом моделировании и проведении вычислительных экспериментов.
9
1. ЗАДАЧИ АППРОКСИМАЦИИ
Аппроксимация заключается в приближенной замене некоторой функции, заданной аналитически или таблично, другой функцией, близкой исходной, но более простой и удобной для вычисления. Такая замена не­обходима, например, при работе на ЭВМ для вычисления значения функ­ции. Другая задача аппроксимации – восстановление функции на некото­ром отрезке по заданным дискретным точкам. В общем случае задача ап­проксимации включает два этапа:
1) выбор вида аппроксимирующей функции, который зависит от не-
которых параметров (коэффициентов);
2) определение числовых значений неизвестных коэффициентов ап­проксимирующей функции на основе выбранного критерия ее близости к исходной функции.
В качестве аппроксимирующей функции в вычислительной матема­тике часто рассматривают алгебраические, тригонометрические многочле­ны, специальные функции, сплайны и т.д., в качестве критерия близости – точное совпадение исходной и аппроксимирующей функции в заданных узловых точках (интерполяция), среднеквадратичное приближение (сгла­живание по методу наименьших квадратов), равномерное приближение (минимизация наибольшего значения абсолютной погрешности между ис­ходной и аппроксимирующей функцией либо на всем отрезке, либо на дискретном множестве точек).
Выбор вида аппроксимирующей функции и критерия близости зави­сят от самих исходных данных (как они получены) и поставленной задачи: для каких целей строится аппроксимирующая функция.
1.1. Методы аппроксимации
Для аппроксимации зависимости, заданной в виде таблицы, необхо­димо выбрать вид аппроксимирующей функции, включающий неизвест­ные параметры. Для определения неизвестных параметров их называют коэффициентами аппроксимирующей функции; существует два основных подхода – интерполяция и сглаживание. Если табличная зависимость пред­ставляет собой опытные данные, результаты некоторых наблюдений, то, как правило, коэффициенты аппроксимирующей функции определяют сглаживанием по методу наименьших квадратов. Результаты опытов, на­блюдений, когда данные получены из физических измерений, принято считать случайными, обладающими некоторым уровнем «шума». Требо­вать прохождение аппроксимирующей функции через опытные точки было бы повторением ошибок, поэтому применяют сглаживание. В тех случаях, когда функция задана табличным способом, или когда табличная зависи­мость является результатом некоторых расчетов, полученных на ЭВМ, на­пример, при решении дифференциальных уравнений и других математиче-
10
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]