Построение параболы по методу наименьших квадратов

Вид материалаДокументы
Подобный материал:
Построение параболы по методу наименьших квадратов


Чтобы понять, как одна величина зависит от другой величины , например
  • длина стержня от температуры
  • показатель преломления данного материала от длины световой волны
  • высота свободно падающего тела от времени ,

мы проводим ряд экспериментов и получаем следующую таблицу.



























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


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

Запишем уравнение параболы в виде

(1) .

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

.

Эти расстояния, измеренные по вертикали, называются невязками. Для параболы, наилучшим образом описывающей наши данные, каждая из этих невязок должна быть мала – они должны быть малы в совокупности. В качестве критерия малости невязок выбирается сумма их квадратов

.

Чем меньше эта сумма, тем лучше парабола описывает наши данные. Таким образом, наилучшая парабола – это та, для которой сумма минимальна.

Минимум суммы ищем, приравнивая ее частные производные нулю

.

Или

.

Это система из трех уравнений с тремя неизвестными , после решения которой получаем параболу , наилучшим образом описывающую наши данные по методу наименьших квадратов.

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

.

В Mathcad имеются собственные средства, реализующие метод наименьших квадратов – это функции regress и interp.