Радиотехника и электроника, 2022, T. 67, № 4, стр. 345-352

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

Г. А. Бакаева *

Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина
394064 Воронеж, ул. Старых Большевиков, 54 “А”, Российская Федерация

* E-mail: gbakaeva@rambler.ru

Поступила в редакцию 21.07.2021
После доработки 21.08.2021
Принята к публикации 17.09.2021

Полный текст (PDF)

Аннотация

Рассмотрено обнаружение ограниченного по спектру фазоманипулированного широкополосного сигнала (ФМШПС) со случайными амплитудой и фазой с учетом искажений в одноконтурной входной цепи приемника. Рассчитаны вероятности обнаружения ФМШПС в зависимости от ширины его спектра, полосы пропускания одноконтурной входной цепи приемника, характеристик сигнала и шума. Определены условия, при которых вероятность обнаружения ФМШПС максимальна.

ВВЕДЕНИЕ

В системах радиосвязи и управления применяют фазоманипулированные широкополосные сигналы (ФМШПС), расширение спектра в которых достигается за счет псевдослучайной перестройки рабочей фазы передаваемого сигнала [1, 2]. Мгновенная полоса частот ФМШПС определяется длительностью одного элемента кода псевдослучайной последовательности (ПСП) и совпадает с шириной полосы частот системы радиосвязи (СРС).

При обнаружении ФМШПС с учетом случайного изменения амплитуды и фазы на фоне шума необходимо определить начало модулирующей ПСП. Для этого максимальное значение выходного сигнала корреляционного приемника сравнивают с “порогом”, при превышении которого выносится решение. На определение времени начала ПСП ФМШПС влияют искажения, вызванные ограничением его спектра [3, 5].

В СРС с ФМШПС возможны также искажения из-за неравномерности амплитудно-частотной и нелинейности фазочастотной характеристик (АЧХ, ФЧХ) входных цепей приемника, которые проявляются тем сильнее, чем шире его спектр. Наиболее распространены одноконтурные входные цепи. Представляет интерес определить параметры одноконтурной входной цепи приемника, при которых вероятность обнаружения ФМШПС с ограниченной шириной спектра с учетом случайного изменения амплитуды и фазы на фоне шумов максимальна.

1. ВЫХОДНОЙ СИГНАЛ КОРРЕЛЯЦИОННОГО ПРИЕМНИКА

Формируемый генератором передатчика ФМШПС без учета информационной составляющей с двоичной фазовой манипуляцией можно представить как [1]

(1)
$s(t) = \sum\limits_{k = 1}^L {{{p}_{k}}rect\left[ {t - (k - 1){{\tau }_{{\text{и}}}}} \right]\cos \left( {{{w}_{0}}t + {{\varphi }_{0}}} \right)} ,$

где ${{w}_{0}} = 2\pi {{f}_{0}}$ (${{f}_{0}}$ – несущая частота ФМШПС); ${{\tau }_{{\text{и}}}}$ – длительность элемента кода ПСП; $L$ – число элементов кода ПСП на длительности бита информации, $L = {T \mathord{\left/ {\vphantom {T {{{\tau }_{{\text{и}}}}}}} \right. \kern-0em} {{{\tau }_{{\text{и}}}}}}$ ($Т$ – длительность одного бита информации); ${{\varphi }_{0}}$ – случайная начальная фаза, $\left| {{{\varphi }_{0}}} \right| \leqslant \pi $;

