Писок рекомендованной литературы 53

ерчь, 2016 г.

Содержание

ведение 4

аздел 1 Математическая обработка результатов наблюдений 5

1.1 Расчеты количественных характеристик погрешностей навигационных

параметров 5

1.2 Расчеты погрешностей навигационных параметров

изменяющихся во времени и определение доверительных оценок 10

1.3 Обработка статистических данных методами линейной корреляции 14

Раздел 2. Способы определение места судна 18

2.1 Определение места судна (ОМС) методом изолиний 18

2.2 Графоаналитический метод 20

2.3 Влияние случайных ошибок измерений на точность места,

определяемого по двум ЛП 21 2.4 Графические методы отыскание вероятнейшего места судна

при избыточном числе линий положения 23

2.5 Аналитическое определение места судна и оценка точности 28

Задания для выполнения задач раздела 1 31

Задания для выполнения задач раздела 2 39

писок рекомендованной литературы 53

Введение

Методические указания по выполнению курсовой работы составлены для учебной дисциплины «Математические основы судовождения»,котораяпредназначена для уровня основной образовательной программы – специалист. Специальность - 26.05.05 «Судовождение». Профиль – «Промысловое судовождение». Статус дисциплины – базовая. Она включена в учебный план набора, начиная с 2015 года

Методические указания составлены в соответствии с рабочей программой отвечающей требованиям ФГОС ВО, компетенций МК ПДМНВ-78 с поправками и рабочего учебного плана с учетом требований ООП и рекомендаций модельных курсов ИМО 3.08, 7.03.

Выполнение курсовой работы направлено на закрепление знаний:

- фундаментальных понятий судовождения и правил использования Мореходных таблиц;

- правил и методов рациональных расчетов и вычислений с заданной точностью;

- теоретических основ морской картографии и построения картографических проекций;

- причин и видов погрешностей измерений;

- свойств и особенностей различных ошибок измерений навигационных параметров с учетом их влияния на безопасность мореплавания;

Теоретических основ получения и использования навигационных изолиний и линий положения (ЛП);

- влияния точности измерений параметров на точность изолиний или ЛП с характеристиками их точности;

- методов получения вероятнейшего места судна и оценки его точности.

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

Уметь:

- вести расчеты искомых величин различными методами с точностью не ниже точности исходных данных;

- выполнять измерения навигационных параметров и определять точность полученных результатов при равноточных и разновесных измерениях;

- на основе анализа точности измерений определять корреляционную связь между измеренными величинами;

- применять знания по ОМС и оценке его точности при решении всех специальных задач судовождения;

- вычислять значения навигационных параметров и их градиенты;

- выполнять расчеты для определения обсервованных координат места судна с оценкой их точности, в том числе и при использовании избыточной навигационной информации;

Мореходных таблиц и с применением вычислительной техники);

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

РАЗДЕЛ 1

МАТЕМАТИЧЕСКАЯ ОБРАБОТКА СТАТИСТИЧЕСКИХ ДАННЫХ.

1.1 Расчеты количественных характеристик

погрешностей навигационных параметров

1.1.1 Обработка результатов равноточных измерений навигационных параметров

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

Эти ошибки по своим свойствам и характеру можно разбить на три основные группы:

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

· Случайные ошибки это ошибки, которые вызваны многообразными и противоречивыми причинами, не поддающимися учёту и существенно не связанные с производством наблюдений, а их величина и знак для каждого измерения свои.

· Промахи определяются как неверные наблюдения или просчёты выходящие за пределы точности данного ряда измерений. Промахи из дальнейшей обработки исключаются.

Обработка наблюдений содержащих систематические ошибки будет рассмотрена в дальнейших работах.

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

1) среднее значение случайных ошибок близко к нулю;

2) вероятность появления ошибок равных по значению, но противоположных по знаку, одинакова;

3) численно небольшие ошибки более вероятны, чем численно большие;

4) случайные ошибки не могут превзойти некоторых границ, связанных с точностью производимых измерений;

5) внутри этих границ случайные ошибки могут принимать любые значения в соответствии с подпунктами 2 и 3.

Одной из важнейших характеристик случайных величин является вероятнейшее значение измеренной величины и оценка его точности.

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

некоторое количество n измерений навигационного параметра a.

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

В результате таких наблюдений мы получим n близких к истинному a_ист значений а₁, a₂...a_n, каждое из которых содержит некоторую ошибку. Так как a_ист нам не известно, то и точное значение этих ошибок нам так же неизвестно.

