Радіотехнічні кола та сигнали
УДК 621.391.17
ПРОЦЕДУРА КОГЕРЕНТНО-НЕКОГЕРЕНТНОЇ ДЕМОДУЛЯЦІЇ ВЗАЄМНОЗАВАЖАЮЧИХ ЦИФРОВИХ СИГНАЛІВ З ДВІЙКОВОЮ
ЧАСТОТНОЮ МОДУЛЯЦІЄЮ
Єрохін В. Ф., д.т.н., професор; Пелешок Є. В.
Інститут спеціального зв ’язку та захисту інформації Національного технічного університету України «Київський політехнічний інститут»,
Київ, Україна
PROCEDURE OF COHERENT-INCOHERENT DEMODULATION OF DIGITAL SIGNALS WITH BINARY FREQUENCY MODULATION
Yerokhin V., Doc. Of Sci (Technics), Professor; Peleshok Y.
Institute of special communications and information protection National Technical University of Ukraine «Kyiv Polytechnic Institute», Kyiv, Ukraine
Вступ
В сучасних умовах прийом радіосигналів здійснюється, як правило, в апріорно невизначеній сигнально-завадовій обстановці, яка обумовлена обмеженістю радіочастотного ресурсу та зростанням кількості та потужності структурних випромінювань різноманітного походження. В загальному випадку на вхід радіоприймального пристрою, крім корисних сигналів, потрапляють завади, які разом з внутрішніми шумами приймача можуть суттєво впливати на якість прийняття сигналу. Тому проблема демодуляції сигналів в умовах впливу завад була і залишається актуальною, а на її вирішення спрямована велика кількість робіт [1].
Окремим важливим класом задач завадозахищеного прийому сигналів є задачі, що пов’язані з прийомом корисного сигналу в умовах впливу потужних синхронних та подібних за структурою завад, що є одними із найне-безпечніших при прийомі [1-3].
У даній роботі пропонується для підвищення завадозахищеності прийому корисного сигналу, що спостерігається на фоні подібної потужної завади, використовувати в демодуляторах приймальних пристроїв компенсаційні процедури із застосуванням когерентної (квазікогерентної) обробки завади та некогерентної обробки корисного сигналу [2,4].
Метою і основним змістом статті є вирішення задачі синтезу процедури когерентно-некогерентної демодуляції взаємозаважаючих цифрових сигналів двостанової частотної модуляції (ЧМ-2) та визначення її потенційної завадозахищеності. Для досягнення поставленої мети доопрацюємо і використаємо методику, що приведена в [2, 4]. *
Вісник Національного технічного університету України «КПІ» 23
Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53
Радіотехнічні кола та сигнали
Методика синтезу процедури когерентно-некогерентного детектування корисного ЧМ-2 сигналу, що спостерігається на фоні потужної та подібної за структурою ЧМ-2 завади
Вважатимемо, що частотні позиції і тактові точки сигналу та завади співпадають, а маніпуляція завади на кожній із двох частотних позицій здійснюється без розриву фази. Остання умова дає можливість використовувати когерентну (квазікогерентну) обробку завади, а корисний сигнал будемо обробляти некогерентно (квадратурно).
Модель спостереження на тривалості тактового інтервалу представимо наступним чином:
У (t ) = s1 ( r1> Фіс >Ф2с > t) + s2 ( r2 >фіз ^2з > t) + n (t ) = r1
^j^cos ( w1t + ф1з )-
-ЛІsin(w1t + ф1з) +(1 -r1) A2cos(w2t + ф2з)-Л^їп(w2t + ф2з)
+
+r2 A21COS (Ю/ + ф1з ) + (1 - Г2 ) A22COS
(®2t + ф2з ) + n (t) >
(1)
де Л1,2 — амплітуди синфазних та квадратурних складових сигналу; Л21, Л22 — амплітуда завади на частотах ш1 та ю2 відповідно; n(t) — адитивний білий гаусівський шум (АБГШ).
В свою чергу
Л1 = Aocos (ф1с- фіз); Лk = Aosin (ф1с- фіз );
Л2 = ACOS (ф2с - ф2з ) ; Л2 = Л0sin ( ф2с - ф2з ) > (2)
де Л0 — амплітуда корисного сигналу, незмінна за частотою.
Будемо також вважати, що стани дискретних параметрів r1 та r2 рівно-імовірні та взаємно незалежні, а початкові фази ф12 с з рівномірно розподілені на інтервалі [0,2ті].
Запишемо умовний функціонал правдоподібності для спостереження (1):
Л(г1,г2,/ф1с> ф2с) = exP
N
1
tk
2 j У (*)• S1 (r1> фіс> ф2с> t)dt
tk-1
+
Ik tk tIk
+2 j У (t )• s2 ( r2 5 ф1з ’ф2з >1) dt - j s12 (r1^1c ’ф2е j ) dt - j s2 ( r2 >фіз ^2з d ) dt -
tk-1 tk-1 tk-1
h
-2 j s1 ( r15 ф1с ’ф2с >1)■ s2 ( r2 5ф1з 5фіз >t) dt
h-1 24
(3)
24 Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53
Радіотехнічні кола та сигнали
З урахуванням моделі спостереження (1) і квадратурного представлення сигналу (2) введемо позначення для запису функціоналу (3) у явному вигляді:
tk
bi = 7Г І у(t)A)cos(фіс-Фіз)■cosЦ*+Ф1з)dt = b1Sos(фю-Фь);
N
0 tk-і
2 tk
b =-— І У(t) A)sin(Фіс - Фіз)■ sin(®1t + Ф1з)dt = -bf0sin(Фк - Ф1з);
N tk-1 2 Ь
b2 =— І У (t) A)c0s ( Ф2с - Ф2з )■ c0s ( ®2t + Ф2з ) dt = b2°c0s ( Ф2о - Ф2з ) ;
No tk-i
2 tk
b2 = -~— І У (t) A)sin (Ф2с - Ф2з )■ sin ( ®2t + Ф2з ) dt = -b20sin (Ф2о - Ф2з );
N t
0 tk-i
2 tk 2 tk
bi = 3^ І У (t) A2ic0s ( ® 2* + Фіз ) dt; b2 = ^ І У (t) A22c0s ( ®2t + Ф2з ) dt •
N) tk-i N) tk-i
i i Ik
h1 =— j 4cos2 ( Фіс - Фіз )■ cos2 (®^ + Ф1з ) dt + — І І12 (Фю - Ф1з )
N 0 tk.і N t
xsin2 (®it + ф1з) dt = his + h^k;
x
1 tk 1 1
h2 =— j A2c0s2 (Ф2с - Ф2з )■ COs2 (®2t + Ф2з ) dt + — І І1 І - Ф2з )
N0 1-і to t
x
xsin2 (®2it + Ф2з ) dt = h^s + h2k;
J tk
h1222)=— І A221,22c0s2 (®1,2t + Фіз,2з ) dt;
N a +
0 tk-i
i tk
R1 =— І A2iA)C0s ( Фіс - Фіз )■ cos2 ( ®1t + Ф1з ) dt = R10 C0S ( Ф1е - Ф1з ) ;
*
0 tk-i
tk
1 lk
R2 =— І A22 A0C0s (Ф2с - Ф2з )■ COs2 ( ®2t + Ф2з ) dt = R^C0S (ф2о - Ф2з ) (4)
N f\ ,
0 tk-1
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53
25
Радіотехнічні кола та сигнали
Зазначимо, що величини h22, h2(2) інваріантні до значень початкових фаз ф1с 2с та ф1з 2з. Крім того, при запропонованій вже відмові від оцінювання амплітуди А0 корисного сигналу очевидне рівняння h1 = h2. Однак припущення про відмову від оцінювання амплітуди корисного сигналу та заміни її величиною A0 << A21, A0 << A22 не дозволяє знехтувати тим, що в загальному випадку A21 ф A22, тому що можливий випадок, коли різниця IA21 A221 спів вимірна з Ao.
З урахуванням позначень (4), функціонал правдоподібності (3) може бути представлений наступним чином
Л(1 r2 1 Ф1с ’Ф2с ) = exP [r1 (Ь0cos (Ф1с - Ф1з ) - b1 0sin (Фіо - Ф1з )) + (! - r1) Х
0cOS(ф2с - Ф2з )) - b2k0sin(Ф2о - Ф2з ) + r2b1 + - r2 )b2 - r1h12 -(2 - r1)h2 -
-r2h12(2) - (1 - r2)h2(2) - 2r1r2Rl cOS(Ф1с - Ф1з ) - 2(1 - r1)(! - r2 )R2cos(Ф2о - Ф2з )] (5)
2
При квазікогерентній обробці завади та за умови h12 >> 1, похибками оцінок початкових фаз ф1з 2з можна знехтувати, тобто ф1з = ф2з = 0. Усереднивши (5) по ф1с 2с на інтервалі [0, 2л] отримаємо безумовний функціонал правдоподібності:
Л(г1 , r2 ) = exp( r'h' -2(1 r1 )h2 ) • exp(r2 (b1 - h12(2) )+ (1 - r2 )(b2 - h22(2) ))x
4л
r2 2n 2n
Х
j j exP r1 (^ГО^ - b1k^Шф^ ) +(1 - r1 )(Ь2°с^Ф2с - b2 ^Іпф2с ) -
0 0
2r1r2^с0^1с - 2(1 - r1 )(1 - r2 )^С0^2с ]^ф1с^ф2с •
(6)
У подальшому множником exp(- r1h12 - (1 - r2 )h12 )/ 2л , що не залежить від r1 та r2 знехтуємо. З урахуванням раніше введених позначень (4) запишемо:
h12(2) A0
R 0
R
A
21
h2(2) A0
A
a1h 2(2); = а^2(2).
(7)
22
Введемо позначення, що полегшують процедуру інтегрування (6): 26
26 Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53
Радіотехнічні кола та сигнали
Bi=$^Wf; в=$bFf^Wf
b
2
k 0
V1 = arctg ^ V 2 = arctg
b
k 0
b
b
s 0
(8)
звідки:
bis0 = B1cosv1; bik0 = B1sinv1; b‘2° = B2cosv 2; b20 = B2sinv 2.
(9)
Тоді (6) з урахуванням (7)-(9) перепишеться у вигляді:
Л(Г1> Г2 ) = ^ЄХР Г2 (Ь1 - Л2<2,) + (1 - Г2 )(b2 - *22
X
2n 2n
j j exp[r1B1 (c0s^1c0s91c -sin^1 sin91c ) + (1 - r1)B2 (cos^2 COS92c - (10)
0 0
- sin^2 Sin92c ) - 2r1r2a1h^2(2)c°s91c - 2(! - r1)(! - r2 )a2^22(2)c°s92c ]d91cd92c .
В загальному випадку правило прийняття рішення (ППР) для рівноімо-вірного дискретного параметру корисного сигналу r1 має вигляд:
r1* = rec і[Л(1, 0) + Л(1, 1) - Л(0, 1) - Л(0, 0)]; rect(x > 0)= 1; rect(x<0) = 0.
Або з урахуванням (10):
(ч 2п
b - ^2(2) ) j exp(B (cos^cos^ - sin^sin^ ))dtylc +
L 0
+exp(b1 - h2(2))- j exp(B1 (cos^1cos91c - sin^sinq^ )- 2a1h2(2)cosq>1c) dq>1c -
2n
exp(b2 - ^22(2)) j exp(B2 (cos^2c0s92c - sin¥2sin92c )- ^Ь^0^ )d^2c
. , 2n
-exp (b1 - h12(2}) j exp (B2 (cos^2COs92c - sin^2sin92c )) d92c
(11)
Еквівалентні та спрощуючі перетворення процедури когерентно-некогерентної демодуляції корисного ЧМ-2 сигналу в умовах адитивного впливу подібної потужної завади
Після заміни змінних інтегрування на £ = щ2 + ф1с 2с отримаємо [5]:
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53
27
Радіотехнічні кола та сигнали
Гу = rect
ьг Л 2)
2 п 2п
j e(Byc°*»di; - j e
2п ^b20-2«2h22|jcos92c-b20 sin92c
dp
2c
-e
bi-h12( 2))
2f e( B*o* )dg - 2fett111 °^ ^ ^
0 0 Введемо позначення аналогічно (8), (9):
(12)
b1 = b!° - 2<х,Л,2(2) ; Bie =т[ЇІ!УЦЇ0? ;
*20 = b|° - 2« 2 h2(2); B2e =^b20fJb0f ;
bk bk0
Лі = arctg-^-; Л2 = arctg 2
s0
le
s 0 2e
звідки
bse = Bye cosrii; b10 = В^іплі; b2e0 = B2ec0sr2 ; b20 = B2eSinr2 .
Після заміни змінних інтегрування у другому та четвертому інтегралі ППР (12) на Лі,2 + 9lc,2c отримаємо:
Л = rect[ exp(b2 - ^22(2))[/0(B1)-І0(B2e )]+ exP(bi - hi2(2))[/0(Bie)-I0(B2 )] ],
де 10 (• • ) — модифікована функція Бесселя нульового порядку. Можна побачити, що при h2(2) »i, h2(2) >> h22
(13)
exp(bi - h2(2)) exp(bi - h2(2))
r2 =i >>i; exp(b2 - h2(2)) r2 =0 = 0 ; exp(b2 - h2(2))
r2 =i
0;
r2 =0 >> i.
Тоді ППР (13) можна замінити асимптотично еквівалентним:
i = rect[ rect(b2 - h22(2))[l0(Bi )-І0(B2e )] + rect(bi - hi2(2))[l0(Bie )-I0(B2)] ]. (14)
*
Таким чином, наближена процедура (14) прийняття рішення r виявляється двоетапною, де на першому етапі приймається рішення про те, на якій з частот випромінюється завада s2 (r2, фь ,ф2з, t ). Якщо енергія завади
суттєво перевищує енергію корисного сигналу si (r, ф^ ,ф2c, t ), то пару правила прийняття рішення rect(bi2 - h2^), через малий вплив похибок на 28
28 Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53
Радіотехнічні кола та сигнали
. *
загальне рішення ri , слід замінити одним правилом прийняття рішення при когерентному (квазікогерентному) прийомі сигналу ЧМ-2:
r2* = rectfo - b2 ).
У результаті (14) перетворюється до вигляду:
r = rect[rect(bi - b2 M^i - В2е )+ rect(b2 - bi )-(віе - В2)], (15)
де враховано, що функція 10 (х) монотонна при x > 0.
Використовуючи наведені вище позначення, ППР (15) можна переписати наступним чином:
*
ri
= rect
(bf0 - 2rect(b1
b2К+ (bf0) -
- (b20 - 2rect(b2 - bi)a2Л22(2))2 - (b*0 f ]. (16)
У разі відсутності завади s2 (r2,фь,ф2з,t), тобто при h2(2) = 0, ППР
(13)—(16) вироджуються в класичні правила некогерентного прийому двійкового ЧМ сигналу.
Виконаємо якісне оцінювання завадостійкості отриманих ППР для асимптотичного випадку необмеженого збільшення середньої потужності завади s2(r2,фь,ф2з,t). Припускаючи, що похибки оцінювання неперервних (A2i, A22, фь ,ф2з) та дискретного параметрів завади будуть наближатися до нуля, отримаємо наступні вирази для bis 0, b2*0 (див. позначення (2), (4)):
bi
s0
Г =r2 =i
2_
N
t*
І Г(Af + A2i) cos (ait + ф1з) + n (t) • Af cos (a1t + ф1з)dt =
0 fk-i
= 2hiS + 2aihi2(2^ + пші;
b
s0
ri=r2 =0
2 tk Г / \ “I
= — I |_( A2 + A22 ) c0S ( ®2t + ф2з ) + n (t )]• A2 c0s ( ®2t + ф2з ) dt =
0 fk-i
= 2h2s + 2a2h2^ + пш2 . (17)
Із співставлення (16) та (17) видно, що при вищезазначених припущеннях (про відсутність похибок оцінювання параметрів завади) складові в
кореляційних інтегралах bi2, що породжуються її наявністю, повністю компенсуються. Шумові складові пш1 пш2, залишаються такими ж, як і для класичного випадку некогерентного прийому ЧМ-2 сигналу. Таким чином, потенційна завадостійкість алгоритму некогерентної демодуляції (13) ЧМ- *
Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53
29
Радіотехнічні кола та сигнали
2 сигналу за умови суттєвого перевищення середньої потужності подібної ЧМ-2 завади над потужністю корисного сигналу і відсутності похибок в оцінюванні її параметрів є такою ж, як і за її відсутності.
Висновки
Процедура когерентно-некогерентної демодуляції взаємозаважаючих цифрових сигналів двостанової частотної модуляції має ряд переваг:
- за умови суттєвого перевищення середньої потужності подібної ЧМ-2 завади над потужністю корисного ЧМ-2 сигналу та відсутності похибок в оцінці параметрів завади потенційна (гранична) завадостійкість процедури некогерентної демодуляції (13) є такою ж, як і за відсутності завади;
- компенсація здійснюється на виходах кореляційної згортки корисного сигналу, що є зручним з точки зору технічної реалізації;
- за відсутності завади дана процедура вироджується у класичну неко-герентну демодуляцію ЧМ-2 сигналу;
- дана процедура може використовуватися при реалізації програм повторного використання частотного ресурсу та при розробці перспективних завадозахищених засобів радіозв’язку.
Література
1. Бобровский В. И. Многопользовательское детектирование / В. И. Бобровский. — Ульяновск. : Вектор — С, 2007. — 348 с.
2. Єрохін В. Ф. Алгоритм демодуляції, що забезпечує повторне використання частот цифрового радіомовлення / В. Ф. Єрохін, І. М. Крутофіст // Захист інформації. — 2005. — № 25. — С. 42—47.
3. Бураченко Д. Л. Потенциальная помехоустойчивость разделения цифровых сигналов. Методика, программы, результаты расчетов / Д. Л. Бураченко, В. Ф. Ерохин, В.
О. Рашич — Л., 1987. — 122 с. — Деп. в ЦСИФ МО 04.03.87, № В-523.
4. Ерохин В. Ф. Демодуляция конфликтующих цифровых сигналов / В. Ф. Ерохин. — К. : КВИУС — ИК им. В. М. Глушкова АН Украины, 1993. — 132 с.
5. Прудников А. П. Интегралы и ряды: Элементарные функции / А. П. Прудников, Ю. А. Брычков, О. И. Маричев. — М. : Наука, 1981. — 797 с.
References
1. Bobrovskyi V. I. Mnohopolzovatelskoie detektirovanie. Ulianovsk, Vector - S, 2007, 348 p.
2. Yerokhin V. F., Krutofist I. M. Alhorytm demoduliatsuyi shcho zabezpechyie povtorne vykorystannia chastot radiomovlennia. Zakhyst informatsyi, 2005. No. 25, pp. 42-47.
3. Burachenko D. L., Yerokhin V. F., Rashych V. O. Potentsyalnaia pomek-houstoychivost razdelenyia tsyfrovykh sihnalov. Metodika, prohramy, rezultaty raschetov. L., 1987, 122 p. Dep. TSSIF MO 04.03.87 # V-523.
4. Yerokhin V. F. Demoduliatsyia konfliktuiushchikh tsyfrovykh sihnalov. K. : KVIUS -IK. im. V. M. Hlushkova AN Ukrainy, 1993, 132 p.
5. Prudnikov A. P., Brychkov Yu. A., Marychev O. A. Intehraly i riady: Elementarnyie funktsyi. M. : Nauka, 1981, 797 p. 30
30 Вісник Національного технічного університету України «КПІ» Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53
Радіотехнічні кола та сигнали
Єрохін В. Ф., Пелешок Є. В. Процедура когерентно-некогерентної демодуляції взаємнозаважаючих цифрових сигналів з двійковою частотною модуляцією. В
статті представлений синтез компенсаційної процедури когерентно-некогерентної демодуляції взаємозаважаючих цифрових сигналів з двійковою частотною модуляцією (ЧМ-2). Показано, що за умови квазікогерентного прийому подібної ЧМ-2 завади та суттєвого перевищення її середньої потужності над потужністю корисного ЧМ-2 сигналу, завадозахищеність прийому останнього наближається до завадозахищеності прийому в каналі без завади з адитивним білим гаусівським шумом.
Ключові слова: радіозв’язок, когерентно-некогерентна демодуляція, потенційна завадозахищеність, двійкова частотна модуляція.
Ерохин В. Ф., Пелешок Е. В. Процедура когерентно-некогерентной демодуляции взаимномешающих цифровых сигналов с двоичной частотной модуляцией. В статье представленный синтез компенсационной процедуры когерентно-некогерентной демодуляции взаимномешающих цифровых сигналов с двоичной частотной модуляцией (ЧМ-2). Показано, что при условии квазикогерентного приема подобной ЧМ-2 помехи и существенного превышения ее средней мощности над мощностью полезного ЧМ-2 сигнала, помехозащищенность приема последнего приближается к помехозащищенности приема в канале без помехи с аддитивным белым гауссовским шумом.
Ключевые слова: радиосвязь, когерентно-некогерентная демодуляция, потенциальная помехозащищенность, двоичная частотная модуляция.
Yerokhin V., Peleshok Y. Procedure of coherent-incoherent demodulation of digital signals with binary frequency modulation.
Introduction. The synthesis procedure of coherent-incoherent demodulation of digital signals with binary frequency modulation (F3E-2) is the basic and the purpose of this article.
Synthesis procedure method of coherent-incoherent demodulation of useful signal with F3E-2 which is observed with a powerful and alike hindrance with F3E-2. Coherent reception of similar hindrance with F3E-2 and incoherent reception of useful digital signal with F3E-2 in a channel with additive white noise are showed.
Simplifying transformations of coherent-incoherent demodulation procedure of useful digital signal with F3E-2 in the conditions of similar powerful hindrance. Simplification of procedure of coherent-incoherent demodulations of digital signals with F3E-2 is accomplished.
Conclusions. In absence of hindrance this procedure degenerates in classic incoherent demodulation of digital signal with F3E-2. This procedure can be used for programs realization of the repeated use offrequency resource and in development of perspective protection from hindrances of radio contact facilities.
Keywords: radio contact, coherent-incoherent demodulations, protecting from hindrances, binary frequency modulation.
Вісник Національного технічного університету України «КПІ» 31
Серія — Радіотехніка. Радіоапаратобудування. — 2013. — №53