$rect\left[ {t - (k - 1){{\tau }_{{\text{и}}}}} \right] = \left\{ \begin{gathered} 1,\,\,\,\,(k - 1){{\tau }_{{\text{и}}}} \leqslant t \leqslant k{{\tau }_{{\text{и}}}} \hfill \\ 0,\,\,\,\,t < (k - 1){{\tau }_{{\text{и}}}},\,\,\,\,t > k{{\tau }_{{\text{и}}}} \hfill \\ \end{gathered} \right.$

– функция единичного скачка; ${{p}_{k}}$ – значения кода ПСП, $\left\{ { + 1, - 1} \right\}$, причем ${{p}_{k}} = {{p}_{{k \pm L}}}$.

На вход одноконтурной цепи приемника поступает смесь ограниченного по спектру выходными цепями передатчика ФМШПС и шума

(2)
$x(t) = a{{s}_{1}}(t - {{\varepsilon }_{0}}) + n(t),$
где $a$ – амплитуда ФМШПС на входе приемника, содержащая постоянную и флуктуирующую составляющие; ${{\varepsilon }_{0}}$ – временная задержка ФМШПС на входе приемника относительно выбранного начала отсчета;
(3)
$\begin{gathered} {{s}_{1}}(t) = \int\limits_0^t {s(\tau )h(t - \tau )d\tau } = \\ = \left\{ \begin{gathered} {{p}_{1}}\int\limits_0^t {\cos ({{w}_{0}}\tau + {{\varphi }_{0}})h(t - \tau )d\tau ,} \,\,\,0 < t < {{\tau }_{{\text{и}}}} \hfill \\ {{p}_{v}}\int\limits_{(v - 1){{\tau }_{{\text{и}}}}}^t {\cos ({{w}_{0}}\tau + {{\varphi }_{0}})h(t - \tau )d\tau } + \hfill \\ + \,\,\sum\limits_{k = 1}^{v - 1} {{{p}_{k}}} \int\limits_{(k - 1){{\tau }_{{\text{и}}}}}^{k{{\tau }_{{\text{и}}}}} {\cos ({{w}_{0}}\tau + {{\varphi }_{0}})h(t - \tau )d\tau ,} \hfill \\ (v - 1){{\tau }_{{\text{и}}}} < t < v{{\tau }_{{\text{и}}}},\,\,\,\,2 \leqslant v \leqslant L{\text{ }} \hfill \\ \end{gathered} \right. \\ \end{gathered} $
– ФМШПС, ограниченный по спектру выходными цепями передатчика;
(4)
$h(t) = {{K}_{0}}\frac{{\Delta w}}{\pi }sinc\frac{{\Delta wt}}{2}\cos {{w}_{0}}t$
– импульсная характеристика выходной цепи передатчика;

$\Delta f = {{\Delta w} \mathord{\left/ {\vphantom {{\Delta w} {2\pi }}} \right. \kern-0em} {2\pi }}$ – рабочая полоса частот выходной цепи передатчика, соизмеримая с шириной спектра элемента кода ПСП ФМШПС;

$\operatorname{sinc} x = {{\sin x} \mathord{\left/ {\vphantom {{\sin x} x}} \right. \kern-0em} x}$; ${{K}_{0}}$ – коэффициент пропорциональности;

(5)
$w(a,{{\varphi }_{0}}) = \frac{a}{{2\pi \sigma _{a}^{2}}}\exp \left[ { - \frac{{{{a}^{2}} + a_{0}^{2} - 2a{{a}_{0}}\cos {{\varphi }_{0}}}}{{2\sigma _{a}^{2}}}} \right]$
– плотность распределения амплитуды и фазы, ${{a}_{0}}$ – постоянная составляющая амплитуды ФМШПС, $\sigma _{a}^{2}$ – дисперсия флуктуирующей составляющей амплитуды, $a \geqslant 0$; $n(t)$ – шум аппаратуры, функция корреляции которого равна
(6)
$\left\langle {n({{t}_{1}})n({{t}_{2}})} \right\rangle = \frac{{{{N}_{0}}}}{2}\delta ({{t}_{1}} - {{t}_{2}}),$
${{N}_{0}}$ – односторонняя спектральная плотность шума.

Для одноконтурной входной цепи приемника, настроенной на несущую частоту ФМШПС ${{f}_{0}}$ с полосой пропускания $\Omega = 2\pi \Delta {{f}_{{{\text{пр}}}}} \ll 2\pi {{f}_{0}}$, соизмеримой с шириной спектра элемента кода ПСП, квадрат модуля частотной характеристики можно представить [4] в виде

(7)
${{\left| {{{k}_{1}}(jw)} \right|}^{2}} = \frac{{K_{1}^{2}{{\Omega }^{2}}}}{{{{w}^{2}} + {{\Omega }^{2}}}},$
где ${{K}_{1}}$ – коэффициент передачи входной цепи на несущей частоте.

На выходе одноконтурной цепи приемника выражения для сигнальной и помеховой составляющих имеют соответственно вид

(8)
$\begin{gathered} {{s}_{2}}(t) = a{{{\hat {s}}}_{2}}(t) = a\int\limits_0^t {{{s}_{1}}(x){{h}_{1}}(t - x)dx} = \\ = a\left\{ \begin{gathered} {{p}_{1}}\int\limits_0^t {dx} \int\limits_0^x {\cos ({{w}_{0}}\tau + {{\varphi }_{0}})h(x - \tau ){{h}_{1}}(t - x)d\tau } , \hfill \\ 0 < t < {{\tau }_{{\text{и}}}} \hfill \\ \left[ {{{p}_{v}}\int\limits_{(v - 1){{\tau }_{{\text{и}}}}}^t {dx} \int\limits_{(v - 1){{\tau }_{{\text{и}}}}}^x {\cos ({{w}_{0}}\tau \, + \,{{\varphi }_{0}})h(x\, - \,\tau ){{h}_{1}}(t\, - \,x)d\tau } \, + } \right. \hfill \\ + \,\,\sum\limits_{n = 1}^{v - 1} {{{p}_{n}}} \int\limits_{(n - 1){{\tau }_{{\text{и}}}}}^{n{{\tau }_{{\text{и}}}}} {dx} \int\limits_{(n - 1){{\tau }_{{\text{и}}}}}^x {\cos ({{w}_{0}}\tau + {{\varphi }_{0}})} \times \hfill \\ \times \,\,h(x - \tau ){{h}_{1}}(t - x)d\tau + {{p}_{{v - 1}}}\int\limits_{(v - 2){{\tau }_{{\text{и}}}}}^t {dx} \times \hfill \\ \times \,\,\int\limits_{(v - 2){{\tau }_{{\text{и}}}}}^{(v - 1){{\tau }_{{\text{и}}}}} {\cos ({{w}_{0}}\tau + {{\varphi }_{0}})h(x - \tau ){{h}_{1}}(t - x)d\tau } + \hfill \\ + \,\,\sum\limits_{n = 1}^{v - 2} {{{p}_{n}}} \int\limits_{(n - 1){{\tau }_{{\text{и}}}}}^{n{{\tau }_{{\text{и}}}}} {dx} \times \hfill \\ \left. { \times \,\,\int\limits_{(n - 2){{\tau }_{{\text{и}}}}}^{(n - 1){{\tau }_{{\text{и}}}}} {\cos ({{w}_{0}}\tau + {{\varphi }_{0}})} h(x - \tau ){{h}_{1}}(t - x)d\tau } \right], \hfill \\ (v - 1){{\tau }_{{\text{и}}}} < t < v{{\tau }_{{\text{и}}}},\,\,\,\,2 \leqslant v \leqslant L \hfill \\ \end{gathered} \right. \\ \end{gathered} $
и
(9)
${{n}_{2}}(t) = \int\limits_0^t {n(\tau ){{h}_{1}}(t - \tau )d\tau } ,$
где
(10)
${{h}_{1}}(t) = {{K}_{1}}\Omega \exp \left[ { - \Omega \left| t \right|} \right]\cos {{w}_{0}}t$
– импульсная характеристика входной цепи приемника.

Выходной сигнал корреляционного приемника с квадратурными каналами определяется по выражению [5]

(11)
$\begin{gathered} M(\varepsilon ) = {{\left[ {\int\limits_{{{\varepsilon }_{{{\text{оп}}}}}}^{T + {{\varepsilon }_{{{\text{оп}}}}}} {\hat {x}(t){{s}_{c}}(t - {{\varepsilon }_{{{\text{оп}}}}})dt} } \right]}^{2}} + \\ + \,\,{{\left[ {\int\limits_{{{\varepsilon }_{{{\text{оп}}}}}}^{T + {{\varepsilon }_{{{\text{оп}}}}}} {\hat {x}(t){{s}_{s}}(t - {{\varepsilon }_{{{\text{оп}}}}})dt} } \right]}^{2}}, \\ \end{gathered} $
где
(12)
$\hat {x}(t) = a{{\hat {s}}_{2}}(t - {{\varepsilon }_{0}}) + {{n}_{2}}(t);$
${{\varepsilon }_{{{\text{оп}}}}}$ – время задержки опорного сигнала квадратурного приемника относительно выбранного начала отсчета;
(13)
$\left. \begin{gathered} {{s}_{c}}(t) \hfill \\ {{s}_{s}}(t) \hfill \\ \end{gathered} \right\} = \sum\limits_{v = 1}^L {{{p}_{v}}} rect\left[ {\frac{{t - (v - 1){{\tau }_{{\text{и}}}}}}{{{{\tau }_{{\text{и}}}}}}} \right]\left\{ \begin{gathered} \cos {{w}_{0}}t \hfill \\ \sin {{w}_{0}}t \hfill \\ \end{gathered} \right\}$
– опорные сигналы квадратурных каналов.

Аддитивная смесь ограниченного по спектру ФМШПС и шума проходит одноконтурную цепь приемника, поступает в квадратурные каналы, перемножается с опорными сигналами, интегрируется, возводится в квадрат и суммируется. С учетом (12), (13) выражение (11) принимает вид

(14)
$M(\varepsilon ) = {{a}^{2}}{{S}^{2}}(\varepsilon ) + 2aS(\varepsilon ){{G}_{1}} + G_{c}^{2} + G_{s}^{2},$
где $\varepsilon = {{\varepsilon }_{0}} - {{\varepsilon }_{{{\text{оп}}}}}$ – время задержки принимаемого ФМШПС относительно начала ПСП опорного сигнала;
(15)
$S(\varepsilon ) = {{\left[ {S_{c}^{2}(\varepsilon ) + S_{s}^{2}(\varepsilon )} \right]}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}$
– огибающая корреляционной функции принимаемого ФМШПС и опорного сигналов [3];

(16)
$\left. \begin{gathered} {{S}_{c}}(\varepsilon ) \hfill \\ {{S}_{s}}(\varepsilon ) \hfill \\ \end{gathered} \right\} = \int\limits_0^Т {{{{\hat {s}}}_{2}}} (t - \varepsilon )\left\{ \begin{gathered} {{s}_{c}}(t) \hfill \\ {{s}_{s}}(t) \hfill \\ \end{gathered} \right\}dt;$
(17)
${{G}_{1}} = {{G}_{c}}\cos [{{\varphi }^{/}} + \chi (\varepsilon )] - {{G}_{s}}\sin [{{\varphi }^{/}} + \chi (\varepsilon )];$

$\chi (\varepsilon ) = - arctg\left[ {\frac{{{{S}_{s}}(\varepsilon )}}{{{{S}_{c}}(\varepsilon )}}} \right]$; ${{\varphi }^{/}} = {{\varphi }_{0}} - {{w}_{0}}\varepsilon $; ${{G}_{c}}$ и ${{G}_{s}}$ – нормальные случайные некоррелированные величины

(18)
$\begin{gathered} \left. \begin{gathered} {{G}_{c}} \hfill \\ {{G}_{s}} \hfill \\ \end{gathered} \right\} = {{K}_{1}}\Omega \sum\limits_{v = 1}^L {{{p}_{v}}} \int\limits_{(v - 1){{\tau }_{{\text{и}}}}}^{v{{\tau }_{{\text{и}}}}} {dt} \left\{ \begin{gathered} \cos {{w}_{0}}t \hfill \\ \sin {{w}_{0}}t \hfill \\ \end{gathered} \right\} \times \\ \times \,\,\int\limits_0^t {n(\tau )} \exp \left[ { - \Omega (t - \tau )} \right]\cos \left[ {{{w}_{0}}(t - \tau )} \right]d\tau , \\ \end{gathered} $
среднее значение которых равно нулю, а дисперсия определяется как