Очевидно, что наиболее близким т.е. вероятнейшим значением a_вер к истинному значению a_ист будет среднее арифметическое ряда равновесных измерений свободных от систематических ошибок:

( 1 )

Необходимо понимать, что формула (8) соответствует фактической мере точности измерений h. Поэтому, если, например, измерения выполняют линейкой разбитой на миллиметры, то ни увеличением числа измерений, ни увеличением числа знаков при вычислениях, ни колдовством нельзя получить результат с точностью, скажем, до 0,01 мм.

Мы уже говорили о том, что каждое из измерений содержало ошибку. Критерием точности единичного измерения на практике чаще всего применяют среднюю квадратичную погрешность m, вычисляемую по формуле

( 2 )

где, v_i уклонения отдельных измерений a_i от a_вер: v_i= a_вер – a_i, соответственно [vv] – сумма квадратов отклонений.

В случае, когда истинное значение параметра a_ист известно, СКП рассчитывается по формуле:

( 3 )

где, x_i ошибки отдельных измерений a_i: x_i= a_ист – a_i, соответственно [xx] – сумма квадратов ошибок.

Критерий точности единичного измерения СКП m даёт нам вероятность случайной ошибки x в границах от –m до m примерно 68.3%.

( 4 )

На практике так же часто применяется величина:

m_пред=3m

называемая предельной ошибкой одного измерения, накрывающая ошибки измерений с вероятностью 99.7%. Измерение у которых уклонения v превышают m_пред, объявляются промахами и из обработки исключаются.

Способ определения СКП одного (любого) измерения по формулам (9), (10) называется методом внутренней сходимости.

( 5 )

Другой способ, называемый методом размаха, позволяет находить СКП значительно проще, хотя и с меньшей точностью:

m = k_nR,

где R – размах, т.е. разность между максимальным и минимальным значением измеряемой величины

R=a_max – a_min

k_n – коэффициент, который зависит от числа измерений в серии, его значения приводятся в таблице 2.

Таблица 2.

n
kn	0.430	0.395	0.370	0.351	0.337	0.325	0.315

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

( 6 )

СКП вероятнейшего значения измеряемой величины m₀ рассчитывается по формуле:

( 7 )

m_0пред=3m₀

Из формулы (13) видно, что с увеличением числа измерений в n раз значение СКП вероятнейшего значения m₀ уменьшается лишь в раз, поэтому нет смысла для повышения точности вероятнейшего значения чрезмерно увеличивать число наблюдений. Следует повышать качество самих измерений. Оптимальное число измерений 7-9, хотя это и не всегда удаётся.

Разберём всё вышесказанное на примере.

Пример

Получена серия из 9 измерений секстаном горизонтального угла.

ОС 36°18,7′ 36°18,3′ 36°16,8′ 36°19,4′ 36°17,6′ 36°18,2′ 36°19,7′ 36°16,0′ 36°16,9′

Рассчитать:

· вероятнейшее значение горизонтального угла;

· СКП одного (любого) измерения двумя способами;

· предельную погрешность одного измерения;

· СКП вероятнейшего значения;

· предельную погрешность вероятнейшего значения

1.Составляем расчётную таблицу, в первой колонке которой серия измерений, во второй уклонения, в третьей квадраты уклонений.

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

OC	v	v v
36°	18.7′	0.7	0.49
	18.3	0.3	0.09
	16.8	-1.2	1.34
	19.4	1.4	1.96
	17.6	-0.4	0.16
	18.2	0.2	0.04
	19.7	1.7	2.89
	16.0	-2.0	4.00
	16.9	-1.1	1.21
OCср 36°18.0′	0,0	[vv]12.18

2. Рассчитываем СКП и предельную погрешность одного измерения:

методом внутренней сходимости по формуле (9):

m_пред = 3m = 3.6′

методом размаха:

ОС_max =36°19,7′

ОС_min =36°16,8′

R = 2,9′

k_n = 0.337

m = k_nR = ±1,1′

3. Рассчитываем СКП и предельную погрешность вероятнейшего значения по формулам (13), (14):

m_0пред= 3*0,4′ = ±1,2′

Контролем результата служит равенство нулю суммы отклонений v.

Все задания выполнены.

1.1.2 Обработка результатов неравноточных измерений навигационных параметров

