Метод наименьших квадратов для аппроксимации экспериментальных данных
Контрольная работа, 10 Октября 2011, автор: пользователь скрыл имя
Краткое описание
Рассчет относительных погрешностей методом наименьших квадратов для аппроксимации экспериментальных данных по заданным формулам.
Содержание работы
Теоретическая выкладка
Задание
Решение
1) Полином 3-й степени Y(x)
2) Тригонометрическая функция
3) Логарифмическая функция
4) Степенная функция
5) Показательная функция
Вывод
Содержимое работы - 1 файл
Оглавление.docx
— 150.30 Кб (Скачать файл)Оглавление
Теоретическая выкладка 2
Задание 4
Решение 5
1) Полином 3-й степени Y(x) 5
2) Тригонометрическая функция 10
3) Логарифмическая функция 13
4) Степенная функция 16
5) Показательная функция 19
Вывод 24
Теоретическая выкладка
Очень часто, особенно при анализе эмпирических данных возникает необходимость найти в явном виде функциональную зависимость между величинами x и y , которые получены в результате измерений.
При
аналитическом исследовании взаимосвязи
между двумя величинами x и y производят
ряд наблюдений и в результате,
получается таблица значений:
| x | ¼ | ¼ | ||||
| y | ¼ | ¼ |
Эта таблица обычно получается как итог каких-либо экспериментов, в которых (независимая величина) задается экспериментатором, а получается в результате опыта. Поэтому эти значения будем называть эмпирическими или опытными значениями.
Между величинами x и y существует функциональная зависимость, но ее аналитический вид обычно неизвестен, поэтому возникает практически важная задача - найти эмпирическую формулу
(где - параметры), значения которой при возможно мало отличались бы от опытных значений .
Обычно указывают класс функций (например, множество линейных, степенных, показательных и т.п.) из которого выбирается функция , и далее определяются наилучшие значения параметров.
Если в эмпирическую формулу (2.1.1) подставить исходные , то получим теоретические значения , где .
Разности называются отклонениями и представляют собой расстояния по вертикали от точек до графика эмпирической функции.
Согласно методу наименьших квадратов наилучшими коэффициентами считаются те, для которых сумма квадратов отклонений найденной эмпирической функции от заданных значений функции будет минимальной.
Построение эмпирической формулы состоит из двух этапов: выяснение общего вида этой формулы и определение ее наилучших параметров.
Если неизвестен характер зависимости между данными величинами x и y , то вид эмпирической зависимости является произвольным.
Большое значение имеет
Для
того, чтобы найти набор
(2.1.3)
Таким образом, нахождение коэффициентов сводится к решению системы (2.1.3). Эта система может быть решена любым известным методом (методом Гаусса, простых итераций, формулами Крамера).
Задание
Получить функциональную зависимость по экспериментальным данным по методу наименьших квадратов.
Функцию представить в виде:
- Полинома 3-й степени Y(x)=
- Тригонометрической
- Логарифмической
- Степенной
- Показательной
Экспериментальные данные:
| X | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 16 |
| Y | 0 | 2 | 5 | 9 | 8 | 3 | 1 | 0 | 1 | 6 | 13 | 10 | 8 | 5 | 2 | 0 |
Рассчитать относительные погрешности по каждому типу функций.
Решение
Для расчетов будем использовать среду Microsoft Office Excel.
- Полином 3-й степени Y(x)=
Полином выглядит, как .
Определим коэффициенты , входящих в эмпирическую формулу.
Так как
известно, что они должны доставлять
минимум функции S, составим уравнение:
Воспользуемся
необходимым условием экстремума функции
нескольких переменных - равенством нулю
частных производных.
Исходя из полученных результатов, можем составить нормальную систему.
Примечание: n=16 – это количество экспериментов.
Теперь с помощью Excel произведем расчет сумм. Для этого:
- занесем в столбик экспериментальные значения x и y;
- в соседних столбцах посчитаем значения i-х элементов сумм, с помощью формул используя ссылки на значения x и y в построчном соответствии;
- с помощью функции Excel “СУММ” найдем суммы;
В общем виде это будет выглядеть так:
| №\пер. | X | Y | X^2 | X^3 | X^4 | X^5 | X^6 | Y*X | Y*X^2 | Y*X^3 |
| 1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
| 2 | 1 | 2 | 1 | 1 | 1 | 1 | 1 | 2 | 2 | 2 |
| 3 | 2 | 5 | 4 | 8 | 16 | 32 | 64 | 10 | 20 | 40 |
| 4 | 3 | 9 | 9 | 27 | 81 | 243 | 729 | 27 | 81 | 243 |
| 5 | 4 | 8 | 16 | 64 | 256 | 1024 | 4096 | 32 | 128 | 512 |
| 6 | 5 | 3 | 25 | 125 | 625 | 3125 | 15625 | 15 | 75 | 375 |
| 7 | 6 | 1 | 36 | 216 | 1296 | 7776 | 46656 | 6 | 36 | 216 |
| 8 | 7 | 0 | 49 | 343 | 2401 | 16807 | 117649 | 0 | 0 | 0 |
| 9 | 8 | 1 | 64 | 512 | 4096 | 32768 | 262144 | 8 | 64 | 512 |
| 10 | 9 | 6 | 81 | 729 | 6561 | 59049 | 531441 | 54 | 486 | 4374 |
| 11 | 10 | 13 | 100 | 1000 | 10000 | 100000 | 1000000 | 130 | 1300 | 13000 |
| 12 | 11 | 10 | 121 | 1331 | 14641 | 161051 | 1771561 | 110 | 1210 | 13310 |
| 13 | 12 | 8 | 144 | 1728 | 20736 | 248832 | 2985984 | 96 | 1152 | 13824 |
| 14 | 13 | 5 | 169 | 2197 | 28561 | 371293 | 4826809 | 65 | 845 | 10985 |
| 15 | 14 | 2 | 196 | 2744 | 38416 | 537824 | 7529536 | 28 | 392 | 5488 |
| 16 | 16 | 0 | 256 | 4096 | 65536 | 1048576 | 16777216 | 0 | 0 | 0 |
| СУММА | 121 | 73 | 1271 | 15121 | 193223 | 2588401 | 35869511 | 583 | 5791 | 62881 |
Посчитав
значение сумм, подставим их в нормальную
систему:
На данном
этапе нахождение коэффициентов сводится
к решению системы.
Для решения используем метод Крамера:
- Отдельно создадим таблицы всех определителей. Значения элементов возьмем из основной таблицы, делая ссылки на значения сумм. Для этого, с помощью оператора “=”(присваивания) сделаем ссылки на ячейки основной таблицы.
- Для того чтобы посчитать значения определителей надо использовать математическую функцию Excel “МОПРЕД”, которая возвращает значение определителя матрицы.
D= 5,46045E+12, D=1,6032E+13, D=-1,6625E+12, D=9,19101E+11, D=-5,4554E+10.
- Имея значения определителей, можем посчитать значения искомых переменных системы. Для этого в ячейках надо записать: = ссылка на ячейку со значением определителя переменной / ссылка на ячейку со значением основного определителя. Получим:
= 2,936017, = -0,30447, = 0,16832, = -0,00999 .
Найдя коэффициенты подставим их в формулу полинома и вычислим значение аппроксимированных Y-ов от X. Результаты запишем в таблицу и создадим графики Yэксп(Х) и Yапрокс(Х).
| № | X | Yэксп | Yапрокс |
| 1 | 0 | 0 | 2,936017 |
| 2 | 1 | 2 | 2,78988 |
| 3 | 2 | 5 | 2,920439 |
| 4 | 3 | 9 | 3,267747 |
| 5 | 4 | 8 | 3,771861 |
| 6 | 5 | 3 | 4,372836 |
| 7 | 6 | 1 | 5,010728 |
| 8 | 7 | 0 | 5,625592 |
| 9 | 8 | 1 | 6,157483 |
| 10 | 9 | 6 | 6,546458 |
| 11 | 10 | 13 | 6,73257 |
| 12 | 11 | 10 | 6,655877 |
| 13 | 12 | 8 | 6,256434 |
| 14 | 13 | 5 | 5,474295 |
| 15 | 14 | 2 | 4,249517 |
| 16 | 16 | 0 | 0,232265 |