(19)
$\begin{gathered} \sigma _{G}^{2} = {{\left( {\frac{{{{K}_{1}}\Omega }}{4}} \right)}^{2}}\sum\limits_{k,v = 1}^L {{{p}_{k}}{{p}_{v}}} \int\limits_{(k - 1){{\tau }_{{\text{и}}}}}^{k{{\tau }_{{\text{и}}}}} {d{{t}_{1}}\exp \left[ { - \Omega {{t}_{1}}} \right]} \times \\ \times \,\,\int\limits_{(v - 1){{\tau }_{{\text{и}}}}}^{v{{\tau }_{{\text{и}}}}} {d{{t}_{2}}} \exp \left[ { - \Omega {{t}_{2}}} \right]\int\limits_0^{{{t}_{1}}} {d{{\tau }_{1}}} \int\limits_0^{{{t}_{2}}} {d{{\tau }_{2}}} \left\langle {n({{\tau }_{1}})n({{\tau }_{2}})} \right\rangle \times \\ \times \,\,\exp \left[ {\Omega \left( {{{\tau }_{1}} + {{\tau }_{2}}} \right)} \right]\cos {{w}_{0}}{{\tau }_{1}}\cos {{w}_{0}}{{\tau }_{2}}. \\ \end{gathered} $

С учетом [3] корреляционную функцию $S(\varepsilon )$ в области высокой корреляции можно записать в виде