Ходовой мостик современного судна оснащён различными по принципу действия техническими средствами. Этим обеспечивается возможность измерения одного и того же навигационного параметра с помощью нескольких приборов. Нередко один и тот же параметр измеряется хоть и одним прибором, но в различных условиях или различными людьми. Во всех таких случаях точность результатов будет неодинаковой, такие измерения называют неравноточными или неравновесными. Для сравнения неравноточных измерений пользуются величинами называемыми весами.

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

Вес p, взятый отдельно, без сопоставления с весами других измерений не несёт никакой информации. Он является сравнительной характеристикой качества наблюдений, и его следует рассматривать совместно с весами других наблюдений.

Принципы, по которым отдельным измерениям приписывают веса, могут быть различными. Например если результаты измерений a_i характеризуются своими СКП m_i, то этим измерениям присваивают веса p_i, обратно пропорциональные квадратам их СКП:

( 8 )

Веса, вычисленные по формуле (16) называют абсолютными.

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

( 9 )

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

хотя и нередко наоборот, что бы наибольший вес был равен единице (в этом случае все веса делят на наибольший вес)

( 10 )

.

( 11 )

В случае если каждое измерение есть результат серии равноточных измерений, то каждой серии могут приписывают веса p_i пропорциональные числу измерений в серии.

p_i= n_i,

( 12 )

и для взвешивания серии применяют величину относительного веса серии

хотя этот принцип менее обоснован, чем предыдущий.

Предположим, что в результате измерений некоторого навигационного параметра получен ряд значений а₁, a₂...a_n с весами p₁, p₂… p_n. За вероятнейшее значение ряда неравновесных измерений принимается весовое среднее или весовая арифметическая середина:

( 13 )

При неравноточных измерениях критерием точности служит СКП m, того измерения, вес которого принят за единицу. Эта погрешность называется СКП единицы веса. Вычисляется СКП единицы веса по формуле подобной (9) и является её обобщением:

( 14 )

В случае если a_ист известно:

( 15 )

СКП mp_i отдельного измерения a_i имеющего вес p_i, вычисляется по формуле:

( 16 )

В конце концов наибольший интерес представляет формула для СКП весового среднего a_вер, являющаяся венцом обработки ряда неравновесных измерений:

( 17 )

Предельная СКП mp_пред весового среднего вычисляется по формуле:

( 18 )

mp_пред = 3m_{0 p}

( 19 )

для случая большого числа измерений с надёжностью a=0.997;

и mp_пред = t m₀ p,

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

Пример 1

Дано: серия из шести пеленгов ИП измеренных с разной точностью и соответствующие им СКП.

ИП	315.6	315.4	316.2	316.1	316.2	315.6
mип	0.4°	0.7°	1.2°	0.8°	0.6°	0.9°

Найти:

1. Веса каждого измерения.

2. Вероятнейшее значение измеренного параметра.

3. СКП единицы веса.

4. СКП вероятнейшего значения и оценить его доверительной оценкой с надёжностью 0.99.

Решение:

1. Составляем расчётную таблицу:

· в первую колонку вносим серию измерений;

· во второй колонке рассчитываем веса измерений по формуле (18), за коэффициент k принимаем квадрат значения наибольшего СКП mип_max=1.2° (третье измерение) k=(1.2)²=1.44

и находим сумму весов [p]=21

· в третьей колонке рассчитываем произведения квадратов уклонений v² на соответствующие веса и находим сумму этих произведений [pvv]=1.93;

ИП	p	ИП*p	v	pv2
315.6	9.0	2840.4	-0.2	0.27
315.4	2.9	926.9	-0.4	0.41
316.2	1.0	316.2	0.4	0.19
316.1	2.3	711.2	0.3	0.25
316.2	4.0	1264.8	0.4	0.76
315.6	1.8	561.1	-0.2	0.05
	21.0	6620.6	[pvv] 1.93

· Рассчитываем вероятнейшее значение по формуле (20).

1. Рассчитываем СКП единицы веса по формуле (21):

1. Рассчитываем СКП и предельную погрешность вероятнейшего значения с заданной надёжностью по формулам (24) и (25):

при a = 0,99, по таблице 3 t=4.03

m_0пред= 4.03*0.01′ = ±0.04°

Истинное значение ИП_ист находится в интервале от 315,76° до 315,84°.

1.2 Расчеты погрешностей навигационных параметров

изменяющихся во времени и определение доверительных оценок

1.2.1 Обработка серии измерений навигационных параметров изменяющих своё значение с течением времени

