- •Численные методы расчеты в microsoft Excel при исследовании автомобильных конструкций
- •Содержание
- •Введение
- •1. Вычисление определенных интегралов Справочная информация
- •Формула трапеций
- •Контрольные задания
- •2. Решение нелинейных уравнений Справочная информация
- •Контрольные задания
- •3. Решение систем линейных алгебраических уравнений Справочная информация
- •Метод Гаусса с выбором главного элемента
- •Программное обеспечение Excel’а
- •Контрольные задания
- •4. Интерполяция таблично заданных функций Справочная информация
- •Интерполяция многочленами
- •Контрольные задания
- •5. Аппроксимация таблично заданных функций Справочная информация
- •Контрольные задания
- •6. Решение задачи коши для обыкновенных дифференциальных уравнений 1-го порядка Справочная информация
- •Усовершенствованный метод Эйлера
- •Оценка погрешностей методов
- •Контрольные задания
- •7. Решение задачи коши для нормальных систем обыкновенных дифференциальных уравнений и уравнений высших порядков Справочная информация
- •Метод Эйлера
- •Оценка погрешности метода
- •Контрольные задания
- •Список литературы
- •Численные методы расчеты в microsoft Excel при исследовании автомобильных конструкций
- •107023, Г.Москва, б.Семёновская ул., 38.
Контрольные задания
С помощью одного из методов интерполяции сформировать образ кривой, визуально совпадающий с графиком заданной функции. В качестве такой функции взять левую часть алгебраического уравнения из лабораторной работы №2, которая соответствует номеру выполняемого варианта. В качестве отрезка интерполяции использовать указанный там же отрезок поиска корней. Точки интерполяции (в количестве 5–7 штук) распределить на заданном отрезке по своему усмотрению наилучшим образом. На рисунке показать график исходной функции, выбранные узловые точки и график интерполирующей функции.
Рис.5.
Рис.6.
5. Аппроксимация таблично заданных функций Справочная информация
И
Рис.1.
Через множество узловых точек таблично заданной функции можно провести бесконечное количество аппроксимирующих кривых. Задача выбора единственной из них делится на две основные подзадачи:
-
выбор аналитических зависимостей, отражающих физику взаимосвязи аргумента и реальной функции, когда должен быть определен общий вид приближающей функции;
-
выбор критерия достоверности описания реальной функции с помощью выбранных зависимостей.
Существует множество подходов к построению вида приближающей функции, как функции с параметрами. Одним из них является выбор в качестве аппроксимирующей зависимости линейной комбинации некоторых известных аналитических функций. Вместе они должны отражать суть физического процесса, описываемого исходной функцией, и быть линейно независимыми на отрезке аппроксимации [x1, xn]
.
Функции φk(x) часто выбираются в виде полиномов, частным случаем которых являются степенные функции
φ1(x) = 1, φ2(x) = x, φ3(x) = x2, φ4(x) = x3,…,
в виде тригонометрических косинусов
,
или в любом другом удобном для исследователя виде.
Другим подходом к построению приближающей функции является её представление сплайнами. Это избавляет исследователя от необходимости подбирать аналитические функции для аппроксимирующей зависимости и часто даёт результат, отвечающий всем требованиям, которые предъявляются к процессу аппроксимации.
В качестве критерия достоверности описания реальной функции Гауссом (1794) и Лежандром (A.M.Legendre, 1805) было предложено использовать сумму квадратов отклонений значений аппроксимирующей функции от ординат узлов таблично заданной
,
где отклонение от каждой узловой точки Δi, показанное на рис.2, вычисляется как
.
С
Рис.2.
Необходимым условием экстремума квадратичной функции многих переменных F является равенство нулю всех её частных производных по параметрам c1, c2, ..., cm
.
Можно показать, что для функции F, являющейся суммой квадратов отклонений, достаточные условия существования её минимума в стационарной точке выполняются тождественно. Поэтому условиями существования экстремума функции F можно пользоваться как условиями её минимума, что позволяет привести задачу аппроксимации n значений табличной функции к задаче решения системы из m линейных алгебраических уравнений с симметричной матрицей относительно этих параметров
где
,
.
Параметры аппроксимации c1, c2,..., cm, определяемые как решение вышеприведённой системы линейных уравнений, которая сформирована для заранее выбранных функций φ1(x), φ2(x), …, φm(x), дают наименьшее значение целевой функции F. Для каждого набора таких функций будет получаться своё наименьшее значение F. Поэтому для выбора наилучшей аппроксимации ориентируются на наименьшее значение погрешности аппроксимации, которая рассчитывается следующим способом для каждого набора функций φ1(x), φ2(x), …, φm(x)
, ,
где под нормой таблично заданной функции понимается евклидова, квадратичная норма
.
Работа метода может быть проиллюстрирована на примере аппроксимации функции, заданной 8-ю узловыми точками, показанными на рис.3, и поиска её значения при х = 1.5.
Для построения
зависимости, аппроксимирующей таблично
или графически заданную функцию,
исследователь должен подобрать
аналитические
функций, которые своей комбинацией
отражают описываемый процесс. В данном
случае исходная функция может быть
описана как комбинация двух функций.
Первая из них должна отражать обратно
пропорциональную зависимость функции
от аргумент
Рис.3.
.
Необходимо заметить, что данное представление аппроксимирующей зависимости не является единственным. Можно подобрать и другие комбинации элементарных функций, которые отражают общий характер рассматриваемой табличной функции.
В соответствии с приведённым выше алгоритмом коэффициенты системы линейных алгебраических уравнений, определяющей параметры аппроксимации c1 и c2, будут вычисляться следующим образом
,
,
,
.
Вычисления по этим формулам удобнее выполнять, сняв с графика координаты узловых точек и сформировав из них следующую таблицу
i |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
x |
0.5 |
0.5 |
1 |
1.5 |
1.5 |
2.5 |
3.5 |
4 |
y |
0.9 |
0.8 |
0.5 |
0.5 |
0.4 |
0.4 |
0.5 |
0.5 |
В этом случае требуемые коэффициенты будут
,
,
,
,
.
Далее решение задачи приводится к решению системы из двух линейных алгебраических уравнений
Эта система имеет следующее решение c1= 0.105, c2= 0.401. Таким образом, аппроксимирующая функция имеет вид
.
Её значения при табличных значения аргумента приведены ниже
x |
0.5 |
1 |
1.5 |
2.5 |
3.5 |
4 |
y |
0.855 |
0.506 |
0.425 |
0.422 |
0.481 |
0.519 |
С их помощью может быть вычислено значение целевой функции F
,
которое определяет погрешность аппроксимации
, ,
где норма таблично заданной функции была вычислена следующим образом
.
Решение этой задачи на ПЭВМ в программе Excel выглядит следующим образом. Оно начинается с ввода исходных данных в виде таблицы значений сглаживаемой функции (см. рис.4) и расчёта базисных функций. Затем выполняется вычисление коэффициентов системы уравнений, определяющей коэффициенты c1 и c2 аппроксимирующей зависимости. Это делается на основе результатов расчёта базисных функций с помощью стандартной функции СУММПРОИЗВ, которая позволяет находить сумму произведений элементов массивов. Система линейных алгебраических уравнений решается с использованием встроенной функции МОБР построения обратной матрицы системы. В заключении рассчитывается аппроксимирующая сглаживающая функция и строится её график, на котором точками отображаются значения исходной табличной функции.
Рис.4.