Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Болдырев, В. С. Методы математической статистики в гидрографии и кораблевождении учеб. пособие

.pdf
Скачиваний:
7
Добавлен:
19.10.2023
Размер:
6.33 Mб
Скачать

§ І б . Уравнивание с условны м и уравнениями

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

вие

х т я + л,

= 180°,

где

х ;

-

истинное

значение

угла.

 

 

 

 

L2

 

L3

 

 

Если обозначить

через

^

,

,

результаты из­

мерения,

то

 

 

 

 

 

 

 

 

 

L, + t 2 *

 

480*+ -и* 9

 

 

 

где и /

- разность

между

суммой углов и

ІЭ0°,

образо­

вавшаяся за счет ошибок измерений. В задачу уравнивания

входит

отыскание таких поправок v-t ,

-с£ ,

к ре­

зультатам

измерений, которые обратили

бы величину си-

в ноль,

то

есть

 

 

 

 

і/ ' + </2 т tf3 + си-= О .

 

 

Это

уравнение называется условным уравнением попра­

вок. В

данном случае одно избыточное

измерение

привело

к образованию одного условного уравнения. Таким образом,

если

сделано п

измерений, а число независимых

условий

равно

т

, то

это

значит, что было произведено

г

из­

быточных измерений.

Обозначим через x t

, х 2

,

 

 

уравненные

значения

неизвестных,

через

£ ,

Л2

 

 

 

результаты измерений.

 

 

 

 

 

Пусть из теоретических соображений известно

т

со­

отношений вида

 

 

 

 

 

 

 

і №! •>л% ? • • • ?

^—^ 7

і у2

 

 

 

.

7

Q7

Подставим вместо х -

£ • и тогда получим

 

СГ

&

 

7^2 V

• »

 

fi (£itC/t 1^2f

°n)~0 ‘

(2.22)

Разложим (2.22) в ряд Тейлора и граничимся первыми членами разложения. Вычтем из полученного ряда соот­ ношение (2 .22):

%

м .

+ M l

и£- = О

 

Ч т

dLn

г

 

 

 

 

 

і жI j 2j 3 , . . . j Т

Таких

равенств

будет г .

Обозначим:

дА

 

ЭІ,

_

s

dl,

6/ ’

іи

Гч

э/.

 

*

вА,

, .

д12 ~ а^

діп ~

* ’

dl2 -

6f

 

^А2 _

£ .

 

" ’

ЭЛ _

r

 

' ’

■df*. = 1п

dlZ ~

2 ’ * '

UBfыг-*

 

С учетом этих обозначений

а, ^

а2 ^

* • • + <*пЧг + “ ) =0

1

5

^

• • • + * п * я +

 

 

*

 

98

Умножим систему на

—2 к{

-2 к г

-2 к 7

Выполнив умножение, получим:

- г к, (a,ts, + а2о2 +

• •

- + а п ^ t u A ^ O ;

- 2 k 2(é l t/l i-

*

•+ én Vn i-u/2)= 0 j

- Z k z ( i , ^ - i 2 ts2 + ••

- + W

« £ ) « 0 ;

 

 

 

 

T C/

 

= m in .

 

 

 

 

гг

 

 

Сложим все полученные уравнения и найдем минимум от

суммы по V-, • ■■<sn :

 

 

 

 

 

 

Ж ,

■■2(/l - Z alkr

2 é l k2----------2 к 27 = о ;

Q.J.

2tS2 2 аг к, ~ 2 ё 2 кг

‘ - 2 к 2 7г О ;

а г

= 2<sn - 2 a n к -2

---------- 2 А 7 7 ^ = 0 .

-

діл

 

 

 

 

 

 

 

Итак мы получили

п

уравнений

с

п неизвестными.

Последнюю систему

перепишем таким образом:

 

 

 

 

«/

 

 

 

 

 

+ 7г кг

аг

 

 

 

* ~*пк2

ап

 

€я

 

99

Произведем умножение, как показано,

и сложим по

столбцам. В результате получим:

 

[аіл]=А( [aa]f aé] +-

• * + кг [<**]

 

[&/] =к, [а£\ + кг [ Щ +

• • + к7 [ fa ] ;

(? .23)

Н = kt [а?] + к2 [£?] + ■ * * + к7 [7 7 ] •

^

Подставим полученные равенства в уравнения поправок:

к, \jaofj +

кг

[а^] + • • •+ [az] + tA = £7 ;

kf

 

_aé']+ /ср (/»] + • ' • + кг [/?] + а/г= О )

к\оп}+

кг [£*]і~ • • + к7[ 77]+ о/?= 0 .

Эта система уравнений называется системой нормаль­ ных уравнений коррелят, а система (2.23) - уравнением коррелят. Неопределенные множители kt $ к г , . . . , к 7

называются коррелятами.

Реная систему нормальных уравнений, найдем корреляты,

подставив которые

в (2 .23),

получим искомые

поправки.

Поправками

, . . . ,

необходимо

исправить

результаты измерений и получить уравненные значения не­

известных.

 

 

 

 

£ =

Ппчеп.

В

треугольнике

измерены три угла

43° 15»

17%

1 2 = 78? 43*

50%

13 =58° 01»

05".Най­

дем л,

, х г

,

х 3 . Составив условное уравнение:

L+ L2 fL 3 ~ і80°= u s ,

u s

= i Z ' .

 

100

Таким образом, уелекнее уравнение справок

р данном случае yew ine одно, значит, будет одно нор­ мальна« уинміеіик и одна коррелята. Уравнения коррелят:

cj к , OlI OCg3 Ol2 I

t/} =k a2

f} ^ k o s f

Тогда

3k + i2 " -0 ; к ~ - ч " ;

o-3 = -ч/ 'l.

Таким образом,

jct = l f - с/ = 't3°I5* 13" f

* i ~ L2 ~ fJ2 * 7 & W 4 6 " i

x3 =L3 ~ b% = 58°0I»0I-.

■Раиаер.точные измерения

При этом, как и в предыдущем случае,

£ р v 1/J = m i n .

Уравнения коррелят получатся такими же, за исключе­ нием того, что г. правой части будет величие веса:

т< = ^ ( аЛ r ^ V • • - t k 2 * > ) ;

ют

^ у г (к*аг ^ 2 кг + " ’ t k z \ ) i

' ' ****/») *

Нормальные уравнения:

L о J 21 f \

Г— ~LLr

1—

1—

7 IP .J

+ ал

=0 \

'

 

j +£ ^

= 0:

1+ us * 0.

LP -I

7

 

Нахождение ошибок функции урине иных элементов при уравнивании с условными уравнениями производится так же, как при уравнивании с уравнениями поправок.

о '17. Сглаживание экспериментальных зависимсстсЛ по методу наимень­

ших квалсатов

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

всего определяется характером экспериментальной кривой или теоретическими соображениями.

ІО?

Пусть имеется таблица экспериментальных данных и выбран общий вид функции

с . . . )