Большинство измерений навигационных параметров судоводитель производит с движущегося судна в течение некоторого промежутка времени. За это время, судно успевает изменить своё положение относительно измеряемого объекта. То есть каждое последующее измерение, должно отличаться от предыдущего. Если характер изменения навигационного параметра близок к линейному, а промежуток времени между измерениями небольшой, то на вероятнейшем значении измеряемой величины это практически не скажется (значения большие и меньшие средней величины взаимокомпенсируются), а вот на значении СКП измеряемой величины это движение отразится в значительно большей степени.

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

Для приведения измеренных высот к одному моменту служит поправка Dh_T. Рассчитывается поправка следующим образом:

· в таблице 17 МТ-75 (таблица 15-б МТ-63) выбирается поправка Dh_T¹⁰ – изменение высоты светила за 10 секунд или в таблицах МТ-63 таблица 15-а поправка Dh_T¹ – изменение высоты светила за 1 минуту. Вход в таблицы по широте j и азимуту на светило А.

· если светило находится в восточной части горизонта - поправка положительна, в западной – отрицательна.

· Умножив Dh_T¹ на прошедшее время в минутах DT^м получим Dh_T

Dh_T=D T^м Dh_T¹

Для приведения измеренных высот к одному зениту служит поправка Dh_z. Рассчитывается поправка следующим образом:

· В таблице 16 МТ-75 (МТ-63) выбирается поправка Dh_z – изменение высоты светила за 1 минуту плавания. Вход в таблицу по скорости V и курсовому углу светила (КУ = А – ПУ). Если КУ<90° поправка положительна, КУ>90° отрицательна. Если нет дрейфа и течения, вместо ПУ берётся ИК.

· Умножив Dh_z¹ на прошедшее время в минутах DT^м получим Dh_z

Dh_z=D T^м Dh_z¹

Совместная поправка Dh = Dh_T + Dh_z или Dh = DT^м(Dh_T¹ + Dh_z¹).

При этом Dh_T¹ = 6*Dh_T¹⁰.

Работа с поправкой Dh_T¹⁰ возможна и в секундах. В этом случае, получаем временной коэффициент делением времени выраженного в секундах на десять. Поправка, выраженная в секундах, получается умножением Dh_T¹⁰ на полученный коэффициент:

1.2.2 Определение с заданной надёжностью доверительных оценок измеряемой величины и её СКП

Хотим подчеркнуть, что методика обработки измерений описанная в предыдущих параграфах, не совсем обоснована, т.к. формулы (9)-(14) выведены в предположении, что число измерений n очень велико. В частности при расчёте предельного значения мы априори брали коэффициент равный 3 и получали вероятность P=0.997, что верно только при n→ ∞. На практике же это число конечно и весьма ограничено. Поэтому возникает важная задача оценки точности величин a_вер, m, m₀ при конечном и ограниченном числе n.

Другими словами требуется установить вероятность P=a того, что при заданном конечном числе измерений n неизвестное значение a_ист заключено в границах a_вер – m_0пред до a_вер + m_0пред. Эти границы называют доверительными, а интервал a_вер ± m_0пред доверительным интервалом с заданной надёжностью a = P. Таким образом, нам известны n, a_вер и задано m_0пред. Найти: a.

Обратная задача формулируется следующим образом: каким при заданном числе измерений n должен быть доверительный интервал, для того, чтобы с заданной надёжностью a можно было утверждать, что неизвестная величина a_ист не выйдет за границы доверительного интервала:

a_вер – m_0пред < a_ист < a_вер + m_0пред

То есть, нам известны n, a_вер и задано a. Найти: m_0пред

Решение этой задачи в обоих постановках достигается при помощи распределения Стьюдента. Оно учитывает ограниченность числа измерений n и при условии n→ ∞ переходит в нормальное распределение. На его основе составлена Таблица 3. связывающая число наблюдений n, надёжность a и коэффициент , выносимый в формулу (11). В этом случае формула (11) приобретает вид: m_0пред = t m₀

Таблица 3. Значения коэффициента t по заданным n и a.

a n	0.5	0.7	0.8	0.9	0.95	0.99	0.999
	0.82	1.39	1.89	2.92	4.30	9.92	31.60
	0.76	1.25	1.64	2.35	3.18	5.84	12.92
	0.74	1.19	1.53	2.13	2.78	4.60	8.61
	0.73	1.16	1.48	2.02	2.57	4.03	6.87
	0.72	1.13	1.44	1.94	2.45	3.71	5.96
	0.71	1.12	1.41	1.89	2.36	3.50	5.41
	0.71	1.11	1.40	1.86	2.31	3.36	5.04
	0.70	1.10	1.38	1.83	2.26	3.25	4.78
	0.70	1.09	1.36	1.80	2.20	3.11	4.44
	0.69	1.08	1.35	1.77	2.16	3.01	4.22
	0.69	1.07	1.34	1.75	2.13	2.95	4.07
	0.69	1.07	1.33	1.73	2.09	2.86	3.88
	0.68	1.05	1.30	1.67	2.00	2.66	3.46
	0.68	1.04	1.29	1.66	1.98	2.62	3.37
