Метод наименьших квадратов для аппроксимации экспериментальных данных

Автор работы: Пользователь скрыл имя, 10 Октября 2011 в 00:21, контрольная работа

Краткое описание

Рассчет относительных погрешностей методом наименьших квадратов для аппроксимации экспериментальных данных по заданным формулам.

Содержание работы

Теоретическая выкладка
Задание
Решение
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)

  (где  - параметры), значения которой при возможно мало отличались бы от опытных значений .

  Обычно  указывают класс функций (например, множество линейных, степенных, показательных  и т.п.) из которого выбирается функция  , и далее определяются наилучшие значения параметров.

  Если  в эмпирическую формулу (2.1.1) подставить исходные , то получим теоретические значения , где .

  Разности  называются отклонениями и представляют собой расстояния по вертикали от точек   до графика эмпирической функции.

  Согласно  методу наименьших квадратов наилучшими коэффициентами   считаются те, для которых сумма квадратов отклонений найденной эмпирической функции от заданных значений функции будет минимальной.

                                    (2.1.2)

  Построение  эмпирической формулы состоит из двух этапов: выяснение общего вида этой формулы и определение ее наилучших параметров.

  Если  неизвестен характер зависимости между  данными величинами x и y , то вид эмпирической зависимости является произвольным.

    Большое значение имеет изображение  полученных данных в декартовых  или в специальных системах  координат (полулогарифмической,  логарифмической и т.д.). По положению  точек можно примерно угадать  общий вид зависимости путем  установления сходства между  построенным графиком и образцами  известных кривых.

  Для того, чтобы найти набор коэффициентов  , которые доставляют минимум функции S , определяемой формулой (2.1.2), используется необходимое условие экстремума функции нескольких переменных - равенство нулю частных производных.  В результате получим нормальную систему для определения коэффициентов :

                          (2.1.3)

  Таким образом, нахождение коэффициентов  сводится к решению системы (2.1.3). Эта система может быть решена любым известным методом (методом Гаусса, простых итераций, формулами Крамера).

 

Задание

Получить  функциональную зависимость по экспериментальным данным по методу наименьших квадратов.

Функцию представить в виде:

  1. Полинома 3-й степени Y(x)=
  2. Тригонометрической 
  3. Логарифмической
  4. Степенной
  5. Показательной
 

Экспериментальные данные:

    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.

  1. Полином 3-й степени Y(x)=

Полином выглядит, как .

Определим коэффициенты , входящих в эмпирическую формулу.

Так как  известно, что они должны доставлять минимум функции S, составим уравнение: 
 

Воспользуемся необходимым условием экстремума функции  нескольких переменных - равенством нулю частных производных. 
 
 
 

Исходя  из полученных результатов, можем составить  нормальную систему.

Примечание: n=16 – это количество экспериментов.

Теперь  с помощью Excel произведем расчет сумм. Для этого:

  1. занесем в столбик экспериментальные значения x и y;

  1. в соседних столбцах посчитаем значения i-х элементов сумм, с помощью формул используя ссылки на значения x и y в построчном соответствии;

  1. с помощью  функции 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
 

Посчитав  значение сумм, подставим их в нормальную систему: 
 

На данном этапе нахождение коэффициентов сводится к решению системы. 

Для решения  используем метод Крамера:

  1. Отдельно создадим таблицы всех определителей. Значения элементов возьмем из основной таблицы, делая ссылки на значения сумм. Для этого, с помощью оператора “=”(присваивания) сделаем ссылки на ячейки основной таблицы.

  1. Для того чтобы посчитать значения определителей  надо использовать математическую функцию  Excel “МОПРЕД”, которая возвращает значение определителя матрицы.

D= 5,46045E+12, D=1,6032E+13, D=-1,6625E+12, D=9,19101E+11, D=-5,4554E+10.

  1. Имея значения определителей, можем посчитать значения искомых переменных системы. Для этого в ячейках надо записать: = ссылка на ячейку со значением определителя переменной / ссылка на ячейку со значением основного определителя. Получим:

= 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

Информация о работе Метод наименьших квадратов для аппроксимации экспериментальных данных