,

 

зависящий от нескольких параметров

а , é

, . . .

Эти параметры и надо подобрать по

способу

наимены к

кіадратов.

 

 

Согласие принципу наименьших квадратов меяго записать условие

Найдем значения

а

, 6

, с

.

обращающие левую

часть равенства

в т і п

,

для

чего возьмем частные про­

изводные по

а

, ё

, с

, . . . ,

приравняем их к нулю и в

результате

получим

систему уравнений

Решая эту систему, найдем неизвестные параметры. Пример. Пусть

» "Э о“ Х

д ё

'

Таким образом,

гг

ю з

L t r ( Q* V ^ ] = ö -

 

Решая данную систему, найдем а

и ё .

Для зависимых случайных величин,

как известно,коэф­

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

Представим случайную величину

У в виде

суммы не­

которой линейной функции от X

и остатиа

Z

У = * ( Х )+ і = ( А Х + б )+ 2 .

 

 

Будем рассматривать

 

 

 

2 = У - * ( Х )

 

 

 

как ошибку приблиаения случайной величины

У

линейной

функцией

 

 

 

* ( X ) = A X + ß .

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

1Л\%г] = т г п .

Запишем искомую линейную Функцию в виде

t ( X ) = A ( X - a ) + ë +C .

Тогда

і-С2- 2 А ( Х - л )(У- Ъ) - 2 С К ß - ë ) - 2 A C M (X - л ) ,

1 0 4

H O

» l ( X - a ) - M ( y - é ) = 0 -

Так как

а и - центры распределения соответствующих слу­ чайных величин,

М ( У - 6 ) ( Х - а ) = ?Х}у Зх ^

и тогда

м [ г г] . б ^ А гв * + с 2- 2 А ^ е . а % ■

 

?

G

2

:

Прибавим и ожииыеы

? .

?

М I*'2] =

 

 

 

^ )

#*«Ма м ш и ,

 

 

деетигает наименьшего зна­

чения яр«

 

 

 

 

A e ^ ^ , f 3 f

°

 

 

С = д

 

 

 

Таким образом, лкяейная функции наилучшего средне-

кмдратмеаиегѳ чриблкяеенил есть

S ( X ) . ? Xb. / £ ( X - * ) ' i ■

График этой .функции - так маммемея прямая среднекеадратичесяой регрессии У на X ■

Уравнение ярлыоК ередиеквадратической регрессіи можно записать и в виде

~&) -

Аналогично может бить зависала п уравнение средне- к^адратической регрессия X на У •

105

Графики этих уравнений показаны на р и с.15.

§18. Применение метода т м м т квадратов для обработки за­

висимых измерений

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

т .п .) . Вернемся еще раз

к

рассмотрению основных теоре­

тических положений метода

наименьших

квадратов.

Обозначим через

,

х г

, . . . ,

неизвестные. Так

же, как и раньше,

можно

записать

 

» Х2 т ” “>х п)~У-г~^

IOfi

Соседние файлы в папке книги из ГПНТБ