∞	0.67	1.04	1.28	1.64	1.96	2.58	3.29

Примерно также решается вопрос и об оценке точности определения СКП единичного измерения, если число измерений конечно и ограничено.

На основе распределения Пирсона с заданной надёжностью a утверждается, что истинное значение СКП единичного измерения m_ист отличается от вычисленного не более чем на величину e.

На основе этого распределения составлена Таблица., связывающая число наблюдений n, надёжность a и коэффициент . То есть выбрав из таблицы 4 коэффициент t и рассчитав величину e=t∙m c надёжностью a можно утверждать, что m – e < m_ист< m + e.

Или наоборот, зная величину e можно определить надёжность a.

Таблица.4. Значения коэффициента t по заданным n и a.

a n	0.5	0.7	0.8	0.9	0.95	0.99	0.999
	0.25	0.40	0.56	0.95	1.37	2.67	5.64
	0.22	0.35	0.48	0.79	1.09	2.01	3.88
	0.20	0.32	0.41	0.65	0.92	1.62	2.98
	0.18	0.29	0.38	0.58	0.80	1.38	2.42
	0.17	0.27	0.35	0.52	0.71	1.20	2.06
	0.16	0.25	0.33	0.49	0.65	1.08	1.80
	0.15	0.23	0.30	0.43	0.55	0.90	1.45
	0.14	0.21	0.28	0.38	0.48	0.78	1.23
	0.13	0.19	0.26	0.35	0.44	0.70	1.07
	0.12	0.18	0.24	0.32	0.40	0.63	0.96
	0.11	0.17	0.22	0.29	0.37	0.58	0.88
	0.05	0.07	0.09	0.12	0.14	0.20	0.27

Пример 2

В широте j = 36.5°N судно следует ИК = 114° со скоростью V = 16 уз. Измерена серия 7 высот светила, азимут которого А=60,7° SO^st.

Моменты наблюдения T_i и отсчёты секстана ОС_i приведены в таблице.

Ti	10ч57м13с	10 57 58	10 58 48	10 59 35	11 00 27	11 01 10	11 02 02
ОСi	29°46,7′	55.5′	30°04,8′	12,0′	20,8′	29,0′	39,5′

Задание:

1. привести ОС_i к одному зениту и моменту времени.

Рассчитать:

2. вероятнейшее значение высоты светила;

3. СКП единичного измерения двумя способами;

4. предельную погрешность единичного измерения;

5. доверительный интервал накрывающий истинное значение СКП единичного значения с надёжностью (вероятностью) 0,90.

6. надёжность a определения СКП для доверительного интервала в 0,5 единицу измеряемого параметра (±0,5′)

7. СКП вероятнейшего значения;

8. предельную погрешность вероятнейшего значения измеренной величины и доверительный интервал накрывающий истинное значение измеряемой величины с надежностью (вероятностью) 0,95.

9. надёжность a для доверительного интервала в одну единицу измеряемого параметра (±1′).

Решение:

· Составляем расчётную таблицу, в первую колонку вносим моменты наблюдений, в пятую отсчёты секстанов.

Ti	DTi		DTi	DhT+Dhz	OСi	OCпр	vi	vi2
10ч	57м	13с	2м	22с	2.37м	25.5′	29°	46.7′	30°	12.2′	-0.2	0.04
					1.62	17.4		55.5		12.9	0.5	0.25
					0.78	8.4		4.8		13.2	0.8	0.64
					0.00	0.0		12.0		12.0	-0.4	0.16
				-52	-0.87	-9.3		20.8		11.5	-0.9	0.81
			-1	-35	-1.58	-17.1		29.0		11.9	-0.5	0.25
			-2	-27	-2.45	-26.4		39.5		13.1	0.7	0.49
									ОС=30°12.4′	0.0	Sv22.64

· Выбираем момент времени T₀ к которому будем приводить все измерения серии, чаще всего это или средний, или последний момент. Возьмём средний момент времени, в нашем случае четвёртый, т.е. T₀ = T₄.

