способ связи с многопараметрической адаптацией
Классы МПК: | H04B7/22 системы, основанные на распространении радиоволн со вторичным излучением при отражении (например тропосферное распространение радиоволн) |
Автор(ы): | Бузин Андрей Юрьевич (RU), Зарубин Андрей Юрьевич (RU), Беседин Александр Борисович (RU), Жуков Валентин Михайлович (RU), Игнатенко Михаил Владимирович (RU) |
Патентообладатель(и): | Тамбовское высшее военное авиационное инженерное училище радиоэлектроники (военный институт) (RU) |
Приоритеты: |
подача заявки:
2006-04-03 публикация патента:
27.04.2008 |
Изобретение относится к радиотехнике и может быть использовано в каналах радиосвязи диапазонов коротких (KB) и ультракоротких волн (УКВ) со случайными параметрами. Техническим результатом является повышение эффективности и максимальной относительной пропускной способности канала связи со случайными параметрами. Для этого в способе последовательность смены параметра адаптации выбирают с учетом рассчитанных вероятностных функций полезности всех возможных решений по управлению радиолинией и требований по достоверности приема информации. 2 ил.
Формула изобретения
Способ связи с многопараметрической адаптацией в ионосферном или тропосферном радиоканале связи, заключающийся в том, что с передающей стороны передают дискретною информацию, а на приемной стороне по принятой дискретной информации принимают решение о многопараметрическом адаптивном управлении процессом связи и посредством передатчика обратного канала передают в приемник обратного канала сигнал управления процессом связи, отличающийся тем, что для реализации алгоритма многопараметрической адаптации определяют параметры сигнально-помеховой обстановки, рассчитывают вероятностные ожидаемые функции полезности при принятии решений на адаптацию по одному из параметров от каждого решения в конкретных условиях связи и средние временные затраты на смену параметра, выбирают интервал То вероятности работы радиолинии без нарушения связи при различных параметрах с учетом интервала квазистационарности, в конце текущего временного интервала квазистационарности на основе рассчитанных функций полезности выбирают параметр адаптации по критерию Гурвица - критерий минимума среднеквадратического значения шумов для анализа качества канала связи, обеспечивающего максимизацию функции полезности на данном интервале квазистационарности канала радиосвязи
где [0, 1] - степень "оптимизма" принятия решения, Fij - функции полезности от принятия решений на адаптацию радиолинии, рассчитанные на основе анализа сигнальной и помеховой обстановки, при этом последовательность смены параметра адаптации выбирают с учетом рассчитанных вероятностных функций полезности всех возможных решений по управлению радиолинией и требований по достоверности приема информации.
Описание изобретения к патенту
Изобретение относится к радиотехнике и может быть использовано в каналах радиосвязи диапазонов коротких (KB) и ультракоротких волн (УКВ) со случайными параметрами.
Известны способы многопараметрической адаптивной радиосвязи в ионосферных и тропосферных радиоканалах связи и передачи дискретной информации, в которых повышение достоверности приема осуществляется путем предварительного зондирования радиоканала с последующей коррекцией характеристики канала связи, непрерывной оценки помеховой обстановки на резервных каналах связи и перестройки таких характеристик канала связи, как несущей частоты, мощности излучения, скорости передачи информационных символов, вида модуляции и способа кодирования.
Наиболее близким по технической сущности к заявляемому изобретению является изобретение SU №1585902 А1, кл. Н04В 7/00, 1990 г., в котором с целью повышения достоверности приема информации в системах радиосвязи KB и УКВ диапазонов по результатам оценки качества приема и в зависимости от помеховой обстановки в канале связи перестраивают следующие параметры: несущую частоту радиостанции, мощности излучения, скорость передачи информационных символов, а также вид модуляции и способ кодирования.
Существенным недостатком данного способа (а также и других известных способов многопараметрической адаптации, например, "Способ адаптивной связи" по AC H04B 7/22, SU 1758882 A1, 1992 г.) является жесткость (неизменяемость) алгоритма выбора параметра адаптации, что фактически заключается в переборе параметров адаптации, и соответственно, время адаптации может значительно превышать интервал квазистационарности канала связи. Это в реальных условиях априорной параметрической неопределенности канала связи может привести к значительному снижению пропускной способности канала.
Техническим результатом предлагаемого способа многопараметрической адаптации является повышение максимальной относительной пропускной способности канала связи (отношение времени передачи информации к времени сеанса связи на выделенном интервале квазистационарности канала связи) со случайными параметрами.
Сущность предлагаемого способа многопараметрической адаптации заключается в том, что в системе передачи дискретной информации, содержащей на передающей стороне приемник обратного канала, а на приемной стороне передатчик обратного канала, для реализации алгоритма многопараметрической адаптации определяют параметры сигнально-помеховой обстановки, рассчитывают вероятностные доходы при принятии решений на адаптацию по одному из параметров и средние временные затраты на смену параметра, в конце текущего временного интервала квазистационарности на основе рассчитанных доходов выбирают параметр адаптации по критерию Гурвица, максимизирующий пропускную способность на данном интервале квазистационарности канала радиосвязи, при этом последовательность смены параметра адаптации выбирают с учетом рассчитанных вероятностных доходов всех возможных решений по управлению радиолинией и требований по достоверности приема информации.
Способ связи с многопараметрической адаптацией по критерию Гурвица поясняется следующими расчетами.
В основе построения алгоритмов адаптивного управления процессом связи лежит модель исследуемого процесса. Представление процесса функционирования адаптивной радиолинии графоаналитической моделью позволяет создать ряд эвристических алгоритмов, обеспечивающих управление в строгом соответствии с установленным порядком применения правил выбора и смены. Однако неоптимальность таких алгоритмов не позволяет в полной мере выполнять требования по достоверной и своевременной передаче информации.
Так как управление состоит из целого ряда последовательно принимаемых решений, процесс выбора оптимального параметра адаптации можно описать с помощью теории принятия решений (Борзенко И.М. Адаптация, прогнозирование и выбор решений в алгоритмах управления технологическими объектами. - М.: Энергоатомиздат, 1984, - 144 с; Вентцель Е.С. Исследование операций: задачи, принципы, методология. - М.: Наука, 1988, - 208 с.).
Задачу принятия решений для начала процесса адаптации целесообразно представить в виде ориентированного графа, называемого деревом решений. При этом принимается, что множества решений, исходов и состояний внешнего фактора дискретны и конечны:
где di - возможное решение; sj - состояние внешнего фактора; х k - исход операции.
При детерминированном управлении, когда отсутствует влияние на исход операции каких-либо случайных воздействий, каждому решению соответствует единственный определенный исход операции, а следовательно, и вполне определена полезность этого решения. В такой ситуации принятие решений происходит в условиях полной определенности, а оптимальное решение принимается по наиболее предпочтительному исходу, соответствующему наибольшей полезности (наименьшим потерям). В задачах принятия решений в условиях определенности аргументом целевой функции является принятое решение и поиск оптимального решения сводится к отысканию экстремума целевой функции при варьировании решения d D.
При стохастическом управлении принимаемое решение не является единственным аргументом целевой функции, так как исход операции в данном случае зависит не только от него, но и от случайных факторов как внутренних, так и внешних. Принятие решений здесь происходит в условиях неопределенности, то есть исход принятия того или иного решения точно не определен.
Дерево решений, соответствующее задаче принятия решений в условиях неопределенности, приведено на фиг.1. Кружками на этом графе показаны узлы, в которых проявляется случайность, обусловленная внутренними и внешними случайными факторами. Из узлов А i выходят дуги, соответствующие различным состояниям внешних факторов при принятии решения di. Из узлов Вij исходят дуги, ведущие к исходам, которые могут быть получены при принятии решения di и si сохранение внешнего фактора. Для выбора решения, оптимального в смысле некоторого критерия, необходимо в каждом узле Аi определить значение целевой функции, соответствующей этому критерию.
Переход от узлов Аi к узлам Вij будет происходить с вероятностями p(sj/d i), причем:
Таким образом, каждый узел Аi характеризуется распределением вероятностей . Если состояние внешнего фактора не зависит от принимаемого решения, то p(sj/di )=p(sj) распределения вероятностей у всех узлов Аi одинаковы.
Если же решение di однозначно определяет состояние внешнего фактора, то неопределенность снимается и из узла А i исходит единственная дуга. Переход от узлов В ij к исходам xk происходит с вероятностями p(xk/di, s j), причем:
Таким образом, каждый узел Вij характеризуется распределением вероятностей . Эти распределения для каждого решения d i можно свести в стохастическую матрицу размера J×K, строками каждой являются распределения вероятностей, приведенные выше.
Условная вероятность событий sj , xk при решении di определяется:
p(sj, x k/di)=p(sj/d i)p(xk/di, sj).
Тогда вероятность исхода x k при решении di будет:
Важной особенностью теории принятия решений является необходимость введения на множестве исходов Х функции полезности (или потерь), которая каждому исходу х Х ставит в однозначное соответствие вещественное число, характеризующее положительный (или отрицательный) эффект данного исхода операции.
Принятие решений в условиях неопределенности основывается на различных критериях оптимальности, важнейшим из которых является байесовский. В соответствии с ним оптимальным является решение, максимизирующее ожидаемую полезность (доход):
где - ожидаемая полезность решения di при состоянии sj внешнего фактора.
Аналогично средний риск решения di будет:
,
где .
Таким образом, байесовское решение d i* находится из уравнения:
.
В соответствии с выбранной теорией принятия решения определяются множества состояний, решений и исходов для радиолинии с многопараметрической адаптацией.
Множество решений будет определяться количеством параметров, по которым возможна адаптация. Пусть это множество состоит из следующих элементов:
d 1 - адаптация по скорости передачи информации (по помехоустойчивому коду);
d2 - адаптация по мощности передающего устройства;
d3 - адаптация по рабочей частоте;
d4 - параметры оставить без изменения.
Состояние внешнего фактора зависит от стратегии радиоэлектронного воздействия противника на данную радиолинию:
s1 - противник ставит преднамеренные помехи;
s2 - преднамеренные помехи в радиолинии отсутствуют.
Воздействие противника и постоянно присутствующие в радиоканале случайные помехи при выборе любого решения приводят радиолинию к следующим исходам:
x1 - радиолиния находится в состоянии ведения связи;
x2 - радиолиния находится в состоянии восстановления или вхождения в связь.
С учетом приведенных выше выражений для выбранных решений d i, i=1, 4 имеем следующее.
Для решения d 1:
p(s1/d1 )+p(s2/d1)=1;
p(x1/d1, s 1)+p(x2/d1, s1)=1;
p(x1 /d1, s2)+p(x 2/d1, s2)=1;
p(x1/d1)=p(s 1/d1)p(x1/d 1, s1)+p(s2 /d1)p(x1/d 1, s2);
p(x 2/d1)=p(s1/d 1)p(x2/d1, s1)+p(s2/d 1)p(x2/d1, s2),
Для решения d 2:
p(s1/d2 )+p(s2/d2)=1;
p(x1/d2, s 1)+p(x2/d2, s1)=1;
p(x1 /d2, s2)+p(x 2/d2, s2)=1;
p(x1/d2)=p(s 1/d2)p(x1/d 2, s1)+p(s2 /d2)p(x1/d 2, s2);
p(x 2/d2)=p(s1/d 2)p(x2/d2, s1)+p(s2/d 2)p(x2/d2, s2),
Для решения d 3:
p(s1/d3 )+p(s2/d3)=1;
p(x1/d3, s 1)+p(x2/d3, s1)=1;
p(x1 /d3, s2)+p(x 2/d3, s2)=1;
p(x1/d3)=p(s 1/d3)p(x1/d 3, s1)+p(s2 /d3)p(x1/d 3, s2);
p(x 2/d3)=p(s1/d 3)p(x2/d3, s1)+p(s2/d 3)p(x2/d3, s2),
Для решения d 4:
p(s1/d4 )+p(s2/d4)=1;
p(x1/d4, s 1)+p(x2/d4, s1)=1;
p(x1 /d4, s2)+p(x 2/d4, s2)=1;
p(x1/d4)=p(s 1/d4)p(x1/d 4, s1)+p(s2 /d4)p(x1/d 4, s2);
p(x 2/d4)=p(s1/d 4)p(x2/d4, s1)+p(s2/d 4)p(x2/d4, s2),
Ожидаемые полезности решений d1-d4 соответственно будут определяться по формулам:
Тогда в соответствии с выбранным байесовским критерием оптимальное решение в выборе параметра адаптации будет определяться максимальным значением функции полезности работы автоматизированной радиолинии при выбранном параметре.
Таким образом, применение теории принятия решений позволяет решить задачу синтеза алгоритма многопараметрической адаптации в условиях сложной сигнально-помеховой обстановки.
Процедуру синтеза алгоритма адаптации можно проводить поэтапно. На начальном этапе необходимо задать множества решений, состояний внешнего фактора и вычислить ожидаемые полезности (доходы) от каждого решения в конкретных условиях ведения связи. В последующем формируется правило выбора решения, обеспечивающего достижения максимума дохода на интервале принятия решения.
Доход, получаемый радиолинией за каждое принятое решение, определяется выбранным критерием оптимизации. Как известно, основным требованием, предъявляемым к связи, как процессу передачи информации, является своевременность передачи сообщения с заданной достоверностью, которая является функцией пропускной способности линий связи. Для линий прерывистой связи, к которым относится радиолиния с многопараметрической адаптацией целесообразно в качестве дохода использовать относительную пропускную способность, определяемую отношением времени передачи сообщений с допустимой потерей достоверности к длительности сеанса связи. В соответствии с этим доход за одно решение можно оценивать как относительную пропускную способность радиолинии С0 на одном интервале принятия решения Т0. При этом для выполнения каждого решения di D требуются временные затраты, связанные с необходимостью перестройки радиолинии, обмена командами управления и синхронизации:
(d1)=2 к+ c,
(d2)=2 к,
(d3)= пер+2 к+ c,
(d4)=0,
где пер - время перестройки радиолинии, 2 к - временные затраты на обмен управляющей и квитирующей командами, c - временные затраты на синхронизацию кодеков.
Так как параметры каналов имеют вероятностный характер, то возможно снижение достоверности принимаемой информации ниже допустимого значения, то есть нарушение связи. Факт отказа рабочего канала устанавливается в течение случайного времени, затрачиваемого на анализ с помощью контрольно-решающего устройства (КРУ). Поскольку момент "срабатывания" КРУ случаен по времени, то его можно полагать равномерно распределенным на интервале Т0:
Т0 - (di)- кру,
где кру - время "срабатывания" КРУ.
Тогда с учетом приведенных выражений доходы решений di будут определяться:
,
,
где (n, k) - параметры помехоустойчивого кода.
Условные вероятности, определяющие состояние внешнего фактора и исход при выборе конкретного управления, имеют следующий физический смысл:
p(s1/di )=Pпп(di)P вк(di), p(s2 /di)=1-Pпп(d i)Рвк(di),
p(x1/di, s 1)=1-Рэп(di ), p(x2/di, s 1)=Pэп(di),
p(x1/di, s 2)=Pсв(di), p(x2/di, s 2)=1-Pэп(di ),
где Pпп(di ) - вероятность принятия противником решения на постановку помех;
Pвк(di) - вероятность временного контакта работы радиолинии с излучением преднамеренной помехи;
Рэп(di ) - вероятность энергетического подавления радиолинии преднамеренной помехой.
В условиях воздействия преднамеренных помех, то есть Ррп(di) 0 выражение примет вид:
,
где - время реакции комплекса радиоподавления.
Тогда ожидаемые доходы от решений будут определяться:
С учетом выражений для конкретных решений d 1-d4 имеем следующее.
Для решения d1 (изменить скорость передачи):
Для решения d2 (изменить мощность передающего устройства):
Для решения d3 (изменить рабочую частоту):
Для решения d4 (параметры оставить без изменения):
,
,
При расчете дохода необходимо определить аналитическую связь между мощностями сигнала на выходе передатчика и на входе приемника.
Мощность сигнала на входе приемника при его оптимальном согласовании с антенной и КПД фидера, примерно равным единице, равна:
,
где Еc - напряженность поля в точке приема; - длина волны (м); G2 - коэффициент усиления приемной антенны; К - поправочный коэффициент, равный 5765.
Напряженность поля в точке приема, выраженная в милливольтах на метр:
,
где rв - путь, проходимый волной от точки передачи до точки приема, км; R - модуль коэффициента отражения от Земли; n - число отражений от ионосферы; Р э - эквивалентная излученная мощность, кВт; Г и - полный интегральный коэффициент поглощения в ионосфере.
Эквивалентная излученная мощность Рэ определяется из следующих предположений:
- мощность, излученная антенной, делится пополам между обыкновенной и необыкновенной составляющими волны, но необыкновенная составляющая сильно поглощается и для приема оказывается полезной только половина излученной мощности;
- в средних и высоких широтах поляризация волны близка к круговой, а прием ведется на антенну с линейной поляризацией, поэтому используемая мощность уменьшается в 2 раза.
В результате принимается:
,
где P1 - мощность, подводимая к передающей антенне; G1 - коэффициент усиления передающей антенны.
Множитель (1+R)/2 учитывает уменьшение коэффициента усиления антенны за счет влияния реальной Земли, множитель Rn-1 - потери при отражении от Земли в случае многоскачкового распространения. При практических расчетах обычно принимают усредненное значение R=0,8. Тогда для односкачковой (n=1) трассы имеем:
.
Длину трассы по лучу rв можно определить по углу наклона траектории луча , для расчета которого необходимо знать критическую частоту для отражающего слоя F2 - f0F2 и f мпч:
,
где l - расстояние между передатчиком и приемником.
Значения f0F2 и fмпч для различных фаз солнечной активности, сезона, времени суток могут быть определены либо по ионосферным картам прогноза, либо по графикам суточного хода МПЧ.
В случае отражения от слоя F2 полный коэффициент поглощения равен
,
где A - суммарный коэффициент поглощения в слоях D, Е, F1; f - рабочая частота; f 1 - продольная составляющая гиромагнитной частоты в средних широтах; ВF2 - поглощение слоя F 2, от которого происходит отражение радиоволны. В средних широтах обычно принимают f1=0,7-0,8f. Ориентировочные данные о критической частоте f0F 2 можно определить A для различных фаз солнечной активности, сезона, времени суток, длины трассы. На частотах, обычно применяемых на среднеширотных радиолиниях значения hд для слоя F2 могут быть приняты следующими: зима, день - 250 км; зима, ночь - 350 км; лето, день - 400 км; лето, ночь - 250 км.
Обозначив через - затухание сигнала в радиоканале, получим мощность сигнала на входе приемника:
.
В реальных условиях ведения связи вероятностные характеристики, описывающие стратегию поведения противника (P пп(di), Рвк (di)), чаще всего неизвестны. В этом случае выбор оптимального параметра адаптации решается с использованием теории принятия решений в условиях неопределенности. При этом учет неопределенных факторов базируется на формировании специальных критериев, на основе которых принимаются решения (Вальд А. Последовательный анализ. - М.: Физматгиз, 1960).
Критерий Вальда (минимаксный) используется, когда лицо, принимающее решение, рассчитывает на реализацию худшего из возможных исходов.
Одна из модификаций минимаксного критерия - критерий Сэвиджа. Он применяется в случаях, когда худшие значения полезности или потерь отличаются незначительно, а лучшие значительно.
Критерий Гурвица является обобщением минимаксного критерия и отражает уровень "оптимизма" лица, принимающего решение, то есть степень его уверенности в получении наилучших результатов при принятии любого решения. Степень "оптимизма" характеризуется величиной [0, 1].
При =0 критерий Гурвица превращается в критерий Вальда, а при =1 - в критерий азартного игрока. В соответствии с критерием Гурвица выбирается решение, обеспечивающее при заданной функции полезности:
.
В этом случае ожидаемые доходы от решений будут определяться:
Для решения d1 (изменить скорость передачи):
.
Для решения d2 (изменить мощность передающего устройства):
.
Для решения d3 (изменить рабочую частоту):
.
Для решения d4 (параметры оставить без изменения):
,
,
.
Числовые значения параметров 1- 4, лежащие в интервале от 0 до 1, будут определяться техническими возможностями противника по обнаружению и идентификации работы радиолинии при данном параметре адаптации, а также текущим этапом боевой операции. В качестве средней точки зрения коэффициенты 1, 2, 4 можно принять равными 0.5, а 3 - 0.7.
В соответствии с байесовским критерием оптимальную стратегию выбора параметра адаптации можно представить в виде системы неравенств:
Решение данной системы неравенств путем сравнения доходов всех возможных решений позволяет определить управление, которое обеспечивает максимум пропускной способности на очередном интервале Т0.
Если требование по достоверности приема информации в одном из направлений дуплексной радиолинии не выполняется, то произойдет ее переход в режим восстановления связи. Этот переход осуществляется по команде КРУ, которое на установление факта нарушения связи затрачивает некоторое случайное время, составляющее наибольшую часть от всех суммарных временных потерь.
Верхней оценкой вероятности связи радиолинии с разработанным алгоритмом при достаточно малом Т 0 и Рэп 0 с достаточной для практики точностью может служить вероятность связи радиолинии, работающей на экстремальных частотах, откуда следует, что даже при относительно плохих условиях связи характеризуемых низкими значениями h2 0 и ограниченным частотным ресурсом вероятность срабатывания КРУ оказывается крайне низкой и не превышает сотые доли единиц. За счет этого и достигается высокая эффективность радиолинии в режиме ведения связи. Однако в условиях радиоподавления вероятность срабатывания КРУ существенно повышается.
Решающее правило в условиях восстановления связи можно синтезировать так же, как и для процесса ведения связи. Подобный оптимальный алгоритм позволит обеспечить перевод радиолинии в состояние ведения связи с наибольшей вероятностью и наименьшими временными затратами.
Множество решений в процессе восстановления связи будет определяться:
d1 - адаптация по скорости передачи информации (по помехоустойчивому коду);
d 2 - адаптация по мощности передающего устройства;
d3 - адаптация по рабочей частоте.
Отсутствие решения d4 объясняется тем, что оно не предполагает применение параметрической адаптации, и поэтому его принятие не обеспечит переход радиолинии в состояние ведения связи.
Решения d1 и d 2 будут элементами множества решений только в случае, если у них есть возможность увеличить выходную мощность или выбрать код с более высокой помехоустойчивостью.
Тогда в соответствии с байесовским критерием оптимальную стратегию выбора параметра адаптации можно представить в виде системы неравенств:
При расчете доходов коэффициенты 1, 2 можно принять равными 0.1-0.05, а 3 - 0.7.
Последовательность принятия решений можно представить алгоритмом многопараметрической адаптации (фиг.2).
Для реализации алгоритма многопараметрической адаптации по критерию Гурвица, максимизирующего пропускную способность, необходимо задать интервал дискретизации, определить средние временные затраты на смену параметров, средние времена реакции комплекса радиоподавления и вероятности работы радиолинии без нарушения связи на интервале Т0 при различных параметрах.
Выбор интервала Т0 должен производиться с учетом интервала квазистационарности процессов изменения во времени параметров сигналов и помех в ДКМ канале и времени, необходимого для оценки уровней сигналов и помех и принятия решения. Исходя из этого длительность интервала Т 0 может быть выбрана в пределах 30-60 сек. Временные затраты на смену параметров нетрудно определить по известным значениям длительности команд управления, времени перестройки и синхронизации для конкретного типа аппаратуры.
Работа радиолинии с описанным алгоритмом многопараметрической адаптации на длительности сеанса связи происходит следующим образом. На каждом интервале Т 0 производятся измерения мгновенных значений амплитуды принимаемого сигнала и определяются параметры распределения Накагами. Одновременно с этим оценивается соотношение сигнал-шум на входе приемника и рассчитываются значения Рсв и Рэп от решений d1 -d4.
При этом величина Р эп(di) определяется:
где h10=Р c/Pп - соотношение сигнал - преднамеренная помеха на входе приемника;
;
Р1п - мощность, подводимая к передающей антенне комплекса радиоподавления;
G 1п - коэффициент усиления передающей антенны комплекса радиоподавления.
В конце текущего интервала Т 0 на основе рассчитанных доходов выбирается оптимальное управление, обеспечивающее максимальную пропускную способность на следующем шаге.
Таким образом, в зависимости от выбранного параметра адаптации происходит заблаговременный (до срабатывания КРУ) переход без потерь времени на анализ отказавшего канала, при этом достаточно высокая интенсивность таких переходов позволяет обеспечить дополнительное повышение максимальной относительной пропускной способности радиолинии и ее помехозащищенности в условиях радиоподавления за счет увеличения времени на доразведку.
Класс H04B7/22 системы, основанные на распространении радиоволн со вторичным излучением при отражении (например тропосферное распространение радиоволн)