Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Downloads / pdf_in_archive (2) / 1_obrabotka.pdf
Скачиваний:
38
Добавлен:
02.04.2015
Размер:
1.05 Mб
Скачать

Число степеней свободы равно здесь n–2, так как для определения прямой регрессии необходимо выполнение двух дополнительных условий. Дисперсии величин a и b равны:

S a2 =

S 2

n

 

 

 

n

 

 

 

,

(2.104)

n

 

n

2

 

 

 

 

nx i2

 

x i

 

 

 

1

 

1

 

 

 

n

S n2 x i2

Sb2 =

 

1

 

 

.

(2.105)

n

 

n

2

 

 

 

 

nx i2

x i

 

 

1

 

1

 

 

 

Для наклона прямой регрессии a и отсекаемого ею отрезка b в принципе справедливы все соображения изложенные в разделе 2.2.3., если бы было известно распределение вероятности для двумерной генеральной совокупности. Здесь вновь возникает вопрос о доверительном интервале, который показывает, с какой статистической достоверностью эти величины можно определить по данной выборке. Если x и y распределены нормально, то доверительный интервал определяется с использованием распределения Стьюдента. Если обе переменные равноправны или между ними нет функциональной зависимости, то для обработки результатов измерений используется корреляционный анализ. Задача о нелинейной регрессии решается аналогично, причем в качестве кривых регрессии используют полиномы разной степени. Ниже приведены соответствующие результаты.

2.3.2. Нелинейная регрессия

Пусть зависимость между величинами x и y дана в виде полинома:

y = a

+ a x +…+a

n

x n . Требуется определить неизвестные параметры

0

1

 

a0 , a1,, an . Алгоритм решения этой задачи строится следующим образом:

1.Проводится N совместных измерений величин x и y (N > n +1).

2.Составляется система условных уравнений:

y

i

= a + a x

i

+…+a

n

x n + ε

i

, i =1,2,, N

 

0 1

 

 

i

 

гдеε i= yi y(xi )

– как и выше отклонение измеренного значения yi от

истинного y(xi).

3.В предположении, что результаты измерений распределены нормально, взаимнонезависимы и ошибкой измерения xi можно пренебречь, оценки

параметров могут быть получены минимизацией суммы квадратов отклонений (невязок):

S = N ε 2i = N [yi (a0 + a1x i +…+an x in ]2 min , i=1 i=1

где параметры ak рассматриваются как неизвестные. Приравнивая к 0 первые производные от S по каждому параметру получаем систему уравнений:

S

= 0;

S

= 0; ;

S

= 0 .

a

a

a

n

 

 

 

0

 

1

 

 

 

После несложных преобразований система нормальных уравнений записывается в виде:

mmax a j [X j+k ]= [YX k ] ,

(2.106)

i=mmin

где mmin – наименьшая степень полинома; mmах – наибольшая степень полинома; k = mmin ,, mmax ; [·] – скобки Гаусса:

[X j+k ]= N

xij+k ;[YX k ]= N

yi xik .

 

 

 

 

 

i=1

 

 

 

 

i=1

 

 

4. Оценки параметров получаются решением системы (2.106):

 

a j

 

 

 

 

 

,

 

 

 

 

 

(2.107)

 

 

 

 

 

 

 

 

= j

 

 

 

 

 

 

 

 

где

 

– главный определитель системы; j – определитель, получаемый

из главного заменой j-го столбца столбцом правых частей.

 

5. Оценки дисперсий параметров даются выражением:

 

S 2

=

 

 

 

jj

S 2 ,

 

 

 

(2.108)

 

 

 

 

 

 

 

 

aj

 

 

 

 

 

 

0

 

 

 

 

где

 

 

 

jj – определитель, получаемый из главного вычеркиванием j-го

столбца и j-й строки;

 

 

S02 =

 

 

 

1

 

 

 

N [yi

(a0 + a1 xi +…+ an xin ]2 .

(2.109)

N

 

 

 

 

 

 

 

 

 

(n +1) i=1

 

 

 

 

6. Оценка дисперсии функции равна:

 

S y2(x) = Sa2+ x 2 Sa2

 

+…+ x 2n Sa2

 

 

 

 

 

 

 

 

 

0

 

1

n

 

7. Доверительные интервалы для параметров и функции:

 

 

=

t(P,k)Saj ; y(x)

=

t(P,k)S y( x) ,

 

 

a j

 

 

 

 

где t(P, k) – коэффициент, определяемый из таблиц распределения Стьюдента для доверительной вероятности P и k=n–1 – степеней свободы.