Введение В настоящее время для автоматического управления движением судна по курсу используются различные виды авторулевых [1, 2]. Несмотря на их различие, основным законом управления при автоматической стабилизации курса является пропорционально-интегрально-дифференциальный (ПИД). Основную информацию для формирования управляющих воздействий авторулевой получает от компаса и датчика угловой скорости. Авторулевые обеспечивают автоматическое ведение судна по заданному курсу в течение продолжительного времени, а также возможность изменения курса [2].Реализация алгоритма управления движением судна по заданной траектории требует использования, кроме курса и угловой скорости вращения корпуса [3], дополнительных навигационных параметров. Для оценки величины отклонения от заданной траектории могут быть использованы глобальные навигационные спутниковые системы (ГНСС).Цель данной статьи – изучить влияние погрешностей измерения координат судна с использованием ГНСС на качество процесса управления движением судна по заданной траектории. Экспериментальные исследования погрешностей измерения координат с использованием ГНСС В работах [4, 5] приведены результаты натурных экспериментов по исследованию характеристик точности позиционирования навигационных приемников при их работе с системами ГЛОНАСС, GPS, а также в совмещенном ГЛОНАСС/GPS режиме. Измерения проводились на открытой местности с дискретностью 1 с в течении 8 суток [4] и 7 суток [5]. В качестве иллюстрации на рис. 1 представлены результаты измерений погрешности позиционирования с использованием системы ГЛОНАСС, полученных авторами [4]. На рис. 1, а – ошибки местоопределения в системе ГЛОНАСС на плоскости (где   σr = 4,7 м – среднеквадратичная погрешность), на рис. 1, б показаны гистограммы распределения ошибок позиционирования соответственно по широте и долготе за 8 суток [4]. Аналогичные результаты (рис. 3, в), полученные за 7 суток, приведены в [5].        а)                                                                               б) в)Рис. 1. Ошибки местоопределения в системе ГЛОНАСС (по результатам [3, 4])Fig. 1. Location errors in the GLONASS system (based on the results of [3,4])Полученные результаты позволяют авторам [4, 5] сделать вывод о возможности аппроксимации погрешностей позиционирования навигационного приемника при моделировании работы навигационной системы случайным процессом с гауссовским законом распределения (следует отметить, что результат получен на интервале измерения в 7…8 суток).Однако при управлении движением судна по заданной траектории требуется генерации управляющих воздействий на движительно-рулевой комплекс с интервалами от нескольких до десятков секунд. Поэтому вопрос о возможности аппроксимации погрешностей позиционирования навигационного приемника случайным процессом с гауссовским законом распределения остается открытым.Для экспериментальной проверки был использован судовой приемник «ТРАНЗАС Т-701», который может работать в системах GPS и ГЛОНАСС, в совмещенном режиме GPS/ГЛОНАСС, а также в указанных режимах с дифференциальными поправками (ДП) [6]. Антенна приемника была установлены на крыше 10-этажного здания (в радиусе 500 м все постройки меньшей этажности), т.е. погрешность многолучевого приема отсутствовала. Измерения в разных системах проводились одновременно в течении одного часа с дискретностью в 1 с.Результаты статистической обработки измерений приведены в табл. 1.Таблица 1Результаты статистической обработки измеренийTable 1The results of statistical processing of measurementsПараметр\СистемаGPSГЛОНАССGPS + ГЛОНАССGPS + ГЛОНАСС + ДПGPS + ДПГЛОНАСС + ДПСреднее значениеШирота56º 19,63305’56º 19,63110’56º 19,63306’56º 19,63304’56º 19,63328’56º 19,63327’Долгота044º 01,00160’044º 01,00071’044º 01,00290’044º 01,00180’044º 01,00097’044º 00,99987’HDOP1,0201,8900,7500,9001,2201,030Среднеквадратическое отклонение, σ, мПо широте1,1491,9450,5921,0981,3160,957По долготе0,7365,1750,6201,8900,8212,243Как и ожидалось, наилучшие результаты получены при совместном использовании GPS и ГЛОНАСС (среднеквадратичное отклонение σ не превышает 0,620 м). При столь малых значениях σ дифференциальные поправки не дают снижения погрешности местоопределения. Эффект использования дифференциальных поправок заметен только для системы ГЛОНАСС (среднеквадратичная погрешность по каждой из координат снизилась почти вдвое).На рис. 2 (а, в, д, ж, и, л) приведены результаты измерений координат приемником «ТРАНЗАС Т‑701» в различных режимах (на каждом графике 3600 точек, снятых с интервалом 1 с). На осях графиков отложены отклонения от усредненной точки установки приемника (56º 19’ с.ш. и 44º 00’ в.д.) в долях минут.Для проверки гипотезы о законе распределения случайных величин измерений на рис. 2 (б, г, е, з, к, м) показаны распределения относительных частот зафиксированных результатов измерений для различных режимов работы навигационного приемника. а)                                                                           б) в)                                                                            г) д)                                                                            е) ж)                                                                          з) и)                                                                                 к) л)                                                                                  м)Рис. 2. Результаты измерения местоположения при работе приемника в различных режимах (а, в, д, ж, и, л) и интервальные распределения (б, г, е, з, к, м), вычисленные на основе измеренийFig. 2. Location measurement results during receiver operation in various modes (а, в, д, ж, и, л) and interval distributions (б, г, е, з, к, м) calculated based on measurementsНа графиках n = ni/n0 – нормированная относительная частота появления измеренной величины ni в заданном интервале, n0 = 3600 – объем выборки. Отклонения от усредненной точки установки приемника (56º 19’ с.ш. и 44º 00’ в.д.) для широты и долготы на осях приведены в долях минут.Вид кривых распределения, приведенных на рис. 2 (б, г, е, з, к, м), позволяет сделать вывод, что погрешности позиционирования навигационного приемника нельзя представить случайным процессом с гауссовским законом распределения, во всяком случае, на приведенной выборке в 3600 измерений в течении 1 часа.  Анализ характера погрешностей измерения координат с использованием ГНСС На рис. 4 а, б представлены две произвольные выборки из 200 точек (с интервалом в 1 с) из 3600 результатов измерений, полученных в режиме GPS+ ГЛОНАСС+ ДП (см. рис. 2, ж). Ошибка по долготе (ДР) меняется в пределах от -1 до 2,5 м, по широте (ШР) – в пределах -0,7 до 0,7 м.На эти графики нанесены изменения широты (ШН) и долготы (ДН), полученные с использованием генератора случайных чисел (ГСЧ) с гауссовским законом распределения. Распределение измеренных величин отклонений в выборке из 3600 измерений в течении часа приближается к виду нормального распределения (рис. 5, г). Характеристики генератора случайных чисел и реальные данные, полученные при измерениях в режиме GPS+ ГЛОНАСС+ ДП, показаны на рис. 5 (среднее и среднеквадратичное отклонения генератора случайных чисел совпадает с соответствующими характеристиками массива данных, полученного в режиме GPS+ ГЛОНАСС+ ДП). Здесь отклонения по широте и долготе пересчитаны в метры.В реальной выборке наблюдается сильная временная корреляция координат. В качестве иллюстрации на рис. 4, в представлены взаимокорреляционные функции трех 30-ти элементных выборок (номер 1, 2 и 3) по 30 с, сдвинутых относительно друг друга на 10 элементов (рис. 3), и исходной экспериментальной выборки из 3600 элементов, полученные при измерениях в режиме GPS+ГЛОНАСС+ДП (см. рис. 5 а, б). Рис. 3. Порядок вычисления взаимокорреляционных функцийFig. 3. The order of calculation of the intercorrelation functions а)                                                                                 б) в)                                                                                 г)Рис. 4. Сравнение реальных зависимостей погрешностей местоопределения и теоретических зависимостей при гауссовском законе распределения измеренных величинFig. 4. Comparison of real dependences of location errors and theoretical dependences under the Gaussian distribution law of measured quantitiesКак видно из рис. 4, в, коэффициент взаимной корреляции в основном превышает ±0,4 (сильная корреляция на малых интервалах времени, характерных для систем управления). На рис. 4, г показаны взаимокорреляционные функции трех 30-ти элементных выборок (номер 1, 2 и 3) по 30 с, сдвинутых относительно друг друга на 10 элементов, и выборки в 3600 элементов, полученной для генератора случайных чисел с гауссовским законом распределения. Коэффициент взаимной корреляции в этом случае не превышает ±0,3. а)                                                                   б) в)                                                                    г)Рис. 5. Сравнительные характеристики ГСЧ и реальных результатов при приёме в режиме ГЛОНАСС+GPS+ДП Fig. 5. Comparative characteristics of GPS and real results when receiving in GLONASS+GPS+DС mode Таким образом, при моделировании динамики судна с учетом погрешностей измерения навигационных приборов некорректно представлять погрешности позиционирования навигационного приемника случайным процессом с гауссовским законом распределения.Экспериментальные данные сняты для фиксированной географической точки. При движении судна координаты меняются и, строго говоря, в каждой новой точки траектории мы имеем свой закон распределения погрешностей определения местоположения. Однако, учитывая, что при решении задач управления рассматриваются весьма ограниченные перемещения судна (масштаб в пространстве десятки – сотни метров, масштаб во времени секунды – минуты), приемники ГНСС работают с одним созвездием спутников, при одинаковых параметрах тропосферы, ионосферы и других факторов, можно полагать, что распределения погрешностей определения местоположения будет одинаковы в разных точках траектории движения судна [7].Далее, при моделировании влияния погрешностей навигационных приемников на динамику судна, будут использованы экспериментальные значения погрешностей датчика позиционирования, приведенные во втором разделе. Влияние характера погрешностей местоопределения судна на динамику процесса управления Для исследования влияния характера погрешности определения местоположения на динамику судна воспользуемся математической моделью инновационного колесного судна «Золотое кольцо», приведенную в [8]. В [9, 10] подробно исследованы динамические характеристики судна при реализации алгоритма удержания судна на заданной траектории при ветровом воздействии. Для формирования управляющих воздействий на гребные колеса используется вспомогательная функция [11]:Rul=kα∙α-αz+kω∙ω+kγ∙y-yz, где α , αz  – текущий и заданный угол курса; ω  – угловая скорость судна; (y-yz)  – отклонение от заданной траектории движения; kα , kω , kγ  – коэффициенты пропорциональности.Моделирование проводилось при следующих условиях:– для измерения значений курса и угловой скорости используются идеальные датчики (погрешность измерения равна 0);– судно движется по траектории, совпадающей с осью х с максимальной скоростью;– боковой ветер (скорость 5 м/с) появляется при t = 300 с, при t = 600 с ветер исчезает.Результаты работы алгоритма удержания судна на прямолинейной траектории (при использовании «идеальных» датчиков измерения курса и угловой скорости, kα=10 , kω=40 , kγ=0,16 ) приведены на рис. 6 (здесь Δy – отклонение от заданной прямолинейной траектории). Использование навигационного датчика с систематической составляющей погрешности в 0,5 и 1,0 м (такого же порядка погрешности получены при экспериментальных исследованиях, табл. 1) не приводит к неприемлемому снижению качественных показателей процесса управления (кривые 2 и 3 на рис. 6). Система управления является статической по отношению к возмущающему воздействию (т.е.  при постоянном воздействии ошибка так же стремится к постоянному значению, зависящему от величины воздействия). Траектория смещается на величину погрешности измерения отклонения Δy, при этом угол дрейфа составляет 1,5 градуса.Результаты моделирования процесса удержания судна на заданной траектории при аппроксимации погрешностей позиционирования навигационного приемника случайным процессом с гауссовским законом распределения (см. рис. 5, в, г), представлены на рис. 7. При этом алгоритм управления не в состоянии обеспечить движение судна по заданной траектории. Для случая с гауссовским законом распределения легко реализовать алгоритм, обеспечивающий устойчивое движение судна с использованием текущего среднего значения отклонения на выборке определенной длины. В работе он не рассматривается, поскольку реальные навигационные датчики имеют другую картину распределения погрешностей позиционирования.  Рис. 6. Базовый вариант динамических характеристик суднаFig. 6. The basic version of the dynamic characteristics of the vessel Рис. 7. Динамика судна при аппроксимации погрешности датчика случайным процессом с гауссовским распределениемFig. 7. The dynamics of the vessel when approximating the error of the positioning sensor by a random process with a Gaussian distribution На рис. 8 показаны результаты работы алгоритма удержания судна на заданной траектории при использовании для моделирования погрешности датчика позиционирования экспериментальных данных (см. рис. 5, а, б) без ветрового воздействия, на рис. 9 – при ветровом воздействии (с 300 по 600 секунду).  Рис. 8. Динамика судна без ветрового воздействия при использовании для моделирования погрешности датчика позиционирования экспериментальных данныхFig. 8. The dynamics of the vessel without wind action when using experimental data to simulate the error of the positioning sensorРис. 9. Динамика судна при ветровом воздействии при использовании для моделирования погрешности датчика позиционирования экспериментальных данныхFig. 9. The dynamics of the vessel under wind influence when using experimental data for modeling the error of the positioning sensorВ результате решения системы дифференциальных уравнений на каждом шаге получается набор параметров, характеризующих поведение судна: скорость (продольная и поперечная), частоты вращения гребных колес, курс судна и угол дрейфа, угловая скорость вращения относительно центра масс, пройденный путь и отклонение от траектории. Каждую секунду к вычисленному значению отклонения от траектории добавляется очередное, экспериментально полученное значение погрешности навигационного приемника (из массива на рис. 5, а). В одном массиве 3600 точек, для моделирования использовались произвольно вырезанные массивы из 900 точек. Как показали результаты моделирования, выбор массива в 900 точек в начале, середины или в конце исходной последовательности не меняет качественные показатели процесса управления. Такой же результат дают и выборки из других массивов (см. рис. 2).Предложенный алгоритм обеспечивает удержание судна на заданной траектории с приемлемой точностью. Отклонения судна от заданной траектории не превышает ±2 м, что вполне приемлемо для судна длиной 83 м (см. рис. 8). При появлении ветра (см. рис. 9) траектория смещается на величину, пропорциональную внешнему воздействию (статическая система управления). При этом характер изменения поперечной координаты определяется характером изменений погрешности навигационного прибора.На рис. 10 приведены параметры процесса управления при изменении силы и направления ветра на 360º – скорость судна V, частоты вращения гребных колес n1 и n2, угол дрейфа αдр, отклонение от траектории Δy. Ось времени разбита на интервалы по 300 с, в каждом интервале меняются параметры внешнего воздействия. В качестве модели распределения погрешностей навигационного приемника использовалась экспериментальная зависимость (см. рис. 5, а).Судно «Золотое кольцо» в силу своих конструктивных особенностей имеет сложную зависимость крутящего момента, вызванного ветровым воздействием, от направления ветра [9]. Алгоритм обеспечивает удержание судна на заданной траектории в пределах от 3,5 м до 4 м при изменении направления ветра на 360º и среднеквадратичной погрешности измерения местоположения порядка 1 м. Это приемлемые показатели для судна длиной 83 м и шириной 14 м. 1,50  Рис. 10. Динамика судна при изменении скорости и направления ветра на 360º при использовании для моделирования погрешности датчика позиционирования экспериментальных данныхFig. 10. The dynamics of the vessel when the wind speed and direction change by 360º when using experimental data for modeling the error of the positioning sensor Заключение При реализации сложных алгоритмов управления судном (удержание судна на заданной траектории, переход судна на заданную циркуляцию, расхождение судов и пр.) требуются датчики контроля местоположения судна. В качестве таковых могут использоваться приемники ГНСС. При моделировании динамических процессов управления судном необходимо использовать адекватные модели характера поведения погрешностей позиционирования навигационных приемников.Проведенные исследования позволяют утверждать, что аппроксимировать погрешности позиционирования навигационного приемника случайным процессом с гауссовским законом распределения можно лишь на значительных интервалах времени (часы и более). На более коротких интервалах времени, характерных для функционирования судовых систем управления (секунды, минуты), модель нормального распределении результатов навигационных определений неприемлема. При моделировании работы алгоритмов управления судном целесообразно использовать экспериментальные зависимости погрешностей измерения местоположения.Как показывают результаты исследования, возможно реализовать алгоритмы, обеспечивающие удержание судна на заданной траектории с приемлемой точностью. Отклонения от траектории, вызванные погрешностями измерения координат навигационными приборами, соизмеримы с величиной отклонений, определяемых конструктивными особенностями корпуса судна и его движителей.