Метод наименьших квадратов. В различных исследованиях приходится пользоваться формулами, составленными на основании экспериментаВ различных исследованиях приходится пользоваться формулами, составленными на основании эксперимента. Одним из лучших способов получения формул является метод наименьших квадратов. Пусть на основании эксперимента необходимо установить функциональную зависимость между переменными величинами x и y, например, между температурой и удлинением прямолинейного металлического стержня. Производим n измерений, по результатам составляем таблицу
При этом вид функции устанавливается из теоретических исследований, или по характеру положения на координатной плоскости экспериментальных точек. Пусть, например, точки, взятые из таблицы, расположены так, как показано на рис. 45. В данном случае естественно предположить, что между и существует линейная зависимость, выражающаяся формулой (1) Мы ограничимся рассмотрением случая линейной зависимости. Так как точки (x1;y1),(x2;y2),…,(xn;yn) приблизительно лежат на одной прямой, то формула (1) является приближенной. Поэтому, подставляя их координаты в формулу (1) вместо и получим следующие равенства: , где некоторые числа, которые назовем погрешностями. Возникает задача – подобрать коэффициенты таким образом, чтобы эти погрешности были возможно, меньше по абсолютной величине. Методом решения этой задачи и является метод наименьших квадратов. Согласно этому методу рассмотрим сумму квадратов погрешностей: где и - заданные числа, а коэффициенты – неизвестные величины, подлежащие определению, т.е. можно рассматривать как функцию двух переменных и исследовать ее на экстремум. Имеем Приравнивая эти частные производные нулю, получаем линейную систему двух уравнений с двумя неизвестными : Система (2) называется нормальной системой метода наименьших квадратов. Из нее мы находим числа , подставляя их в уравнение (1), получаем форму искомой прямой. Во – первых, для разрешимости системы (2) потребуется условие = (3) Лемма. Величина в правой части (3) равна и, следовательно, больше 0. Доказательство. Правая часть этого равенства равна Эту сумму легко сгруппировать и получить . Итак, такие a, b, чтобы выполнялась система (2), существуют. Чтобы проверить, что в этих точках функция S(a,b) действительно имеет минимум, вычислим . Следовательно, определитель
, по лемме имеет положительные главные миноры, поэтому найденная точка – точка минимума. Пример. Пусть в результате эксперимента получены пять значений искомой функции у при пяти значениях аргумента х (n=5), которые записаны в таблице: дем искать функциональную зависимость между x и y в виде линейной функции y=ax+b. При составлении нормальной системы (2) для определения коэффициентов a и b, вычисляем Система (2) принимает вид 25a+5b=16,5, 5a+5b=8. Решая эту систему, находим: a=0,425,b=1,175. Отсюда формула искомой прямой есть y=0,425x+1,175.
|