- •Численные методы решения задач строительства
- •Часть 1
- •Предисловие
- •В ведение
- •Общие сведения о вычислительном эксперименте и математическом моделировании
- •Численные методы.
- •Погрешности вычислений
- •Понятия точности, устойчивости и сходимости при численном решении
- •1.1. Матрицы и векторы.
- •1.2. Матрицы специального вида
- •Глава 1 Основные понятия матричного исчисления
- •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
2.1. Первый этап. Отделение корней
Отделить (локализовать) корни – это значит выделить из области допустимых значений функции y=f(x) отрезки, в каждом из которых содержится единственный корень.
Для функции общего вида не существует универсальных методов решения задачи локализации корней. Отделить корни уравнения f(x)=0 можно разными способами. Рассмотрим некоторые из них.
Табулирование функции. Строится таблица значений функции у=f(x) на некотором отрезке x[a, b]. Если окажется, что для соседних значений аргумента значения функции имеют разные знаки, то хотя бы один корень уравнения 2.1 может находиться между ними.
Графический метод. Если удается построить график функции y=f(x), то можно определить количество и расположение нулей функции, выделяя те промежутки оси Х, где график y=f(x) пересекает ось Х.
Если построение графика y=f(x) затруднительно, то исходное уравнение f(x)=0 заменяется эквивалентным ему уравнением (x)=(x) и строятся графики функций у1=(x) и у2=(x). Искомый корень является абсциссой точки пересечения графиков этих функций.
Исходя из физического смысла задачи.
Убедиться в том, что на данном отрезке x[a, b] (например, грубо определенном графическим способом) действительно имеется единственный корень уравнения (2.1), можно аналитическим способом, в основе которого лежит известная теорема математического анализа [3]:
Теорема 2.1. Если непрерывная на отрезке [a, b] функция y=f(x) принимает на концах его противоположные знаки, т.е. f(a)f(b)<0, то внутри этого отрезка содержится, по меньшей мере один корень уравнения f(x)=0. Корень заведомо будет единственным, если производная функции y=f '(x) существует и сохраняет постоянный знак внутри интервала (a, b), т.е. функция монотонна на отрезке [a, b]. |
Пример 2.1. Отделить (локализовать) корни уравнения
4 – ех – 2х2 = 0, (2.2)
т.е. выяснить, сколько корней имеет это уравнение, и найти интервалы, в которых находятся по единственному корню.
Способ 1. Составим таблицу значений функции f(x)=4 – ех – 2х2 на промежутке [–3, 1] (табулирование функции).
Таблица 2.1
x |
-3 |
-2 |
-1 |
0 |
1 |
f(x)=4-ех-2х2 |
–14,049 |
–4,135 |
1,632 |
3,012 |
–0,718 |
Из таблицы видно, что на отрезках [–2,–1] и [0,1] существуют, по крайней мере, по одному корню уравнения (2.2).
Способ 2. Уравнение (2.2) эквивалентно уравнению (x)=(x), где. (x)=4–2х2 ; 3) (x)= ех.
Построим графики функций у=(x) и у=(x) (рис. 2.2). Они пересекаются в двух точках, абсциссы которых х*1 и х*2 являются решениями уравнения (x)= (x), т.е. решением уравнения (2.2).
Рис. 2.2. Графический способ отделения корней
Пример 2.2. Локализовать корни уравнения
х – tg(x)=0, (2.3)
которое получается
при решении задачи устойчивости стержня
(рис.2.3), где
При решении задач устойчивости нас обычно интересует наименьшее значение критической нагрузки, т.е. надо найти наименьший положительный корень уравнения (2.3).
Т
аблица
2.2
|
f(x)=х-tg(x) |
0 |
0,000 |
0,5 |
–0,046 |
1 |
–0,557 |
1,5 |
–12,601 |
2 |
4,185 |
2,5 |
3,247 |
3 |
3,143 |
3,5 |
3,125 |
4 |
2,842 |
4,5 |
–0,137 |
5 Рис.2.3. К зада-че устойчи-вости стержня |
8,381 |
Рис.2.4. Локализация корней уравнения х–tg(x)=0
Если при решении данной задачи отделение корней производить на основании таблицы табулирования (табл.2.2), то можно допустить ошибку, предположив, что корень уравнения находится на отрезке [1.5, 2], где функция меняет знак.
В действительности, на этом участке функция f(x)=х–tg(x) терпит разрыв (т.е. не выполняются условия теоремы 2.1) и это хорошо видно на рис.2.4.
Таким образом, искомый корень уравнения находится на отрезке [4, 4.5], где выполняются условия теоремы 2.1.

х