· Рассчитываем промежутки времени DT_i между T₀ и текущим моментом T_i, DT_i = T₀ – T_i и вносим результаты во вторую колонку. В третью колонку внесём те же промежутки DT_i, но секунды выразим в десятых долях минуты.

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

из таблицы 17 МТ-75 по широте j и азимуту на светило А

Dh_T¹⁰ = +1,75′; Dh_T¹ = 6*1,75′=10,50′.

Из таблицы 16 МТ-75(63) по скорости V = 16 уз и КУ=119,3°– 114° =5,3°

Dh_z¹ = +0,27.

Рассчитываем совместную поправку за 1 минуту:

Dh¹ = Dh_z¹+ Dh_T¹=10,77

· Рассчитываем произведения Dh_i= DT_i Dh¹, результаты вносим в четвёртую колонку.

· В шестой колонке рассчитываем приведенные отсчёты секстана

ОС_{пр i}=Dh_i+OC_i

· Находим вероятнейшее значение ОС на T₀ (четвёртый момент времени) – внизу колонки среднее арифметическое от ОС_{пр I}.

· Рассчитываем СКП:

(методом внутренней сходимости)

В последние две колонки вносим уклонения v_i и квадраты уклонений v_i² соответственно. Внизу колонок находим сумму уклонений и сумму квадратов уклонений.

(методом размаха)

ОС_max =30°13.2′

ОС_min =30°11.5′

R = 1.7′

k_n = 0.370

m = k_nR = ±0.63′

Предельная погрешность единичного измерения m_пред = 3m = 2.0′

· Находим доверительный интервал единичного измерения,

по известному количеству наблюдений n = 7 и заданному в условии a = 0,90 из Таблица. находим значение коэффициента t = 0,65.

Рассчитываем e= t ∙m = 0.65∙0.66′ = 0.44′

Находим m-e=0,22 и m+e=1,1 следовательно 0,22≤ m_ист≤1.1.

· Рассчитываем надёжность a для заданного доверительного интервала m±e в примере e= 0.5′.

находим t=e/m=0.75

обратным входом в Таблица. по n=7 измерениям и по известному t=0.75 получаем надёжность a > 0.9.

· Рассчитываем СКП вероятнейшего значения:

· Находим предельную погрешность вероятнейшего значения и доверительный интервал:

по известному количеству наблюдений n = 7 и заданному в условии a = 0,95 из Таблицы 3 находим значение коэффициента t = 2.45.

m_0пред=tm₀=2,45∙0.25′ 0.6′ следовательно 30°11,8′≤ ОС_ист≤30°13,0′.

· Рассчитываем надёжность a для заданного m_0пред=±1′

рассчитываем t= m_0пред/m₀=4

для чего, обратным входом в Таблицу 3 по n=7 измерениям и по найденному t получаем надёжность a >0.99

1.3 Обработка статистических данных методами

линейной корреляции

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

Предположим, что какая-либо величина z является функцией некоторого количества различных навигационных параметров x₁, x₂,…x_n,

z=f(x₁, x₂,…x_n);

каждый из которых имеет среднюю квадратичную погрешность m_xi, тогда формула для расчёта m_z будет иметь вид:

( 20 )

В общем случае по формуле (15) рекомендуется следующий порядок расчётов:

1. Рассчитывают по исходной формуле значение определяемой величины.

2. Рассчитывают частные производные по переменным (измеренным) величинам.

3. Преобразовывают общую формулу СКП для конкретного случая.

4. В рабочую формулу подставляют исходные значения и рассчитывают СКП m_z.

Предположим, что мы произвели одновременные измерения двух серий величин x и y. Разумеется, эти измерения были подвержены влиянию случайных ошибок. Как уже говорилось выше, причины возникновения случайных ошибок многочисленны и практически не поддаются учёту. Среди них могут быть ошибки, выражающие некоторую неизвестную нам закономерность, или ошибки, одновременно влияющие на обе величины. То есть эти серии измерений уже нельзя считать независимыми друг от друга. С другой стороны, сказать, что они связаны жёсткой функциональной зависимостью вида y=ax+b, тоже нельзя. В нашем случае взаимосвязь является статистической, т.е. при изменении величины x величина y имеет тенденцию также изменяться.

Эта тенденция соблюдается лишь в большей или меньшей степени, а в каждом конкретном случае возможны отклонения от неё. Зависимость в таких случаях называется корреляционной и характеризуется коэффициентом корреляции r (-1 ≤ r ≤ 1). Если r>0, то большим значениям x в среднем соответствуют большие значения y и наоборот. Чаще всего эта тенденция носит линейный характер. Если r=±1, то зависимость становится функциональной. При r<|0.2| можно считать, что величины не коррелированны.

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

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

