Научная статья на тему 'Определение весовых коэффициентов двухступенчатого автокомпенсатора на основе LMS-алгоритма и корреляционных характеристик помех'

Определение весовых коэффициентов двухступенчатого автокомпенсатора на основе LMS-алгоритма и корреляционных характеристик помех Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
407
101
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
AUTOMATIC INTERFERENCE COMPENSATORS / METHOD OF LEAST SQUARES / LMS -ALGORITHM / АВТОКОМПЕНСАТОРЫ ПОМЕХ / МЕТОД НАИМЕНЬШИХ КВАДРАТОВ / LMS-АЛГОРИТМ / АВТОКОМПЕНСАТОРИ ЗАВАД / МЕТОД НАЙМЕНШИХ КВАДРАТіВ

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Жук С. Я., Семибаламут К. М.

Рассмотрены задачи определения весовых коэффициентов двухступенчатого автокомпенсатора c ортогонализацией сигналов части компенсационных каналов на основе LMS-алгоритма, а также корреляционных характеристик помех. Анализ эффективности рекуррентного двухступенчатого цифрового автокомпенсатора на основе LMS-алгоритма с одновременным функционированием ступеней выполнен с помощью статистического моделирования при разной обусловленности корреляционной матрицы помех и разном количестве источников помех.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Жук С. Я., Семибаламут К. М.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Definition of weight coefficient in two-stage automatic compensation based on LMS-algorithm and interference correlation characteristics

Introduction. Two-stage digital automatic compensators of active noise interference based on LS and RLS-algorithms provide parallel processing of signals and capability to connect/disconnect additional compensation channel block. LMS-algorithm is widely used for synthesis of digital automatic compensators of active noise interference. Its advantage is low computational cost. When solving certain practical problems correlation characteristics of interference signals are either given or their estimates are obtained. Problem statement. The first stage of the two-stage digital automatic compensator provides compensation of interference in the main channel using compensation channel signals of one of the blocks and orthogonalization of signals in the other block. The second stage provides compensation of the residue interference using the obtained orthogonalized signals. Modules of the first and the second stage are single type weight totalizers. We need to determine weight coeeficient of the two-stage digital automatic compensator with orthogonalization of signals of part of compensation channels based on LMS-algorithms and interference correlation characteristics. Using LMS-algorithm there was obtained two-stage digital automatic compensator with orthogonalization of part of compensation channels that is quasi-optimal regarding the criterion of minimum mean square error and it’s feature is lower computational complexity compared to two-stage digital automatic compensators based on LS and RLS-algorithms. It can operate in modes when preocessing of the complete training package of the input data is carried out in the first and then in the second stage as well as in the mode of simultaneous operation of the stages. There have been acheived analytical expressions for defining weight coefficients of modules in the first and in the second stage of the automatic compensator determined on the basis of selective and precisely known correlation characteristics of input interference signals. LMS-method based two-stage digital automatic compensator efficiency analysis with simultaneously operating stages and its comparison with a known one-stage automatic compensator with the same number of compensation channels were conducted with the help of static modeling at various conditions of interference correlation matrix and various numbers of interference sources. Conclusions. The obtained two-stage digital automatic compensator with orthogonalization of part of compensation channels is quasi-optimal regarding the criterion of minimum mean square error and it has lower computational complexity compared to two-stage digital automatic compensators based on LS and RLS-algorithms. Weight coefficients of the automatic compensators determined on the basis of selective and precisely known correlation characteristics of input interference signals are optimal regarding the criterion of least squares and minimum mean square error respectively. The obtained two-stage digital automatic compensators provide capability to connect/disconnect additional compensation channel block and appropriate computational stage modules proportionately to the interference situation.

Текст научной работы на тему «Определение весовых коэффициентов двухступенчатого автокомпенсатора на основе LMS-алгоритма и корреляционных характеристик помех»

УДК 621.396.963:621.391.26

ОПРЕДЕЛЕНИЕ ВЕСОВЫХ КОЭФФИЦИЕНТОВ

ДВУХСТУПЕНЧАТОГО АВТОКОМПЕНСАТОРА НА ОСНОВЕ LMS-АЛГОРИТМА и корреляционных характеристик

ПОМЕХ1

1 2 Жук С. Я.д.т.н., профессор , Семибаламут К. М.

1 Национальный технический университет Украины

«Киевский политехнический институт», г. Киев, Украин,

[email protected]

2Военно-дипломатическая академия имени Евгения Березняка,

г. Киев, Украина

