Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
УДК 621.396.946
Повышение энергетической скрытности узкополосных сигналов систем спутниковой связи при трансионосферном распространении радиоволн через неоднородную ионосферу за счет понижения несущих частот и пространственно-разнесенного приема
Часть 1. Постановка задачи повышения энергетической скрытности
и ее решение
Ляхов А. В.
Актуальность статьи обусловлена тем, что решение задачи проектирования перспективных спутниковых систем связи (ССС) должно базироваться на комплексном подходе, который предполагает учет всех параметров и ограничений, обеспечивающих как требования по эффективности ее функционирования, включая требования, как по помехоустойчивости, так и по ее скрытности от средств радио- и радиотехнической разведки (РРТР). Известен способ повышения энергетической скрытности ССС при близком размещении приемника радиообнаружения за счет понижения несущей частоты и применения пространственно-разнесенного приема сигналов. В качестве количественного критерия достигаемого уровня энергетической скрытности ССС при этом используется такой общепризнанный показатель, как коэффициент энергетической скрытности (КЭС) ССС. Однако при этом известная методика определения КЭС является графоаналитической. Поэтому она не позволяет автоматизировать процесс определения этого показателя с помощью современных средств вычислительной техники, что затрудняет решение задачи проектирования перспективных ССС. Целью работы является анализ возможности повышения энергетической скрытности узкополосных сигналов систем спутниковой связи с пониженными несущими частотами на основе разработки методики, позволяющей сократить временные затраты и уменьшить погрешность оценки КЭС спутниковых систем связи, за счет применения вместо графических численных методов уточнения корней нелинейных уравнений обнаружения и помехоустойчивости и их реализации в MATLAB. Результаты и их новизна. Разработанная методика оценки КЭС спутниковых систем связи отличается от известных тем, что она позволяет сократить время оценки этого коэффициента с часов до единиц секунд и при этом уменьшить погрешность его оценки не менее чем на три порядка. Кроме того, впервые получены оригинальные аналитические выражения для вычисления первой и второй производной от Q-функции Маркума по параметру энергетическое отношение сигнал-шум, необходимые для программной реализации численных методов Ньютона и Ньюто-на-хорд уточнения корней нелинейного уравнения обнаружения. Практическая значимость заключается в том, что разработанная методика оценки КЭС спутниковых систем связи будет полезна техническим специалистам для обоснования исходных данные, необходимых для проведения расчетов по уточнению параметров перспективных спутниковых систем связи с учетом одновременного обеспечения ими требований, как по требуемой помехоустойчивости, так и по скрытности связи.
Библиографическая ссылка на статью:
Ляхов А. В. Повышение энергетической скрытности узкополосных сигналов систем спутниковой связи при трансионосферном распространении радиоволн через неоднородную ионосферу за счет понижения несущих частот и пространственно-разнесенного приема. Часть 1. Постановка задачи повышения энергетической скрытности и ее решение // Системы управления, связи и безопасности. 2022. № 4. С. 256-312. DOI: 10.24412/2410-9916-2022-4-256-312 Reference for citation:
Lyakhov A. V. Increasing the energy secrecy of narrow-band signals of satellite communication systems during transionospheric propagation of radio waves through an inhomogeneous ionosphere due to a decrease in carrier frequencies and space-diversified reception. Part 1. Statement of the problem of increasing energy secrecy and its solution. Systems of Control, Communication and Security, 2022, no. 4, pp. 256-312 (in Russian). DOI: 10.24412/2410-9916-2022-4-256-312
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
Кроме того, данная методика может представлять интерес научным работникам и соискателям, ведущим научные исследования в области спутниковой связи.
Ключевые слова: низкочастотная система спутниковой связи, райсовские замирания, трансионосферное распространение радиоволн, численный метод, коэффициент энергетической скрытности, спутниковая система связи, Q-функция Маркума.
Введение
Современные спутниковые системы связи (ССС) специального назначения являются системообразующим компонентом стратегического звена государственного управления. Именно эти ССС обеспечивают связь с удаленными пунктами органов государственной власти, дипломатическими представительствами Российской Федерации (РФ) в других странах, передачу команд боевого управления подразделениям силовых и специальных ведомств, дислоцированных за рубежом. При организации связи через ССС с вышеуказанными абонентами к ней предъявляются повышенные требования по безопасности, помехоустойчивости и скрытности.
Известно, что скрытность радиоэлектронного устройства - это его способность противостоять мерам радиоразведки, при этом, энергетическая скрытность характеризует способность противостоять мерам, направленным на обнаружение сигнала разведывательным приемным устройством [1].
И именно энергетическая скрытность ССС на линиях «вниз» «искусственный спутник Земли (ИСЗ) - земная станция спутниковой связи (ЗсСС)» при близком (на удалении не более 10 км) размещении приемника (ПРМ) радиообнаружения средства радио- и радиотехнической разведки (РРТР) стала объектом рассмотрения в настоящей статье.
Актуальность
Для обоснования актуальности, а также формальной постановки и решения задачи в работе введены обозначения, сведенные в таблицу 1 .
Таблица 1 - Обозначения
Обозначение Физический смысл обозначения
Р Р по з - заданная вероятность правильного обнаружения сигнала разведываемой системы связи
Р 1 лт з - заданная вероятность ложной тревоги при обнаружении сигнала разведываемой системы связи
hp2 - энергетическое ОСШ на входе ПРМ радиообнаружения средства РРТР
h 2 hp з - энергетическое ОСШ на входе ПРМ радиообнаружения, необходимое для достижения заданной вероятности правильного обнаружения Рпо з (например, Рпо з = 0,99) сигнала ССС при заданной вероятности ложной тревоги Рлт з (например, Р лт з = 10 "4)
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
Обозначение Физический смысл обозначения
R р - удаление ПРМ радиообнаружения средства РРТР от ПРМ ЗсСС
P 1 ош тр - требуемая вероятность ошибки Рош тр приема бита информации
h2 птр - ОСШ на входе ПРМ ЗсСС, определяемое из условия обеспечения требуемого качества связи (помехоустойчивости), оцениваемого таким показателем, как требуемая вероятность ошибки Рош тр приема бита информации (например, р _ р _ 10"5 ) 1 ош 1 ош тр 10 )
h 2 h тр 4 - ОСШ на входе ПРМ ЗсСС, определяемое из условия обеспечения требуемого качества связи (помехоустойчивости), при пространственно разнесенном приеме на четыре антенны
Y эс - коэффициент энергетической скрытности разведываемой системы связи
Уэс тр - требуемое значение коэффициента энергетической скрытности разведываемой системы связи
B - база сигнала разведываемой системы связи
Ni - электронная концентрация в /-той неоднородности нормальной ионосферы
N - среднее (фоновое) значение электронной концентрации в нормальной ионосфере
ANi - флуктуации электронной концентрации относительно среднего (фонового) значения N в /-той неоднородности нормальной ионосферы
Аф i - относительный фазовый сдвига луча, приходящего в точку приема, рассеянного /-той неоднородностью электронной концентрации нормальной ионосферы (Д^)
/о - несущая частота сигнала разведываемой системы связи
f0 гр - граничная пониженная несущая частота, на которой коэффициент энергетической скрытности ССС достигает требуемого значения у эс тр
рр - мощности регулярной (не случайной) составляющей принимаемого сигнала
Рфл - мощности флуктуационной (случайной) составляющей принимаемого сигнала
Y 2 - параметра Райса принимаемого сигнала, который определяется как отношение мощностей его регулярной и флуктуационной составляющей
( Y 2 )• v ' min - минимальное значение параметра Райса принимаемого сигнала
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
Обозначение Физический смысл обозначения
R - коэффициент пространственной корреляции быстрых замираний в элементах приемной антенной системы
AR - диапазон изменения коэффициента пространственной корреляции быстрых замираний в элементах приемной антенной системы
a - угол а трансионосферного распространения радиоволн относительно вертикали (зенитный угол)
Q(gb h01) - функция Маркума от аргументов gi и hoi , которые должны быть неотрицательными вещественными числами (gl, hoi ^ 0)
10(gl* ) - модифицированная функция Бесселя первого рода нулевого порядка
- среднеквадратическое отклонение относительного фазового сдвига лучей, приходящих в точку приема
( аф Lax - максимальное значение среднеквадратическое отклонение фазового фронта принимаемой волны
c - скорость света
ls - характерный размер ионосферных неоднородностей
h 'm - высота максимума ионизации в слое F ионосферы
N 1 v m - максимальная средняя электронная концентрация в слое F ионосферы на высоте максимума ионизации
¿э - эквивалентная толщина ионосферы со средней электронной концентрацией, равной ее максимальному значению
° AN - среднеквадратическое отклонение флуктуаций электронной концентрации в неоднородностях ионосферы
ви - интенсивность ионосферных неоднородностей, которая определяется как отношение среднеквадратическое отклонение флуктуаций электронной концентрации в неоднородно-стях ионосферы к максимальной средней электронной концентрации в слое F ионосферы на высоте максимума ионизации
î - вектор параметров ионосферы с элементами hэ, од^, ls или кэ, N m, в и, ls в зависимости от контекста
n - количество элементов приемной антенной системы
Ap - расстояние между элементами приемной антенной системы в случае их линейного эквидистантного разнесения
Ap к - интервал пространственной корреляции быстрых замираний принимаемых сигналов
^ош - вероятности ошибки приема бита информации
^k, ^ р - собственные значений корреляционной матрицы быстрых замираний в элементах приемной антенной системы {к,р = \,2,...,п)
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
Обозначение Физический смысл обозначения
X - вектор собственных значений корреляционной матрицы быстрых замираний в элементах приемной антенной системы
I - единичная диагональная матрица
h 2 - среднее значение ОСШ в ветвях пространственного разнесения приемных антенн
h 2 - в зависимости от контекста ОСШ на входе ПРМ радиообнаружения или среднее значение ОСШ в ветвях пространственного разнесения приемных антенн
Is a 1треб - требуемая процентная погрешность уточнения корня нелинейного уравнения обнаружения или помехоустойчивости
N - номер итерации итерационного процесса уточнения корня нелинейного уравнения обнаружения или помехоустойчивости
s 1° a\N - процентная погрешность уточнения корня нелинейного уравнения обнаружения или помехоустойчивости, достигаемая на итерации N
hi, hU - нижняя и верхняя граница отрезка локализации корня нелинейного уравнения обнаружения или помехоустойчивости
hrN - середина отрезка локализации корня нелинейного уравнения обнаружения или помехоустойчивости на итерации с номером N при использовании метода половинного деления (бисекций)
hU - итерационный процесс уточнения корня нелинейного уравнения обнаружения или помехоустойчивости численным методом половинного деления (бисекций)
h] - итерационный процесс уточнения корня нелинейного уравнения обнаружения или помехоустойчивости численным методом хорд
i 2 hN - итерационный процесс уточнения корня нелинейного уравнения обнаружения методом Ньютона
h 2 'чтив - итерационный процесс уточнения корня нелинейного уравнения обнаружения методом Ньютона-хорд
h0 - начальное приближение корня нелинейного уравнения обнаружения для численного метода Ньютона
f '( h2 ) - первая производная от нелинейного уравнения обнаружения по параметру ОСШ к
f '( h2 ) - вторая производная от нелинейного уравнения обнаружения по параметру ОСШ к
h 2 'чтив - итерационный процесс уточнения корня нелинейного уравнения обнаружения или помехоустойчивости методом хорд
M эс - абсолютная погрешность определения коэффициент энергетической скрытности ССС у эс
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
Обозначение Физический смысл обозначения
Мр з г\ - абсолютная погрешность определения параметра Нр з
длТр - абсолютная погрешность определения параметра й-2р
SY i эс - относительная погрешность определения коэффициент энергетической скрытности ССС у эс
$ M 2 MJ з 'Л - относительная погрешность определения параметра Нр з
$ m 2 Mip 'Л - относительная погрешность определения параметра Нтр
f $Ahp з J V р з y числ. - относительная погрешность определения параметра Нр з, найденная на основе численных методов уточнения корней нелинейного уравнения обнаружения
f $мТр J f тр У числ. - относительная погрешность определения параметра 8д/г 2 ,
найденная на основе численных методов уточнения корней нелинейного уравнения помехоустойчивости
( $Ahp2 з У f р з у гр. аналит. - относительная погрешность определения параметра Нр з, найденная графоаналитическим методом уточнения корней нелинейного уравнения обнаружения
($AhTp У f тр у гр. аналит. - относительная погрешность определения параметра 8. ,2 , ^лйтр найденная графоаналитическим методом уточнения корней нелинейного уравнения помехоустойчивости
П - число, показывающее во сколько раз найденная графоаналитическим методом относительная погрешность определения коэффициент энергетической скрытности ССС у эс превосходит аналогичную погрешность, но найденную численными методами решения нелинейных уравнений обнаружения и помехоустойчивости
Как известно, обнаружение сигнала происходит в условиях, когда на ПРМ радиообнаружения средства РРТР действуют помехи (шумы), и может сопровождаться ошибками двух видов: пропуск сигнала при его наличии на входе ПРМ радиообнаружения и ложное обнаружение (ложная тревога) при отсутствии сигнала. Эти ошибки носят вероятностный характер. Общепризнанной количественной мерой энергетической скрытности, от обнаружения сигнала разведываемой системы связи, является заданная вероятность правильного обнаружения Рпо з при заданных ограничениях на вероятность ложной тревоги Рлт з, которые в свою очередь зависят от энергетического отношения сигнал-шум (ОСШ) в рассматриваемой радиолинии и правила принятия решения на обнаружение сигнала [1-3].
Обеспечиваемая ПРМ радиообнаружения вероятность правильного обнаружения Рпо сигнала разведываемой системы связи при заданной (фиксиро-
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
• (1)
Рлт =const
ванной) вероятности ложной тревоги Рлт = const в общем случае связана некоторой функциональной зависимостью (Z), именуемой функцией обнаружения
[4], с энергетическим ОСШ кр на входе этого ПРМ (здесь и далее индекс «р» означает радиообнаружение):
Рпо = Z ( Рлт , кр )
В этом случае условием обеспечения энергетической скрытности системы
О 1
связи является выполнение неравенства к2 з > кр, т. е. снижение фактического
о
энергетического ОСШ на входе ПРМ радиообнаружения кр по сравнению со
значением кр з, необходимым для достижения заданной вероятности правильного обнаружения Рпо з (например, Рпо з = 0,99) сигнала ССС при заданной вероятности ложной тревоги Р лт з (например, Р лт з = 10 -4) [1-4]
кр з = Z (Рпо, Рлт )> кр , (2)
где Z 1 (Рпо, Рлт) - функция, обратная к функции Z (Рпо, к р).
г\ г\
Условие (2) кр з > кр принято выражать через коэффициент энергетической скрытности (КЭС) системы связи следующим образом [2]
Yэс = кр з/кр > 1 или (Уэс)дБ =(к) з)д£-(к)) > (3)
Из соотношения (3) следует, что КЭС системы связи показывает во сколько раз (на сколько дБ) значение энергетического ОСШ кр) з на входе ПРМ
радиообнаружения, необходимое для достижения заданного качества обнаружения, т. е. обеспечения заданной вероятности правильного обнаружения Рпо з сигнала системы связи при заданной допустимой вероятности ложной тревоги
-л
Рлт з, превосходит фактическое значение энергетического ОСШ кр) на входе этого ПРМ.
Примечание - С практической точки зрения, необходимо, чтобы КЭС системы связи был не просто больше единицы, а был бы не менее некоторого требуемого значения уэс тр. Например, если уэс тр = 20, то тогда должны выполняться неравенства уэс > уэстр = 20 или (уэс)дБ > (уэстр) Б = 13 дБ. Традиционно этого достигают за счет [1-3, 5]:
1) снижения уровня сигнала в ПРМ радиообнаружения при сохранении достаточного для нормальной работы уровня сигнала в ПРМ системы связи (амплитудный метод). Этот метод реализуется путем повышения чувствительности ПРМ системы связи и одновременного снижения мощности передатчика системы связи до минимально необходимого уровня (работа на пределе чувствительности). Признано, что такой метод практически нецелесообразен, поскольку при этом ухудшается помехоустойчивость разведываемой системы связи;
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
2) применения узконаправленных передающих антенн с малым уровнем боковых лепестков диаграммы направленности (ДН), применяемых в передатчике системы связи (пространственный метод).
Очевидно, что для ССС выполнение условия обеспечения помехоустойчивости при одновременном обеспечении ее энергетической скрытности, от обнаружения сигнала, сложнее всего реализовать в случае, когда ПРМ радиообнаружения расположен близко к ПРМ ЗсСС, т. к. в этом случае пространственный метод повышения энергетической скрытности оказывается неэффективным.
Актуальность такого допущения подтверждается высокой вероятностью ведения противоборствующей стороной (например, специальными службами недружественных стран и (или) незаконными террористическими формированиями) РРТР сигналов ССС, передаваемых на радиолинии «ИСЗ-ЗсСС», в случае если ЗсСС расположена за рубежом, а передаваемые сигналы содержат конфиденциальную информацию.
Ведение РРТР в таких случаях может производиться путем размещения ПРМ радиообнаружения непосредственно вблизи ПРМ ЗсСС, таким образом, чтобы они оба попадали в луч вещания ИСЗ. Например, в соседнем здании от дипломатического представительства РФ за рубежом или вблизи подразделений силовых и специальных ведомств РФ, дислоцированных за рубежом, обеспечивая при этом энергетическую и пространственную доступность средства РРТР к сигналам в радиолинии «ИСЗ-ЗС».
Поэтому, в дополнение к традиционным способам повышения энергетической скрытности ССС, был предложен способ повышения энергетической скрытности радиолинии «ИСЗ-ЗсСС», основанный на использовании в ССС широкополосных сигналов (ШПС), обладающих высокой помехоустойчивостью и скрытностью за счет большой базы (B >> 1) сигнала и энергетического уровня сигнала на уровне шумов на входе ПРМ радиообнаружения [1, 2, 5].
Однако при этом необходимо учитывать, что в УКВ диапазоне ССС специального назначения преимущественно используются не ШПС, обладающих высокой помехоустойчивостью и скрытностью за счет большой базы сигнала и энергетического уровня сигнала на уровне шумов, а узкополосные, т. е. простые сигналы (B = 1) с шириной спектра порядка 25 кГц. Примером такой ССС может служить отечественная узкополосная ССС «Гонец», которая осуществляет связь в УКВ диапазоне в интересах обслуживания органов государственной власти РФ, расположенных, в том числе, и за рубежом. Другим примером современной узкополосной военной ССС является система Mobile User Objective System (MUOS) [6].
В связи с этим, представляется актуальным известный способ обеспечения энергетической скрытности ССС при близком размещении ПРМ радиообнаружения от ПРМ ЗсСС (на удалении не более 10 км) и передаче в радиолинии «ИСЗ-ЗсСС» узкополосных, простых сигналов (B = 1) [7].
Этот способ учитывает известные [8] особенности распространения радиоволн (РРВ) диапазона частот от 30 до 100 МГц в ионосфере, заключающиеся в том, что в нормальной ионосфере электронная концентрация (ЭК)
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
N = N + АЫг- содержит флуктуации АЫг- относительно среднего (фонового) значения N (рис. 1).
Пониженные несущие частоты (ПНЧ): /о = 30...100 МГц
Геостационарный ИСЗ
В общем случае многолучевое РРВ приводит к райсовским замираниям принимаемых сигналов произвольной глубины
В частном случае, когда параметр Райса равен нулю, многолучевое РРВ приводит к релеевским
замираниям (РлЗ) принимаемых сигналов
В ПРМ ЗсСС используется прием сигналов на четыре антенны (п = 4)
Рис. 1. Модель размещения передатчика и приемника ЗсСС с ПНЧ при близком размещении приемника радиообнаружения средства РРТР и пространственно-разнесенном приеме на четыре антенны
Вследствие этого, при трансионосферном РРВ с пониженными несущими частотами (ПНЧ) до /о = 30...100 МГц через неоднородную ионосферу к ПРМ ЗсСС приходит множество (/ = 1 ...М) рассеянных неоднородностями ЭК (ДА^ ) лучей с относительными фазовыми сдвигами Аф г ~ АЫг //о . Значения фазовых сдвигов Аф I определяют условия возникновения интерференционных быстрых замираний райсовского (0 < Аф ^ < 2п) или рэлеевского (Аф г >> 2п) типа, причем райсовские замирания могут быть произвольной глубины. При традиционно используемых в ССС несущих частотах /о =1...10 ГГц относительные фазовые сдвиги лучей отсутствуют (т. к. Аф г ~ АЫг//о = 0), и поэтому интерференционные быстрые замирания сигналов не возникают.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Глубину райсовских замираний сигналов принято количественно оценивать при помощи параметра Райса у2, который определяется как отношение мощности регулярной (не случайной) составляющей принимаемого Рр сигнала
к мощности флуктуационной (случайной) Рфл составляющей у2 = Рр/Рфл [9]. Таким образом, параметр Райса может принимать произвольные значения от нуля (у2 = 0) при рэлеевских замираниях, когда мощность регулярной составляющей принимаемого сигнала равна нулю (Рр = 0) и вплоть до бесконечности
(у2 = да) при отсутствии замираний, когда мощность флуктуационной (случайной) составляющей принимаемого сигнала равна нулю (Рфл = 0), т. е.
0 < у2 < да.
Для того чтобы скомпенсировать отрицательные эффекты, возникающие при приеме сигналов с быстрыми замираниями в ЗсСС, в рамках известного способа повышения энергетической скрытности ССС [7], предлагается использовать пространственно-разнесённый приём сигналов на четыре антенны. В то время как в ПРМ радиообнаружения скомпенсировать отрицательные эффекты, возникающие при приеме сигналов с замираниями, невозможно, т. к. этот приемник, в связи с ограничениями на массогабаритные характеристики средства РРТР, вынужден использовать традиционные методы одиночного приема сигнала (т. е. на одну антенну).
При этом для ССС с ПНЧ (здесь термин «ССС с ПНЧ» указывает на то, что ССС работает в диапазоне, ниже минимального значения /о =150...200 МГц диапазонов частот, отведенных Регламентом радиосвязи МККР для ССС) формула (3) для вычисления КЭС у эс была в [7] представлена в виде:
У эс = Н р з/ Нтр или (У эс ) дБ =(Н р з ) -(Н ^ ) , (4)
где йтр - ОСШ на входе ПРМ ЗсСС, определяемое из условия обеспечения
требуемого качества связи (помехоустойчивости), оцениваемого таким показателем, как требуемая вероятность ошибки Рош тр приема бита информации
(например, РоШ = Рош тр = 10-5).
О О
Примечание - Так как параметры Нр з и Н^,, входящие в формулу (4),
представляют собой соответствующие энергетические ОСШ (альтернативное название ОСШ), то для упрощения далее в тексте термин «энергетические» будет упускаться.
Как видно, при переходе от формулы (3) к (4) сделано допущение о том,
-Л -Л
что Н2 = Н тр. Справедливость такого допущения обосновывается в [7] именно
пространственной близостью ПРМ радиообнаружения средства РРТР к ПРМ ЗсСС.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Полагая, что в ПРМ радиообнаружения используется схема оптимальной некогерентной обработки узкополосных, т. е. простых (В = 1) сигналов со случайной начальной фазой и амплитудой, распределенной по рэлеевскому закону,
параметр кр з, входящий в формулу (4) в качестве числителя дроби, определялся в [7] на основе анализа уравнения обнаружения [4]
р =Р (1+кр2 Г =К)^ (5)
1 по 1 лт V1 лт • К-'/
Уравнение (5) является нелинейным, поэтому найти его аналитическое решение в виде кр з =£-1 (Рпо з, Рлт з) невозможно. Поэтому в [7] параметр
г\
кр з определялся графическим методом:
1) строился график функции обнаружения (5) при фиксированном значении вероятности ложной тревоги Р лт з = 10 -4;
2) находилась точка пересечения этого графика с линией постоянного уровня вероятности правильного обнаружения Рпо з = 0,99, которая далее проецировалась на ось абсцисс. Таким образом, в [7] было найдено значение параметра кр; з =30 дБ.
В свою очередь, параметр к-2р, входящий в формулу (4) в качестве знаменателя дроби, определяется на основе анализа уравнения помехоустойчивости Рпо = £ (к2), связывающего между собой вероятность ошибки Рош приема бита
информации и ОСШ на входе ПРМ ЗсСС. Известно, что в случае рэлеевских быстрых замираний и пространственно-разнесенном приеме простых сигналов (например, типа частотной телеграфии) на п антенн (п > 2) в ПРМ ЗсСС, а также при отсутствии корреляции (коэффициент пространственной корреляции замираний Я = 0) этих замираний в элементах приемной антенной системы (АС), определяется соотношением [10]:
Рош =С21п_х1 ( к 2 ) П , (6)
где СШ = ш!/п!(п -т)! - число сочетаний из т по п.
'Л
Параметр ктр в [7] также определялся графическим методом:
1) строился график функции помехоустойчивости (6), при фиксированном значении количества антенн, используемом при пространственно-разнесенном приеме в ЗсСС, а именно при приеме на четыре антенны (п = 4);
2) находилась точка пересечения этого графика с линией постоянного уровня вероятности ошибки Рош тр приема бита информации Р)ш = р>ш тр =
= 10-5, которая далее проецировалась на ось абсцисс. В результате было найде-
-л
но значение параметра к-2р «16 дБ.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
После подстановки найденных значений Нр з = 30 дБ и Н-2р ~ 16 дБ в формулу (4) было получено значение КЭС ССС с ПНЧ у эс = 30 -16 = 14 дБ.
Это позволило сделать вывод о том, что, если в ССС используются ПНЧ (до /о =30... 100 МГц), обуславливающие рэлеевские замирания принимаемых сигналов, и пространственно-разнесенный прием на четыре антенны, то можно обеспечить достижение высокой энергетической скрытности ССС [7] на уровне КЭС равного 16 дБ. Причем такой эффект достигается при близком размещении ПРМ радиообнаружения к ПРМ ЗсСС и при передаче в радиолинии «ИСЗ-ЗсСС» простых сигналов.
Недостатками методики определения КЭС ССС с ПНЧ, использованной в [7], являются:
1) не учитывается возможность возникновения в канале радиообнаружения райсовских замираний принимаемых сигналов произвольной глубины
(0 < у2 < да) вследствие эффектов, сопровождающих трансионосферное РРВ в диапазоне несущих частот /о = 30.. .100 МГц, [11];
2) не учитывается наличия пространственной корреляции быстрых замираний сигналов в элементах приемной АС, используемых для приема в ЗсСС [10, 12];
3) использование графических методов определения параметров Н^р и
Нр з, что приводит к значительным затратам времени и методическим погрешностям вычисления КЭС по формуле (4), т. к. исходные данные для таких вычислений определяются достаточно грубо, «на глаз».
В настоящее время известен ряд решений задачи повышения энергетической скрытности ССС (например, представленные в работах [13-16]).
Анализ работ [13-16] показал, что все они посвящены анализу возможностей повышения энергетической скрытности ССС от радиоперехвата (РПХ) сигналов, а не от их обнаружения. Однако более актуальной является задача обеспечения энергетической скрытности ССС при близком размещении не ПРМ РПХ, а ПРМ обнаружения сигналов средства РРТР. Это обусловлено тем, что подсистема радиоперехвата сигналов средства РРТР задействуется только после обнаружения факта прихода сигнала на вход ПРМ радиоразведки этого средства. К тому же, традиционно энергетическая скрытность характеризует способность собственно противостоять именно мерам, направленным на обнаружение сигнала ПРМ радиоразведки средства РРТР, а не его перехват.
Предметом рассмотрения в настоящей статье является научно-методический аппарат (НМА) определения энергетической скрытности ССС с ПНЧ от обнаружения сигналов.
Целью работы является анализ возможности повышения энергетической скрытности узкополосных сигналов систем спутниковой связи с пониженными несущими частотами на основе разработки методики, позволяющей уменьшить временные затраты и повысить точность оценки КЭС ССС, за счет применения
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
вместо графических численных методов уточнения корней нелинейных уравнений обнаружения и помехоустойчивости и их реализации в МЛТЬЛВ.
Авторский материал по оценке КЭС ССС на основе численных методов, ввиду его большого объема, был разделен на две относительно независимые части. Первая часть, представленная в данной статье, посвящена постановке задачи повышения энергетической скрытности узкополосных сигналов систем спутниковой связи с пониженными несущими частотами и анализу возможных подходов к ее решению на основе численных методов, которые могут быть применены для оценки КЭС ССС, а также сравнительной оценке погрешности определения этого параметра на основе численных и графоаналитических методов. Вторая часть будет посвящена рассмотрению методики автоматизированной оценки КЭС ССС и её реализации в независимом приложении МЛТЬЛВ, а также результатам вычислительного эксперимента по оценке КЭС ССС и их анализу.
Материал данной статьи, ввиду своей объемности, был декомпозирован на ряд логических подразделов.
1. Постановка общей задачи повышения энергетической скрытности узкополосных сигналов систем спутниковой связи с пониженными несущими частотами.
2. Постановка и решение частной задачи определения значения отношения сигнал-шум, необходимого для обнаружения с заданным качеством сигнала спутниковой системы связи, на основе численных методов.
2.1. Постановка частной задачи и анализ возможных подходов к ее решению.
2.2. Численный метод половинного деления.
2.3. Численный метод хорд.
2.4. Численный метод Ньютона.
2.5. Численный метод Ньютона-хорд.
3. Постановка и решение частной задачи определения значения отношения сигнал-шум, необходимого для обеспечения требуемой помехоустойчивости при пространственно-разнесенном приеме сигналов земной станцией спутниковой связи, на основе численных методов.
3.1. Постановка частной задачи и анализ возможных подходов к ее решению.
3.2. Численный метод половинного деления.
3.3. Численный метод хорд.
4. Сравнительная оценка погрешности определения коэффициента энергетической скрытности спутниковой системы связи на основе численных и графоаналитических методов.
Результаты, представленные в статье, основаны на исследованиях автора, которые были опубликованы в журнальных статьях только частично [17].
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
1. Постановка общей задачи повышения энергетической скрытности узкополосных сигналов систем спутниковой связи с пониженными
несущими частотами
Дано: С геостационарного ИСЗ, излучается узкополосный простой сигнал с базой близкой к единице (В = 1) на несущей частоте, принадлежащей диапазону ПНЧ /о е30...100 МГц (рис. 1).
Угол а трансионосферного РРВ относительно вертикали (зенитный угол) составляет 1,1344640138 рад или 65°.
ПРМ радиообнаружения средства РРТР расположен на одной высоте при малом удалении Яр < 10 км от приемника ЗсСС.
Противоборствующая сторона обладает полным объемом информации о разведываемом сигнале (известны несущая частота /0 и ширина спектра Ар0 разведываемого сигнала), модели распространения радиоволн применительно к радиолинии «ИСЗ-ЗсСС» и актуальных параметрах ионосферы. Кроме того, предполагается, что средством РРТР поиск сигнала разведываемой ССС с ПНЧ по несущей частоте /0 и в пространстве (по направлению на ИСЗ) уже осуществлен.
В ПРМ радиообнаружения средства РРТР реализована схема оптимальной некогерентной (НК) обработки сигналов с флуктуирующей начальной фазой. Заданные значения вероятности правильного обнаружения Рпо сигнала
ССС и ложной тревоги Р лт равны: Рпо з =0,99, Рлт з =10-4. При этом вероятность правильного обнаружения сигнала ССС определяется уравнением обнаружения [18]:
Рпо=е
gl h2, у2 ), hoi (h2,y2, Рлт)], (7)
где: Q(gi,hoi) - функция Маркума от аргументов gi и hoi, которые должны
bi
быть неотрицательными вещественными числами (gi, hoi ^ 0); у2 - параметр
о
Райса; h2 - ОСШ на входе ПРМ радиообнаружения средства РРТР. Функция Маркума определяется формулой:
ö(gbhoi) = J *exp
¿01
2 2 ^ *2 + gi2
Io(gix)dx, (8)
где Io (gix) обозначает модифицированную функцию Бесселя первого рода нулевого порядка.
Аргументы gi и hoi функции Маркума в (8) при условии, что в ПРМ обнаружения используется схема оптимальной НК обработки сигналов, зависят от
9 9
параметров h , у2, Рлт следующим образом:
i/2
, (9)
g1( h 2, y2 )= 2h2 y2 (l + y 2 + h2 )
-1
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
hOl (h2, Y2,Рлт) =
2ln I Рлт-l )(l + Y2 )(l + Y2 + h2
12
(10)
Зависимость между параметром Райса у2 и среднеквадратическим отклонением (СКО) фазовых сдвигов (Дф ~ AN /f0) приходящих лучей оф в радиоканале с трансионосферным распространением радиоволн описывается выражением вида [12]:
Y2 = Рр/Рфл = [ехр(аф) -1]"1. (11)
В (11) параметр оф, рад, определяется соотношением:
Оф «(80,8^с)Vls'hseca(одм/fo), (12)
где: c = 3-108 - скорость света, м/с; ls я 400 - характерный размер ионосферных
неоднородностей, м; Нэ я 2,5 • 105 - эквивалентная толщина ионосферы, м, со средней ЭК, равной ее максимальному значению N(hm) = Nm =(2...4)4012, эл/м3; N(hm) = Nm - максимальная средняя ЭК в слое F ионосферы (на высоте h = hm максимума ионизации), эл/м3; a - угол трансионосферного распространения радиоволн относительно вертикали (зенитный угол), рад; f0 - несущая частота, Гц; о д^ - СКО флуктуаций ЭК в неоднородностях ионосферы, эл/м3, вычисляемое по формуле:
OAN = ßиNm = 10-3 х(2...4)-1012 =(2...4)409; (13)
В формуле (13) ри = oan / Nm = 10-3 - интенсивность ионосферных неоднородностей [12].
Требуемая вероятность ошибки приема бита информации в приемнике
ЗсСС
составляет />ш тр = 10 .
В приемнике ЗсСС используется пространственно-разнесенный некогерентный (НК) прием сигналов на несколько (n > 2) антенн. Для упрощения вычислений предполагается, что используется вариант линейного эквидистантного разнесения элементов приемной АС на расстояние Др = 400 м, причем количество таких элементов и е 2,3,..., 10.
Функциональная зависимость вероятности ошибки приема бита информации в ССС при пространственно-разнесенном НК приеме сигналов на несколько (n > 2) антенн с квадратичным сложением в условиях релеевских общих замираний сигналов определяется уравнением помехоустойчивости [10, 19]:
n ХП-1) n Хр + 1
ош ^ Хк + 2 П (Хк - Хр )(Хр + Хр + ХkХр ), ( )
р ф k
где Х к, Хр - собственные значения корреляционной матрицы KQ, определяемые из решения уравнения:
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
det (KQ - XI ) = 0. (15)
Здесь: X =(Xj,X2,---,Xn) - вектор собственных значений матрицы KQ; I - единичная диагональная матрица. В случае линейного эквидистантного разнесения элементов приемной АС на фиксированное расстояние Ар = const корреляционная матрица KQ является квадратной, размером n х n:
KQ = h
1 R R 4
R 1 R
Ra R 1
R(n-1)2 R(n-2)2 R(n-3)2
R(n-1)
2 Л
R R
(п-2У {п-ЪУ
1
(16)
— 7
где: h 2 - среднее значение ОСШ в ветвях пространственного разнесения антенн; R - коэффициент пространственной корреляции быстрых замираний в линейно разнесенных на эквидистантный интервал Ар элементах приемной АС, вычисляемый по формуле [12]:
R = expI-(Ар/Арк)2
exp
(Араф/ ls )2
= exp|-[Ар^(80,8л/сlsh3 (gа^/f 0)>lseca/1s
(17)
где Ар к = Оф - интервал пространственной корреляции быстрых замираний
принимаемых сигналов, м.
Требуется:
2
1. Определить значения ОСШ на входе ПРМ радиообнаружения кр з, не-
= 0,99
= 10 -4 )
обходимое для обнаружения с заданным качеством (Рпо з = 0,99 , р лт з сигнала ССС с на основе численных методов уточнения корня уравнения обнаружения (7) с относительной (процентной) погрешностью не более 10-4 %.
2. Определить значение ОСШ на входе ПРМ ЗсСС ктр, необходимое для
обеспечения требуемой помехоустойчивости (/^ш тр = Ю 5) в радиолинии
«ИСЗ-ЗсСС», на основе численных методов уточнения корня уравнения помехоустойчивости (14) с относительной (процентной) погрешностью не более 10-4 %.
3. По результатам решения задачи по пунктам 1 и 2, провести сравнительную оценку погрешности определения КЭС ССС на основе численных и графоаналитических методов.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
2. Постановка и решение частной задачи определения значения отношения сигнал-шум, необходимого для обнаружения c заданным качеством сигнала систем спутниковой связи с пониженной несущей частотой, на основе численных методов
2.1 Постановка частной задачи и анализ возможных подходов к ее
решению
Так как в соответствии с формулировкой общей задачи несущая частота сигнала /о может, принимать произвольное значение в интервале от 30 до
о
100 МГц, то для дальнейших расчетов примем /о = 108 Гц.
Дано: Рпо = Рпо з = 0,99, = Рлт з = 10"4, 13 « 400 м, Нэ я 2,5 • 105 м,
аду = 2 • 10 9 эл/м3, а = 650, /0 = 108 Гц.
Требуется найти все или некоторые корни уравнения (7) относительно ар-
2 2
гумента Н = Нр з - ОСШ на входе ПРМ радиообнаружения.
Решение. Анализ возможных подходов к решению задачи показывает, что по условиям задачи соотношением (7) задана непрерывная функция обнаружения Рпо = ( Н 2 ,у 2, Рлт | аналитической зависимости вероятности правильного
обнаружения сигналов (Рпо) от
ОСШ (Н 2) на входе ПРМ радиообнаружения средства РРТР, параметра Райса (у ) и вероятности ложной тревоги (Р1Т):
Рио = ^2,12(/оХа),Рлт) =
= 0{81^2,У2(/оД,а)],/7О1[/72,у2(/оД,а),Рлт]}, (18)
где: О(ё1,/го1) - функция Маркума, определяемая соотношением (8) от аргументов и йоъ вычисляемых по формулам (9) и (10); Х = (/гэ,аддг,15) или и~~ вектор параметров ионосферы.
При этом требуется определить такое значение ОСШ Нр з, при котором вероятность правильного обнаружения Рпо сигнала ССС в ПРМ обнаружения средства РРТР примет заданное значение Рпо з = 0,99, если Рлт з = 10 "4:
Ф2)|Н2 =Н2 = Рпо з . (19)
P
\h 2=крз
Очевидно, что уравнение (19) является нелинейным. Поскольку нет воз-
можности получить обратную функцию Ç 1 РПоз>^лтз>У2(/оД>а)
т. е. зна-
о
чение Нр з нельзя вычислить аналитически, поэтому для его определения необходимо воспользоваться численными методами. Наиболее близкими являются численные методы поиска всех или некоторых корней нелинейного уравнения с одним неизвестным [20-22] вида:
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
2 ,2 -Puos=0. (23)
h2 =hp з
/ (* ) = 0, (20) где / (х) - некоторая непрерывная функция. Действительно, если принять:
/(к2 ) = ф2 )-/по з, (21)
то уравнение (19) можно представить в виде:
/(к2)=0' (22)
или:
/(ь2) = Рио-Риоз = ^з^2,у2(/оХа)
Вычисление функции / (к2) по формуле (23) реализовано на языке
МЛТЬАВ в виде программного компонента (см. приложение ко второй части статьи).
К численным методам поиска корня кр з нелинейного уравнения (7),
представленного в виде (23) относятся метод половинного деления, метод хорд, метод Ньютона, комбинированный метод Ньютона-хорд и др. [20-22].
В рамках указанных методов задача поиска корней уравнения (23) относительно аргумента к распадается на несколько подзадач [20]. Во-первых, надо исследовать количество, характер и расположение корней (локализовать корни). Во-вторых, найти приближенные значения корней (получить начальное приближение). В-третьих, выбрать из них интересующие корни и вычислить их с требуемой точностью на основе численных методов.
Первая и вторая подзадачи решаются аналитическими и графическими методами.
В рассматриваемом случае в рамках первой подзадачи на основе аналитических соотношений (7) - (13), (23) и оговоренных ранее исходных данных
были рассчитан и построен график функции / (к2) (рис. 2) при изменении от-
2
ношения ОСШ к в интервале от 45 до 75.
На рис. 2 видно, что график исследуемой функции / (к2) пересекает ось
абсцисс один раз в промежутке [55,60]. На концах этого интервала функция / (к2) имеет разные знаки. Поэтому уравнение (23) имеет один действительный корень приближенно равный кр з « 57,5. В рамках третьей подзадачи выбираем интересующий нас корень кр з «57,5, локализованный на отрезке [55,60], и уточняем его с требуемой точностью на основе численных методов.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
x10"3
гхиремь - и 1 .и
45 50 55 60 65 70 75
2
Значение ОСШ h
Рис. 2. График функции f |h2 j при Рлт з = 10"4, ls » 400 м, кэ » 2,5 • 105 м,
aAN = 2 • 109 эл/м3, а = 65°, f0 = 108 Гц
При использовании численных методов для решения сформулированной задачи учтем ее особенности:
1) нужно найти всего один корень;
2) относительная (процентная) погрешность уточнения корня |еа\ должна
быть не более 10-4 %;
3) временные затраты на поиск корня не критичны, но точность и достоверность результата должны быть обеспечены.
Эффективность численных методов будем оценивать по таким критериям, как количество итераций, требуемых для уточнения корня с заданной точностью, фактически достигаемая относительная (процентная) погрешность и усредненные временные затраты уточнения корня.
Для этого рассматриваемые численные методы реализованы на языке MATLAB в виде программных компонентов (приведены в приложении ко второй части статьи). Для их тестирования и контрольных прогонов, а также проведения сравнительной оценки усредненных временных затрат уточнения корня использовался персональный компьютер следующей аппаратно-программной конфигурации: процессор Intel Core i5-4570 3,2 ГГц, объем оперативной памяти 16,00 Гбайт и 64-разрядная операционная система Windows 10.
2.2. Численный метод половинного деления
Постановка частной задачи: в п. 2.1 был отделен единственный корень
-л
уравнения (23) приближенно равный hp з « 57,5, который заключен в проме-
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
жутке [55,60]. Требуется уточнить этот корень методом половинного деления с процентной погрешность не более |еа| б = 104 %.
Решение. Итерационный процесс половинного деления (его называют также методом бисекций) заключается в следующем [20-22]. До начала итерационного процесса предполагается, что искомый корень локализован на отрезке
hi, hf
Начиная с отрезка
hi, hl
на концах которого функции / ( к2) име-
ет разные знаки, алгоритм вычисляет среднюю точку этого отрезка: к.2=( к2+к 2)/2.
(24)
Если / ( кГ | = 0, корень найден, и алгоритм завершает работу. В против-
hf, hl
либо на от-
ном случае он продолжает поиск корня либо на отрезке
2 2_
резке кг , ки - в зависимости от того, на какой из двух половин отрезка функция / ( к21 имеет разные знаки на концах нового отрезка изменения аргумента. Если требуется найти корень с относительной точностью |еа| б, %, то
деление отрезка поиска корня пополам продолжается до тех пор, пока на итерации N выполнится следующее неравенство:
'a\N
(hfN - hf( N-1))
/ h
rN
100 % < s,
11треб '
2 2 hi( N-1), hu( N-1)
(25)
на предыду-
где через -1) обозначена середина отрезка
щей итерации с номером N -1. Тогда середина к^ последнего отрезка на итерации N даст значение корня уравнения (23) с требуемой точ-
hlN, huN
ностью
I а1треб
к^ = ( кN + к2ш) /2. (26)
Ход итерационного процесса определения допустимого значения ОСШ методом половинного деления иллюстрирует таблица 2.
Из данных таблицы 2 видно, что итерационный процесс поиска корня уравнения (23) методом половинного деления завершается на восемнадцатой итерации (N=18), когда искомый корень принимает значение 2 2
кГ18 = кр з = 57,3887 с относительной процентной погрешностью вычисления
корня |ва|18 = 6,6471 х 10 5 %, удовлетворяющей неравенству (25).
2
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Таблица 2 - Пошаговое выполнение метода половинного деления
для уравнения ( v23)
N Крайняя слева точка, h2 Средняя точка, h2 Крайняя справа точка, h2 Значение ошибки, , % Значение функции f ( h} j
1 55 57,5000 60 - 4,9406x10-5
2 55 56,2500 57,5000 4,3478 4,9406x10-5
3 56,2500 56,8750 57,5000 2,2222 -5,2609x10-4
4 56,8750 57,1875 57,5000 1,0989 -2,3260x10-4
5 57,1875 57,3438 57,5000 0,5464 -9,0198x10-5
6 57,3438 57,4219 57,5000 0,2725 -2,0051x10-5
7 57,3438 57,3828 57,4219 0,1361 1,4763x10-5
В 57,3828 57,4023 57,4219 0,0681 -2,6224x10-6
9 57,3828 57,3926 57,4023 0,0340 6,0758x10-6
10 57,3828 57,3877 57,3926 0,0170 1,7281x10-6
11 57,3877 57,3901 57,3926 0,0085 -4,4681x10-7
12 57,3877 57,3889 57,3901 0,0043 6,4071x10-7
13 57,3877 57,3883 57,3889 0,0021 9,6974x10-B
14 57,3883 57,3886 57,3889 0,0011 -1,7491x10-7
15 57,3886 57,3888 57,3889 5,3177х10-4 -3,8968x10-B
16 57,3886 57,3887 57,3888 2,6588х10-4 2,9003x10-B
17 57,3887 57,3887 57,3888 1,3294х10-4 -4,9822x10-9
18 57,3887 57,3887 57,3887 6,6471х10-5 1,2010x10-B
Как следует из таблицы 2 скорость сходимости итерационного процесса поиска корня уравнения (23) невелика: за одну итерацию точность увеличивается примерно вдвое, т. е. уточнение трех цифр требует 10 итераций.
Зато точность результата уточнения корня гарантируется. Таким образом, приходим к выводу, что метод половинного деления определения искомого
значения ОСШ hр з прост и очень надежен, при этом он устойчив к ошибкам
округления.
К недостаткам рассмотренного метода можно отнести то, что до начала расчета надо найти отрезок, на котором исследуемая функция меняет знак. Метод половинного деления реализован на языке MATLAB в виде программного компонента bisect.m (см. приложение ко второй части статьи). Усредненные временные затраты уточнения корня составили 0,022971 с.
2.3. Численный метод хорд
Постановка частной задачи: в п. 2.1 был отделен единственный корень
-л
уравнения (23) приближенно равный h^ з « 57,5, который заключен в промежутке [55,60]. Требуется уточнить этот корень методом хорд с относительной процентной погрешностью не более |еа| б = 104 %•
Решение. Метод хорд, который так же называют метод ложного положения, или по латыни метод - regula falsi [20], стали использовать потому, что,
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
как было показано в п. 2.2, метод половинного деления сходится весьма медленно (за одну итерацию точность корня увеличивается примерно вдвое).
Как и алгоритм метода половинного деления, на каждой N итерации ме-
кт, киы
тод хорд получает некоторый отрезок
содержащий корень непре-
рывной функции / (к2), которая имеет значения противоположных знаков в
9 9
точках кN и к^. В отличие от метода деления пополам очередное приближе-
hN, huN
а как точка
ние на итерации N вычисляется не как средина отрезка пересечения оси абсцисс с прямой линией, проведенной через точки
(кн, / (к,N)) и (к2ш, / (к2ш)).
Формула для точки пересечения выглядит следующим образом:
7 2
hsN =
hNf (h
'\iN )- huNf ( hlN ) / f ( huN )- f ( hlN )
(27)
Из (27) следует три возможности:
1) если / (кN) и / (к^и) имеют различные знаки, то нуль принадлежит интервалу к2ы,к2ы ;
2) если / (к^п) и / (кип) имеют различные знаки, нуль лежит на интер-
вале
2 2 hsN, huN
lsN-
3) если / (к2м) = 0, значит, нулем является к'&
Так как требуется найти корень с точностью |еа| б, то итерационный
процесс поиска корня продолжается до тех пор, пока на итерации N не выполниться следующее неравенство:
' ~ '100% < Ы.. (28)
'a In
(hfN - hfN )/
h
sN
'а1треб'
Применительно к сформулированной задаче ход указанного итерационного процесса определения искомого значения ОСШ методом хорд при
= 10-4 % иллюстрирует таблица 3.
1треб
Как видно из данных таблицы 3 вычислительный процесс поиска корня
уравнения (23) методом хорд завершается на шестой итерации ^ = 6), когда
? ?
искомый корень принимает значение к^ = к^, з = 57,3887 с относительной процентной погрешностью вычисления корня |в тельно меньше |е „I = 10-4 %.
I а 1треб
al 6
9,9086 х 10-6 %, что значи-
s
a
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Таблица 3 - Пошаговое выполнение метода хорд для уравнения (23)
N Крайняя слева точка Ит Точка пересечения, И^ Крайняя справа точка, 7 2 Значение ошибки, еа\, % Значение функции, f (hSN )
1 55 57,5861 60 - 8,737117х10"5
2 55 57,4035 57,5861 0,3181 6,568805х10"6
3 55 57,3898 57,4035 0,0238 4,908097х10"7
4 55 57,3888 57,3898 0,0018 3,665542х10"8
5 55 57,3887 57,3888 1,3268х10-4 2,737462х10"9
6 55 57,3887 57,3887 9,9086х10-6 2,044358х10"10
По данным таблиц 2 и 3 на рис. 3 приведены графики динамики сходимости итерационного процесса метода половинного деления и метода хорд к корню уравнения (23).
Проанализировав рис. 3, можно сделать вывод, что сходимость итераций методов половинного деления и хорд не монотонна. Скорость сходимости итерационного процесса половинного деления замедляется от итерации к итерации. Что касается метода хорд, то он сходится со скоростью геометрической прогрессии (при этом говорят, что он имеет линейную скорость сходимости) [20, 21].
Сравнивая между собой данные таблиц 2 и 3, приходим к выводу, что корни уравнения (23), уточненные методом половинного деления и методом
хорд, совпали: Ир з = 57,3887. Используя метод хорд, удалось уточнить этот корень за шесть итераций (N = 6), что в три раза меньше чем 18 итераций, потребовавшихся методу половинного деления. При этом относительная процентная
погрешность |е а|6 = 9,9086 х10-6 % корня, найденного методом хорд, на порядок меньше, чем методом половинного деления |е а|18 = 6,6471 х 10-5 %.
К недостаткам рассмотренного метода, как и у метода половинного деления, можно отнести то, что до начала расчета надо найти отрезок, на котором исследуемая функция меняет знак. Кроме того, следует учитывать, что в знаменателе формулы (27) для точки пересечения стоит разность значений функции.
Вдали от корня это несущественно, но вблизи корня значения функции /1И |
малы и очень близки (таблица 3). При этом может возникнуть потеря значащих цифр, приводящая к «разболтке» счета (поочередному приближению и удалению от корня) [20, 21]. Например, такое может произойти в случае, если функция / ^ И2 ^ в области пересечения с осью абсцисс достаточно пологая. В этом
случае очередная хорда может практически лечь на ось абсцисс, то есть полностью попасть в полосу «шумов». В этой ситуации произойдет сильное увеличение ошибки вычислений, так как уже отмечалось ранее, в (27) разность двух
близких величин /1И/1ИN | стоит в знаменателе. В этом смысле метод
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
половинного деления значительно устойчивее. Это ограничивает точность, с которой можно уточнить искомый корень, но для простых корней это ограничение невелико. Метод хорд реализован на языке МЛ^ЛВ в виде программного компонента regula_falsi.m (см. приложение ко второй части статьи). Усредненные временные затраты уточнения корня составили 0,013206 с, что в 1,74 раза меньше, чем у метода половинного деления 0,022971 с.
Номер итерации n
Рис. 3. Динамика сходимости итерационного процесса методов половинного деления (И2 ) и хорд (И5 ) к корню уравнения (23)
2.4. Численный метод Ньютона
Постановка частной задачи: в п. 2.1 был отделен единственный корень
о
уравнения (23) приближенно равный Ир з « 57,5, который заключен в промежутке [55,60]. Требуется уточнить этот корень методом Ньютона с относительной процентной погрешностью не более |еа| б = 10 4 %.
Решение. Известно [20-22], что метод Ньютона один из наиболее распространенных методов уточнения корней нелинейных уравнений. Его называют также методом Ньютона-Рафсона или методом касательных.
Если И02 начальное приближение корня уравнения (23), то последовательные приближения методом Ньютона применительно к рассматриваемой задаче находят по формуле:
' " 4 ' ' " 4 (29)
h2N=hN _i - f ( hN _i )/ f '( hN _i )
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
где: N - номер итерации (Я = 1,2, ...); /( ИN-11 - значение функции /(И21 при И2 = И_2"-1; /'(И^-11 = ё/|И2^ё|И2- значение первой производной
функции /1И21 по параметру ОСШ И2 при И2 = И\
.2 N •
о
В методе Ньютона начальное приближение корня И0 уравнения (23) вы-
бирают на интервале
ht, К
руководствуясь следующим правилом: в каче-
' 1
стве исходной точки И02 выбирается тот конец интервала
h/, ht
, которому
отвечает ордината того же знака, что и знак второй производной / И 21 [20]:
h2 И
h/2, если f ( h
f ( h2 ) f '( h2 )> 0;
h 3, если f ( h2 hf ) > 0.
(30)
Из (29) следует, что для реализации метода Ньютона требуется существование первой непрерывной производной функции /(И21, т. е. должно выполняться условие: /И^-10. Чтобы не было потери точности, отношение
/1ИN -1/ИN-11 надо вычислять достаточно аккуратно. К остальным погрешностям расчета метод Ньютона хорошо устойчив [18].
о
Так как для вычисления очередного приближения И^ по формуле (29)
нужно знать только одно предыдущее значение И N-1, то алгоритм поиска корня методом Ньютона называют одношаговыми. Геометрически итерационный процесс (29) означает замену на каждой итерации графика функции /(И21 касательной к нему [20-22].
Основным фактором, ограничивающим прикладное использование метода Ньютона, является необходимость вычислять первую и вторую производную функции, что не всегда удобно. Так, применительно к рассматриваемой задаче дифференцированию подлежит функция (23), которую учетом соотношений (7) - (10), можно представить в виде:
2И2 у2 (1 + у 2 + И21 1
' | = | х ехр <
2ln( ^T_1)(l+Y2 )(l+Y2 +h2 )_1
У x
X Ir
^2h 2 y2 (l + y2 + h 2 ) 1 x
dx - P
по з •
(31)
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Из (31) следует, что функция, подлежащая дифференцированию, задана
интегралом, содержащим параметр ОСШ И . Причем не только подынтегральное выражение в формуле (31) содержит это параметр, но и нижний предел интеграла зависит от него. Таким образом, вычисление даже первой производной
/ И 21 = /1И 2 ёИ 2 предполагает дифференцирование под знаком интеграла в (31), что в рассматриваемом случае представляет достаточно сложную задачу, не говоря уже о второй производной /И 21 = ё 2/(И 2 |Д ёИ 21 , не только на
этапе получения аналитического выражения для производной, но и на этапе реализации численного алгоритма Ньютона средствами вычислительной техники. Известный НМА оценки скрытности ССС с ПНЧ не содержит известных аналитических решений для вычисления первой / И 21 и второй / И 21 производной от вероятности правильного обнаружения сигналов (Рпо), выражаемой че-
2
рез функцию Маркума, по параметру ОСШ (И ).
Известно, что дифференцирование под знаком интеграла, содержащим параметр, упрощается, если интеграл содержит параметр только в подынтегральном выражении [23-25]. В связи с этим воспользуемся альтернативным представлением функции Маркума, приведенным в [26]:
Q( а,р ) = 1 -1 j oßcosx + ß
Па 2 +2aßcosx + ß2
1 - exp
f а2 + 2aßcosT + ß2A
dx
(32)
Для того чтобы записать выражение для функции, подлежащей дифференцированию согласно методу Ньютона, учтем, что каждый из аргументов а и
в функции Маркума (8) зависит от параметра ОСШ И на основании соотношений (3) и (4) следующим образом:
а ( И2 ) = В! ( И2, у2 ) = 2И2 у2 (1 + у2 + И2 ^
12
в (h2hoi (h2 ,у2, Рлт )= 2ln (Рлт"1 )(l + Y2 )(l + Y2 + h2 )_1 Тогда выражение (32) можно переписать следующим образом:
/ \ / \ г / чп2
п
12
(33)
(34)
1
f (h 2 1 - П {г
01 а
х<
1 - exp < -0,5
а ( h 2 ) ß ( h2 ) cosr + ^ß ( h2 ) ( h 2 )] 2 + 2а ( h2 ) ß ( h2 ) cosx + ß ( h2 )
( h2 )2 + 2a ( h 2 ) ß ( h2 ) cosx + ß ( h2 )2
■ х
dx - Pno з .(35)
Очевидно, теперь дифференцирование под знаком интеграла, содержа-
2
щим параметр ОСШ И2 , упрощается, т. к. интеграл в формуле (35) содержит
2
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
указанный параметр только в подынтегральном выражении. Введем обозначение:
>2
Ф ( а,Р,т ) =
aPcosx + в2
a2 + 2aPcosx + в2
1 - exp
Г a2 + 2aecosT + в
V
Также, введем следующие дополнительные обозначения:
f 2 2 Л n2 a + 2aBcosx + В ц = aBcosT + в ; V = exp--11-—
2 2 u = a + 2aвcosт + в ; к = a + Bcost
(36)
(37)
Для вычисления производной f '( h ) воспользуемся формулой для диф
(38)
ференцирования под знаком интеграла [24, 25]: /(х,у) ¿г = 4х
а а
и правилом вычисления производной от сложной функции [23]:
¿и _ ди ёх + ди ёу + ди ¿г дх ду д2
где и = / (х, у, 2) - функция, определенная в открытой области О, причем каждая из переменных х, у, г в свою очередь является функцией от переменной ? в некотором промежутке:
х = ф(*), у = V('), 2 = х(*). (40)
Тогда, взяв производную от функции / (И2), представленную выражени-
2
ем (35), по параметру ОСШ И получим:
df ( h
d ( h2 )
1 п
=-1 /I
ni
д[ф(i,B,t)] da(h2) дф(i,B,t) в(h2)
д a
d
( h2 )
дв d
( h2 )
dT.
(41)
С учетом (33), (34) и (36)-(38) формулы для расчета производных в подынтегральном выражении (41) принимают следующий вид:
где:
д[ ф ( а,в,т )] = a( a^,T ) + в ( |ДТ ) + С ( I,B,T ),
A( a,P,T ) = 2ц ( v -1) ки-2 ;
B (a^,T) = ц ( v -1) ки-1;
С(a^,T) = -bcost(v -1)u-1. Далее имеем:
^ Ф ( g|f,T )] = D ( ^в,т ) + E ( a,в,т ) + F ( a,в,т ):
(42)
(43)
(44)
(45)
(46)
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
где:
D(a,ß,x) = 2ц(v-l)(ß + acosx)и 2;
E(a,ß,x) = -2(v -1)(2ß + acosx)и 1 ;
F (a,ß,x) = vp, (ß + acosx) и 1.
В свою очередь:
da ( h
d ( h2 )
= V2
Y
,2 2 h y
1 + y 2 + h2
1 + y2 + h2 )2
22 h y
Л-1
1 + y2 + h2
(47)
(48)
(49)
(50)
dß ( h2
r 1 Л
d ( h 2 )
= -V2ln - (1 + y2 ) х
P
V лт
х
r 1л
P
V лт y
1 + Y2
1 + y 2 + h2
(1+y 2 + h
-1
(51)
Расчет первой производной функции й/(И2 ^й (И21 по формулам (41)-
(51) реализован на языке МЛТЬАВ в виде программного компонента derivаtive.m (см. приложение ко второй части статьи). После того, как получены
соотношения, позволяющие определить значение первой производной /И21
функции / (И21 для расчетов очередного приближения корня по формуле (29)
необходимо определить начальное приближение корня И02 и критерий завершения итерационного процесса по методу Ньютона. Если требуется уточнить корень с процентной погрешностью |еа | б %, то итерации метода Ньютона продолжаются до тех пор, пока на итерации N не выполниться следующее неравен-
ство:
-а N
(hN - hN-1 )/hN
100% < k
1треб '
(52)
2
Для выбора же начального приближение корня И02 согласно (30) необходимо определить вторую производную /"(И 21 функции / (И21. Для того, чтобы упростить формулы для вычисления второй производной /И21 в дальнейшем условимся использовать следующие обозначения для производных ар-
2
гументов аир функции Маркума по параметру ОСШ И : йа (И2 У й (И2 )= йа/ й (И2) и йв (И2й (И2 )= й (И2).
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
В соответствии с определением второй производной [24, 25], на основании соотношения (41) вторую производную /"(к2) = d2/(к^d(к2) найдем по формуле:
d2f (h ) d
d
(h2 f~ d (h2 )
1 n 1J
ni
Эф ( a,ß,T ) da Эф ( a,ß,T ) dß
Э a + öß djh2)
dx
1 n d
J
П 0 d ( h2 )
1 fi d = "~J
Эф ( a,ß,T ) da Эф ( a,ß,T ) dß
àa + öß dp)
dx =
п о I d ( h
Производные
Эф ( a,ß,T ) da öa d(h2
+ ■
d
d (h
Эф ( a,ß,T ) dß
Э ß dÜ2
dx.(53)
( h2 )
d
[•] в подынтегральном выражении (53) определим,
воспользовавшись правилом взятия производной от произведения двух функций [24]:
(и • V) = и • V + и • V'. (54)
Тогда, с учетом (54) получим:
d Эф ( a,ß,T ) da
( h2 ) öa d ( h 2 )J
d Э ф ( a,ß,T ) da Э ф ( a,ß,T ) d da
d ( h 2 ) Э a d (h 2 ) ' Э a d (h 2 ) d (h 2 )
d Э ф ( a,ß,T ) da 1 Э ф ( a,ß,T ) d 2a
d ( h 2 ) Э a d ( h 2 ) Э a d ( h 2 )2 1 '
d Эф ( a,ß,T ) dß
( h2 ) Эß d(h 2 )
d Э ф ( a,ß,T ) dß ( Э ф ( a,ß,T ) d dß
d ( h 2 ) Э ß d ( h 2 ) Э ß d ( h 2 ) d (h 2 )
d Э ф ( a,ß,T ) dß ( Э ф ( a,ß,T ) d 2ß
d ( h 2 ) Э ß d ( h 2 ) Э ß d ( h 2 ) 2.
(55)
DOI: 10.24412/2410-9916-2022-4-256-312
(56)
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
Определим производную-
d
d
( *2 )
0ф ( a,p,x )
ôa
, входящую в (55) пользуясь
правилом вычисления производной от сложной функции:
d
d ( h 2 )
0ф ( a,p,x )
ôa
ô 2ф (a,p,x) da ô 2ф (a,p,x) dp ôa2 dp)+ ôaôp dp)
(57)
С учетом (57) выражение (55) перепишем в виде:
d
d ( h 2 )
ôф ( a,P,x ) da ôa d(h 2)
ô 2ф ( a,P,x ) da ô 2ф ( a,P,x ) dp
ôa2 d(h2) ôaôp d(h2)
+
+
ôф (a,p,x) d 2a _ ô 2ф (a,p,x ) ôa d (h2 )2 " " ôф ( a,P,x ) d 2a
ôa2
da
d(h2)
da
d ( h 2 )
+
ôa
d
( h2 )2
ô 2ф ( a,P,x ) da dp + ôaôp dp) dp)+
(58)
Определим производную
d
d
( h2 )
ôф ( a,e,T ) ôp
входящую в (56) пользуясь
правилом вычисления производной от сложной функции:
d
d ( h 2 )
дф ( a,P,x ) ôp
ô ф(a,p,x) da ô ф(a,p,x) dp
ôPôa d(h2)
ôp2
d ( h 2 )
(59)
С учетом (46) выражение (43) перепишем в виде:
d
d ( h 2 )
ôф ( a,P,x ) dp ôp dp]
ô 2ф ( a,P,x ) da ô 2ф ( a,P,x ) dp
ôpôa d(h2) ôp2 d(h2)
dp
+
ôф ( a,P,x ) d 2P ô 2ф ( a,P,x ) da dp ô 2ф ( a,P,x )
+ ■
ôP d (h 2 )2 ôpôa d ( h 2 ) d (h 2 ) ôp2 ôф ( a,P,x ) d 2P
d ( h 2 )
2
+
dp
d ( h 2 ) (60)
+
d(h2)2
Подставив (59) и (60) в (53) получим следующее выражение для второй d2f ( h ) .
производной
dh
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
d 2f ( h ) = _ ! ï TT "
d(h2)2 П
1 71 1 {
TT J
d
d
(h2)
d ф ( a,ß,x ) da da d(h2)
d
d
(h2)
дф ( a,ß,x ) dß dß djh2)
dr =
n •
d2 ф ( а, ß,x ) da
da2 d (h2 )
d2 ф ( a,ß,x ) da dß дф ( a,ß,x ) d2
а
dadß d (h2 ) d (h2
da d(h2
dr -
1 ж -1 \
Я Î
d 2ф ( a,ß,x )
( \ dß
dß2
1"2 )
d
V v U
d 2ф ( a,ß,x ) da dß 5ф ( a,ß,x ) d 2ß
dßda d ( h2 ) d ( h2
dß
d
dr. (61)
В (61) подлежат определению следующие производные:
(h2 )2
d 2ф ( a,ß,x )
da2
d 2ф ( a,ß,x ) d 2a d 2ф(аДт) d 2ф(a,ß,т) d 2ß
-, -т, -, -;;- И -
dadß
а
d ( h 2 )
2
для второй производной
dßda
d 2ф ( a,ß,x )
dß2
d ( h 2 )
том (42):
_d_ da
5ф ( a,ß,x )
da
da
d
-—-, взяв производную —
2 da
j. Получим выражение d Г5ф(a,ß,x)
d a
с уче-
—Г A ( a,ß,x ) + B ( a,ß,x ) + С ( a,ß,x )]
da
dA (a,ß,x) dB (a,ß,x) dC (a,ß,x)
da da
Из (43) - (45) следует: dA ( a,ß,x )
da
da
A1 ( a,ß,x ) + A2 ( a,ß,x ) + A3 ( a,ß,x ) + A4 ( a,ß,x ) ;
где, с учетом (37), можем записать: A1 ( a,ß,x ) = 2ц ( v -1) и
-2.
A2 (a,ß,x) = -2vp,(a + ßcosx)2 u 2 ; A3(a,ß,x) =-8ц(а + ßcosx) (v - 1)u A4 ( a,ß, x ) = 2ß ( v -1) KU-2 cos x ;
dB ( a,ß,x ) = B1 ( a,ß,x ) + B2 ( a,ß,x ) + B3 ( a,ß,x ) + B4 ( a,ß,x ),
-3.
da
где:
B1 (a,ß,x) = pvu- ; B2 ( a,ß,x ) = -цук2 u-1 ; B3 ( a,ß,x ) = -2^k 2 u-2 ;
(62)
(63)
(64)
(65)
(66)
(67)
(68)
(69)
(70)
(71)
0
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
B4 ( а,в,т ) = Pvku 1 cost .
ÔC ( а,в,т )
Эа
= Ci( а,в,т ) + C2( а,в,т ):
где:
Ci(а,в,т) = -2p(v-1)ки 2cost
c2 ( а,р,т ) = -pvku 1 cost .
Ô 2ф ( а,р,т )
Получим выражение для второй производной 2
(72)
(73)
(74)
(75)
, взяв производ-
Ô
ную — ÔP
Эф( а,в,т)
Ô р д_ ÔP
с учетом (46):
Эф ( а,в,т ) ÔP
Ô
= — Г D ( а,Р,т ) + E ( а,Р,т ) + F ( а,Р,т )] = Ôp
где:
_ ÔD(а,р,т) ÔE(а,р,т) ÔF(а,в,т)
= ÔP ÔP ÔP .
Из (46)-(48) следует:
ÔD () = D1 ( а,в,т ) + D2 ( а,Р,т ) + D3 ( а,в,т ) + D4 (а,Р,т ), D ( а,Р,т) = 2ц( v -1) и-2 ;
2 9
D2(а,Р,т) = -2pv(в + аcosт) и- ; D3 ( а,в,т ) = 4 ( v -1) и-2 (в + а cost )2 ;
Я 2
D4(а,в,т) = -8ц(v -1)и-3 (в + а^т) .
ÔE ( !:в,Т ) = Ei ( а,в,т ) + E2 ( а,в,т ) + E3 ( а,в,т ), Ôp
Ei (а,в,т) = 2(v -1)и-2 (2в + аcosт)(в + аcosг) ;
E2 (а,в,т) = -2 ( v -1) и-1 ;
E3 ( а,в,т ) = vu-1 (в + аcosг)( 2в + аcosт ).
ÔF ( !:в,Т ) = F1 ( а,в,т ) + F2 ( а,в,т ) + F3 ( а,в,т ) + F4 (а,в,т ),
ôp
где, с учетом обозначений (37), имеем: F1 ( а,в,т ) = v^u-1 ;
F2 ( а,в,т ) = -v^ (в + аcosт )2 и-1 ;
2 1
F3(а,в,т) = v(аcosт + в)(в + аcosт) и- ;
где:
(76)
(77)
(78)
(79)
(80) (81) (82)
(83)
(84)
(85)
(86)
(87)
(88) (89)
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
2 , -2
F(а,в,т) = -2vp,(в + acosT) u В рассматриваемом случае смешанные производные
(90)
д 2ф ( а,в,т ) дадв
и
д 2ф ( а,в,т ) дрда
отличающиеся лишь порядком дифференцирования, равны между
собой, т. к. выполняется условие их непрерывности: а2ф(а,в,т)_а2ф(а,в,т)_ а ЭаЭр араа аа Соответственно, с учетом (42), можно записать:
дф ( а,в,т ) _ д дф (а,в,т )
дв дв да
д_ да
дф ( а,в,т ) дв
= f[ D ( а,в,т ) + Е ( а,в,т ) + F ( а,в,т )] = да
дЭ ( а,в,т ) дЕ ( а,в,т ) дF ( а,в,т )
да да
Из (47)-(49) следует:
дЭ ( а,в,т )
да
да
= G1 ( а,в,т ) + G2 ( а,в,т ) + G3 ( а,в,т ) + G4 ( а,в,т )
где:
Gi (а,в,т) = 2ц(v -1)и 2 cost ;
где:
где:
G2 (а,в,т) = -8ц(v -1)ки- (в + acosT) ;
G3 ( а,в,т ) = -2^vkd-2 (в + a cost ) ;
G4(а,в,т) = 2в(v -1)и-2(в + ассБт)cost .
дЕ ( а,в,Т ) = K1 ( а,в,т ) + K2 ( а,в,т ) + K3 ( а,в,т ), да
K1 (а,в,т) = 2(v -1)и-2 (2в + a cost)к ; K2 (а,в,т) = (v -1)и-1 cost ; K3 ( а,в,т ) = vu-1 ( 2в + а cost ) к.
dF( а,в,Т ) = L1 ( а, в, т ) + L1 ( а, в, т ) + L3 ( а, в, т ) + L3 ( а, в, т ), да
Li ( а,в,т ) = vu-Vcost ; L2 ( а,в,т ) = 2v^u-2k (в + aссsт ) ; L3 (а,в,т Hvu-1 (в + aссsт)cost ; L4 ( а,в,т ) = -v^u-1 (в + acosT ) к.
(91)
(92)
(93)
(94)
(95)
(96)
(97)
(98)
(99) (100) (101)
(102)
(103)
(104)
(105)
(106)
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
d2a(h 2
Получим выражение для второй производной--—, взяв производную
d ( h 2 ) '
d
d ( h ■
da ( h '
d ( h ■
с учетом (50):
d
d ( h '
da
d (h
= H ! ( y 2, h 2 ) + H 2 ( y 2, h 2 ),
(107)
где:
H ! ( y 2, h 2 | =
y
j 2 2
h y
(l + y2 + h2 ) (l + y2 + h H2 ( y 2," h 2 ) =
f ,22 V1/2
h y 1 + y 2 + h 2
(108)
-0,25>/2
Y
,2 2 h y
1 + Y2 + h2 (l + y2 + h2) Выражение для второй производной
f 7 2 2 Y 2 h y 2
1 + y 2 + h2
(109)
d2в(h 2
d ( h 2 )
Y получим, взяв производ-
ную
d
( a
d ( h
de ( h d ( h 2
d
d ( h 2 )
d|
с учетом (51) в ( * 2 )
d ( h 2 )
= 11 ( Y 2, * 2, Рлт ) +12 ( Y 2, * 2, Рлт ),
(110)
где:
' 1 Л
11 ( Y 2, * 2, Рлт )
V2ln — (1 + y 2 )
P
V лт
С Л \
V
ln
P V лт
1 + y 2
1 + y 2 + h 2
(1 + y 2 + h 2 ) '
(111)
2
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
12 ( Y 2, h 2, Рлт )
V2in
M^2 p
V лт у
(1 + y 2 )2
(112)
4
'■n ( 1 Л
p
v V p лт у
1 + Y '
1 + y 2 + h 2
(l + y 2 + h 2 )4
Теперь, воспользовавшись критерием (30) и полученными соотношениями (53)-(112) для второй производной /"(к2), определим начальное прибли-
о
жение к корня уравнения (23) для метода Ньютона. При этом учтем, что как следует из рис. 2, / (к} ) = / (55) < 0 и (к2 ) = / (60 )> 0.
На рис. 4 приведен график изменения второй производной /"(к2) при
= [55,60] уточнения кор-
2
изменении аргумента ОСШ к в промежутке ня уравнения (23).
h/2, hl
Как видно на рис. 4 в промежутке
h2, h2
= [55,60] вторая производная
/"(к2) принимает отрицательные значения. Следовательно, на основании кри-
терия (30) принимаем f (h/2)f " (h/2) > 0.
/ /
начальное приближение \ = Щ = 55, т. к.
3
2
4x16 5
6x16 51---
55 56 57 58 59 60
2
Значение ОСШ п
Рис. 4. График изменения второй производной /"(к2
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
Ход итерационного процесса уточнения корня уравнения (23) методом
22 *р з
2
Ньютона, т. е. значения ОСШ кр з при выборе начального приближения корня
2
kf = 55 , иллюстрирует таблица 4.
Таблица 4 - Пошаговое выполнение метода Ньютона
N Приближение корня, к2 Значение ошибки е a , % Значение функции, f (k2N ) Значение производной, f'(kN)
1 55 — —0,0011 5,1977х10-4
2 57,2124 3,8670 —7,8983х10-5 4,5046х10-4
3 57,3877 0,3055 —4,3751х10-7 4,4548х10-4
4 57,3887 0,0017 —1,3598х10-11 4,4545х10-4
5 57,3887 5,3193х10-8 1,1102х 10-16 4,4545х10-4
Сравнивая данные таблиц 2, 3 и 4, приходим к выводу, что корни уравнения (23), уточненные методом половинного деления, методом хорд и методом Ньютона совпали с точностью до четвертого знака после запятой:
к2 з = 57,3887 (метод Ньютона), к^ з = 57,3887 (метод бисекций) и
^ з = 57,3887 (метод хорд). Тем самым подтверждается корректность полученных аналитических выражений (41)-(51), (53), (55)-(58), (60)-(112) для вычисле-
9 9
ния первой /'(к ) и второй /"(к ) производной от вероятности правильного обнаружения сигналов (Рпо), выражаемой через ^-функцию Маркума, по параметру ОСШ (к ), т. к. сами алгоритмы, реализующие эти методы, существенно отличаются.
Используя метод Ньютона, удалось уточнить этот корень за пять итераций (Ы = 5), что на одну итерацию меньше чем при использовании метода хорд (Ы = 6) и более чем в три раза меньше чем 18 итераций (Ы = 18), потребовавшихся методу половинного деления.
Относительная процентная погрешность корня уравнения (23), уточнен-
I I —8
ного методом Ньютона, |еа|5 = 5,3193 х 10 8 % на два порядка меньше, чем относительная ошибка |еа|6 = 9,9086 х 10—6 %, корня уточненного методом хорд и на три порядка меньше, чем относительная ошибка метода половинного деления |еа|18 = 6,6471 х10—5 %. Таким образом, приходим к выводу, что метод
Ньютона при уточнении корня уравнения (23) предпочтительнее метода половинного деления и метода хорд по критериям количества итераций и точности. К недостаткам рассмотренного метода можно отнести то, что до начала расчета необходимо найти начальное приближение корня, а для этого требуется рас-
считать вторую производную /"(к ) на интервале к/,ки =[55,60] поиска уточненного значения корня.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Усредненные временные затраты уточнения корня составили 1,1439 с, что в 49,79 раза больше, чем у метода половинного деления (0,022971 с) и в 86,62 раза больше, чем у метода хорд (0,013206 с).
Численный метод Ньютона реализован на языке МЛ^ЛВ в виде программного компонента Newton.m (см. приложение ко второй части статьи).
2.5. Численный метод Ньютона-хорд
Сравнивая между собой данные таблиц 3 и 4, иллюстрирующих пошаговое выполнение метода хорд и метода Ньютона для уравнения (23), приходим к выводу, что метод хорд формирует последовательные приближения к искомому
значению корня с избытком км, а метод Ньютона с недостатком км (N = 1,2,...). Соединяя способ хорд и метод Ньютона, получаем комбинированный метод Ньютона-хорд или метода касательных-хорд [21], на каждом этане которого находим значения по недостатку и значения по избытку уточненного
2 2 корня к1гив уравнения
/(к 2) = 0 . Отсюда, в частности, вытекает, что цифры,
2 —2 2 общие для км и км обязательно принадлежат уточненному корню к(гие [21].
Постановка частной задачи: требуется уточнить корень уравнения (23)
2 2 2 приближенно равный кр з ~57,5 и заключенный в промежутке [к/,ки] = [55,60]
методом Ньютона-хорд с относительной (процентной) погрешностью не более
|е1 = 10-4 %.
I а 1треб
9 9 9 9
Решение. Как было показано ранее /(к/) • /(ки) < 0, a /'(к ) и /"(к )
9 9
непрерывны и сохраняют постоянные знаки / '(к ) > 0, /"(к ) < 0 на рассмат-
риваемом отрезке
h2, hl
= [55,60] локализации корня уравнения (23).
Итерационный процесс метода Ньютона-хорд реализуем на основе выражений:
N
h 2 =
^/(Км) - Км/Ш )\\1_Г(кш) - /»)\ , (113)
к% = к2ы-1 - /(к2 -1V/'(к2 -1), м = 1,2,., (114)
полагая, в соответствии с критерием (30), начальное приближение корня для метода Ньютона к0 = к2 = 55.
Как показано в [21], особенностью метода Ньютона-хорд является то, что
2 —2
на каждом шаге метод хорд применяется к новому отрезку [к^,к-^].
2 2 —2 Из сказанного выше следует, что к^ < кГие < км и
0 < к1ие - к2м < Ц - к2,. (115)
Так как допустимая относительная погрешность приближенного корня
2 I | -4
км задана заранее и равна |ва| б = 10 %, то процесс сближения заканчивается на той итерации, когда будет обнаружено, что
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
a\N - (4 " 4)/41100% < |ee|Треб . (116)
По окончании итерационного процесса в соответствии с рекомендациями [21], за истинное (англ. true) значение корня htrue примем среднее арифметическое полученных последних значений:
htrue N - (hN + hNo/^. (117)
Ход итерационного процесса уточнения корня уравнения (23) методом Ньютона-хорд иллюстрирует таблица 5.
Таблица 5 - Пошаговое выполнение метода Ньютона-хорд
N Приближение корня с недостатком, 4 Усредненное приближение корня .2 htrue N Приближение корня с избытком, Значение ошибки, N а |, % Значение функции, f ( htrue N )
1 57,2124 57,3992 57,5861 — 4,6859x10—6
2 57,3877 57,3888 57,3898 0,0036 2,6562x10—8
3 57,3887 57,3887 57,3887 1,1284x10-'7 8,2545x10—13
Из данных таблицы 5 следует, что вычислительный процесс уточнения корня уравнения (23) завершается за три итерации (N = 3). Значение уточненного корня Н^ие N в соответствии с (117) принимает значение 3 = кр з =
= 57,3887 с относительной процентной погрешностью |еа| = 1,128425 х 10_7 %,
что на три порядка меньше, чем |еа| б = 10"4 %.
Сравнивая между собой данные таблиц 2, 3, 4 и 5, приходим к выводу, что корни уравнения (23), уточненные методом половинного деления, методом
2
хорд, методом Ньютона и методом Ньютона-хорд совпали: Нр з = 57,3887.
По данным таблиц 4 и 5 на рис. 5 построены графики динамики сходимости итерационного процесса методом Ньютона и Ньютона-хорд к корню уравнения (23).
Как следует из рис. 5, сходимость итераций метода Ньютона монотонна вблизи корня и немонотонна вдали от него, а вот у метода Ньютона-хорд сходимость монотонна, как вблизи корня, так и вдали от него. При этом скорость сходимости итерационного процесса для простого корня достаточно велика.
Используя метод Ньютона-хорд, удалось уточнить этот корень за три итерации (N = 3), что на две итерации меньше чем при использовании метода Ньютона (N = 5), в два раза меньше чем при использовании метода хорд (N = 6) и в шесть раз меньше чем 18 итераций (N = 18), потребовавшихся методу половинного деления.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
ГЧ 4
■с
,_,
G; s см -S
I ■С"
ф
CL го т ч Q_ 0 X 1
Ю го
О I
G; о
s I Q
ф т -О X
m ГО s
Q. > см > -г
9
I CL го т
(J ü о I-
ü
0) л
^ т
i ш о о 1-
ф * s ц
io ф
s о. s
1=
57.4
57.3975
57.395
57.3925
57.39
57.3875
57.385
\ - \ /
h2 157.3887
\ / \/ \ / N
/ \ H*
j \ У tme J57.38871
/ / /
1
Номер итерации N
58
57
56
55
Рис. 5. Динамика сходимости итерационного процесса методов Ньютона (к^) и Ньютона-хорд (к{гие) к корню уравнения (23)
Относительная процентная погрешность корня, уточненного методом Ньютона-хорд, составляет |еа| = 1,1284 х 10-7 %, что сопоставимо с относительной погрешностью корня, уточненного методом Ньютона,
I I -8
|еа| = 5,3193 х 10 %. Следовательно, наиболее более эффективным по критериям количества итераций и точности из рассмотренных численных методов уточнения корня уравнения (23) оказался метод Ньютона-хорд. Основным недостатком метода Ньютона-хорд является то, что до начала расчета необходимо найти начальное приближение корня, а для этого требуется рассчитать вторую
производную /"(к ) на интервале его локализации.
Усредненные временные затраты уточнения корня уравнения (23) методом Ньютона-хорд составили 0,17462 с, что в 7,60 раза больше, чем у метода половинного деления (0,022971 с), в 13,22 раза больше, чем у метода хорд (0,013206 с) и в 6,55 раза меньше, чем у метода Ньютона (1,1439 с). Таким образом, приходим к выводу, что из рассмотренных методов самым быстрым методом уточнения корня уравнения (23) является численный метод хорд, а далее, по мере возрастания времени счета, располагаются следующие численные методы: метод половинного деления, Ньютона-хорд и Ньютона.
Численный метод Ньютона-хорд реализован на языке МА^АВ в виде программного компонента еотЫ.т (см. приложение ко второй части статьи).
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
3. Постановка и решение частной задачи определения значения отношения сигнал-шум, необходимого для обеспечения требуемой помехоустойчивости при пространственно-разнесенном приеме сигналов земной станцией спутниковой связи, на основе
численных методов
3.1. Постановка частной задачи и анализ возможных подходов к ее решению
В соответствии с постановкой общей задачи функциональная зависимость вероятности ошибки приема бита информации P()m в ССС при пространственно-разнесенном НК приеме сигналов на несколько (n > 2) антенн с квадратичным сложением в условиях релеевских общих замираний сигналов определяется уравнением помехоустойчивости (14) и может быть записана в общем виде, как:
Рош=^Й2,/7,р[Ар,Арк(/оД,а)]}. (118)
В дальнейшем будем полагать, что в выражении (16) выполняется равен— 9 9
ство h 2 = h 2, т. е. среднее значение ОСШ в ветвях пространственного разнесе-
— 9 9
ния антенн h 2 равно искомому ОСШ h2.
С учетом того, что коэффициент пространственной корреляции быстрых замираний в линейно разнесенных на эквидистантный интервал Ар элементах
приемной АС Р^Ар,Арк(/оД,а) , вычисляется по формуле (17), то при заданных исходных данных и фиксированном значении несущей частоты f, принадлежащей диапазону ПНЧ /о е30...100 МГц, он принимает вполне определенное числовое значение, и в дальнейшем может рассматриваться в качестве
константы R = Р^Ар,Арк(/оД,а)
= const.
Таким образом, в правой части формулы (118) можно не учитывать зависимость вероятности ошибки приема бита информации Pош от коэффициента R пространственной корреляции быстрых замираний в линейно разнесенных
-л
на эквидистантный интервал Ар элементах приемной АС, т. е. ^^ = С,(Н 2, п) .
Если, кроме того, по условию общей задачи можно зафиксировать количество ветвей (кратность) пространственного разноса элементов приемной АС (например п = 4), то вместо ОСШ на входе ПРМ (выходе приемной АС) ЗсСС
( Н 2 ), можно рассматривать тот же параметр, но с учетом кратности простран-
9 9
ственного разноса элементов приемной АС, т. е. Н = Н2 .
Тогда, требуется определить такое значение ОСШ Н2 , при котором вероятность ошибки приема бита информации Pош в ЗсСС при пространственно-разнесенном НК приеме сигналов на п антенн (например п = 4) с квадратичным сложением в условиях релеевских общих замираний сигналов, опреде-
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
лимое уравнением помехоустойчивости (14), примет требуемое значение
P = P = 10-5 •
1 ош 1 ош тр 10 •
Рош = ^,2 ,2 = Рош тр = 10"5 . (119)
кп =к тр п
Анализ уточненной задачи показывает, что уравнение (118) является нелинейным. В настоящее время аналитическое решение уравнения (118) не известно, а, следовательно, для поиска его решения необходимо применить численные методы. При этом наиболее близкими являются численные методы поиска всех или некоторых корней нелинейного уравнения с одним неизвестным [20-22]. Действительно, если принять:
/ (к2) = а к2) - Рош тр, (120)
то уравнение (118) можно представить в виде:
/(к2) = 0. (121) В развернутом виде уравнение (121) можно записать так:
) — Рош _ Рош тр —
>-5
= ф2,/0,я,я[Ар,АрК(/0 Д,а)
2 2 Роштр — (122)
к =к тр п
о
Вычисление функции /(к ) по формуле (122) реализовано на языке МЛТЬАВ в виде программного компонента 1Ъг.т (см. приложение ко второй части статьи).
Таким образом, задачу поиска ОСШ, необходимого для обеспечения требуемой помехоустойчивости при пространственно-разнесенном приеме сигналов в ЗсСС можно конкретизировать следующим образом:
Дано: Рош = Рош тр = 10-5, ¡5 - 400 м, кэ ~ 2,5 • 105 м, ад^ = 2 • 109 эл/м3,
о
а = 650, /0 = 10 8 Гц, п = 4 и Ар = 400 м. Требуется, используя численные методы, найти все или некоторые корни уравнения (122) с относительной процентной погрешностью не более к! ~ = 10-4. г I а 1треб
Решение. Как было показано в п. 2, для поиска корня (корней) нелинейного уравнения с одним неизвестным разработан ряд численных методов, не требующих вычисления производной. При этом наиболее известными из них являются метод половинного деления, метод хорд и др.
Как отмечалось в п. 2.1, в рамках указанных методов задача поиска корней нелинейного уравнения с одним неизвестным распадается на несколько подзадач: во-первых, надо исследовать количество, характер и расположение корней (локализовать корни); во-вторых, найти приближенные значения корней (получить начальное приближение); в-третьих, выбрать из них интересующие корни и уточнить их с требуемой точностью на основе численных методов.
Как было показано в п. 2.1, первая и вторая подзадачи решаются аналитическими и графическими методами. В рассматриваемом случае, в рамках первой подзадачи на основе аналитических соотношений (14)-(17) и оговоренных
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
ранее исходных данных, при п = 4, был построен график функции (122)
9 9
/(кпри изменении ОСШ к2р4 в интервале от 154 до 159 (рис. 6).
Рис. 6. График функции f (hпри /0ш тр = 10 5, h ~ 400 м, кэ ~ 2,5 • 105 м, аду = 2 • 109 эл/м3, а = 1,1344640138 рад, f0 = 108 Гц, Лр = 400 м и n = 4
На рис. 6 видно, что график исследуемой функции /(к2р4) пересекает
-у
ось абсцисс один раз в интервале значений ОСШ на входе ПРМ ЗС ССС (к 2р4)
9 9 9
от к2 = 154 до к2 = 159. На концах этого интервала функция /(к^4) имеет разные знаки. Поэтому уравнение (122) имеет один действительный корень приближенно равный к«156,4.
В рамках третьей подзадачи выбираем интересующий нас корень к тр 4 ~ 156,4, локализованный на отрезке [156,157], и уточним его с требуемой точностью на основе численных методов.
3.2. Численный метод половинного деления
Постановка частной задачи: в п. 3.1 был отделен единственный корень
уравнения (122) приближенно равный к2р4 «156,4, который заключен в про-
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
межутке [156,157]. Требуется уточнить этот корень методом половинного деления с относительной процентной погрешностью не более |еа| б = 104 %.
Решение. Итерационный процесс метода деления пополам (бисекций)
был подробно рассмотрен в п. 2.2. Поэтому перейдем непосредственно к расо
смотрению хода итерационного процесса определения значения ОСШ кметодом половинного деления, который иллюстрирует таблица 6.
Таблица 6 - Пошаговое выполнение метода половинного деления
для уравнения (122)
Крайняя слева Средняя Крайняя справа Значение ошиб- Значение функ-
N точка, к1 точка, к2 точка, к2 ки, ед , % ции, f ( hT j
1 156 156,5000 157 - -3,5937х10"8
2 156 156,2500 156,5000 0,3195 -3,5937х10"8
3 156,2500 156,3750 156,5000 0,1600 1,0403х10"8
4 156,2500 156,3125 156,3750 0,0799 -1,2802х10"8
5 156,2500 156,2813 156,3125 0,0400 -1,2080х10"9
6 156,2813 156,2969 156,3125 0,0200 4,5955х10"9
7 156,2969 156,3047 156,3125 0,0100 1,6932х10"9
8 156,3047 156,3086 156,3125 0,0050 2,4243х10"10
9 156,3047 156,3066 156,3086 0,0025 -4,8284х10"10
10 156,3047 156,3057 156,3066 0,0012 -1,2021х10"10
11 156,3057 156,3062 156,3066 6,2478х10-4 6,1106х10"11
12 156,3057 156,3059 156,3062 3,1239х10-4 -2,9554х10"11
13 156,3059 156,3060 156,3062 1,5619х10-4 1,5776х10"11
14 156,3059 156,3060 156,3060 7,8097х10-5 -6,8892х10"12
Из данных таблицы 6 видно, что итерационный процесс уточнения корня уравнения (122) методом половинного деления завершается на четырнадцатой
итерации (Я = 14), когда искомый корень принимает значение кг14 = ктр 4 =
= 156,3060с относительной процентной погрешностью |еа|14 = 7,8097х 10-5%,
удовлетворяющей неравенству (25).
При этом следует учесть, что значение коэффициент пространственной корреляции быстрых замираний в линейно разнесенных на эквидистантный интервал Ар = 400 м элементах приемной АС равен Я = 0,88683, т. е. является значительным [10], что и приводит к такому большому значению требуемого ОСШ к2р 4.
Таким образом, метод половинного деления может быть рекомендован для уточнения корня уравнения (122), т. к. он обеспечивает высокую надежность и точность определения значения ОСШ ктр п, а скорость сходимости вы-
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
числительного итерационного процесса в этом случае не велика, но, как отмечалось ранее, этот показатель не является критичным.
Усредненные временные затраты уточнения корня уравнения (122) методом хорд составили 0,025815 с.
3.3. Численный метод хорд
Постановка частной задачи: в п. 3.1 был отделен единственный корень
о
уравнения (122) приближенно равный к~ 156,4, который заключен в промежутке [156,157]. Требуется уточнить этот корень методом хорд с относительной процентной погрешностью не более |е а| б = 0,0001 %.
Решение. Итерационный процесс метода хорд был подробно рассмотрен в п. 2.3. Поэтому перейдем непосредственно к рассмотрению хода итерационного
процесса определения значения ОСШ к методом хорд, который иллюстрирует таблица 7.
Таблица 7 - Пошаговое выполнение метода хорд для уравнения (122)
N Крайняя слева точка, hlN Точка пересечения, h 2 nsN Крайняя справа точка, h^ Значение ошибки, |е a 1, % Значение функции, f ( h2 )
1 156 156,3086 157 — -4,7955х10"10
2 156 156,3060 156,3086 0,0016 -1,7822х10"12
3 156 156,3060 156,3060 6,1181 х10_6 -6,6240х10"15
Из данных таблицы 7 следует, что в рассматриваемом случае итерационный процесс завершается на третьей итерации (Ы = 3), когда искомый корень
2 2
принимает значение к^ = к^ = 156,3060 с относительной процентной погрешностью |еа = 6,1181 х 10-6 %, что меньше требуемой еа I „ = 10-4 %.
По данным таблиц 6 и 7 на рис. 7 построен график динамики сходимости итерационного процесса методов метода половинного деления и хорд к корню уравнения (122).
Как следует из рис. 7, скорость сходимости итерационного процесса метода половинного деления уточнения корня уравнения (122), как и в случае поиска корня уравнения (23), невелика: за одну итерацию точность увеличивается примерно вдвое, т. е. уточнение трех цифр требует 10 итераций. При этом точность результата уточнения корня уравнения (122) гарантируется.
На рис. 7 видно, что сходимость итераций метода хорд не монотонна. Скорость сходимости итерационного процесса замедляется от итерации к итерации.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Рис.
156.3-1-1-1-1-1-1-1-1-1-1-1-1-1-1-155
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Номер итерации N
7. Динамика сходимости итерационного процесса методов половинного
2 2 деления (к.) и хорд (к.) к корню уравнения (122)
Сравнивая между собой данные таблиц 6 и 7, приходим к выводу, что корни уравнения (122), уточненные методом половинного деления и методом
7 7
хорд, совпали: кр 4 = 156,3060 (метод половинного деления) и к. 4 = 156,3060
(метод хорд), что и показано на рис. 7. Используя метод хорд, удалось уточнить этот корень за четыре итерации (Ы = 3), что в более, чем в четыре раза меньше чем 14 итераций, потребовавшихся методу половинного деления. При этом относительная процентная погрешность |еа|3 = 6,1181 х10_6 % корня, найденного методом хорд, в 12,76 раза меньше, чем у метода половинного деления Ын = 7,8097 х10"5 %.
Таким образом, приходим к выводу, что метод хорд при поиске корня уравнения (122) превосходит метод половинного деления по таким критериям эффективности, как количество итераций и точность. Усредненные временные затраты уточнения корня уравнения (122) методом хорд составили 0,074557 с, что в 2,89 раза меньше, чем у метода половинного деления 0,025815 с. Следовательно, в данном случае самым быстрым также является численный метод хорд.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
4. Сравнительная оценка погрешности определения коэффициента энергетической скрытности спутниковой системы связи на основе численных и графоаналитических методов
Для оценки погрешности определения КЭС ССС (у эс) по формуле (4) введем следующие обозначения:
Ду эс - абсолютная погрешность определения КЭС ССС у эс;
Г\ Г\
Дкр з - абсолютная погрешность определения параметра кр з;
Г\ Г\
Дк2р - абсолютная погрешность определения параметра к 2р.
При расчетах по формуле (4) погрешность определения КЭС ССС Ду эс
определяется на основе теоремы о погрешности частного, согласно которой относительная погрешность частного не превышает сумму относительных погрешностей делимого и делителя [21]. На основании следствия из этой теоремы, можно записать:
(123)
8Y = 8 . , 2 Y эс A h р2 з
8 9
Ah Тр
где: 8
Y эс
AY эс / Y
эс
относительная погрешность определения КЭС у
эс 5
8д/г 2 = Акр з / к2 з - относительная погрешность определения параметра кр2 з;
8
9 9 9
д/г2 = Ак2 / ктр - относительная погрешность определения параметра к 2
тр
Таким образом, для получения сравнительной оценки погрешности определения КЭС ССС Ду эс на основе численных и графоаналитических методов
необходимо вначале определить относительные погрешности 8 ^ 2 и 8
р з
Ah
тр
достигаемые при определении КЭС Ду эс графоаналитическим методом и на основе применения численных методов решения нелинейных уравнений (23) и (122). Сравнительная оценка погрешности определения КЭС ССС Ду эс на основе численных и графоаналитических методов может быть получена на основе сопоставления относительной погрешности определения значения параметров найденных, соответственно, на основе численных методов
8 Л7 2 и Ahp з
Г \
8 л, 2 V Ahp з у
8AhTp
числ.
Г \
8 м 2
V AhT2p у
и графоаналитическим методом
числ.
Г \
8 м 2
V Ahp2 з у
гр. аналит.
/
Л
8 м 2
V Ahip у
по следующей формуле:
гр. аналит.
Г \
V ^ р2 з у
п = -
г
+
Л
гр. аналит.
V ^ т2р у
гр. аналит.
/ Л
V ^ Р з у
+
числ.
^ т2р у
(124)
числ.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
где п _ число, показывающее во сколько раз найденная графоаналитическим методом относительная погрешность определения КЭС ССС у эс превосходит аналогичную погрешность, но найденную численными методами решения нелинейных уравнений (23) и (122).
В п. 2 и 3 были использованы аналитические и графическими методы на этапах локализации и получения начального приближения корней нелинейных уравнений (23) и (122) при исходных данных, сформулированных в п. 1 при постановке общей задачи исследования.
Так, для нелинейного уравнения (23) оказалось, что единственный корень локализован на отрезке [55,60], а начальное приближение корня приближенно
равно кр з « 57,5 (рис. 2). Аналогично, для нелинейного уравнения (122) установлено, что единственный корень локализован на отрезке [156,157], а началь-
-л
ное приближение корня составляет к. 4 «156,4 (рис. 6).
Кроме того, в п. 2 и 3 были использованы одни и те же численные методы для уточнения корней нелинейных уравнений (23) и (122): метод половинного деления и метод хорд. Оказалось, что уточненное значение корня нелинейного уравнений (23) при использовании численного метода половинного деления и
метода хорд совпали: кр з = 57,3887. При этом относительная погрешность
|еа1б = 9,9086 х10_6 % корня, уточненного методом хорд за шесть итераций
вычислительного процесса, на порядок меньше, чем относительная погрешность метода половинного деления, достигнутая за 18 итераций
|е Д8 = 6,6471 х10 "5 %.
Соответственно уточненное значение корня нелинейного уравнения (122) при использовании численного метода половинного деления и метод хорд также совпали кр 4 = 156,3060. При этом относительная погрешность
|еа|3 = 6,1181 х10_6 % корня, уточненного методом хорд за три итерации вычислительного процесса, на порядок меньше, чем относительная погрешность метода половинного деления, достигнутая за 14 итераций
|еа|14 = 7,8097 х10"5 %.
Учитывая сопоставимость относительных погрешностей уточнения корня нелинейных уравнений (23) и (122) методом половинного деления и методом хорд, сравнительную оценку погрешности определения КЭС ССС на основе численных и графоаналитических методов по формуле (124), проведем для случая уточнения корней этих уравнений, как для метода половинного деления, так для метода хорд. При этом учтем, что числители при вычислениях по формуле (124), как для метода половинного деления, так для метода хорд будут одинаковыми. Тогда, для метода половинного деления имеем:
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
I \ Ah23 |57,5 - 57,38871
( ôAh 2 I 100 % = ^---L100 % =
\ F / гр. аналит. h23 57,3887
*P3
0,1113
57,3887
100% = 0,1939 %; (125)
I \ AhTp n |156,4 -156,30601
(SAh2 I 100 % = J-!-!-1100 % =
\ тр / гр. аналит. hTp n
0,0946
hTP n 156,3060
100% = 6,0138-10"2 %; (126)
156,3060
(5дкр2 з) = 6,6471х 10 -5 % и
^ р з > числ.
(5дк2 ) = 7,8097х 10-5 %. (127)
V тр / числ.
Подставляя (125)-(127) в (124) получим:
(0,1939 + 6,0138-10 -2) % 0,25404
П (6,6471 х 10 -5 + 7,8097 х 10 -5) % 14,4568 х 10 -5
= 0,0175724 х 105 = 1757,24. (128)
Как следует из (128), относительная погрешность определения КЭС ССС графоаналитическим методом превосходит аналогичную погрешность численного метода половинного деления уточнения корней нелинейных уравнений (23) и (122) более чем в тысячу раз или на три порядка. Аналогично, для метода хорд:
(§Дк23) = 9,9086 х 10 -6 % и
' 3 ' числ.
тр ' числ.
6,1181 х 10 -6 %. (129)
Подставляя (125), (126) и (129) в формулу (124) получим: (0,1939 + 6,0138-10 -2) % 0,25404
П
(9,9086 х 10 -6 + 6,1181 х 10 -6) % 16,0267 х 10 -6
= 0,01585105 х 106 = 15851,05. (130)
Как следует из (130), относительная погрешность определения КЭС ССС графоаналитическим методом превосходит аналогичную погрешность численного метода хорд уточнения корней нелинейных уравнений (23) и (122) более чем в 15 тысяч раз или на четыре порядка.
Таким образом, приходим к выводу, что разработанная методика оценки КЭС ССС на основе применения численных методов уточнения корней нелинейных уравнений (23) и (122) позволяет, как минимум, на три порядка уменьшить относительную процентную погрешность определения этого параметра по сравнению с известными графоаналитическими методами.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Заключение
С целью исследования зависимости от значения несущей частоты /о отношения сигнал-шум, необходимого для обнаружения с заданным качеством сигнала ССС средством РРТР (кр з) и для пространственно-разнесенного приема узкополосных простых сигналов с требуемой вероятностью ошибки приема бита информации в ЗсСС на четырехэлементную приемную АС ССС (к 4),
были проведены дополнительные расчеты указанных параметров для частот от 30 и до 90 МГц с шагом 10 МГц. Таким образом, в сводной таблице 8 представ-
1 2
лены значения этих параметров кр з и к^ 4 при изменении частоты /о от 30 и до 100 МГц с шагом 10 МГц.
Таблица 8 - Значения параметров кр з и к2 4 при изменении частоты /о
fo, МГц 30 40 50 60 70 80 90 100
hp з> ДБ 29,393 28,527 26,865 24,57 22,134 20,112 18,641 17,588
к^р 4, дБ 16,318 16,956 17,858 18,829 19,749 20,572 21,296 21,94
По данным таблицы 8 на рис. 8 построены графики, иллюстрирующие,
? 2
как изменяются параметры к^ з и ктр 4 при изменении несущей частоты /о.
На рис. 8 видно, что указанные графики являются нелинейными. При этом график зависимости кр2 з с ростом частоты /о монотонно убывает от
29,393 до 17,588 дБ, а график зависимости А^р 4, наоборот, с ростом частоты /о
монотонно возрастает от 16,318 до 21,94 дБ. Кроме того, из анализа рис. 8 следует, что эти графики пересекаются в точке, соответствующей значению частоты, примерно равной 78 МГц. Следовательно, на этой частоте значение числителя и знаменателя дроби, представленной выражением (4) для вычисления КЭС ССС, становятся равными друг другу, а сам КЭС принимает значение равное нулю в логарифмических единицах измерения (дБ). Поэтому, для исходных данных, соответствующих формулировке общей задачи, (см. п. 1) энергетическая скрытность рассматриваемой ССС обеспечивается в интервале частот от 30 до 78 МГц, т. к. в этом случае КЭС (у эс)дБ > 0, а начиная с частоты 78 МГц
и выше, энергетическая скрытность этой ССС уже не обеспечивается, т. к. для этих частот (у эс) < 0.
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
30
Ш 28
26
24
n
Csj Q.
3 22 О
О
Ф
20
® 18
CO
16
1 -1- г\ - - —1-
У эс J ^р23-/4>0 Уэс = 0 Г 'эс <0
- н ' <rh2 Л "тр 4 1 N. 1
30
40
50
60
70 78 SO
90
100
Несущая частота Л, МГц
V 7
Рис. 8. Графики, иллюстрирующие изменение параметров кр з и к^ 4
при изменении частоты /о
Таким образом, анализ полученных результатов позволяет сделать следующие выводы:
1. Задача определения с высокой точностью, когда относительная (процентная) погрешность не превышает 10-4 %, значения ОСШ, необходимого для обнаружения сигнала ССС ПРМ радиообнаружения средства РРТР с заданным
качеством (Рпо з = 0,99 и Р лт з = 10 ~4) и значения ОСШ, гарантирующего пространственно-разнесенный прием сигналов на несколько антенн (п > Г) антенн ССС с требуемой вероятностью ошибки приема бита информации (Рош =
= Рош тр = 10 5) была успешно решена на основе численных методов поиска
корня нелинейных уравнений (23) и (122) с одним неизвестным: метод половинного деления, метод хорд, метод Ньютона и методом Ньютона-хорд (для уравнения (23)) и методами половинного деления и метода хорд (для уравнения (122)).
2. В рамках проведенных исследований впервые были получены оригинальные аналитические выражения (41)-(51) для вычисления первой / кГ | и
(53), (55)-(58), (60)-(112) для вычисления второй кГ | производной от вероятности правильного обнаружения сигналов (Рпо), выражаемую через
^-функции Маркума, по параметру ОСШ (к ) на входе ПРМ обнаружения средства РРТР, что позволило использовать численные методы Ньютона и
DOI: 10.24412/2410-9916-2022-4-256-312
Systems of Control, Communication and Security
ISSN 2410-9916
Ньютона-хорд применительно к рассматриваемой задаче уточнения корней нелинейного уравнения (23).
3. Хронометраж временных затрат численных методов при определении корня нелинейных уравнений (23) и (122) показал, что самым быстрым является численный метод хорд, а далее по мере возрастания времени счета располагаются: метод половинного деления, метод Ньютона-хорд и метод Ньютона.
4. Разработанный метод оценки КЭС ССС на основе применения численных методов уточнения корней нелинейных уравнений (23) и (122) позволяет, как минимум, на три порядка уменьшить относительную процентную погрешность определения этого параметра по сравнению с известными графоаналитическим методами.
5. Для исходных данных, соответствующих формулировке общей задачи (см. п. 1), энергетическая скрытность рассматриваемой ССС обеспечивается в интервале частот от 30 до 78 МГц, т. к. в этом случае КЭС (у эс) Б > 0, а начиная с частоты 78 МГц и выше, энергетическая скрытность ССС уже не обеспечивается, т. к. на этих частотах (у эс) < 0.
Полученные результаты и их новизна. В настоящей статье предложен метод оценки КЭС ССС с ПНЧ, который отличается от известных тем, что он:
1) впервые, на этапе формирования исходных данных для непосредственного вычисления КЭС по формуле (4), использует численные методы половинного деления, хорд, Ньютона и метода Ньютона-хорд, применяемые традиционно для уточнения корня нелинейного уравнения с одним неизвестным, что в результате позволило повысить точность определения значений ОСШ на входе
ПРМ ЗсСС ктр и на входе ПРМ радиообнаружения станции РТР кр з, а также
повысить точность определения КЭС не менее чем на три порядка и сократить время на проведение таких расчетов с часов до единиц секунд;
2) позволяет оценить этот коэффициент с позиций возможного радиообнаружения сигнала ССС при одновременном обеспечении требуемой помехоустойчивости связи, причем, с учетом того, что в канале радиообнаружения могут возникать райсовские замирания принимаемых сигналов произвольной глубины;
3) учитывает факт наличия пространственной корреляции быстрых замираний сигналов в элементах приемной АС, а также зависимость коэффициента пространственной корреляции этих замираний Я от значения несущей частоты /о, параметров ионосферы, угла а трансионосферного распространения радиоволн и расстояния Ар между элементами приемной АС.
Литература
1. Тузов Г. И, Сивов В. А., Прытков В. И., Урядников Ю. Ф., Дергачев Ю. А., Сулиманов А. А. Помехозащищенность радиосистем со сложными сигналами. - М.: Радио и связь, 1985. - 264 с.
2. Борисов В. И., Зинчук В. М., Лимарев А. Е. Помехозащищенность систем радиосвязи с расширением спектра сигналов методом псевдослучайной
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
перестройки рабочей частоты / под ред. В.И. Борисова. - М.: РадиоСофт, 2008. - 512 с.
3. Буга Н. Н. Основы теории связи и передачи данных. Часть 1. - Л.: ЛКВВИА, 1968. - 548 с.
4. Коростелев А. А., Клюев Н. Ф., Мельник Ю. А., Веретягин А. А., Губин В. А., Дулевич В. Е., Зиновьев Ю. С., Петров А. В. Теоретические основы радиолокации / под ред. В.Е. Дулевича. - М.: Советское радио, 1978. -608 с.
5. Макаренко С. И., Иванов М. С., Попов С. А. Помехозащищенность систем связи с псевдослучайной перестройкой рабочей частоты: монография. -СПб.: Свое издательство, 2013. - 166 с.
6. Макаренко С. И. Описательная модель системы спутниковой связи MUOS // Системы управления, связи и безопасности. 2019. № 3. С. 89-116. doi: 10.24411/2410-9916-2019-10306.
7. Чипига А. Ф. Анализ энергетической скрытности низкочастотных систем спутниковой связи от обнаружения сигналов // Известия ЮФУ. Технические науки. 2014. № 2. С. 209-216.
8. Чипига А. Ф., Сенокосова А. В. Защита информации в системах космической связи за счет изменения условий распространения радиоволн // Космические исследования. 2007. Т. 45. № 1. С. 59-66.
9. Кириллов Н. Е. Помехоустойчивая передача сообщений по линейным каналам со случайно изменяющимися параметрами. - М.: Связь, 1971. - 256 с.
10. Андронов И. С., Финк Л. М. Передача дискретных сообщений по параллельным каналам. - М.: Советское радио, 1971. - 408 с.
11. Калинин А. И. Распространение радиоволн на трассах наземных и космических радиолиний. - М.: Связь, 1979, 1971. - 296 с.
12. Пашинцев В. П., Солчатов М. Э., Гахов Р. П. Влияние ионосферы на характеристики космических систем передачи информации. - М.: Издательство физико-математической литературы, 2006. - 191 с.
13. Chipiga А.Б., Pashintsev V.P. Procedure for Calculating the Dependence of the Energy Concealment Factor on Carrier Frequency Selection for Low-Frequency Satellite Communications System // Automatic Control and Computer Sciences. 2016. Vol. 50. № 6. P. 408-414.
14. Раshintsev V. P., Chipiga A. F., Tsymbal V. A., Zelenevskiy V. V. Low-Frequency Satellite Communication System Technical Means' Parameters Synthesis by the Requirements for Energetic Concealment and Noise Immunity // Automatic Control and Computer Sciences. 2018. Vol. 52. № 3, P. 243-249.
15. Pashintsev V. P., Chipiga A. F., Anzin I. V. Energetic concealment of low-frequency sattelite communication system with arbitrary recession of radiointercepting receiver // REMS 2018. Multidisciplinary Symposium on Computer Science and ICT. Proceedings of the 2018 Multidisciplinary Symposium on Computer Science and ICT, Russia, (October 15, 2018, Stavropol). - Aachen, 2018. CEUR Workshop Proceedings Vol. 2254. P. 237-244.
16. Pashintsev V. P., Chipiga A. F., Kostyuk D. V., Mikhajlichenko R. M. Estimation of energy secrecy of low-frequency satellite communication system with arbitrary number of spatial diversity reception of signals // The 2nd International Scientific and Practical Conference on Innovations in Engineering and Technology
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
(Veliky Novgorod, 27-28 June 2019). - Bristol, 2019. IOP Conference Series: Materials Science and Engineering Vol. 656. 012041. doi: 10.1088/1757-899X/656/1/012041.
17. Ляхов А. В. Определение зависимости коэффициента энергетической скрытности систем спутниковой связи при обнаружении сигналов от выбора рабочей частоты на основе численных методов решения нелинейных уравнений // Современная наука и инновации. 2017. № 1 (17). С. 42-49.
18. Ляхов А. В., Пашинцев В. П., Белов А. Д., Бессмертный М. Ю. Методика оценки вероятности правильного обнаружения сигналов с райсовскими замираниями // Радиотехнические и телекоммуникационные системы. 2012. № 3 (7). С. 21-26.
19. Коржик В. И., Финк Л. М., Щелкунов К. Н. Расчет помехоустойчивости систем передачи дискретных сообщений: справочник / под ред. Л.М. Финка - М.: Радио и связь, 1981. - 232 с.
20. Калиткин Н. Н., Альшина Е. А. Численные методы: в 2 кн. Кн. 1. Численный анализ. - М.: Издательский центр «Академия», 2013. - 304 с.
21. Демидович Б. П., Марон И. А. Основы вычислительной математики: Учеб. пособие для вузов. 3-е изд., испр. - М.: Наука. Главная редакция физико-математической литературы, 1966. - 664 с.
22. Chapra S. С. Applied Numerical Methods with MATLAB for Engineers and Scientists. - Mcgraw Hill Book Co, Boston MA, 2006. - 584 p.
23. Фихтенгольц Г. М. Курс дифференциального и интегрального исчисления. В 3 т. Т. I / пред. и прим. А.А. Флоринского. 8-е изд. - М.: ФИЗМАТЛИТ, 2003. - 680 с.
24. Фихтенгольц Г. М. Курс дифференциального и интегрального исчисления. В 3 т. Т. II / пред. и прим. А.А. Флоринского. 8-е изд. - М.: ФИЗМАТЛИТ, 2003. - 864 с.
25. Смирнов В. И. Курс высшей математики, Т. 2. - Издательство «Наука», 1974. - 479 с.
26. Чучин Е. В. Свойство специальных функций в моделях качества информационных систем // Auditorium: электронный научный журнал Курского государственного университета. 2014. № 1. С. 80-89 - URL: https://api-mag.kursksu.ru/api/v1/get_pdf/1408/ (дата обращения: 15.01.2022).
Reference
1. Tuzov G. I., Sivov V. A., Prytkov V. I., Uryadnikov U. F., Dergachev U. A., Sulimanov A. A. Pomekhozashchishchennost' radiosistem so slozhnymi signalami [Noise immunity of radio systems with complex signals]. Moscow, Radio i sviaz' Publ., 1985. 264 p. (in Russian).
2. Borisov V. I., Zinchuk V. M., Limarev A. E. Pomekhozashchishchennost' sistem radiosviazi s rasshireniem spektra signalov metodom psevdosluchainoi perestroiki rabochei chastity [Noise immunity of radio communication systems with the spread of the spectrum of signals by the method of pseudo-random tuning of the operating frequency]. Edited by Borisov V. I. Moscow, RadioSoft Publ., 2008. 512 p. (in Russian).
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
3. Buga N. N. Osnovy teorii sviazi i peredachi dannykh. Chast' 1. [Fundamentals of the theory of communication and data transmission. Part 1]. Leningrad, LRMAEA, 1968. 548 p. (in Russian).
4. Korostelev A. A., Klyuyev N. F., Melnik U. A., Veretiagin A. A., Gubin V. A., Dulevich V. E., Zinov'ev Y. S., Petrov A. V. Teoreticheskie osnovy radiolokatsii [Theoretical foundations of radiolocation]. Edited by Dulevich V. E. Moscow, Sov. radio Publ., 1978. 608 p. (in Russian).
5. Makarenko S. I., Ivanov M. S., Popov S. A. Pomekhozashchishchennost' sistem sviazi s psevdosluchainoi perestroikoi rabochei chastoty [Noise immunity of communication systems with pseudo-random tuning of the operating frequency]. Saint Petersburg, Svoe izdatel'stvo, 2013. 166 p. (in Russian).
6. Makarenko S. I. Descriptive Model of the MUOS Satellite Communication System. Systems of Control, Communication and Security, 2019, no. 3, pp. 89-116. doi: 10.24411/2410-9916-2019-10306 (in Russian).
7. Chipiga A. F. Analiz energeticheskoi skrytnosti nizkochastotnykh sistem sputnikovoi sviazi ot obnaruzheniia signalov [Analysis of energy secrecy of low-frequency satellite communication systems from signal detection]. Izvestiya SFedU. Engineering Sciences, 2014, no. 2, pp. 209-216. (in Russian).
8. Chipiga A. F., Senokosova A. V. Zashchita informatsii v sistemakh kosmicheskoi sviazi za schet izmeneniia uslovii rasprostraneniia radiovoln [Protection of information in space communication systems by changing the conditions of propagation of radio waves]. Cosmic Research, 2007, vol. 45, no. 1, pp. 52-59 (in Russian).
9. Kirillov N. E. Pomekhoustoichivaia peredacha soobshchenii po lineinym kanalam so sluchaino izmeniaiushchimisia parametrami [Noise-immune messaging over linear channels with randomly changing parameters]. Moscow, Connection, 1971. 256 p. (in Russian).
10. Andronov I. S., Fink L. M. Peredacha diskretnykh soobshchenii po parallel'nym kanalam [Transmission of discrete messages over parallel channels]. Moscow, Soviet radio Publ., 1971. 408 p. (in Russian).
11. Kalinin A. I. Rasprostranenie radiovoln na trassakh nazemnykh i kosmicheskikh radiolinii [Propagation of radio waves on the paths of surface and space radio links]. Moscow, Sviaz' Publ., 1979, 1971. 296 p. (in Russian).
12. Pashintsev V. P., Solchatov M. E., Gahov R. P. Vliianie ionosfery na kharakteristiki kosmicheskikh sistem peredachi informatsii [Influence of the ionosphere on the characteristics of space information transmission systems]. Moscow, Publishing house of physical and mathematical literature, 2006. 191 p. (in Russian).
13. Chipiga A. F., Pashintsev V. P. Procedure for Calculating the Dependence of the Energy Concealment Factor on Carrier Frequency Selection for Low-Frequency Satellite Communications System. Automatic Control and Computer Sciences, 2016, vol. 50, no. 6, pp. 408-414 (in English).
14. Pashintsev V. P., Chipiga A. F., Tsymbal V. A., Zelenevskiy V. V. Low-Frequency Satellite Communication System Technical Means' Parameters Synthesis by the Requirements for Energetic Concealment and Noise Immunity. Automatic Control and Computer Sciences, 2018, vol. 52, no. 3, pp. 243-249 (in English).
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
15. Pashintsev V. P., Chipiga A. F., Anzin I. V. Energetic concealment of low-frequency sattelite communication system with arbitrary recession of radiointercepting receiver. REMS 2018. Multidisciplinary Symposium on Computer Science and ICT, Stavropol, 2018, CEUR Workshop Proceedings, vol. 2254, pp. 237244. (in English).
16. Pashintsev V. P., Chipiga A. F., Kostyuk D. V., Mikhajlichenko R. M. Estimation of energy secrecy of low-frequency satellite communication system with arbitrary number of spatial diversity reception of signals. The 2nd International Scientific and Practical Conference on Innovations in Engineering and Technology, Veliky Novgorod, 2019, vol. 656, 012041. doi: 10.1088/1757-899X/656/1/012041. (in English).
17. Lyakhov A. V. Opredelenie zavisimosti koeffitsienta energeticheskoi skrytnosti sistem sputnikovoi sviazi pri obnaruzhenii signalov ot vybora rabochei chastoty na osnove chislennykh metodov resheniia nelineinykh uravnenii [Determination of the dependence of the coefficient of energy secrecy of satellite communication systems upon detection of signals on the choice of operating frequency based on numerical methods for solving nonlinear equations]. Modern science and innovations, 2017, no. 1 (17), pp 42-49 (in Russian).
18. Lyakhov A. V., Pashintsev V. P., Belov A. D., Bessmertni M. Y. Metodika otsenki veroiatnosti pravil'nogo obnaruzheniia signalov s raisovskimi zamiraniiami [Methodology for estimating the probability of correct detection of signals with Rician fading]. Radio and telecommunication systems, 2012, no. 3 (7), p. 21-26 (in Russian).
19. Korzhik V. I., Fink L. M., Shchelkunov K. N. Raschet pomekhoustoichivosti sistem peredachi diskretnykh soobshchenii: Spravochnik [Calculation of noise immunity of discrete message transmission systems: a handbook]. Edited by Finn L. M. Moscow, Radio i sviaz' Publ., 1981. 232 p. (in Russian).
20. Kalitkin N. N., Al'shina E. A. Chislennye metody: v 2 kn. Kn. 1. Chislennyi analiz [Numerical methods: in 2 volumes. Vol. 1. Numerical analysis]. Moscow, the Publishing center "Academy", 2013. 304 p. (in Russian).
21. Demidovich B. P., Maron I. A. Osnovy vychislitel'noi matematiki [Fundamentals of Computational Mathematics]. Moscow, Science Publ., The main editorial office of physical and mathematical literature, 1966. 664 p. (in Russian).
22. Chapra S. C. Applied Numerical Methods with MATLAB for Engineers and Scientists. Mcgraw Hill Book Co, Boston MA, 2006. 584 p. (in English).
23. Fihtengol'c G. M. Kurs differentsial'nogo i integral'nogo ischisleniia. V3 t. T. I [Differential and Integral Calculus Course. In 3 volumes. Vol. 1]. Intr. and notes by. Florinsky A. A., 8th ed. Moscow, FIZMATLIT, 2003. 680 p. (in Russian).
24. Fihtengol'c G. M. Kurs differentsial'nogo i integral'nogo ischisleniia. V3 t. T. II [Differential and Integral Calculus Course. In 3 volumes. Vol. 2]. Intr. and notes by Florinsky A. A., 8th ed. Moscow, FIZMATLIT, 2003. 864 p. (in Russian).
25. Smirnov V. I. Kurs vysshei matematiki [Higher Mathematics Course, Vol. 2]. Moscow, Science Publ., 1974. 479 p. (in Russian).
26. Chuchin E. V. Svoistva spetsial'nykh funktsii v modeliakh kachestva informatsionnykh sistem [Properties of Special Functions in Information Systems Quality Models]. Auditorium: electronic scientific magazine of Kursk state university,
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
2014, no. 1. Available at: http://auditorium.kursksu.ru/pdf/001-013.pdf. (accessed 15 January 2022) (in Russian).
Статья поступила 14 ноября 2022 г.
Информация об авторе
Ляхов Алексей Владимирович - кандидат технических наук, доцент кафедры «Информационная безопасность автоматизированных систем». СевероКавказский федеральный университет. Область научных интересов: влияние возмущений ионосферы на показатели качества систем спутниковой связи; математическое моделирование; численные методы и вычислительные эксперименты в MATLAB. E - mail: lyachov_st@mail.ru
Адрес: 355017, Россия, г. Ставрополь, ул. Пушкина, д. 1.
Increasing the energy secrecy of narrow-band signals of satellite communication systems during transionospheric propagation of radio waves through an inhomogeneous ionosphere due to a decrease in carrier frequencies and space-diversified reception
Part 1. Statement of the problem of increasing energy secrecy
and its solution
A. V. Lyakhov
The relevance is due to the fact that the solution of the problem of designing advanced satellite communication systems (SCS) should be based on an integrated approach, which involves taking into account all the parameters and restrictions that provide both the requirements for the efficiency of its operation, including the requirements for both noise immunity and its secrecy from means of radio and electronic intelligence (REI). A known method of increasing the energy secrecy of the SCS with a close placement of the radio detection receiver by lowering the carrier frequency and the use of space-diversified signal reception. As a quantitative criterion of the achieved level of energy secrecy of the SCS, such a generally recognized indicator as the coefficient of energy secrecy (ECS) of the SCS is used. However, the well-known method for determining the ECS is graphic-analytical. Therefore, it does not allow automating the process of determining this indicator using modern computer technology, which makes it difficult to solve the problem of designing promising SCS. In the first part of the article, the numerical methods of bisection, chords, Newton and the Newton-chord method, which are traditionally used to refine the root of a nonlinear equation with one unknown, were considered and adapted to determine the values of the energy signal-to-noise ratio (SNR) at the input of a satellite earth station receiver communication and at the input of the radio detection receiver of the RRTR means. Thus, the necessary prerequisites have been implemented for the implementation of the methodology for the automated assessment of CES SCS, which allows, even at the stage of pre-contract work, to determine the potential energy secrecy of the designed CCC, taking into account the capabilities of the party implementing the RRTR. The aim of the work is to analyze the possibility of increasing the energy secrecy of narrow-band signals of satellite communication systems with low carrier frequencies based on the development of a methodology that allows to reduce time costs and improve the accuracy of estimating the ECS of the SCS, by using numerical methods instead of graphical methods to refine the roots of nonlinear detection and noise immunity equations and their implementation in MATLAB. Results and their novelty. The developed methodology for evaluating the ECS of satellite communication systems differs from known in that it allows you to reduce the time for assessing this coefficient from hours to units of seconds and at the same time reduce the error of its assessment by at least three orders of magnitude. In addition, for the first
DOI: 10.24412/2410-9916-2022-4-256-312
Системы управления,связи и безопасности №4. 2022
Systems of Control, Communication and Security ISSN 2410-9916
time, original analytical expressions were obtained for calculating the first and second derivatives of the Marcum Q-function with respect to the signal-to-noise energy ratio parameter, which are necessary for the software implementation of the numerical Newton and Newton-chord methods for refining the roots of the nonlinear detection equation. The practical significance is that the developed methodology for evaluating the ECS of satellite communication systems will be useful to technical specialists to substantiate the source data necessary for calculations on clarifying the parameters of promising satellite communication systems, taking into account the simultaneous provision of requirements, both for the required noise resistance and in secretions connections. In addition, this technique may be of interest to scientists and applicants leading scientific research in the field of satellite communication.
Keywords: low-frequency satellite communication system, Rician fading, transionospheric propagation of radio waves, numerical method, energy secrecy coefficient, satellite communication system, Marcum Q-function.
Information about Author
Alexey Vladimirovich Lyakhov - Ph.D. of Engineering Sciences. Associate Professor at the Department of Information Security of Automated Systems. North Caucasus Federal University. Field of research: influence of ionospheric disturbances on the quality indicators of satellite communication systems; mathematical modeling; numerical methods and computational experiments in MATLAB. E-mail: lyachov_st@mail .ru
Address: Russia, 355017, Stavropol, ul. Pushkina, 1.
DOI: 10.24412/2410-9916-2022-4-256-312