Для определения коэффициента корреляции необходимо измерить при более или менее идентичных условиях n парных значений двух навигационных параметров x_i и y_i (например, два пеленга береговых ориентиров). После определённой статистической обработки данных коэффициент корреляции рассчитывается по формуле:

( 21 )

,

где v_xi и v_yi уклонения величин x_i и y_i от вероятнейших значений соответственно.

Сама зависимость будет задаваться уравнением линейной регрессии:

( 22 )

,

где и вероятнейшие значения измеряемых величин.

Смысл уравнения линейной регрессии можно выразить так, каждому значению x_i с наибольшей вероятностью будет соответствовать значение y_xi.

Теперь обратимся к такому важному случаю, когда некоторая искомая величина a определяется, как функция непосредственно измеряемых величин x, y: a=f(x, y).

Предположим, что величины x, y коррелированны, т.е. если налицо взаимная зависимость случайных ошибок в аргументах x и y, в этом случае, СКП величины a определяется следующим образом:

( 23 )

В случае отсутствия корреляционной связи между величинами x и y, формула (30) примет вид (15).

Пример: Дано: две серии пеленгов ИП₁, ИП₂.

ИП1	30.8°	34.9°	36.5°	39.2°	38.9°	40.0°	45.8°	44.1°
ИП2	81.0°	78.8°	78.2°	85.0°	80.6°	84.4°	83.5°	87.7°

Найти:

1. вероятнейшие значения измеренных величин;

2. коэффициент корреляции между пеленгами;

3. составить уравнение линейной регрессии, построить график;

4. СКП каждого навигационного параметра;

5. рассчитать вероятнейшее значение горизонтального угла a;

6. рассчитать СКП вероятнейшего значения горизонтального угла a.

Решение:

1. Составляем расчётную таблицу:

1. во вторую и восьмую колонку вводим ИП2 и ИП1 соответственно и рассчитываем вероятнейшее значение (среднее арифметическое) каждой величины;

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

3. в четвёртой и шестой колонках квадраты уклонений;

4. в пятой их произведения;

5. внизу 4 – 6 колонок находим их суммы.

n	ИП2	vип2	v2ип2	vип2*vип1	v2ип1	vип1	ИП1
	81.0	1.4	2.0	11.2	63.6	8.0	30.8
	78.8	3.6	13.0	14.0	15.0	3.9	34.9
	78.2	4.2	17.6	9.6	5.2	2.3	36.5
	85.0	-2.6	6.8	1.1	0.2	-0.4	39.2
	80.6	1.8	3.2	-0.2	0.0	-0.1	38.9
	84.4	-2.0	4.0	2.4	1.5	-1.2	40.0
	83.5	-1.1	1.2	7.7	49.4	-7.0	45.8
	87.7	-5.3	28.1	28.2	28.4	-5.3	44.1
	82.4	S	75.86	73.95	163.2		38.8

2. По формуле (28) находим коэффициент корреляции:

3. Составляем уравнение линейной регрессии

ИП1
30.8	79.9
34.9	81.2
36.5	81.7
39.2	82.5
38.9	82.4
40.0	82.8
45.8	84.6
44.1	84.0

Строим зависимость от ИП₁

4. По формуле (9) рассчитываем СКП единичных значений

;

и формуле (13) СКП вероятнейших значений пеленгов.

5. Рассчитываем значение горизонтального угла a.

6. Рассчитываем СКП горизонтального угла a

6. без учёта корреляционной зависимости по формуле (15);

Находим частные производные пеленгов:

;

7. и с учётом корреляционной зависимости по формуле (29)

РАЗДЕЛ 2

СПОСОБЫ ОПРЕДЕЛЕНИЕ МЕСТА СУДНА

2.1 Определение места судна (ОМС) методом изолиний

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

Геометрические величины, получаемые из наблюдения внешних объектов для определения места судна, называют навигационными параметрами.

Совокупность точек, в которых навигационный параметр сохраняет своё значение постоянным, называется изолинией U навигационного параметра. Изолиния навигационного параметра является функцией координат:

U = f(j, l) = const.

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

Основными методами определения места судна являются – метод изолиний и обобщённый метод линий положения (ЛП).

Методом изолиний место судна можно определять, как графически, так и аналитически.