(20)
$S(\varepsilon ) = \frac{{{{K}_{0}}{{K}_{1}}L{{\tau }_{{\text{и}}}}}}{{4\pi }}\hat {S}(\varepsilon ),$
$\hat {S}(\varepsilon ) = \left\{ \begin{gathered} \left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}} - \frac{1}{c}} \right){{Y}_{1}} + \frac{{\cos b\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right) - 1}}{b} + \frac{{\exp \left\{ { - с\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} \right\}}}{c}{{Y}_{2}},\,\,\,\,0 < \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}} < 1; \hfill \\ \left( {1 + \frac{{еxp\left\{ { - с\left( {2 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} \right\} - еxp\left\{ { - 2с} \right\} - еxp\left\{ {\frac{{с\varepsilon }}{{{{\tau }_{{\text{и}}}}}}} \right\}}}{c}} \right){{Y}_{3}} - \left( {\frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}} + \frac{1}{c}} \right){{Y}_{4}} + \frac{{еxp\left\{ { - с\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} \right\}}}{c}{{Y}_{5}} - \hfill \\ - \,\,\frac{{еxp\left\{ {\frac{{с\varepsilon }}{{{{\tau }_{{\text{и}}}}}}} \right\}}}{c}{{Y}_{6}} + \left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right){{Y}_{7}} - {{Y}_{8}} + \left( {1 - \frac{1}{c}} \right){{Y}_{9}} + \frac{{еxp\left\{ { - с\left( {2 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} \right\}}}{c}{{Y}_{{10}}} - \frac{{еxp\left\{ { - 2с} \right\}}}{c}{{Y}_{{11}}} - \hfill \\ - \frac{{еxp\left\{ { - с\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} \right\}}}{c}{{Y}_{{12}}} + \frac{{еxp\left\{ { - с} \right\}}}{c}{{Y}_{{13}}} + \frac{1}{c}{{Y}_{{14}}} + \hfill \\ + \frac{{2\cos b - \cos \frac{{b\varepsilon }}{{{{\tau }_{{\text{и}}}}}} + \cos b\left( {2 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right) - \cos 2b - \cos b\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)}}{b},\,\,\,\, - 1 < \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}} < 0; \hfill \\ \end{gathered} \right.$
где $2b = \Delta w{{\tau }_{{\text{и}}}} = 2\pi \Delta f{{\tau }_{{\text{и}}}}$ – произведение длительности элемента кода ПСП на ширину спектра ФМШПС на выходе передатчика; 2c = 2Ωτи = $ = 4\pi \Delta {{f}_{{{\text{пр}}}}}{{\tau }_{{\text{и}}}}$ – произведение длительности элемента кода ПСП на полосу пропускания одноконтурной входной цепи приемника;