DEFINITION OF WEIGHT COEFFICIENT IN TWO-STAGE AUTOMATIC COMPENSATION BASED ON LMS-ALGORITHM AND INTERFERENCE CORRELATION CHARACTERISTICS

Zhuk S. Ya., Doctor of Engineering, Professor; Semibalamut K. M.

National Technical University of Ukraine "Kyiv Polytechnic Institute ", Kyiv, Ukraine Yevhenii BerezniakMilitary-Diplomatic Academy, Kyiv, Ukraine

Вступление

Задача когерентной компенсации помех в цифровой антенной решетке в большинстве случаев сводится к определению весовых коэффициентов адаптивного фильтра [1]. Перспективным является подход, основанный на ортогональных преобразованиях сигналов [2], который позволяет получить алгоритмы, обеспечивающие параллельную обработку сигналов.

В статье [3] на основе теоремы разбиения регрессоров на два блока с использованием LS- и RLS-алгоритмов синтезированы двухступенчатые цифровые автокомпенсаторы активных шумовых помех, которые обеспечивают параллельную обработку сигналов и возможность подключения/отключения блока компенсационных каналов и соответственно части вычислительных модулей полученного устройства.

Широкое распространение для синтеза цифровых автокомпенсаторов активных шумовых помех нашел LMS-алгоритм, известный также как алгоритм Уидроу и Хоффа [4]. Он относится к классу алгоритмов стохастического градиентного поиска и позволяет получить решение задачи адаптивной фильтрации в условиях априорной неопределенности характеристик помех. Достоинством LMS-алгоритма является наименьшая вычислительная сложность.

При решении некоторых практических задач корреляционные характеристики помеховых сигналов являются известными или получены их оценки. Это позволяет выполнить расчет весовых коэффициентов до нача-

1 http://radap.kpi.ua/radiotechnique/article/view/1328

ла функционирования компенсатора. Кроме того, показатели эффективности алгоритмов компенсации помех, синтезированных при известных их корреляционных характеристиках позволяют определить потенциально достижимые границы эффективности для адаптивных алгоритмов.

В статье рассмотрены задачи определения весовых коэффициентов двухступенчатого автокомпенсатора c ортогонализацией сигналов части компенсационных каналов на основе LMS-алгоритма, а также корреляционных характеристик помех.

Постановка задачи

Пусть цифровая антенная решетка имеет основной канал, а также m

компенсационных каналов. Вектор yT =(y(1),...,y(n)) содержит отсчеты комплексной огибающей помехового сигнала y (i), i = 1, n на входе основного канала. Компонентами векторов xTy =( xy. (1),..., Xj (n)) являются отсчеты комплексных огибающих помеховых сигналов x. (i), i = 1, n, поступающие на входы компенсационных каналов j = 1, m . Векторы x , j = 1, m объединены в блоки в виде матриц X = К, ..., xk\ и X2 = |x¿+j, ..., xm|.

В работе [3] на основе теоремы разбиения регрессоров на два блока с использованием LS- и RLS-алгоритмов синтезированы двухступенчатые цифровые автокомпенсаторы активных шумовых помех c ортогонализацией сигналов части компенсационных каналов. Структурная схема двухступенчатого цифрового автокомпенсатора с модулями в виде весового сумматора показана на рис. 1.

Рис. 1. Структурная схема двухступенчатого цифрового автокомпенсатора с модулями в виде весового сумматора

В первой ступени выполняются компенсация помех в основном канале с использованием части сигналов компенсационных каналов х ,у = 1, к и

определяется вектор остатков е1, а также ортогонализация оставшихся векторов х ., у = к +1, т. После их получения, во второй ступени выполняется компенсация остатков помех е1 с использованием полученных орто-гонализированных регрессоров , у = к +1, т и вычисляется вектор остатков е . Как следует из результатов синтеза [3], при применении LS- и RLS-алгоритмов, переход к обработке во второй ступени выполняется после завершения обработки обучающего пакета в первой ступени. Вместе с тем, применение RLS-алгоритма обеспечивает возможность одновременного функционирования первой и второй ступеней автокомпенсатора, что имеет важное практическое значение при поступлении больших объемов данных.

Модули первой ступени представляют собой однотипные по структуре устройства в виде весового сумматора с одним основным и к компенсационными входами. Модуль второй ступени является весовым сумматором с числом компенсационных входов к — т. В каждом модуле первой и второй ступеней по критерию наименьших квадратов с использованием RLS-алгоритма решается задача определения соответствующих весовых коэф-

/V- - /V- /V ^