При помощи графического способа, место судна определяется непосредственно на карте. Этот способ применяется, когда изолинии имеют простую форму, удобную для нанесения на карту. Графический способ подробно рассматривается в курсе навигации, посему, его рассматривать не будем.

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

Задачу опрелеления места судна методом изолиний решают в следующем порядке:

· выполняются измерения навигационных параметров;

· составляются системы уравнений изиолиний измеренных параметров;

· подставляются в систему полученные параметры и система решается относительно j_o и l_o.

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

Если, изолиния показывает постоянное значение навигационного параметра, то градиент показывает направление и величину его изменения, разумеется, градиент всегда направлен перпендикулярно изолинии.

( 24 )

В общем виде, градиент равен изменению навигационного параметра на единицу расстояния, в судовождении милю:

,

то есть показывает, как изменяется навигационный параметр с расстоянием.

Один из простых, но показательных случаев, расчёт градиента глубин. Снимаем разность глубин, на рисунке – DU=10 метров, измеряем расстояние между изолиниями – Dn=0,4 мили, направление градиента идёт по увеличению глубины.

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

Таблица 7.

Навигационный параметр	Изолиния	Градиент	Направление градиента
Локсодромический пеленг с судна на ориентир
Локсодромический пеленг с ориентира на судно
Ортодромический пеленг с ориентира на судно		град/миля
Ортодромический пеленг с судна на ориентир		град/миля
Дистанция на плоскости		миля
Дистанция на сфере		миля
Горизонтальный угол		град/миля	Снимается графически, направлен из счислимой точки к центру описанной окружности.
Вертикальный угол
Высота светила		миля

В таблице:

j_с – счислимая широта

l_с – счислимая долгота

j_А – широта измеряемого объекта

l_А – долгота измеряемого объекта

D_А, D_{В –} дистанция до измеряемого объекта

Dc – счислимая дистанция

d – расстояние между ориентирами

Пс – счислимый пеленг

g - сферическое схождение меридианов

Для ортодромического пеленга с судна на ориентир

2.2 Графоаналитический метод

Этот метод основан, на обобщённом методе линий положения (ЛП). Суть способа состоит в том, что на коротком отрезке, изолинию сколь угодно сложной формы можно заменить, на прямую, касательную к этой изолинии. Впервые этот метод был применён для определения места судна по наблюдениям светил и проще его пояснить именно на этом примере.

Предположим, что наблюдатель находится в точке М с координатами j и l, которые ему не известны и которые предстоит определить с достаточной точностью. В этой точке он измеряет высоту h светила С.

Зная высоту светила можно провести круг равных высот с радиусом z = 90° - h. Разумеется, точка М будет находиться где-то на этой окружности.

При этом наблюдателю известны счислимые (приблизительные) координаты j_С и l_С точки. Из параллактического треугольника можно рассчитать счислимые высоту и азимут светила:

(25)

Эту высоту наблюдатель измерил бы, если б находился в точке М_С. Через эту точку так же можно провести круг равных высот с радиусом z_С = 90° - h_С. Разность n = h – h_C даст нам расстояние в милях между действительным и счислимым кругами равных высот. Проведя азимут А_С на светило и отложив на нём со своим знаком расстояние n, мы найдём определяющую точку К на действительном круге равных высот. Проведя через неё перпендикуляр, мы получим Высотную Линию Положения (ВЛП).

Измерив высоту другого светила и произведя аналогичные расчёты, мы получим вторую ВЛП. Пересечение обоих ВЛП даст нам обсервованное место судна М₀.

Учитывая то, что радиус круга равных высот, как правило, на несколько порядков больше расстояния между точками М_С и М, замена дуги на прямую линию практически не отразится на точности расчётов. То есть мы можем считать, что полученная нами точка М₀ практически совпадёт с действительной точкой М.

Переходя к обобщённому методу линий положения можно сказать, что навигационным параметром U, является истинная высота светила h, градиент g при этом равен единице, направление градиента совпадает с азимутом на светило.

Обобщённый порядок расчётов при графоаналитическом способе выглядит следующим образом:

1. Измеряются навигационные параметры U_о1 и U_о2;

2. рассчитываются счислимые параметры U_с1 и U_с2 (по формулам приведенным в таблице),на моменты измерений U_о1 и U_о2;

3. вычисляются разности

DU₁= U_о1 - U_с1 ; DU₂= U_о2 - U_с2;

4. рассчитываются модули g₁ и g₂ градиентов навигационных параметров и их направления t₁ и t₂;

5. по формулам

;

рассчитываются переносы;