$\begin{gathered} {{Y}_{1}} = \int\limits_0^{b\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} {\operatorname{sinc} xdx} ,\,\,\,\,{{Y}_{2}} = \int\limits_0^{b\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} {\operatorname{sinc} xеxp\left( {\frac{{cx}}{b}} \right)dx} , \\ {{Y}_{3}} = \int\limits_0^b {\operatorname{sinc} xdx} ,\,\,\,{{Y}_{4}} = \int\limits_0^{\frac{{b\varepsilon }}{{{{\tau }_{{\text{и}}}}}}} {\operatorname{sinc} xdx} , \\ {{Y}_{5}} = \int\limits_0^b {\operatorname{sinc} xеxp\left( {\frac{{cx}}{b}} \right)dx} ,\,\,\,{{Y}_{6}} = \int\limits_0^{ - \frac{{b\varepsilon }}{{{{\tau }_{{\text{и}}}}}}} {\operatorname{sinc} xеxp\left( {\frac{{cx}}{b}} \right)dx} , \\ {{Y}_{7}} = \int\limits_{b\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)}^{b\left( {2 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} {\operatorname{sinc} xdx} ,\,\,\,\,{{Y}_{8}} = \int\limits_b^{2b} {\operatorname{sinc} xdx} , \\ {{Y}_{9}} = \int\limits_{2b}^{b\left( {2 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} {\operatorname{sinc} xdx} ,\,\,\,\,{{Y}_{{10}}} = \int\limits_0^{b\left( {2 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} {\operatorname{sinc} xеxp\left( {\frac{{cx}}{b}} \right)dx} , \\ {{Y}_{{11}}} = \int\limits_0^{2b} {\operatorname{sinc} xеxp\left( {\frac{{cx}}{b}} \right)dx} , \\ {{Y}_{{12}}} = \int\limits_{ - b}^{b\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} {\operatorname{sinc} xеxp\left( {\frac{{cx}}{b}} \right)dx} , \\ {{Y}_{{13}}} = \int\limits_{ - b}^b {\operatorname{sinc} xеxp\left( {\frac{{cx}}{b}} \right)dx} ,\,\,\,\,{{Y}_{{14}}} = \int\limits_b^{b\left( {1 - \frac{\varepsilon }{{{{\tau }_{{\text{и}}}}}}} \right)} {\operatorname{sinc} xdx} . \\ \end{gathered} $

После замены порядка интегрирования в (19), вычисления интегралов и пренебрежения слагаемыми, имеющими порядок малости $0\left( {{1 \mathord{\left/ {\vphantom {1 L}} \right. \kern-0em} L}} \right)$, выражение (19) примет вид

(21)
$\sigma _{G}^{2} = \frac{{{{N}_{0}}K_{1}^{2}L{{\tau }_{{\text{и}}}}}}{{16}}\hat {\sigma }_{G}^{2},$
где $\hat {\sigma }_{G}^{2} = 1 - \frac{{1 - еxp\left( { - c} \right)}}{c}$.

2. ВЕРОЯТНОСТЬ ОБНАРУЖЕНИЯ ФМШПС КВАДРАТУРНЫМ ПРИЕМНИКОМ

Определим вероятность обнаружения ФМШПС как

(22)
$D = \int\limits_{{{z}_{{{\text{пор}}}}}}^\infty {W(z)dz} ,$
$W(z)$ – плотность распределения вероятностей выходного сигнала корреляционного приемника $z = M(\varepsilon )$; ${{z}_{{{\text{пор}}}}}$ – уровень “порога”, при превышении которого выносится решение о наличии ФМШПС на входе приемника.

В соответствии с [5] выражение для плотности распределения сигнала со случайной начальной фазой, равномерно распределенной на $\left[ { - \pi ;\pi } \right]$, и амплитудой, содержащей постоянную и флуктуирующую составляющие с совместным распределением

(23)
$w(a,{{\varphi }_{0}}) = \frac{a}{{2\pi \sigma _{a}^{2}}}\exp \left[ { - \frac{{{{a}^{2}} + a_{0}^{2} - 2a{{a}_{0}}\cos {{\varphi }_{0}}}}{{2\sigma _{a}^{2}}}} \right],$
на выходе корреляционного приемника будет иметь вид
(24)
$\begin{gathered} W(z) = \frac{1}{{2\sigma _{G}^{2}\left[ {1 + QH(\gamma ,b,c)} \right]}} \times \\ \times \,\,\exp \left\{ { - \frac{{z + \sigma _{G}^{2}{{\eta }^{2}}QH(\gamma ,b,c)}}{{2\sigma _{G}^{2}\left[ {1 + QH(\gamma ,b,c)} \right]}}} \right\} \times \\ \times \,\,{{I}_{0}}\left[ {\sqrt {\frac{{z{{\eta }^{2}}QH(\gamma ,b,c)}}{{\sigma _{G}^{2}\left\{ {1 + QH(\gamma ,b,c)} \right\}}}} } \right], \\ \end{gathered} $
где ${{I}_{0}}(x)$ – функция Бесселя нулевого порядка от мнимого аргумента; $Q = {{\sigma _{а}^{2}L{{\tau }_{{\text{и}}}}} \mathord{\left/ {\vphantom {{\sigma _{а}^{2}L{{\tau }_{{\text{и}}}}} {{{N}_{0}}}}} \right. \kern-0em} {{{N}_{0}}}}$ – отношение сигнал/шум для флуктуирующей составляющей амплитуды; $\eta = {{{{a}_{0}}} \mathord{\left/ {\vphantom {{{{a}_{0}}} {{{\sigma }_{a}}}}} \right. \kern-0em} {{{\sigma }_{a}}}}$ – отношение постоянной составляющей амплитуды к среднеквадратичному отклонению флуктуирующей составляющей;
(25)
$\begin{gathered} H(\gamma ,b,c) = \frac{1}{{{{\pi }^{2}}\hat {\sigma }_{G}^{2}}} \times \\ \times \left\{ \begin{gathered} \left[ {\left( {1 - \gamma - \frac{1}{c}} \right){{Y}_{1}} + \frac{{\cos b(1 - \gamma ) - 1}}{b}} \right. + {{\left. {\frac{{\exp \left\{ { - c(1 - \gamma )} \right\}}}{c}{{Y}_{2}}} \right]}^{2}},\,\,\,\,{\text{ }}0 \leqslant \gamma < 1; \hfill \\ \left[ {\left( {1 + \frac{{еxp\left\{ { - c(2 - \gamma )} \right\} - еxp\left\{ { - 2c} \right\} - еxp\left\{ {c\gamma } \right\}}}{c}} \right){{Y}_{3}} - \left( {\gamma + \frac{1}{c}} \right){{Y}_{4}} + \frac{{еxp\left\{ { - c(1 - \gamma )} \right\}}}{c}{{Y}_{5}} - \frac{{еxp\left\{ {c\gamma } \right\}}}{c}{{Y}_{6}} + } \right. \hfill \\ + \,\,(1 - \gamma ){{Y}_{7}} - {{Y}_{8}} + \left( {1 - \frac{1}{c}} \right){{Y}_{9}} + \frac{{еxp\left\{ { - c(2 - \gamma )} \right\}}}{c}{{Y}_{{10}}} - \frac{{еxp\left\{ { - 2c} \right\}}}{c}{{Y}_{{11}}} - \frac{{еxp\left\{ { - c(1 - \gamma )} \right\}}}{c}{{Y}_{{12}}} + \hfill \\ + \,\,\frac{{еxp\left\{ { - c} \right\}}}{c}{{Y}_{{13}}} + \frac{{{{Y}_{{14}}}}}{c} + {{\left. {\frac{{2\cos b - \cos b\gamma + \cos b(2 - \gamma ) - \cos 2b - \cos b(1 - \gamma )}}{b}} \right]}^{2}},\,\,\,\, - 1 < \gamma < 0; \hfill \\ \end{gathered} \right. \\ \end{gathered} $
– коэффициент, зависящий от параметров одноконтурной входной цепи приемника и выходной цепи передатчика; $\gamma = {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{\tau }_{{\text{и}}}}}}} \right. \kern-0em} {{{\tau }_{{\text{и}}}}}}$ – время задержки принимаемого ФМШПС относительно опорного, нормированное на длительность элемента кода ПСП; ${{x}_{{{\text{пор}}}}} = {{{{z}_{{{\text{пор}}}}}} \mathord{\left/ {\vphantom {{{{z}_{{{\text{пор}}}}}} {2\sigma _{G}^{2}}}} \right. \kern-0em} {2\sigma _{G}^{2}}}$ – уровень “порога”, нормированный на мощность шумовой составляющей, связанный с вероятностью ложной тревоги как

${{F}_{{{\text{л}}{\text{.т}}{\text{.}}}}} = \exp \left[ { - \frac{{{{z}_{{{\text{пор}}}}}}}{{2\sigma _{G}^{2}}}} \right].$

С учетом (24), (25), выражение (22) для вероятности обнаружения ФМШПС принимает вид

(26)
$\begin{gathered} D(\gamma ,b,c,Q,\eta ,{{F}_{{{\text{л}}{\text{.т}}{\text{.}}}}}) = 1 - \frac{1}{{1 + QH(\gamma ,b,c)}} \times \\ \times \,\,\exp \left\{ { - \frac{{{{\eta }^{2}}QH(\gamma ,b,c)}}{{2\left[ {1 + QH(\gamma ,b,c)} \right]}}} \right\} \times \\ \times \int\limits_0^{ - \ln {{F}_{{{\text{л}}{\text{.т}}{\text{.}}}}}} {\exp \left\{ { - \frac{x}{{1 + QH(\gamma ,b,c)}}} \right\}{{I}_{0}}\left[ {\sqrt {\frac{{2x{{\eta }^{2}}QH(\gamma ,b,c)}}{{1 + QH(\gamma ,b,c)}}} } \right]} \,dx. \\ \end{gathered} $

Для каждого из значений, характеризующих ограниченный спектр ФМШПС на выходе передатчика с параметрами: $b = 1.8\pi $, $2.0\pi $, $2.5\pi $, $3.0\pi $, $3.5\pi $, $4.0\pi $, $4.5\pi $, $5.0\pi $, отношения сигнал/шум $Q = 10$, $30$, 50 для флуктуирующей составляющей амплитуды и отношения $\eta = 0.3$, $1.0$, $3.0$ для постоянной составляющей амплитуды к среднеквадратичному отклонению ее флуктуирующей составляющей, построены семейства вероятностей обнаружения ФМШПС (26) как функции от нормированной на длительность элемента кода ПСП временной задержки между принимаемым и опорным сигналом $\gamma = {\varepsilon \mathord{\left/ {\vphantom {\varepsilon {{{\tau }_{{\text{и}}}}}}} \right. \kern-0em} {{{\tau }_{{\text{и}}}}}}$ для параметров одноконтурной входной цепи приемника при $c = 4\pi $, $6\pi $, $8\pi $, $10\pi $, $12\pi $. В качестве минимального значения параметра входной цепи приемника $c$ принято $4\pi $, при котором искажения, вызванные неравномерностью ее АЧХ, не приводят к существенному искажению передаваемой информации. Для каждого графика из семейства характеристик обнаружения определены его наибольшие значения. Из заданного множества наибольших значений семейства кривых выделены максимальные, характеризуемые временной задержкой между принимаемым и опорным сигналами и полосой пропускания входной цепи приемника. Результаты расчетов представлены в табл. 1 для ${{F}_{{{\text{л}}{\text{.т}}}}} = {{10}^{{ - 2}}}$ и в табл. 2 для ${{F}_{{{\text{л}}{\text{.т}}}}} = {{10}^{{ - 4}}}$.

Таблица 1.

Параметры рабочих частот передатчика, одноконтурной входной цепи приемника и максимальные значения вероятности обнаружения ФМШПС для заданных отношений сигнал/шум регулярной и флуктуирующей составляющих амплитуды и вероятности ложной тревоги ${{F}_{{{\text{л}}{\text{.т}}}}} = {{10}^{{ - 2}}}$

$b = \pi \Delta f{{\tau }_{{\text{и}}}}$ $c = 2\pi \Delta {{f}_{{{\text{пр}}}}}{{\tau }_{{\text{и}}}}$ ${{\gamma }_{{\max }}} = \frac{{{{\varepsilon }_{{\max }}}}}{{{{\tau }_{{\text{и}}}}}}$ ${{D}_{{\max }}}$
η
0.3 1.0 3.0
Q = 10
1.8π 4 π –0.1322 0.238 0.337 0.843
2.0π 4 π –0.123 0.238 0.336 0.842
2.5 π 4 π –0.112 0.24 0.339 0.845
3.0 π 12 π –0.0419 0.244 0.345 0.849
3.5 π 8 π –0.0659 0.256 0.358 0.859
4.0 π 6 π –0.0733 0.257 0.36 0.861
4.5 π 6 π –0.0697 0.258 0.361 0.862
5.0 π 12 π –0.0387 0.259 0.363 0.863
      Q = 30
1.8π 4 π –0.1322 0.547 0.659 0.975
2.0π 4 π –0.123 0.546 0.658 0.975
2.5 π 4 π –0.112 0.549 0.661 0.975
3.0 π 12 π –0.0419 0.554 0.666 0.976
3.5 π 8 π –0.0659 0.568 0.678 0.978
4.0 π 6 π –0.0733 0.57 0.679 0.978
4.5 π 6 π –0.0697 0.571 0.68 0.978
5.0 π 12 π –0.0387 0.572 0.681 0.979
      Q = 50
1.8π 4 π –0.1322 0.683 0.774 0.989
2.0π 4 π –0.123 0.682 0.773 0.989
2.5 π 4 π –0.112 0.684 0.775 0.989
3.0 π 12 π –0.0419 0.689 0.779 0.99
3.5 π 8 π –0.0659 0.7 0.788 0.99
4.0 π 6 π –0.0733 0.701 0.789 0.991
4.5 π 6 π –0.0697 0.702 0.789 0.991
5.0 π 12 π –0.0387 0.703 0.79 0.991
Таблица 2.

Параметры рабочих частот передатчика, одноконтурной входной цепи приемника и максимальные значения вероятности обнаружения ФМШПС для заданных отношений сигнал/шум регулярной и флуктуирующей составляющих амплитуды и вероятности ложной тревоги ${{F}_{{{\text{л}}{\text{.т}}}}} = {{10}^{{ - 4}}}$

$b = \pi \Delta f{{\tau }_{{\text{и}}}}$ $c = 2\pi \Delta {{f}_{{{\text{пр}}}}}{{\tau }_{{\text{и}}}}$ ${{\gamma }_{{\max }}} = \frac{{{{\varepsilon }_{{\max }}}}}{{{{\tau }_{{\text{и}}}}}}$ ${{D}_{{\max }}}$
η
0.3 1.0 3.0
Q = 10
1.8π 4 π –0.1322 0.057 0.108 0.599
2.0π 4 π –0.123 0.056 0.107 0.597
2.5 π 4 π –0.112 0.058 0.109 0.603
3.0 π 12 π –0.0419 0.06 0.113 0.611
3.5 π 8 π –0.0659 0.065 0.122 0.633
4.0 π 6 π –0.0733 0.066 0.124 0.636
4.5 π 6 π –0.0697 0.067 0.124 0.638
5.0 π 12 π –0.0387 0.067 0.125 0.64
      Q = 30
1.8π 4 π –0.1322 0.299 0.428 0.927
2.0π 4 π –0.123 0.298 0.427 0.926
2.5 π 4 π –0.112 0.302 0.43 0.928
3.0 π 12 π –0.0419 0.307 0.436 0.93
3.5 π 8 π –0.0659 0.322 0.453 0.935
4.0 π 6 π –0.0733 0.324 0.455 0.936
4.5 π 6 π –0.0697 0.325 0.456 0.936
5.0 π 12 π –0.0387 0.327 0.457 0.937
      Q = 50
1.8π 4 π –0.1322 0.466 0.595 0.97
2.0π 4 π –0.123 0.465 0.593 0.97
2.5 π 4 π –0.112 0.468 0.597 0.97
3.0 π 12 π –0.0419 0.474 0.602 0.971
3.5 π 8 π –0.0659 0.49 0.616 0.973
4.0 π 6 π –0.0733 0.492 0.618 0.974
4.5 π 6 π –0.0697 0.493 0.619 0.974
5.0 π 12 π –0.0387 0.494 0.62 0.974

В таблицах приведены значения ширины спектра ФМШПС в выходной цепи передатчика $b$, нормированные на длительность элемента кода ПСП; значение полосы пропускания одноконтурной входной цепи приемника $c$, нормированное на длительность элемента кода ПСП, при котором функция вероятности обнаружения $D(\gamma )$ располагается выше других из заданного семейства кривых при b = const; значение нормированной на длительность элемента кода ПСП временной задержки между принимаемым ФМШПС и опорным ${{\gamma }_{{\max }}} = {{{{\varepsilon }_{{\max }}}} \mathord{\left/ {\vphantom {{{{\varepsilon }_{{\max }}}} {{{\tau }_{{\text{и}}}}}}} \right. \kern-0em} {{{\tau }_{{\text{и}}}}}}$, при котором наибольшее значение вероятности обнаружения семейства кривых максимально. Так, например, для рабочей полосы частот передатчика $b = \pi \Delta f{{\tau }_{{\text{и}}}} = 1.8\pi $ вероятность обнаружения достигает максимума ${{D}_{{\max }}}$, если принимаемый сигнал опережает опорный на $\left| \varepsilon \right| = 0.1322{{\tau }_{{\text{и}}}}$ и параметр полосы пропускания приемника $c = 2\pi \Delta {{f}_{{{\text{пр}}}}}{{\tau }_{{\text{и}}}}$ при этом равен $4\pi $.

В табл. 1 и 2 для заданных значений $Q$ и $\eta $ рассчитаны максимальные вероятности обнаружения ${{D}_{{\max }}}$ при указанных параметрах рабочих частот передатчика и одноконтурной цепи приемника. Так, например, если параметр рабочих частот передатчика $b = 1.8\pi $, то параметр одноконтурной входной цепи приемника, при котором вероятность обнаружения максимальна, равен $c = 4\pi $, при этом принимаемый сигнал опережает опорный на $\left| \varepsilon \right| = 0.1322{{\tau }_{{\text{и}}}}$. Для отношения сигнал/шум флуктуирующей составляющей амплитуды $Q = 10$ и постоянной составляющей амплитуды к среднеквадратичному отклонению ее флуктуирующей составляющей $\eta = 0.3$, вероятности ложной тревоги ${{F}_{{{\text{л}}{\text{.т}}}}} = {{10}^{{ - 2}}}$, максимальная величина вероятности обнаружения составит ${{D}_{{\max }}} = 0.238$. А при $Q = 50$, $\eta = 0.3$ и неизменных остальных параметрах вероятность обнаружения составит ${{D}_{{\max }}} = 0.683$. При увеличении отношения регулярной составляющей амплитуды к среднеквадратичному отклонению до $\eta = 3$ и неизменными другими параметрами, максимумы вероятностей правильного обнаружения будут составлять ${{D}_{{\max }}} = 0.843$ при $Q = 10$ и ${{D}_{{\max }}} = 0.989$ при $Q = 50$. При уменьшении вероятности ложной тревоги ${{F}_{{{\text{л}}{\text{.т}}}}} = {{10}^{{ - 4}}}$ (см. табл. 2) максимум вероятности обнаружения при тех же параметрах частот цепей передатчика и приемника составит: ${{D}_{{\max }}} = 0.057$ для $Q = 10$ и $\eta = 0.3$; ${{D}_{{\max }}} = 0.466$ для $Q = 50$ и $\eta = 0.3$; ${{D}_{{\max }}} = 0.599$ для $Q = 10$ и $\eta = 3$; ${{D}_{{\max }}} = 0.97$ для $Q = 50$ и $\eta = 3$.

С расширением полосы частот выходной цепи передатчика наибольшее значение функции вероятности обнаружения ФМШПС увеличивается, также увеличивается максимальное из всех наибольших значений вероятностей обнаружения заданного семейства кривых, абсолютная величина временной задержки между принимаемым и опорным сигналами (для случая опережения принимаемым ФМШПС опорного) при этом уменьшается, параметр избирательной цепи приемника, при котором достигается максимум, увеличивается. Если параметр, пропорциональный полосе рабочих частот передатчика, увеличивается от $b = 1.8\pi $ до $b = 5.0\pi $, то максимальное значение вероятности обнаружения увеличивается на $0.4...17{\text{ }}\% $, в пределах заданных значений $Q$, $\eta $ и ${{F}_{{{\text{л}}{\text{.т}}}}}$. При ${{F}_{{{\text{л}}{\text{.т}}}}} = {{10}^{{ - 2}}}$, $Q = 10$, $\eta = 0.3$ и изменении параметра полосы частот выходной цепи передатчика $b$ от $1.8\pi $ до $5.0\pi $ максимум вероятности обнаружения увеличивается на 8%, а при ${{F}_{{{\text{л}}{\text{.т}}}}} = {{10}^{{ - 4}}}$ – на 17%. При увеличении отношений сигнал/шум для флуктуирующей и постоянной составляющих амплитуды максимум вероятности обнаружения ФМШПС в заданных пределах рабочих частот передатчика изменяется на $0.4...7{\text{ }}\% $.

Параметр одноконтурной входной цепи приемника, при котором наибольшее значение вероятности обнаружения из заданного семейства кривых максимально, с увеличением полосы пропускания рабочих частот передатчика увеличивается на интервалах $b = 1.8\pi ...3.0\pi $ и $b = 4.5\pi ...5.0\pi $.

ЗАКЛЮЧЕНИЕ

Вероятность обнаружения ФМШПС со случайной начальной фазой и амплитудой, содержащей регулярную и флуктуирующую составляющие, определяется характером и параметрами рабочей полосы частот передатчика и одноконтурной входной цепи приемника, искажающими принимаемый сигнал, а ее максимум достигается при конечном времени опережения принимаемого сигнала относительно опорного. Расширение рабочей полосы частот выходной цепи передатчика увеличивает вероятность его правильного обнаружения. Увеличение полосы пропускания входной цепи приемника приводит, с одной стороны, к увеличению энергии принимаемого сигнала, с другой – к увеличению энергии шума, а также к уменьшению абсолютного значения временной задержки между принимаемым и опорным сигналами, соответствующей максимуму вероятности обнаружения.

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

Список литературы

  1. Борисов В.И., Зинчук В.М., Лимарев А.Е. и др. Помехозащищенность систем радиосвязи с расширением спектра сигналов модуляцией несущей псевдослучайной последовательностью. М.: Радио и связь, 2003.

  2. Деев В.В. Методы модуляции и кодирования в современных системах связи. М.: Наука, 2007.

  3. Бакаева Г.А. // Радиолокация, навигация, связь: труды XXII международной научно-технич. конф., 15–17 апреля 2016 г. Воронеж, 2016. Т. 1. С. 594.

  4. Буга Н.Н., Фалько А.И., Чистяков Н.И. Радиоприемные устройства. М.: Радио и связь, 1986.

  5. Нахмансон Г.С., Бакаева Г.А. // Теория и техника радиосвязи. 2008. № 1. С. 26.

Дополнительные материалы отсутствуют.