CC BY
13Транспорт
УДК 004.042
DOI: 10.30987/article 5db95e85b8bb12.44709607
В.А. Ермолаев, Ю.А. Кропотов
ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ МОДЕЛЕЙ АКУСТИЧЕСКИХ СИГНАЛОВ СИСТЕМ ТЕЛЕКОММУНИКАЦИЙ АУДИООБМЕНА В УСЛОВИЯХ ПОМЕХ И АКУСТИЧЕСКОЙ ОБРАТНОЙ СВЯЗИ
Описывается решение задачи проектирования моделей систем обмена информацией с акустической обратной связью с применением методов теории линейных функционально-
дифференциальных уравнений. Исследование устойчивости моделей выполняется частотным методом. Рассматривается подход к оцениванию корреляционных и спектральных функций акустических сигналов, основанный на параметрическом
представлении акустических сигналов и компонент шумовых составляющих.
Ключевые слова: модели акустических сигналов, системы телекоммуникаций аудиообме-на, акустическая обратная связь, эхосигналы, речевая связь, озвучивание, функционально-дифференциальные уравнения, аппроксимация функций распределения, корреляционная функция, адаптивный фильтр.
V.A. Yermolaev, Yu.A. Kropotov
MODEL PARAMETER IDENTIFICATION OF ACOUSTIC SIGNALS IN AUDIO EXCHANGE TELECOMMUNICATION SYSTEMS UNDER CONDITIONS OF INTERFERENCE AND ACOUSTIC FEEDBACK SYSTEM
The paper reports the solution description for the problem of designing the models of the systems of information exchange with the acoustic feedback with the use of the methods of the linear functional-differential equation theory. With the aid of them it is ensured a stability increase in systems of operational-command loud-speaker communication and notice, and also the problem of echo decrease at the inputs of speech transforming devices is solved that conditions on quality increase of their functioning. In the paper there are presented models having a continuous consequence and taking into account a character of acoustic echo-signals in closed spaces more complete which results in authenticity increase of modeling results at echo-signal compensation. To solve the echo identification problem the functions of echo lag distribution
Введение
При построении моделей систем с дискретным и распределенным запаздыванием в самых различных областях - механики и электротехники, автоматического управления и регулирования, синхронизации, нейронных сетей [1], а также акустики [2] - находят применение функционально-дифференциальные уравнения. Становление этого направления выразилось в потоке публикаций, в частности по динамике и стабилизации процессов в сетях, состоя-
are approximated by a number of exponents which simplifies the equations and allows accepting a consequence both at a finite and infinite interval. The investigation of the stability of models of systems with the acoustic feedback in the paper is carried out by a frequency method. At that there is considered an approach to the estimate of correlate and spectral functions of acoustic signals which is based on a parametric presentation of acoustic signals and components of noise constituents.
Key words: models of acoustic signals, audio exchange telecommunication systems, acoustic feedback, echo-signals, voice communication, postsynchronizing, functional-differential equations, approximation of distribution functions, correlation function, adaptive filter.
щих из множества релаксационных нейро-подобных автогенераторов и мемристоров.
Во многих случаях на проблему решения функционально-дифференциальных уравнений накладываются сложности, обусловленные неопределенностью параметров моделируемой системы. Наблюдающийся в последние годы рост интереса к системам с неопределенностью нашел выражение в понятиях робастности, робаст-ной устойчивости, робастного управления
и развитии соответствующих методов анализа. При этом неопределенность, в зависимости от характера задачи, может заключаться в отсутствии точных сведений о параметрах линейных или нелинейных элементов модели; параметров, обусловленных их естественным разбросом и изменением во времени, а в случае систем с запаздыванием (или последействием) - в отсутствии точных сведений о характере и величине запаздывания.
В системах с акустической обратной связью неопределенность обусловливается отсутствием точной информации о характере звукового поля, резонансах и направлениях прихода звуковых волн в точке наблюдения, параметрах отражающих поверхностей, об ориентации и диаграмме направленности приемника звука и, соответственно, о распределении компонентов наблюдаемого сигнала по величине запаздывания и фазе. При всей важности проблемы неопределенности о моделях систем с акустической обратной связью, описываемых функционально-дифференциальными уравнениями, известно немногое [3].
Значимость решения проблемы неопределенности заключается в стремлении повысить качество речевой связи и озвучивания, в частности повысить разборчивость речи, которая, как известно [4], зависит от характеристик акустического эха и акустической обратной связи. За формирование акустического эха отвечает ряд прямых и отраженных звуковых волн, а также
диффузное излучение, обусловленное реверберацией, т.е. рассеянием энергии, заключенной в резонансных колебаниях акустической среды. Особо негативно акустическое эхо влияет на качество сжатия речи, а при наличии акустической обратной связи - на устойчивость акустической системы.
Неопределенность акустической системы может иметь как стационарный, заключающийся в естественном разбросе ее характеристик, так и динамический характер: всякое непроизвольное (порядка десятка сантиметров) смещение или разворот микрофона, как и кратковременные изменения условий распространения звука, обусловливают изменение не только уровня сигнала, но и, что более существенно, его фазы.
Названные причины неопределенности требуют усложнения представляющих их моделей, к числу которых можно отнести модели с переменными параметрами, в частности с переменным запаздыванием, модели, описываемые детерминированными и стохастическими функционально-дифференциальными уравнениями, модели с элементами адаптации и робастные модели. Помимо сложности анализа перечисленных моделей, требующих применения численных методов, встает и проблема идентификации параметров моделей, относящаяся в общем случае к разряду нелинейных.
Механизм формирования акустической обратной связи поясняется рис. 1.
Рис. 1. Механизм формирования акустической обратной связи
Модель акустического эха в простейшем случае представляют линейной комбинацией u(t) запаздывающих копий излучаемого сигнала v(t) :
u(t) = kv(t -хк к=1
(1)
где v(t) = J f (0)uq (t - 0)d0 * F(s)U0 (s).
0
Здесь хк и а - соответственно время распространения (запаздывание) и амплитуда звуковой волны, поступающей на вход приемника по к -пути, f (0) - импульсная функция усилителя входного сигнала Uo(t), значок * является символом эквивалентности оригинала функции действительного аргумента ее изображению, а г - число значимых путей распространения звука.
Модель (1), фактически заданная дифференциально-разностным уравнением [5], справедлива, строго говоря, лишь при дискретном распределении эха по величине задержки и отсутствии реверберации
Модели формирования эха и реверберации
Поставленную цель можно достичь в классе функционально-дифференциальных систем или систем с распределенным последействием, представляя распределение задержек звука непрерывными или обобщенными функциями, а резонансные свойства акустической среды - резонансами систем с сосредоточенными параметрами. Схема соответствующей модели приведена
- ситуации, возникающей, например, при озвучивании открытых площадей с небольшим числом удаленных отражателей.
Построение моделей эха и систем с акустической обратной связью позволяет получить представление о непрерывном распределении эха по величине задержки и о резонансных процессах реверберации. Непрерывное распределение эха по величине задержки моделируется при этом ядрами детерминированных функционально-дифференциальных уравнений [3], а резонансные колебания акустической среды -передаточными функциями соответствующих резонаторов. В рассматриваемом линейном приближении проблема устойчивости систем с акустической обратной связью решается частотными методами, из которых обсуждаются только критерии Найквиста и Понтрягина. Проблема идентификации ядер считается решенной, например методом наименьших квадратов. Приводятся уравнения, использованные при рассмотрении примеров простых систем.
на рис. 2, где функции gk (0) аппроксимируют распределения задержек, а передаточные функции Нк (з) - резонансные свойства акустической среды. При этом (0) - 0 только в области значений 0 е [0, ]. В противном случае они равны нулю.
Рис. 2. Модель формирования эха, обусловленного последействием и реверберацией
Эта модель описывается функционально-дифференциальным уравнением
r 0
u(t) =1 J gk (0)Vk (t-Xk +0)d0. k=1
где Vk (0) * Hk (s)V(s).
r
t
Здесь ук (г) - сигнал на выходе резо-наторного блока Ик (з), [0, 3к ] - интервал значимой величины непрерывной задержки звука, распространяющегося по к -пути. Как и в уравнении (1), функция г
у(г) = | /Ф)ио (г ~ 0)С0 представляет сиг-
о
нал на выходе блока ^(Х), а набор параметров тк, к = 1,..., г, - конечное множество начальных значений задержки.
Очевидно, что если все передаточные функции И к (з) = 1, а функции распределения задержек представлены дельта-функциями вида gk (0) = ак^(0), то уравнения (1) и (2) оказываются эквивалентными. Если же эхо аппроксимируется суперпозицией волн, как непрерывно, так и дискретно распределенных по величине задержки, то уравнение (2) записывается (после очевидных замен) в виде
q г 3к
и(г) = ХаV(г-т/) + X |§к(0)Ук(г-Тк-0)С0.
/=1 к=q+1 0
В области изображений уравнению (3) отвечает выражение
^ 5 ^(з) = £а/И/(з)е
^(з) ^(з^(з) Й 1 1 ( )
Г/ + X е"зТк | gk (0)Ик (з)е-30С0.
к=q +1 0
(3)
(4)
При этом очевидно, что функция Ж(з) является передаточной функцией модели, изображенной на рис. 2. Функции И к (з), считающиеся в рамках рассматриваемой модели независимыми от параметра задержки 0 , выносятся в (4) из-под знаков интеграла. Конечно, учет зависимости функций И к (з) от параметра 0, когда И к (з) = И к (з, 0), позволяет обобщить модель и потенциально повысить ее точность. Однако ниже эта возможность не рассматривается.
Следуя [6], отметим, что первая группа слагаемых (числом q ) в правой части уравнений (3) и (4) моделирует основной звук и ранние отражения, а вторая группа - реверберацию. Ранние отражения - это небольшое число волн, запаздывающих относительно основного звука на величину около 60 мс, а реверберация - волны, запаздывающие на 60.. .300 мс и более. Конечно, эти цифры, зависящие от геометрии помещения и расположения отражающих поверхностей, носят ориентировочный характер.
Аппроксимация функций распределения запаздывания по величине
Чтобы придать передаточной функции Ж (з) явный вид, ее ядра (распределения запаздываний по величине) gk (0) необходимо каким-то способом аппрокси-
мировать, например последовательностями экспонент или ступенчатыми функциями. В первом, простейшем случае их можно представить выражениями вида
gk (0) =
ск (1 - е - Рк 0 )е 0 = = ск (е —к0 - е + рк)0), 0е[о, 3к ] 0,
(5)
При этом максимум функции (5), равный
gk,тах = 0, 5Ск (1 + рк / qk 'Рк рк / (дк + рк ) ,
достигается при значении
0*[0, 3к ]
0тах = (1/Рк )1п(1 + Рк / qk ) > 0 • На рис- 3 изображен график одного такого ядра при значении 0тах = 1/34 с, что соответствует десятиметровому пути распространения
3
к
г
<
звук^ тах = 0,5
и
ck = 1, откуда
pк /чк = 3,4, pk = 50,37 и дк = 14,8.
Подстановка (5) в соответствующее слагаемое уравнения (4) и последующее интегрирование позволяет получить выражение
{ gk (0)e "э6d0 = Ck
e k (эщ k) + Pk (э + qk + Pk)(э + qk)
где щ =е
= е-1Л [е-РА _])и ^ =е-чА-{дк +Рк).
(6)
Рис. 3. Пример функции распределения запаздывания (5)
9
k
0
Поскольку распределения запаздываний (5) и выражения (6) справедливы при всех значениях k (отличаются только па-
ч
раметрами Ck, Pk, qk и ), передаточная функция
Ж(э) = ¥(э)ЪщН№ 1 +т I скНк(э)
I=1
k=ч+1
"Э(Т k )(эЛk + Ук ) + е~эх kPk (э+qk + Pk)(э+qk)
В том случае, когда ^ ^^, выражение (7) принимает вид
Щэ) = ¥(э)!а1Н1 (э)е-П1 +¥(э) I скНк(э) I=1
'■Рк
(7)
(8)
Передаточные функции ¥ (э) и Нk (э), k = 1,..., г, представляют соответственно свойства аппаратных средств передачи звука и свойства акустической среды. Модуль функции ¥(]а) в рабочем диапазоне частот обычно считается равномерным (специально вводимая иногда в целях повышения разборчивости речи неравномерность не меняет существа зада-
k=ч +1 (э+qk + Pk)(э+qk)
чи); фазовая характеристика названной функции значимой при этом не считается. В общем случае можно записать: ¥(э) = Л(э)1 В(э), где А(э) и В(э) - многочлены соответственно степени т и п > т.
В ряде случаев ядра функционально-дифференциальных уравнений проще аппроксимировать ступенчатыми функциями вида
откуда следует, что
м
gk (0) =1 gk (тД - А / 2)(1(0 - (т -1)А) -1(0 - тД)):
т=1
г
е
е
3 1 м
| gk (0)е-звёд = -(1 - е-зА )Х gkme
-з(т-1)А _
Здесь gkт = gk(тА-А /2), А - период дискретизации аппроксимируемого ядра и 1(0) - единичная функция, т.е.
Г1, 0 > 0,
1(0) ч
I 0, 0 < 0.
,1 + Х (^к,т+1 gk,т )е ,ме
-зМ А
Соответственно передаточная функция разомкнутой системы
~зТ/
1
Ж (з) = Р (з) Ха/И/ (з)е /=1
Особенность рассматриваемой задачи заключена в том, что аппроксимации подлежат ядра подынтегральных выражений - по своему смыслу заведомо положительные и убывающие на бесконечности функции. Аппроксимация ядер по системе экспоненциально убывающих функций, примером чему является выражение (5), позволяет решить эту проблему на всей действительной полуоси. Аналогично ре-
+ - Р (з)(1 з
г М
- е"зА) XXИк(з) Zgk,те'з(т-1)А .
к=q+1 т =1
шается, только на конечном интервале, задача аппроксимации ступенчатыми функциями. Конечно, задача может быть решена и по системе степенных функций. В этом случае требуется принять дополнительные меры по выделению области допустимых значений результатов аппроксимации, включая расходимость на бесконечности. Примером такой аппроксимации является функция
gk (0) =
Л2(1 - е-*9 )-Л10, 0е[0, 3к ], 0, 0й[0, 3к ]
Л2 >Л1 > 0, Х = --^1п(1 -Л13к / Л2 )> 0, 3к > 0. 3к
Подстановка этой функции под знак интеграла функционально-
дифференциального уравнения и интегрирование дает выражение
| gk (0)е-з0С0 = [Л1 (з"2 + 0з-1 )е-з0 + Л2 (е-(з+Х)0 (з + X)-1 - е-з0з-1)]
0
к
1
0
2 Ч[(Л2Х-Л1)з-ЛхХ]+[(Л2(ек -1) + Л:3к)з2 + (тЩк +1)-^з +Лх^]е"з3к} з 2(з + Х)
к
з + Ь
2
з + Ь2 з + Ьз -з3к
= а
1 з 2(з + Х)
что усложняет решение и, в силу неестественной формы распределения запазды-
1 + а," : ' "3 е "з3к.
2 з2 (з + Х)
ваний по величине, обусловливает появление сомнений в его адекватности.
Модели акустических систем с запаздывающей обратной связью
Согласно рис. 4, прямая часть системы с акустической обратной связью (т.е. система с разомкнутой обратной связью) эквивалентна модели эхосигнала у(г) в области его смешивания с входным сигналом х(г) и определяется передаточной функцией Ж(з). При этом для передаточ-
ной функции замкнутой системы имеет место известное выражение
ад =
Ж (з)
(9)
1 + Ж (з)
Соответственно изображение сигнала у(г) на выходе источника звука с передаточной функцией Р (з) имеет вид
0
q
3
V (э) = ¥(э) X (э).
1 + Ж (э)
Функция ¥ (э), которая, как правило, может считаться характеристикой полосового фильтра или близкой к нему системы, а также, по определению, и Ж (э) являются
устойчивыми. По этой причине факт устойчивости или неустойчивости модели с акустической обратной связью полностью определяется свойствами характеристической функции 1 + Ж (э) - свойствами знаменателя формулы (9).
Рис. 4. Модель системы с акустической обратной связью
Функции Н1 (э) = Р1 (э)/<2г (э),
I = 1,..., ч, и Нk (э) = Рк (э)/& (э), k = ч +1,..., г, определяют соответственно передачу звуковых волн по прямым акустическим каналам и волн реверберации -переносчиков энергии, накопленной в акустических резонаторах озвучиваемого объекта. (Приближенно можно положить, что
Н1 (э) = 1, а Нк(э) = 1/(э2 + Ъкэ + йк).) Если к тому же обозначить
(э + чк + Рк)(э + чк) = Ок(э), то характеристический многочлен (квазиполином) замкнутой системы с передаточной функцией (9), при условии что разомкнутая система задана выражением (8), можно записать в виде
ч
В(э)П0п(э) П От(э) + Л(э)П0п(э) П От(э)Р(э)е~
п=1 т=ч+1 /=1 п=1 т=ч+1
пф1
г
г
г
г
ч г
+Л(э) I СкРк П 0-п О
к=ч +1 п=1
пф к
Для устойчивости замкнутой системы (при разомкнутой обратной связи система считается устойчивой) все корни характеристического многочлена должны находиться в левой полуплоскости комплексной переменной. Поскольку нахождение корней квазиполинома (10) - сложная задача, то обычно прибегают к частотным методам проверки устойчивости. Используются также и имеющие большую общность критерии устойчивости во временной области, что требует перехода к описанию модели системой дифференциальных уравнений (обычно первого порядка).
I п От (э)Рк (э)е~эТк =0. (10)
т=ч+1 тф к
Из известных частотных критериев -Михайлова, Найквиста, Понтрягина и D-разбиения - ниже обсуждаются только критерии Найквиста и Понтрягина. В рамках критерия Найквиста утверждается, что годограф комплексной функции Ж(уш) устойчивой системы не должен охватывать точки (-1, у'0) комплексной плоскости. Другой подход к анализу устойчивости, критерий Понтрягина, основывается на анализе косвенных условий отсутствия нулей квазиполинома системы в левой полуплоскости комплексной переменной.
Пусть комплексная функция Ж(уш), соответствующая (7) или (8), имеет вид
q
Ж О) =Ха/Ж/ (ш)е
/=1
-Дф/(ш)+шт /)
+
ХХскЖк (ш)е~ /(ф к(ш)+шт к).
к=q+1
(11)
Здесь амплитудные и фазовые характеристики компонентов разомкнутой системы Ж(з) , определяемой выражением
(8), описываются входящими в (11) функциями
Ж/ (ш) = |р(>)И/ С/ш)| и Жк н =|р(/ш)Ик мРк/(-ш2+М^к + Рк)+(qk + Рк ^) I
Ф/(ш) = аг§(р(/ш)И/ (/ш)) и фк (ш) = р(>)Ик О^к /(-о2 + + Рк) + ^к + Рк)qk).
Топологический анализ годографа функции (11) или соответствующего характеристической функции 1 + Ж (/ш) квазиполинома Р( /ш) в принципе решает задачу устойчивости системы с акустической обратной связью методами Найквиста и Понтрягина соответственно. При этом наличие неизбежных погрешностей идентификации параметров модели, в частности погрешностей аппроксимации ядер, удается учесть посредством их соответствующей модификации. Это нашло выражение в появлении области робастного анализа и робастной устойчивости.
В одной из модификаций критерия Найквиста предполагается, что разомкнутая система, которая при охвате ее обратной связью сохраняет свою устойчивость, характеризуется априорно известной (номинальной) передаточной функцией
Ж 0(/ш), а ее неопределенность задается либо аддитивным, либо мультипликативным неравенством. Пусть передаточной функции разомкнутой системы с неопределенными параметрами отвечает выражение (11). Передаточную функцию Ж 0(/ш) номинальной системы можно задать аналогичным выражением, обозначив ее но-
минальные компоненты как Жп (/ш),
Ж/°(ш) и ф°(ш) . При этом аддитивная мера неопределенности определяется неравенством
Ж(/ш) - Ж0 (/ш) <Ху (ш),
а мультипликативная мера, которую можно свести к аддитивной мере, - неравенством
Ж(/ш) - Ж0 (/ш) < Ху (ш) Ж0 (/ш)
где у/(ш) - некоторая заранее заданная положительная функция.
В случае аддитивной неопределенности и устойчивости номинальной разомкнутой системы робастная устойчивость замкнутой системы имеет место, если го-
1 + Ж 0(,/ш)
дограф функции Жу (/ш) =
-1
не охватывает круга радиусом X с центром в точке (-1, /0) . Критерий Понтрягина основывается на полученной им теореме, согласно которой если корни квазиполинома, например вида (10), находятся в левой полуплоскости комплексной переменной, т.е. система устойчива, то: 1) коэффициент при главном слагаемом - слагаемом с максимальными показателями степени квазиполинома - отличен от нуля; 2) выполняется неравенство
Я'(&У(ш) - Я(ш) Р '(ш) > 0, ¥ше Я,
(12)
где Л(ш) и J(ш) - соответственно действительная и мнимая части квазиполинома - характеристического многочлена, записанного как Р(/ш) = Я(ш>) + р (ш), и корни Я(ш) и Р(ш) являются простыми, действительными и перемежающимися.
Обратно, все корни квазиполинома находятся в левой полуплоскости комплексной переменной: 1) когда все корни Л(ш) и Р(ш) являются простыми, действительными, выполняется неравенство (12) и коэффициент главного слагаемого отличен от нуля; 2) когда все корни ^(ш)
и 3(ш) являются действительными, простыми и обеспечивающими выполнение неравенства (12), а коэффициент главного слагаемого отличен от нуля.
В моделях множественных отражений ядра интегральных выражений - непрерывные (возможно, за исключением
конечного числа точек разрыва); распределения интенсивностей эха по величине запаздывания подлежат аппроксимации неотрицательными функциями, графики которых аналогичны по форме графику, приведенному на рис. 3.
Уравнения двухканальной модели прямого и распределенного эха
При единственном канале прямого распространения эха и одном канале рас-
пределенного запаздывания уравнения (7) и (8) принимают соответственно вид
Ж (э) = ¥ (э)аН1 (э)е~ эТ1 + ¥ (э)сН2 (э)
Гэ(т 2 +&)(эЛ + ^) + е~эт 2 ^ (э + ч + р)(э + ч)
Ж (э) = ¥ (э)аН1 (э)е - эТ1 + ¥ (э)сН2 (э)
(э + ч + р)(э + ч)
(13)
(14)
Здесь л = е "ч9(е - р9
-1), ц = е"ч9(че"р9- (ч + р)) и | g(0)е
0
-э0 Л = се + Р
(э + ч + р)(э + ч)
Любое из полученных выше уравнений можно представить известными способами и в форме системы функционально-дифференциальных уравнений первого порядка - форме, обычно принимаемой в работах по данной теме в качестве основной. Однако, при всех достоинствах, ниже эта возможность не рассматривается. Из уравнений (13) и (14) видно, что они имеют самостоятельное значение, заключаю-
щееся, помимо упрощения численного анализа, в установлении более прозрачной зависимости характеристик модели от ее параметров, в том числе в установлении роли каналов распространения прямого и распределенного эха.
Из уравнения (14), в частности, следует, что соответствующий ему квазиполином имеет вид
2
е
Р(э) = Л(э)(э + ч + р)(э + ч)02 (э)аР1 (э)е"эТ1 + А(э)с01 (э)рР2 (э)е"
12 -
= 0.
Оценка параметров модели методом максимального правдоподобия
При оценивании параметров модели систем с акустической обратной связью, например, может быть использован метод, основанный на согласовании корреляционных функций. В этом случае используются методы, основанные на собственных значениях и сингулярных разложениях ковариационных матриц, имеющих тёплице-ву структуру. Достижимая точность приближения наблюдаемых данных функцией регрессии зависит в значительной степени от размеров области ее определения. С увеличением размеров области, в пределах которой наблюдаемые данные не стремятся к нулю, точность снижается. Устранить этот недостаток можно, воспользовавшись методом локальной аппроксимации. При-
ближение наблюдаемых данных в этом случае обеспечивается последовательностью функций регрессии, каждая из которых задана на своем конечном интервале. Одновременно это позволяет аппроксимировать нестационарные сигналы и системы.
Вопрос сопряжения отдельных функций регрессии можно решить, если дополнить ограничения, представленные в задаче минимизации, условиями согласования значений этих функций и, возможно, значений их производных в узлах сопряжения [7].
Поставим задачу определения параметров в наблюдаемом сигнале х(^), являющемся аддитивной суммой оценивае-
мого сигнала з(г ) и акустической помехи л(г), которая также считается комплекс-
ной функцией. В дискретной форме этот сигнал имеет вид
Р (пш к +Фп ) х(к) = з(к) + л(к) = Xапе *0 + л(к) .
п=1
(15)
Здесь /0 - частота дискретизации.
Задачу определения параметров функции (15) можно решить также методом максимального правдоподобия, применение которого осложнено недостаточной надежностью априорной информации о распределениях помех. Метод максимального правдоподобия, как известно, в случае одинаково распределенных независимых гауссовых величин эквивалентен методу наименьших квадратов. В этом случае параметры функции (15) можно, в принципе, найти методами нелинейного
з - х
= (з - х, з - х) =
N
X
к=1
программирования, а именно посредством решения задачи минимизации:
{ш 1, 3, ф}= arg miп ||з - х|
{со1,а,ф}
|2
Если ввести векторы выборок выде-
т
ляемого сигнала 8 = (з(1) з(2) ...з(N)) и наблюдаемого сигнала
X = (х(1) х(2) ...х(N))Т,
векторы
а = (<1,...,<р )Т и ф = ^ф1,...,фр)Т, то норму можно записать в виде
Р /(Пш к +фп) ,9
Xапе *0 - х(к) |2.
п=1
Если к тому же помеха является коррелированной, с корреляционной функцией Я(г) и корреляционной матрицей
Я(к, п) =
Г ^(1,1) ^(1,2) Я(2,1) Я(2,2)
Л(1, К )> Я(2, К )
ч Я(К ,1) Я(К ,2) Я(К, К ) ,
то квадрат нормы в задаче минимизации принимает вид
и 2
.. 2 1 \ Т —1 *
||з - х| = (з - х, з - х = (з - х) Я (з - х)
Можно рассматривать метод максимального правдоподобия как задачу нелинейной регрессии, если принять, что оцениваемый сигнал имеет вид многочлена
а ф(г) по системе линейно независимых функций. Это приводит к параметрической задаче максимизации функции правдоподобия Р(уЗЬ |а), решение которой по наблюдаемым на интервале [а, Ь] данным у аЬ записывается в виде
а = aгgsup р( Уа |а).
(16)
Метод максимального правдоподобия может быть использован при определенных условиях, если наблюдаемые дан-
ные представляют собой последовательность независимых случайных величин с плотностью вероятности Р(Ук |а). При этом функция правдоподобия имеет вид
, п
Р(уа |а) = П Р(ук |а) . Если известны одно-
к=1
мерная плотность вероятности р( у к |а) и условная плотность вероятности Р( Ук+1У к, а), то функция правдоподобия описывается выражением
, п
Р( УЪа |а) = Р( У1 |а) П Р( у к+1У к, а). к=1
В связи с рассмотренными выше методами параметрической оптимизации и
2
а
локальной аппроксимации следует отметить, что по своему содержанию они практически эквивалентны проекционным ме-
тодам решения операторных уравнений, проекционным и интерполяционным методам анализа и расчета систем.
^руктуры абонентских устройств в телекоммуникационных системах связи, согласованных с моделями сигналов
В свободном пространстве эхосигна-лы образуются в результате отражений, обусловленных особенностями рельефа местности, расположения строений и крупногабаритных объектов. Значительные размеры территории приводят к большим задержкам распространения эхосиг-налов и снижению разборчивости речи. Кроме того, в акустическом поле с множеством различных каналов эха часто возни-
кают зоны молчания, что, в частности, необходимо учитывать при проектировании систем оповещения. Импульсные функции в каналах распространения эха можно при этом принять за некоторые постоянные коэффициенты передачи [8].
С учетом этого модель множественных отражений принимает вид
) = Iаки( - тк ) = Iаки^- БкТ) .
(17)
к=1
Задача заключается в определении параметров затухания а к и запаздывания
тк = БкТ для г каналов эхосигнала.
По условию сигнал и(^), многократные отражения которого формируют сигнал ), является известным и нестационарным. Это позволяет применить при оценивании указанных параметров технику корреляционного анализа.
Наличие в наблюдаемом сигнале встречной речи характерно для систем оперативно-командной связи и оповещения. Структура абонентского устройства с компенсационным каналом показана на рис. 5. Входной сигнал х(/) состоит из аддитивной суммы речевого сигнала и0^), сигнала множественных отражений у(^) и акустического шума ) :
г
х(Г) = и0 (I) + I аки(1 - Бк ) + у(1). к=1
В абонентском устройстве (рис. 5) имеются шумоподавляющее устройство, реализующее алгоритм адаптивного подавления акустических помех с формант-ным распределением полос режекции, адаптивный фильтр компенсации эха, блок вычисления долговременных параметров
к=1
эха - а к и Бк, блок формирования опорного сигнала адаптивного фильтра [8; 9] в виде
г
У к(п) = Iа ки0(п - Бк), к=1
сумматор, на выходе которого получают разность эха и опорного сигнала в соответствии с выражением
г
ивых(п) = Iек (п) + и0 (п) + у(п) , к=1
г г г
где I ек(п) = Iу(п - Бк) - I у к(п - Бк). к=1 к =1 к =1
Эхосигнал в замкнутых помещениях, в отличие от открытой территории, следует считать интегральной суммой компонентов, непрерывно распределенных по величине задержки. Достаточно хорошим приближением к действительности может явиться использование модели эхосигнала в виде (2). Входящие в эти выражения импульсные функции /(0) имитируют при этом непрерывное распределение задержек в пределах относительно небольших окрестностей [Тк -в, тк + в ], к = 1,...,г,
выделенных главных каналов
эха.
Рис. 5. Абонентское устройство с дополнительным компенсационным каналом
Заключение
Построение модели систем с акустической обратной связью, модели систем с дискретным и распределенным запаздыванием методами решения функционально-дифференциальных уравнений является одним из решений проблемы устранения неопределенности при повышении качества речевой связи и оповещения.
Рассмотренная выше в представлениях геометрической акустики модель множественных отражений предполагает необходимость нахождения ядер - неотрицательных весовых функций, характеризующих распределение интенсивности отраженных звуковых волн, содержащих соответствующую информацию, от величины их запаздывания. В отличие от этой модели, в модели резонансных мод акустическая среда рассматривается как резонансная система, описываемая в зависимости от геометрии озвучиваемой зоны набором соответствующих резонансных частот, и, соответственно, эхо и реверберация моделируются реакцией резонансной системы на звуковые колебания.
Методы численного решения функционально-дифференциальных уравнений,
характеризуемые в целом большим разнообразием, зависят от типа уравнения, которому отвечает квазиполином (10). Так, если степень первого слагаемого в (10) выше степени многочленов - коэффициентов в остальных слагаемых, то уравнение относится к разностному типу. Если же указанные степени равны, то уравнение относится к нейтральному типу.
Представленные выше модели систем с акустической обратной связью обеспечивают важное, но только первое приближение к действительности; в них, в частности, не учитывается возможность возникновения в трактах распространения звука нелинейных искажений, причиной которых могут стать в первую очередь нелинейность некоторых операторов и нестабильность распределения времени запаздывания по величине. Нелинейность операторов, во многом обусловленную нелинейностью источника звука (громкоговорителя), можно учесть в рамках, например, модели Вольтерра. Нестабильность запаздывания моделируется функционально-дифференциальными уравнениями с
переменными или случайными параметрами.
Таким образом, создание модели системы с акустической обратной связью позволяет решать вопросы повышения эффективности речевой связи и озвучивания методами создания систем обработки в телекоммуникациях. Также для решения
СПИСОК ЛИТЕРАТУРЫ
1. Erneux, T. Applied delay differential equations / T. Erneux. - New York: Springer Sciences, 2009. -199 p.
2. Kahrs, M. Applications of digital signal processing to audio and acoustics / M. Kahrs, K. Brandenburg.
- New York: Kluwer Academic Publisher, 2002. -572 p.
3. Хейл, Д. Теория функционально-дифференциальных уравнений / Д. Хейл. - М.: Мир, 1984. - 421 с.
4. Кропотов, Ю.А. Моделирование и методы исследований акустических сигналов, шумов и помех в системах телекоммуникаций: монография / Ю.А. Кропотов, В.А. Ермолаев. - М.Берлин: Директ-Медиа, 2016. - 256 с.
5. Беллман, Р. Дифференциально-разностные уравнения / Р. Беллман, К.Л. Кук. - М.: Мир, 1967. - 548 с.
6. Kuttruff, H. Room acoustics / H. Kuttruff. - London-New York: Spon Press, 2009. - 374 p.
7. Кропотов, Ю.А. Аппроксимация закона распределения вероятности амплитуд речевого сигнала
1. Erneux, T. Applied delay differential equations / T. Erneux. - New York: Springer Sciences, 2009. -199 p.
2. Kahrs, M. Applications of digital signal processing to audio and acoustics / M. Kahrs, K. Brandenburg.
- New York: Kluwer Academic Publisher, 2002. -572 p.
3. Hail, D. Theory of Functional-Differential Equations / D. Hail. - M.: Mir, 1984. - pp. 421.
4. Kropotov, Yu.A. Modeling and Investigation Methods of Acoustic Signals, Noises and Interferences in Telecommunication Systems: monograph / Yu.A. Kropotov, V.A. Yermolaev. - M.:-Berlin: Direct-Media, 2016. - pp. 256.
5. Bellman, R. Differential-Difference Equations / R. Bellman, K.L. Cook. - - M.: Mir, 1967. - pp. 548.
6. Kuttruff, H. Room acoustics / H. Kuttruff. - London-New York: Spon Press, 2009. - 374 p.
7. Kropotov, Yu.A. Approximation of law of speech signal amplitude probability distribution / Yu.A.
данной задачи можно применять метод минимизации функции потерь и метод максимального правдоподобия. Моделирование внешних помех эхосигналов можно рассматривать как средство преодоления априорной неопределенности, основанное на извлечении данных из результатов наблюдений [10].
/ Ю.А. Кропотов, Я.Ю. Кульков // Радиотехника.
- 2006. - № 11. - С. 63-66.
8. Кропотов, Ю.А. Методы проектирования алгоритмов обработки информации телекоммуникационных систем аудиообмена: монография / Ю.А. Кропотов, А.А. Парамонов. - М.-Берлин: Директ-Медиа, 2015. - 226 с.
9. Kropotov, Y.A. Algorithms for processing acoustic signals in telecommunication systems by local parametric methods of analysis / Y.A. Kropotov, V.A. Ermolaev // SIBCON 2015 - Proceedings: International Siberian Conference on Control and Communications. - 2015. - Р. 7147109.
10. Кропотов, Ю.А. Методы проектирования телекоммуникационных информационно-управляющих систем аудиообмена в сложной помеховой обстановке / Ю.А. Кропотов, А.А. Белов, А.Ю. Проскуряков, А.А. Колпаков // Системы управления, связи и безопасности. - 2015.
- № 2. - С. 165-183.
Kropotov, Ya.Yu. Kulkov // Radio Engineering. -2006. - No.11. - pp. 63-66.
8. Kropotov, Yu.A. Design Methods of Algorithms for Information Processing of Telecommunication Systems of Audio Exchange: monograph / Yu.A. Kropotov, A.A. Paramonov. - M.: Berlin: Direct-Media, 2015, - pp. 226.
9. Kropotov, Y.A. Algorithms for processing acoustic signals in telecommunication systems by local parametric methods of analysis / Y.A. Kropotov, V.A. Ermolaev // SIBCON 2015 - Proceedings: International Siberian Conference on Control and Communications. - 2015. - P. 7147109.
10. Kropotov, Yu.A. Methods for design of telecommunication information-control systems of audio exchange in complex interference situation / Yu.A. Kropotov, A.A. Belov, A.Yu. Proskuryakov, A.A. Kolpakov // Systems of Control, Communication and Safety. - 2015. - No.2. - No.2. - pp. 165-183.
Ссылка для цитирования:
Ермолаев, В.А. Идентификация параметров моделей акустических сигналов систем телекоммуникаций аудио-обмена в условиях помех и акустической обратной связи /В.А. Ермолаев, Ю.А. Кропотов //Вестник Брянского государственного технического университета. - 2019. - № 10. - С. 25 - 38. -DOI: 10.30987/article 5db95e85b8bb12.44709607.
Статья поступила в редакцию 09.09.19 Рецензент: д.т.н., профессор Брянского государственного
технического университета, главный редактор журнала «Вестник БГТУ»
Киричек А.В.
Статья принята к публикации 20. 09. 19.
СВЕДЕНИЯ ОБ АВТОРАХ:
Ермолаев Валерий Андреевич, к.т.н., доцент кафедры «Электроника и вычислительная техника» Муромского института (филиала) Владимирского государственного университета им. А.Г. и Н.Г. Столетовых, е-mail: kaf-eivt@yandex.ru.
Yermolaev Valery Andreevich, Can. Sc. Tech., Assistant Prof. of the Dep. "Electronics and Computer Engineering", Murom Institute (Branch) of Stoletovs State University of Vladimir, е-mail: kaf-eivt@yandex.ru.
Кропотов Юрий Анатольевич, д.т.н., профессор, зав. кафедрой «Электроника и вычислительная техника» Муромского института (филиала) Владимирского государственного университета им. А.Г. и Н.Г. Столетовых, е-mail: kaf-eivt@yandex.ru.
Kropotov Yury Anatolievich, Dr. of Sc. Tech., Prof., Head of the Dep. "Electronics and Computer Engineering", Murom Institute (Branch) of Stoletovs State University of Vladimir, е-mail: kaf-eivt@yandex.ru.
