§ 5 Метод наименьших квадратов (мнк).

Помимо двух вышеописанных способов оценки погрешности результата при косвенных измерениях, иногда применяют еще так называемый «метод наименьших квадратов» или сокращенно МНК. Этот метод можно использовать, если известен вид функциональной зависимостимежду измеряемыми физическими величинами, а требуется определитькоэффициенты, входящие в эту функцию. В наших лабораторных работах предлагается применять этот метод для определения параметров линейной зависимости.

Пусть в эксперименте можно измерить ряд значений некоторой величины xи, соответствующие им значения, величиныy. И пусть при этомизвестно, что между ними справедлива зависимость вида: y = ax + b.Как известно, такая зависимость графически представляется прямой линией (рис.4). Однако измеренные значенияx_i и y_i включают в себя погрешность и, в результате, не лягут идеально на прямую линию.

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

	Для аналитического выражения коэффициентов применяется метод наименьших квадратов. Утверждается, что наилучшей будет та прямая, сумма квадратов расстояний до которой, от всех экспериментальных точек будет минимальной. Расстояние (вдоль оси y) от точки с координатамиxi,yi до искомой прямой определяется выражением: (axi +b yi),
Рис.4

тогда сумма квадратов расстояний будет равна:

Решение задачи на нахождение минимума этого выражения (см. приложение, §7) приводит к следующим выражениям для коэффициентовaиb.

	(25)
	(26)

Дисперсию отклонения экспериментальных точек от прямой S₀²и дисперсию коэффициентовaиb  S_a² иS_b² можно вычислить по формулам:

S02 =	(27)
	(28)
Sb2=	(29)

Доверительные интервалы для коэффициентов aиb определяются как обычно:

a_p = S_a b_p = S_b

Если график исследуемой зависимости проходит через начало координат, то есть b=0, формулы 25, 27 и 28 существенно упрощаются. В этом случае коэффициентaи его дисперсию можно рассчитать по следующим формулам:

	(30)
	(31)
	(32)

Однако следует иметь в виду, что формулы 25-28 включают разности больших величин, мало отличающихся друг от друга, что легко может привести к ошибкам при вычислениях, если их проводить с недостаточным числом значащих цифр. Поэтому все промежуточные вычисления следует выполнять с большим числом значащих цифр, без округления. Если вы проводите вычисления с помощью компьютера, то это условие выполняется. В том случае, если расчеты выполняются «вручную», или с помощью не очень совершенного калькулятора, а результаты измерений имеют более трех верных знаков, то велика вероятность получить неправильный результат вычислений. В этом случае рекомендуется для вычисления коэффициентовaиb,а также доверительных границ их погрешности вместо вышеуказанных формул использовать выражения, преобразованные к другому виду:

	(33)
	(34)
	(35)
	(36)
	(37)

Во всех случаях проведения расчетов по МНК для определения коэффициентов линейной зависимости и их погрешностей необходимо вычислить некоторые суммы.Для этого удобно воспользоваться следующими таблицами – алгоритмами вычислений. Подсчитав суммы в каждом столбце и подставив их в соответствующие формулы, легко определить и значенияa,b,a и b

Таблица 5

xi	yi	xi2	yi2	xi yi
x1	y1	x12	y12	x1 y1
x2	y2	x22	y22	x2 y2
…	…	…	…	…
xn	yn	xn2	yn2	xn yn
xi	yi	xi2	yi2	x1 y1

Или в случае использования формул 33 -37

Таблица 6

xi	yi	xi	(xi )2	(xi )yi	xi2	(axi +b yi)2
x1	y1	x1	(x1 )2	(x1 )y1	x12	(ax1 +b y1)2
x2	y2	x2	(x2 )2	(x2 )y2	x22	(ax2 +b y2)2
…	…			…	…
xn	yn	xn	(xn )2	(xn )yn	xn2	(axn +b yn)2
xi	yi	(xi )	(xi )2	(xi )yi	xi2	(axi +b yi)2