фициентов Рр у = к +1, т, (30, Р0.

Необходимо определить весовые коэффициенты двухступенчатого автокомпенсатора c ортогонализацией сигналов части компенсационных каналов на основе LMS-алгоритма, а также корреляционных характеристик помех. Провести сравнительный анализ эффективности двухступенчатого цифрового автокомпенсатора на основе LMS-алгоритма с одновременным функционированием ступеней и известного одноступенчатого цифрового автокомпенсатора.

Определение весовых коэффициентов двухступенчатого автокомпенсатора на основе LMS-алгоритма

Широко известным в силу своей простоты и универсальности методом адаптивной подстройки весовых коэффициентов в цифровых антенных решетках является LMS-алгоритм. Он обеспечивает квазиоптимальное решение задачи определения весовых коэффициентов весового сумматора по критерию минимума среднего квадрата ошибки (МСКО).

Известно, что векторы весовых коэффициентов сумматора, определяемые с помощью RLS-алгоритма сходятся к винеровскому решению, оптимальному по критерию МСКО [4]. Это обуславливает допустимость применения для вычисления весовых коэффициентов модулей автокомпенсатора активных шумовых помех c ортогонализацией сигналов части компенсационных каналов LMS-алгоритма. При этом структурная схема двух-

ступенчатого цифрового автокомпенсатора с использованием LMS-алгоритма совпадает со схемой, представленной на рис. 1.

Следуя методике синтеза двухступенчатого цифрового автокомпенсатора активных шумовых помех с использованием RLS-алгоритма [3], а также учитывая взаимосвязь между RLS- и LMS-алгоритмами [1, 4], можно показать, что с использованием LMS-алгоритма, вычисление весовых

л 1 -

коэффициентов для модулей первой ступени (3р j = к +1, m , может быть выполнено рекуррентно по формулам

в(0 = Pj0 -1) + 2цх*(0(X (О - Xi (Oft(i -1)), (1)

где Xj (i) = (x(О,---,X(i)) - вектор, содержащий отсчеты помеховых сигналов на входах компенсационных каналов j = 1, к на i -м шаге; ц - шаговый

множитель, определяющий скорость сходимости итерационного процесса и устойчивость процесса адаптации; "*" - знак комплексного сопряжения. Шаговый множитель выбирается из компромисса между скоростью сходимости и точностью настройки адаптивного фильтра. Антиградиент мгновенного значения квадратичной ошибки имеет вид

Vj (i -1) = 2x*(i)(х- (i) - x[ (z)ßj (i -1)). (2)

Вычисление весовых коэффициентов ß1 выполняется с помощью рекуррентного алгоритма аналогичного (1), в котором, вместо х (i) используется i -я компонента y(i) вектора y

р0(0 = ßUi -1) + 2цх*(0(y(i) - < (i)p-1)). (3)

При раздельном функционировании ступеней значения весовых коэффициентов, используемые в модулях первой ступени определяются на последнем n -м шаге (3j = ßj(n), j = к +1,m, (31 = ß0(n), после чего вычисля-

±

ются ортогонализированные регрессоры х , j = к +1, т и вектор остатков е1 по формулам

I -

Ху1 = Ху - Х1Р1, у = к +1, т; е1 = у - . (4)

Вычисление весовых коэффициентов с использованием LМS-алгоритма для модуля второй ступени (3 2 выполняется рекуррентно по формулам

Р2(0 = Р0(/ -1) + 2^(0 (Ю) - ^ (ОРОО -1)), (5)

где х^(г) = (х^+1(/),...,х^(г)) - вектор, включающий г -е компоненты орто-гонализированных регрессоров , у = к +1, т ; е (0 - г -я компонента

вектора остатков е выхода первой ступени.

Значения весовых коэффициентов модуля второй ступени определяются на последнем п -м шаге Р0 = Р0 (п), после чего вычисляется вектор остатков е по формуле

е = е—Х^Р 0, (6)

где х£ =

1 1 хк+1, •••, хт

- матрица, столбцами которой являются вектора

х1, 1 = к +1, т.

Для выполнения LМS-алгоритмов в модулях первой (1,3), и второй (6) ступеней должны быть заданы начальные условия, которые имеют вид

,.,01 = 0" 1 = К + I т п„ (0) = и" и р0'

р!(0) = 0;у = к +1,т, р 1(0) = 0; и р0(0) = 0. Алгоритмы (1,3,6), относятся к

классу итерационных алгоритмов стохастического градиентного поиска, известных как алгоритмы Уидроу и Хоффа [1, 4].

Двухступенчатый автокомпенсатор с использованием LMS-алгоритма также обеспечивает возможность одновременного функционирования первой и второй ступеней. При этом весовые коэффициенты Л1 - Л1 л2

Р°(0, У = к +1, т, Р 1(/), Р°(0 определяются на каждом г -м шаге на основе

рекуррентных алгоритмов (1,3,6), а сигналы на выходах модулей первой и второй ступеней определяются по формулам

х}(0 = Ху (0 — хТ (0Р1 (/ — 1), у = к +1, т; (7)

^(0 = у (г) — х\ (г)Р 0(/ — 1); (8)

е(г) = в1(г) — х1Т (г)Р 0(г — 1). (9)

Определение весовых коэффициентов двухступенчатого автокомпенсатора на основе корреляционных характеристик помеховых сигналов

Синтезированные двухступенчатые цифровые автокомпенсаторы активных шумовых помех с использованием LS-, RLS- и LMS-алгоритмов были получены без конкретизации характеристик помеховых сигналов, и являются адаптивными.

Представляет интерес зависимость между весовыми коэффициентами Л1 - Л1 Л2

Ру = к +1, т, Р0, Р0 и корреляционными характеристиками входных помеховых сигналов. Будем полагать, что комплексные огибающие помеховых сигналов на входах компенсационных х (г), У = 1, к и основного у(г)

каналов являются гауссовскими стационарными последовательностями с нулевыми математическими ожиданиями.

В работе [3] показано, что векторы весовых коэффициентов для модулей первой и второй ступеней с использованием LS-алгоритма определя-

ются по формулам

в 1 = (X? Х )_1 Х?х;, } = шт; (10)

р0 = (Х?Х )-1 Х?у; (11)

в2 =( Х2±н Х^ )_1 Х2±не1. (12)

Выражение (10) можно представить в виде

в1 = (Х?Х )" ХЧ = (2п)^ Щ\2п)аи = ]"Й и, } = к +1, т, (13)

где = 0,5х* (г )х[ (г) - выборочная корреляционная матрица (КМ) входных помеховых сигналов х1 (г) = ( х (г),..., хк (г)) компенсационных каналов I = 1, к с элементами 1, I, я = 1, к, которые определяются по формуле

1 п

К — X х*(г)х(г); (14)

2п1=1

<1^ = 0,5х* (г) х? (г), j = к +1, т - векторы выборочной взаимной корреляции

X (г) и сигналов х] (г) компенсационных каналов j = к +1, т с элементами С1/7,/ = 1, к, которые определяются по формуле

1 п

= х*(г) х;(г) • (15)

2п г=1

Необходимо отметить, что при вычислении выборочных корреляционных характеристик используется определение полной комплексной КМ помеховых сигналов на компенсационных входах в виде Кх = 0,5£[х*(г)х1 (г)], где х1 (г) = ( х^),..., Хт (г)) [1].

Аналогично вектор (31, вычисляемый по формуле (11), можно представить в виде

в0 = ]^ 10 , (16)

где ¿10 = 0,5х* (г)у(г) - вектор выборочной взаимной корреляции х (0,

г = 1, п и сигнала в основном канале у(г) с элементами Сш ,1 = 1, к, которые определяются по формуле

1 п

¿и 0х*(г) у(г) • (17)

2п

Выражения (13), (16) совпадают с известным уравнением для определения весовых коэффициентов оптимального по критерию наименьших квадратов автокомпенсатора по выборочным корреляционным характери-

стикам [1, 4]. Таким образом, ортогонализированные регрессоры х1, у = к +1, т формируются из ошибок компенсации помеховых сигналов

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

х., у = к +1, т, полученных с использованием сигналов компенсационных каналов у = 1, к.

Для анализа коэффициента Р0 выражение (12) преобразуем с учетом свойств проекторов [3] к виду

Р 0=(Хн Х)—1Х* Му=(х" М1Х0)—1 (МХ )н М1У=(Х^М1Х0 )—1XX у, (18) где М - ортогональный проектор, который имеет вид

М =1-Х! (Хн Х )-1 ХН.

При анализе формулы (18) отдельно рассмотрим выражения Х^МХ и ХНМу . Для первого выражения справедливы преобразования

I— Х1 ( хН Х1 ) Х1Н Х0 = хН Х0 - хН Х1 ( хН Х1 ) хН Х0 = (19)

= 0п ( Я2 —Я2°К12) = 0и1^1,

^ * ~Т г-

где Н2 = 0.5х2 (г)х2 (г) - выборочная КМ входных помеховых сигналов хТ (г) = ((г),..., (г)) компенсационных каналов у = к +1, т;

111 = 0,5х12 (г)х!Т (г) - выборочная КМ ортогонализированных помеховых сигналов ^(г) = (х1+1(г),•••,(г)); Ню = 0,5х2(/)х0(г), 1101 = 0,5х0(г>Т(г) -выборочные взаимно корреляционные матрицы входных помеховых сигналов х (г) и х2 (г), при этом 1121 = 11н.

Второе выражение также можно преобразовать к виду

ХН М1У = ХН У — ХН Х (ХН Х)—1 ХНу = 0п (а 00 —110111—Я10) = 0пй ^ (20)

/V

где с120 = 0,5х2 (г)у(г) - вектор выборочной взаимной корреляции входных сигналов х2 (г) компенсационных каналов I = 1, к и сигналом в основном

/V . .

канале у(г); с1 ^ = 0,5х22 (г)е(г) - вектор выборочной взаимной корреляции

выходных сигналов х1 (г) компенсационных каналов I = 1, к и ошибки ех(г)

на выходе первой ступени основного канала.

С учетом (19), (20) выражение (18) можно представить в виде

/V / /V /V /V /V \—1 ^ ^ 1 ^ \

Р0 = (и0 — И01И—11к?1) (а00 —^к—а 10)• (21)

Таким образом, если известны выборочные корреляционные характе-

ристики 1, Я2, Я215 сС10, сС20, векторные коэффициенты

Л1 - Л1 я2

(3] = к +1,т, Р0, (30 могут быть найдены по формулам (13), (16) и (21) и

являются оптимальными по критерию наименьших квадратов.

Следуя методике [4], при неограниченном увеличении тактов наблюдения п выборочные корреляционные характеристики

111, 112, 112112,сС10, сС20,сС^ стремятся к корреляционным характеристикам

1,1, К21, Я2, Со, Со, С 20. При этом, весовые коэффициенты модулей первой и второй ступеней определяются по формулам

р 0ско = rl1d1 ,, ] = к +1, т, р 0ско = Rг1dl0; (22)

Р0СКО = ^С2 = (R2 1 ) (^0 1 R21RГ1d10 ) • (23)

Выражения для весовых коэффициентов модулей первой ступени совпадают с винеровским решением, оптимальным по критерию МСКО [1, 4]. Соответственно дисперсии ст0.ско ортогонализированных сигналов на вы-

ходе первой ступени х. (г), у = к +1, т определяются по формулам

<ско = < - К-Ч у, у = к +1, т, (24)

где а°, у = к +1, т - дисперсии сигналов х (г) на входе компенсационных

каналов у = к +1, т.

Дисперсия а01СКО ошибки е (г) на выходе первой ступени основного канала описывается выражением

^0ЮКО =^0 - СГо1- Ч0, (25)

где а0 - дисперсия сигнала у(г) на входе основного канала.

Поскольку при неограниченном увеличении интервала наблюдения оценка по критерию наименьших квадратов сходится к оценке по критерию МСКО [4], дисперсию ст0ско ошибки в(1) на выходе второй ступени основного канала также можно представить в виде

°еСКО °е!СКО С20 12 С20, (26)

20 2

С учетом (23, 25) выражение (26) можно представить в виде

а°СКО = аУ -

у "10"1 "10

н /_ _ _ тт \-1

0

(27)

-(С00 -1011- Чо ) (10 -1011- ^) (С00 -1011- Чо )•

Для рассмотренной задачи, при известной полной КМ помеховых сигналов 1 и векторе взаимной корреляции сигналов в основном и дополни-

тельных каналах СХу одноступенчатый оптимальный по критерию МСКО

автокомпенсатор может быть получен на основе фильтра Винера [1, 4]. При этом дисперсия ошибки компенсации описывается выражением

^0еКО — сХу^хСу. (28)

Выполнив разбиение полной КМ помеховых сигналов 1 на блоки 1, Я12, Я21,1, на основе формулы Фробениуса для обращения блочной матрицы [5] после несложных преобразований выражение (28) можно привести к виду (27). Таким образом, полученный двухступенчатый автокомпенсатор с частичной ортогонализацией сигналов компенсационных каналов (22, 23) является оптимальным по критерию МСКО. При этом он обеспечивает параллельную обработку сигналов и возможность подключения/отключения блока компенсационных каналов и соответственно части вычислительных модулей полученного устройства.

Анализ эффективности двухступенчатого автокомпенсатора на основе

LМS-алгоритма с одновременным функционированием ступеней

Анализ эффективности двухступенчатого цифрового автокомпенсатора на основе LMS-алгоритма с одновременным функционированием ступеней выполнен с помощью статистического моделирования, в соответствии с соотношениями (1,3,6-9). В качестве показателей эффективности использованы: скорость настройки (сходимости) автокомпенсатора до установившегося режима (количество итераций п), мощность ошибок компенсации помех. Общее количество компенсационных каналов m=6. Компенсационные каналы разбивались на два блока с одинаковым числом каналов, равным трем (&=3). Суммарная мощность помех в 600 раз превышает мощность собственных шумов приемных каналов, которая полагалась равной единице, что соответствует 27,8 дБ. Моделирование проведено для поме-ховых ситуаций, создаваемых L постановщиками помех. Помеховые ситуации характеризуются хорошей (единицы децибел) и плохой (десять и более децибел) обусловленностью КМ Кх [3]. Рассмотрены помеховые ситуации при Ь < к и Ь > к. Сравнительный анализ проведен с цифровым автокомпенсатором по схеме m-канального одноступенчатого весового суммирования, в котором вычисление весовых коэффициентов выполняется на основе LMS-алгоритма [1]. Значение шагового множителя ц выбрано: для одноступенчатого цифрового автокомпенсатора - ц^^Лт^х), для модулей двухступенчатого цифрового автокомпенсатора - ц2=2ц1.

На рис. 2 представлены зависимости мощности ошибок компенсации помеховых сигналов на выходах 14, 10 первой и 20 второй ступеней цифрового автокомпенсатора (первое число соответствует номеру ступени, а второе - номеру канала) от количества итераций. Рассмотренная помехо-вая ситуация соответствует случаю плохой обусловленности КМ помех

и

30

а о о

и 25

и §

® 20

при действии двух постановщиков активных шумовых помех. Для сравнения на рис. 2 представлена зависимость мощности ошибок компенсации помеховых сигналов одноступенчатого цифрового автокомпенсатора помех.

Как известно, эффективная компенсация одновременно действующих L постановщиков активных шумовых помех осуществляется при наличии не менее L вспомогательных каналов, что в данном случае выполняется уже на первой ступени. Поэтому, как следует из результатов моделирования (выходы модулей 10 и 14, рис. 2), обеспечивается практически

■а н о о я 3

о

В 15 Я

а

и 10 я

о

о И

0 10

10 10 10 К о л и ч е с т в о и т е р а ц и й,

10

Рис. 2. Мощность ошибок компенсации помех: т=6; £=3; L =2; л=21 дБ

одинаковая эффективность подавления помехи в основном канале первой и второй ступени. Такая же эффективность сохраняется и при хорошей обусловленности КМ помех.

На рис. 3-5 представлены результаты модели-ро-вания для случаев хорошей и плохой обусловленности КМ помех в ситуациях Ь = к и Ь > к. В установив-шемся режиме двухступен-чатый и одноступенчатый автокомпенсаторы имеют одинаковую эффектив-ность. Скорость настройки (рис. 3) двухступенчатого автокомпенсатора на переходном участке работы опережает скорость нас-

0

10° 101 102

К о л и ч е с т в о и т е р а ц и й, п

Рис. 3. Мощность ошибок компенсации помех: т=6; £=3; L =3; ^=2,5 дБ

10

тройки одноступенчатого автокомпенсатора примерно на 10 итераций (по уровню 5 дБ).

10

10

В установившемся режиме, мощность ошибки компенсации помех на выходе 10 на 1,5 дБ больше, чем на выходе второй ступени 20 (рис. 3, 4), что подтверждает необходи-мость применения двух-

ступенчатой схемы автокомпенсатора при Ь = к как при хорошей, так и при плохой обусловленности КМ помеховых сигналов. При плохой обусловленности КМ помех и количестве постановщиков помех

10 10 10 К о л и ч е с т в о и т е р а ц и й, п

Рис. 4. Мощность ошибок компенсации помех:

т=6; к=3; Ь =3; ^=19 дБ

Ь = 3 (рис. 4), скорость сходимости двухступенчатого автокомпенсатора в

3 раза выше скорости сходимости одноступенчатого автокомпенсатора (по

уровню 5 дБ).

30 ......... ......... ......... ......... Двухступенчатой автокомпенсатор также сохраняет работоспособность, когда количество постановщиков помех больше числа дополнительных каналов в первой ступени рис. 5. Так, при плохой обусловленности КМ помех, скорость сходимости двухступенчатого автокомпенсатора в 5 раз выше скорости сходимости одноступенчатого автокомпенсатора. Мощность ошибок компенсации помех на выходе второй ступени 20 по сравнению с выходом первой ступени 10 меньше на 9 дБ.

Зависимости мощности ошибок компенсации помех выхода 14 первой ступени на рис. 2-5 иллюстрируют работу модулей, в которых выполняется ортогонализа-ция помеховых сигналов дополнительных каналов.

На рис. 5 также показана рассчитанная для рассмот-ренного примера по

о

10о 101 102 103

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

К о л и ч е с т в о и т е р а ц и й, п

Рис. 5. Мощность ошибок компенсации помех:

т=6; к=3; Ь =4; л=21 дБ;

104

формуле (27) мощность а°ско =4,1 дБ ошибки на выходе второй ступени основного канала (20opt). По сравнению с оптимальным решением, мощность ошибки компенсации помех автокомпенсаторами на основе LMS-алгоритма в установившемся режиме больше на 1,5 дБ, что обусловлено шумом градиента [1] и подтверждает их квазиоптимальный характер.

Выводы

Полученный с использованием LMS-алгоритма двухступенчатый цифровой автокомпенсатор с ортогонализацией части компенсационных каналов является квазиоптимальным по критерию МСКО и характеризуется значительно меньшей вычислительной сложностью по сравнению с двухступенчатыми автокомпенсаторами на основе LS- и RLS-алгоритмов [3]. Он может функционировать в режимах, когда обработка всего обучающего пакета входных данных выполняется сначала в первой, а затем второй ступенях, а также при одновременном функционировании ступеней.

Векторные коэффициенты ßj = k +1, m, ß0, ß0 модулей первой и второй ступеней, определяемые по формулам (13), (16) и (21) на основе выборочных корреляционных характеристик Rj, R2, R2рd10,d20 являются оптимальными по критерию наименьших квадратов и равны соответствующим коэффициентам, определяемым на основе LS- алгоритма [3]. Векторные

коэффициенты ß°ско, j = k +1, m, ßjCKO, ß°ско модулей первой и второй ступеней, определяемые по формулам (22, 23) на основе корреляционных характеристик R , R , R ,d ,d являются оптимальными по критерию МСКО.

Для рассмотренных примеров скорость сходимости двухступенчатого автокомпенсатора не ниже чем у одноступенчатого с тем же числом компенсационных каналов, а при числе постановщиков помех L > 3 выигрыш составляет от 3 до 5 раз. При L < 3 мощности ошибок компенсации помех в основном канале на выходах первой и второй ступеней практически одинаковы, что позволяет отключать вторую ступень и переходить к одноступенчатой структуре автокомпенсатора с меньшим числом компенсационных каналов.

Список источников

1. Monzingo R. A. Introduction to adaptive arrays. 2nd ed. / R. A. Monzingo, R. L. Haupt, T. W. Miller. - Scitech publishing, 2011. - 510 p. doi: 10.1049/sbew046e

2. Ратынский М. В. Адаптация и сверхразрешение в антенных решетках / М. В. Ра-тынский. - М. : Радио и связь, 2003. - 200 с.

3. Жук С.Я. Двухступенчатая адаптивная компенсация активных шумовых помех с ортогонализацией сигналов части компенсационных каналов / С. Я. Жук, К. М. Семи-баламут // Вестник НТУУ "КПИ". Серия Радиотехника, Радиоаппаратостроение. -2016. - № 64. - С. 61-74.

4. Джиган В. И. Адаптивная фильтрация сигналов: теория и алгоритмы / В. И. Джи-ган. - М. : Техносфера, 2013. - 528 с.

5. Гантмахер Ф. Р. Теория матриц / Ф. Р. Гантмахер. - М. : Физматлит, 2004. -560 с.

References

1. Monzingo R. A., Haupt R. L. and Miller T. W. (2011) Introduction to adaptive arrays. 2nd ed., Scitech publishing, 510 p. doi: 10.1049/sbew046e

2. Ratynskii M. V. (2003) Adaptatsiya i sverkhrazreshenie v antennykh reshetkakh [Adaptation and superresolution in antenna arrays]. Moskow, Radio i svyaz', 200 p.

3. Zhuk, S. Ya., Semibalamut, K. M. (2016) Two-stage adaptive compensation of active noise interference with signals orthogonalization of a part of compensation channels. Visn. NTUUKPI, Ser. Radioteh. radioaparatobuduv., no. 64, pp. 61-74. (in Russian).

4. Dzhigan V. I. (2013) Adaptivnaya fil'tratsiya signalov: teoriya i algoritmy [Adaptive Signal Filtering: Theory and Algorithms]. Moskow, Tekhnosfera, 528 p.

5. Gantmakher F. R. (2004) Teoriya matrits [The theory of matrices]. Moskow, Fizmatlit, 560 p.

Жук С. Я., Семгбаламут К. М. Визначення вагових коеф1ц1ент1в двохступеневого автокомпенсатору на основ1 LMS-алгоритму i кореляцшних характеристик завад.

Розглянуто задач1 визначення вагових коефщент1в двохступеневого автокомпенсатору з ортогонал1зац1ею сигнал1в частини компенсацтних канал1в на основ1 LMS-алгоритму, а також кореляцшних характеристик завад. Анал1з ефективност1 реку-рентного двохступеневого цифрового автокомпенсатору на основ1 LMS-алгоритму з одночасним функцюнуванням ступетв виконано за допомогою статистичного моде-лювання при р1знт обумовленост1 кореляцтног матриц завад i р1знт к1лькост1 джерел завад.

Ключовi слова: автокомпенсатори завад, метод найменших квадратiв, LMS-алгоритм.

Жук С. Я., Семибаламут К. М. Определение весовых коэффициентов двухступенчатого автокомпенсатора на основе LMS-алгоритма и корреляционных характеристик помех. Рассмотрены задачи определения весовых коэффициентов двухступенчатого автокомпенсатора c ортогонализацией сигналов части компенсационных каналов на основе LMS-алгоритма, а также корреляционных характеристик помех. Анализ эффективности рекуррентного двухступенчатого цифрового автокомпенсатора на основе LMS-алгоритма с одновременным функционированием ступеней выполнен с помощью статистического моделирования при разной обусловленности корреляционной матрицы помех и разном количестве источников помех.

Ключевые слова: автокомпенсаторы помех, метод наименьших квадратов, LMS-алгоритм.

Zhuk S. Ya., Semibalamut K. M. Definition of weight coefficient in two-stage automatic compensation based on LMS-algorithm and interference correlation characteristics.

Introduction. Two-stage digital automatic compensators of active noise interference based on LS and RLS-algorithms provide parallel processing of signals and capability to connect/disconnect additional compensation channel block. LMS-algorithm is widely used for synthesis of digital automatic compensators of active noise interference. Its advantage is low computational cost. When solving certain practical problems correlation characteristics of interference signals are either given or their estimates are obtained.

Problem statement. The first stage of the two-stage digital automatic compensator provides compensation of interference in the main channel using compensation channel signals of one of the blocks and orthogonalization of signals in the other block. The second stage provides compensation of the residue interference using the obtained orthogonalized signals. Modules of the first and the second stage are single type weight totalizers. We need to determine weight coeeficient of the two-stage digital automatic compensator with orthogonalization of signals of part of compensation channels based on LMS-algorithms and interference correlation characteristics.

Using LMS-algorithm there was obtained two-stage digital automatic compensator with orthogonalization of part of compensation channels that is quasi-optimal regarding the criterion of minimum mean square error and it's feature is lower computational complexity compared to two-stage digital automatic compensators based on LS and RLS-algorithms. It can operate in modes when preocessing of the complete training package of the input data is carried out in the first and then in the second stage as well as in the mode of simultaneous operation of the stages.

There have been acheived analytical expressions for defining weight coefficients of modules in the first and in the second stage of the automatic compensator determined on the basis of selective and precisely known correlation characteristics of input interference signals.

LMS-method based two-stage digital automatic compensator efficiency analysis with simultaneously operating stages and its comparison with a known one-stage automatic compensator with the same number of compensation channels were conducted with the help of static modeling at various conditions of interference correlation matrix and various numbers of interference sources.

Conclusions. The obtained two-stage digital automatic compensator with orthogonaliza-tion of part of compensation channels is quasi-optimal regarding the criterion of minimum mean square error and it has lower computational complexity compared to two-stage digital automatic compensators based on LS and RLS-algorithms. Weight coefficients of the automatic compensators determined on the basis of selective and precisely known correlation characteristics of input interference signals are optimal regarding the criterion of least squares and minimum mean square error respectively. The obtained two-stage digital automatic compensators provide capability to connect/disconnect additional compensation channel block and appropriate computational stage modules proportionately to the interference situation.

Keywords: automatic interference compensators, the method of least squares, LMS -algorithm.

i Надоели баннеры? Вы всегда можете отключить рекламу.