- •Численные методы решения задач строительства
- •Часть 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
5.3.4. Коэффициент корреляции
Всегда ли существует функциональная зависимость между экспериментальными данными, заданными, например, таблицей 5.1.
Оценить функциональную близость (в линейном смысле) значений x, y можно с помощью коэффициента корреляции:
где
Если зависимость между значениями x, y линейная (y=ax+b), то:
для а > 0 -1< Rxy < 0,
для а < 0 0 < Rxy < 1
При Rxy =0 cвязь отсутствует.
Принято считать:
Rxy 0.3 – наблюдается слабая линейная связь,
Rxy = 0.3 – 0.7 – средняя,
Rxy 0.7 –сильная,
Rxy = 1 – линейная зависимость, все точки лежат на одной прямой.
5.3.5. Квадратичное (параболическое) приближение
Чаще всего линейная аппроксимация является достаточно грубым приближением. При решении задач строительства возникает необходимость использования более сложной аппроксимирующей функции.
Ограничимся случаем m=2. Уравнение регрессии 2-го порядка в этом случае называется квадратичным (или параболическим) и имеет вид
y=(x,a,b,c)=a+bx+cx2 . (5.21)
Неизвестные параметры a, b, c согласно МНК находим из условия минимизации функции S(a,b,c), суммы квадратов отклонений:
После дифференцирования и соответствующих преобразований получим нормальную систему для определения неизвестных параметров a, b, c.
Решая систему (5.23), получим уравнение регрессии 2-го порядка, степень точности такого приближения для исследуемого процесса оценивается по величине среднеквадратичного отклонения (5.19).
Если точность этого приближения не устраивает, повышают степень аппроксимирующей функции m (но надо помнить, что m<n–1).
5.3.6. Эмпирические формулы с двумя параметрами (метод выравнивания)
Для описания некоторых технологических процессов удобно использовать эмпирические формулы, содержащие два параметра:
y=(x, a, b). (5.24)
Например: y=a xb , y=a+b/x , y=x/(a+bx), y=a bx и т.д.
Пусть заранее известно, что экспериментальные точки M(xi,yi), i=1,…n, заданные табл. 5.1 не лежат на одной прямой. Для нахождения параметров a, b используется метод выравнивания.
Идея метода. Вводятся новые переменные (новая система координат X*OY*):
так, чтобы преобразованные точки M*(xi*, yi*), i=1,…n, в плоскости X*OY* могли быть аппроксимированы линейной зависимостью
y*=A+Bx*. (5.26)
Здесь
(i=1,2,…,n)
(рис.5.4).
Параметры А и В находятся методом наименьших квадратов (см. пункт 5.3.3).
Рис.5.4. Схема метода выравнивания
Пример 5.2. Пусть заранее известно, что экспериментальные точки M(xi, yi), i=1,2,…n, заданы в виде табл.5.1 и не лежат на одной прямой. А эмпирическая формула имеет вид:
y=a xb (5.27)
Для определения новой системы координат прологарифмируем выражение (5.27)
ln y=ln a + b ln x
и введем новые переменные:
y* = ln y; x* =ln x.
Обозначив A=lna; B=b, получим линейную эмпирическую формулу в новой системе координат X*OY* :
y*=A+Bx*
Неизвестные параметры А, В находим, используя МНК, и по аналогии с (5.17) строим нормальную систему
Переходя к старым переменным, получим систему уравнений для определения параметров a, b.
Решив эту систему относительно a, b (например, методом Крамара) и подставив их значения в выражение (5.27), получим нужную эмпирическую формулу. Насколько она хороша, можно оценить приведенным выше способом.
