CC BY
26ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2017, том 27, № 2, c. 83-90
МАТЕМАТИЧЕСКИЕ МЕТОДЫ г _
И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ
УДК 621.391
© В. Н. Якимов, А. В. Машков
ЗНАКОВЫЙ АЛГОРИТМ АНАЛИЗА СПЕКТРА АМПЛИТУД И ВОССТАНОВЛЕНИЯ ГАРМОНИЧЕСКИХ СОСТАВЛЯЮЩИХ СИГНАЛОВ В УСЛОВИЯХ ПРИСУТСТВИЯ НЕКОРРЕЛИРОВАННЫХ
ФОНОВЫХ ШУМОВ
Рассмотрена задача обнаружения и различения информативных составляющих многокомпонентных сигналов на фоне аддитивных некоррелированных шумов. Решение данной задачи основано на гармоническом анализе Фурье с использованием знакового аналого-стохастического квантования реализации исследуемого сигнала. Данный вид квантования позволяет осуществлять предельно грубое двухуровневое (бинарное) квантование без систематической погрешности независимо от статистических свойств анализируемых сигналов. Разработанные алгоритмы гармонического анализа обеспечивают высокую вычислительную эффективность. Практическая реализация этих алгоритмов в основном сводится к выполнению простейших арифметических операций суммирования и вычитания дискретных значений фильтрующих функций, вычисленных в моменты времени, когда результат знакового аналого-стохастического квантования меняет свой знак. Для принятия решения об обнаружении гармонических составляющих используется критерий превышения порогового значения оценками амплитуд гармонических составляющих. Результаты экспериментальных исследований показывают, что использование знакового аналого-стохастического квантования позволило эффективно обеспечить устранение влияния аддитивного некоррелированного шума и получить высокое качество восстановления непрерывных многокомпонентных полезных сигналов.
Кл. сл.: многокомпонентный сигнал, аддитивный шум, гармонический анализ, знаковое аналого-стохастическое квантование, знаковый сигнал
ВВЕДЕНИЕ. ПОСТАНОВКА ЗАДАЧИ
Проблемы обнаружения и выделения информативных составляющих в составе многокомпонентных сигналов на фоне внешних шумов возникают при решении многих практических задач, связанных с проведением экспериментальных исследований физических свойств объектов в самых различных областях науки и техники. В частности, методы обнаружения и анализа состава сигналов играют существенную роль при решении задач в биомедицинской и экологической инженерии, при обработке телекоммуникационных и акустических сигналов и т. п. Как правило, шумы являются результатом различного рода физических флюктуаций, носят случайный характер и являются аддитивными и некоррелированными. Следует отметить, что спектр шумов может занимать широкую полосу частот и перекрывать спектр информационного сигнала по всему его частотному диапазону. В этих условиях одновременная регистрация с основным информационным сигналом мешающих шумов может привести к его существенному искажению во временной области и получению аномальных результатов измерений. По-
этому задача надежного обнаружения и восстановления информационного сигнала без существенного искажения его полезных составляющих является актуальной задачей. При этом в качестве критерия оптимальности может служить степень искажения формы (спектрального состава) восстановленного сигнала.
Под обнаружением сигнала понимают констатацию факта его наличия. В случае работы в условиях присутствия внешних шумов задача обнаружения сигнала сводится к принятию решения: имеет ли место информационный сигнал, маскируемый шумом, или присутствует только шум. При положительном ответе возникает задача выделения информационных составляющих сигнала из смеси с шумом и восстановление сигнала. Решение этой задачи состоит в получении сигнала, который был бы наименее отличающимся по своим характеристикам от передаваемого [1, 2].
Ввиду того, что внешние шумы, а во многих случаях и полезный сигнал, носят вероятностный характер, современные подходы к анализу зашум-ленных непрерывных во времени многокомпонентных сигналов основываются на измерении их частотно-временных характеристик с использова-
нием методов статистической обработки. Одними из наиболее эффективных среди этих методов являются методы, направленные на проведение гармонического анализа. В процессе такого анализа осуществляется выделение гармонических компонент на фоне внешнего широкополосного шума.
С математической точки зрения наиболее адекватной формой представления многокомпонентного сигнала являются полигармонические функции, которые описываются суммой гармонических колебаний. Что же касается внешнего шума e(t), то в пределах времени анализа его обычно рассматривают как фоновый с непрерывным спектром, с нулевым математическим ожиданием и постоянной дисперсией cre2. В соответствии с этим математическая модель центрированного, т. е. имеющего нулевое математическое ожидание, полезно-
o
го сигнала x(t) в аддитивном шуме e(t) имеет вид
o o
Xe (t) = x(t) + e(t) =
K
= ^ (ak cos 2л ft + bk sin 2лk/0t) + e(t), (1)
k=1
где ak и bk — коэффициенты ряда Фурье; / — основная частота гармонических составляющих.
Модель (1) имеет нулевое математическое ожидание. Это является оправданной абстракцией, т. к. постоянная составляющая в наблюдаемом сигнале всегда может быть удалена путем использования операции центрирования. Допущение об отсутствии постоянной составляющей позволяет сосредоточить основное внимание на идентификации непосредственно гармонических составляющих.
С инженерной точки зрения более рациональной формой представления гармонических составляющих является их запись в полярных координатах. В этом случае (1) примет вид
o K
Xe(t) = Y,ak cos(2n/kt + (Pk) + e(t). (2)
k =1
В (2) для k-й гармонической составляющей с частотой /k = / будут соответствовать амплитуда ak и начальная фаза рк:
ak =Vak + bk , Pk = arctg ("bklak). (3)
Амплитуды ak образуют амплитудный спектр, а начальные фазы рк — фазовый спектр.
В соответствии с (1)-(3), уточняя общую задачу гармонического анализа, подчеркнем, что он связан с получением на дискретных частотах /к = k/0
оценок ак и Ък коэффициентов ряда Фурье и вычислением оценок ак амплитудного и фк фазового спектров на фоне широкополосного шума . Отметим, что амплитудный и фазовый спектры полигармонического сигнала являются дискретными.
Ввиду явного преимущества на практике для обработки сигналов широко используется цифровой гармонический анализ. Существующие цифровые алгоритмы гармонического анализа, основу которых составляют процедуры быстрого преобразования Фурье, позволяют обеспечить эффективную обработку цифровых отсчетов анализируемого сигнала. Однако в реальных условиях при низких значениях отношения сигнал/шум становится сложно различить полезный сигнал, и малое число разрядов АЦП оказывается недостаточным для точного представления сигналов в широкополосном шуме. Решение данной проблемы за счет простого увеличения числа разрядов АЦП приводит к усложнению процедур цифрового представления и обработки сигналов. При этом выбор числа уровней квантования может быть ограничен, а также может зависеть от разрешающей способности датчиков, регистрирующих сигналы. Таким образом, простое увеличение числа разрядов не может быть абсолютным решением качественного представления сигналов в цифровом виде. Кроме того, при равномерном шаге квантования помехозащищенность от шумов будет разной для полезных сигналов с малой амплитудой и с большой. Относительная погрешность квантования для сильных сигналов будет меньше, в то время как для слабых сигналов она будет больше. В соответствии с этим следует использовать методы первичного преобразования сигналов, обеспечивающие рациональное соотношение между точностью цифрового представления сигналов и простотой технического решения.
ЦИФРОВОЙ АЛГОРИТМ АНАЛИЗА СПЕКТРАЛЬНОГО СОСТАВА И ВЫДЕЛЕНИЕ ГАРМОНИЧЕСКИХ СОСТАВЛЯЮЩИХ МНОГОКОМПОНЕНТНЫХ СИГНАЛОВ
Для повышения эффективности обнаружения компонент информативных сигналов (особенно слабых составляющих) в аддитивной смеси с широкополосным шумом предлагается осуществлять цифровую обработку сигналов на основе знакового аналого-стохастического квантования, которое позволяет осуществлять предельно грубое двухуровневое (бинарное) квантование без систематической погрешности независимо от статистических свойств преобразуемых сигналов [1, 3-5].
Знаковое аналого-стохастическое квантование представляет собой нелинейную операцию веро-
ЗНАКОВЫЙ АЛГОРИТМ АНАЛИЗА СПЕКТРА АМПЛИТУД.
85
ятностного округления реализации непрерывного сигнала, где в качестве меры квантования выступает вспомогательный случайный сигнал £ (1). В результате такого квантования получаем ограниченный по уровню знаковый сигнал:
z(t) =
(4)
+1, если Хе (1) + ) > 0;
О
-1, если Хе(1) + ) < 0.
Вспомогательный сигнал ) принимает значения в пределах от -£тах до +£тах и имеет
О
равномерное распределение. Здесь £тах >| Хе (1) |тах,
О
где | Хе (1) |тах — максимально возможное по абсолютной величине значение, которое может
О О
принять Хе (1) = х(1) + е(1).
Оценки коэффициентов ряда Фурье ак и Ък в процессе гармонического анализа в пределах интервала времени 10 < 1 < 10 + Т будем искать в следующем виде [6]:
10 +Т
ак = 2£тахТ- | ¿(1)со$(2лк^^,
'0
10 +Т
4 = 2£тхТ- | ¿(1 )ът(2лк^)dt,
(5)
(6)
а,
к = 2£тхТ-1 £ )0С8(2^к/01 )dt,
\ = 2£тхТ-1 £ } ¿(1 )яп(2якГ01 )dt.
Знаковый сигнал 2(1) принимает значения -1 и +1, которые чередуются строго последовательно во времени и остаются постоянными в переделах интервалов ^ < 1 < 1*+1. Поэтому достаточно знать только его начальное значение 2(10). Принимая это во внимание, интегралы в (7) и (8) можно проинтегрировать аналитически, и оценки ак и Ък примут вид:
= 2$тЛккГ *„)£ (-1)'-1 sш(2f),
1=1
К = $тЛк*Г 2(о£(-1)1 я, ^(2жк/1),
(9) (10)
к =
1 , при 1 = 0, 1 = I ;
2 , при 1 < 1 < I -1 .
где к = 1,2,3,...,К ; ^ = 1/Т; со8(2жк/01) и 81п(2^к/01) — фильтрующие функции.
Величина ^ = 1/ Т соответствует величине предельного разрешения по частоте дискретных гармонических составляющих.
Обозначим 12 = 10 и = 10 + Т. Не нарушая общности и строгости математических рассуждений, будем считать, что 10 = 0 и, следовательно,
= Т . Пусть множество |, где 1 < 1 < I -1 и = , определяет те моменты времени, в которые знаковый сигнал 2(1) в пределах интервала 0 < 1 < ^ меняет текущее значение на противоположное значение. Тогда (5) и (6) можно записать следующим образом:
В (9) и (10) перейдем к числовому представлению моментов времени 1*. Тогда будем иметь = ' А1, где А1 — период следования счетных импульсов и 1 = 1,2,3,...,I -1. При этом = 'А1 и ^ = Т = ЫА1. В соответствии с этим интервалу времени < 1 < , в пределах которого осуществляется вычисление оценок ак и Ък, будет соответствовать его дискретный аналог ' А1 < пА1 < г]2 А1, или просто ' < п < . Кроме того, если 10 = 0, то ] = 0 и ] = N .
Подставив в (9) и (10) числовые оценки моментов времени ^, окончательно получаем
ак = 2£тЛк*У1 2(о£(-1Г^1п(2^-'), (11)
Ък = 2£пах(к^Г ¿(О х
(
х
0
1 + (-1) I I -1 ^
—+ £ (-1)" cos(2nkN)
V 2 1=1 /
(12)
(7)
(8)
Из (12) следует, что если I является нечетным, то будем иметь
Ък = 2£шЛкп)- 2(10)£(-1)' cos(2лkN').
1=1
Таким образом, вычисление оценок коэффициентов ряда Фурье ак и Ък в цифровом виде в конечном итоге свелось к формированию двух матриц размером К х (I -1) с дискретными значениями фильтрующих функций sin(2жkN~l'i) и cos(2л kN), обработка которых осуществля-
ется согласно (11) и (12). При этом важной особенностью такой обработки является то, что основными операциями являются простейшие арифметические операции суммирования и вычитания. Практически она не требует осуществлять операции цифрового умножения.
Для принятия решения об обнаружении гармонических составляющих используется критерий, основу которого составляет принцип превышения оценок амплитуд гармонических составляющих порогового значения. Согласно этому критерию устанавливается уровень значимости гармонических составляющих в зависимости от характера решаемой задачи. Если значение полученной оценки амплитуды ак гармонической составляющей с частотой 2nk/N , где к = 1,2,3,...,K, меньше установленного уровня е0, то считается, что данная гармоническая составляющая не является значимой, и она не принимается во внимание.
РЕЗУЛЬТАТЫ ЧИСЛЕННОГО ЭКСПЕРИМЕНТА
На основе соотношений (11) и (12), а также критерия принятия решения с учетом задания уровня значимости гармонических составляющих разработан цифровой алгоритм обнаружения сигналов на фоне внешних широкополосных шумов. Моделирование практической реализации данного алгоритма проводилось с использованием имита-
ционного моделирования процедуры знакового аналого-стохастического квантования.
В процессе формирования модели (2) генерировались гармонические составляющие с нормированными частотами значения которых изменялись от нуля до 0.5. Амплитуды ак принимались от нуля до единицы. Начальные фазы рк задавались с помощью генератора равномерно распределенных величин в пределах —п < рк < п . Некоррелированный аддитивный шум г($) представлял собой белый шум с нулевым математическим ожиданием и единичной дисперсией. В частности модель (2) содержала одиннадцать гармонических составляющих, параметры которых приведены в таблице. В этой же таблице приведены результаты вычисления оценок ак и ррк, а также представлены значения относительных погрешностей 8а и 8р для этих оценок.
о
На рис. 1 показан вид полезного сигнала ), соответствующий набору гармонических составляющих из таблицы. На рис. 2 представлен сигнал
о
хе ((), определяемый моделью (2). На рис. 3 приведен результат гармонического анализа, представляющий собой нормированную оценку спектра амплитуд ак, которые вычислены с разрешением 0.001 единиц нормированной частоты.
o
Параметры гармонических составляющих модели сигнала x(t) и их оценки
k ак «к 5 , % а > Рк, рад 4>к 5 , % ср ' а к /, дБ
1 0.225 0.1 0.1012 1.20 -1.2961 -1.3226 2.04 -20.00
2 0.23 0.2 0.2008 0.40 1.3905 1.3800 -0.76 -13.98
3 0.235 0.3 0.2980 -0.67 -2.3089 -2.3199 0.48 -10.46
4 0.24 0.4 0.4008 0.20 -1.5957 -1.6012 0.34 -7.96
5 0.245 0.5 0.5021 0.42 -0.6312 -0.6331 0.30 -6.02
6 0.25 1 1.0019 0.19 -0.0833 -0.0857 2.88 0.00
7 0.255 0.5 0.4994 -0.12 -0.6962 -0.7080 1.69 -6.02
8 0.26 0.4 0.4006 0.15 1.3449 1.3486 0.28 -7.96
9 0.265 0.3 0.2970 -1.00 2.6846 2.6777 -0.26 -10.46
10 0.27 0.2 0.1986 -0.70 1.0151 1.0315 1.62 -13.98
11 0.275 0.1 0.1016 1.60 -1.7407 -1.7742 1.92 -20.00
знаковый алгоритм анализа спектра амплитуд...
87
Рис. 1. Модель полезного сигнала х(¿)
Рис. 2. Модель сигнала хе )
Рис. 3. Нормированная оценка спектра амплитуд ак
Рис. 4. Нормированная оценка спектра амплитуд ак в пределах частотного диапазона, занимаемого гармоническими составляющими сигнала
Для большей наглядности разрешающей способности разработанного алгоритма на рис. 4 представлена нормированная оценка спектра амплитуд ак в пределах частотного диапазона, занимаемого гармоническими составляющими по-
о
лезного сигнала х^). На рис. 5 приведена оценка спектра амплитуд ак в децибелах. Хорошо видно, что спектральная оценка шума после обработки
о
исходного сигнала хе ) остается на достаточно низком уровне (в среднем около -50 дБ) по отношению к оценкам гармонических составляющих. Это свидетельствует о высокой устойчивости разработанного алгоритма к внешним некоррелированным аддитивным шумам. На рис. 6 пред-
ставлен вид непрерывного восстановленного
о
во временной области сигнала Хв (^) , полученного по оценкам ак.
Приведенные результаты моделирования свидетельствуют о том, что положения оценок гармонических составляющих в спектре точно соответствуют табличным данным. Ложных линий в спектре не наблюдается, спектральные линии четко различимы. Они не маскируют друг друга. Значения оценок амплитуд ак и начальных фаз Ррк получены с относительными погрешностями, которые не превышают соответственно двух и трех процентов, что позволяет говорить о высокой точности их вычисления в статистическом смысле.
0,06 0.09 0.1? 0.15 0.13 0.21 0,24 0.27 0 3 0.33 0 36 0.39 0А2 0 .«5
Нормированная частота
Рис. 5. Оценка спектра амплитуд ак в децибелах
Рис. 6. Восстановленный сигнал Хв (1)
ЗАКЛЮЧЕНИЕ
На основе использования знакового аналого-стохастического квантования разработан цифровой алгоритм, который достаточно эффективно устраняет влияния внешнего некоррелированного аддитивного шума и позволяет получить высокое качество оценивания спектра гармонических составляющих непрерывных многокомпонентных
информационных сигналов. Кроме того, данный алгоритм практически исключает необходимость выполнения многоразрядных операций цифрового умножения, которые свойственны классическому цифровому гармоническому анализу Фурье. Как было отмечено выше, в нашем случае основными операциями являются простейшие арифметические операции суммирования и вычитания дискретных значений фильтрующих функций, соот-
ЗНАКОВЫЙ АЛГОРИТМ АНАЛИЗА СПЕКТРА АМПЛИТУД... 89
ветствующих моментам времени, в которые результат знакового аналого-стохастического квантования меняет свой знак на противоположный. Уменьшение числа операций умножения ведет к сокращению временных затрат, что обеспечивает вычислительную эффективность гармонического анализа.
Работа выполнена при финансовой поддержке РФФИ (проект № 16-08-00269 -А).
дискретного знакового преобразования непрерывных сигналов // Вест. Самар. гос. техн. ун-та. Сер. Технические науки. 2000. Вып. 8. С. 190-192.
6. Якимов В.Н. Цифровой гармонический анализ многокомпонентных случайных процессов // Измерительная техника. 2006. № 4. С. 22-26.
СПИСОК ЛИТЕРАТУРЫ
Богданович В.А., Вострецов А.Г. Теория устойчивого обнаружения, различения и оценивания сигналов. М.: Физматлит, 2004. 320 с.
Белецкий Ю.С. Методы и алгоритмы контрастного обнаружения сигналов на фоне помех с априори неизвестными характеристиками. М.: Радиотехника, 2011. 436 с.
Мирский Г.Я. Характеристики стохастической взаимосвязи и их измерения. М.: Энергоиздат, 1982. 320 с. Макс Ж. Методы и техника обработки сигналов при физических измерениях. Т. 1. М.: Мир, 1983. 312 с. Якимов В.Н. Математическое представление потоков
Самарский государственный технический университет, Россия
Контакты: Машков Андрей Валерьевич, [email protected]
Материал поступил в редакцию: 17.04.2017
ISSN 0868-5886
NAUCHNOE PRIBOROSTROENIE, 2017, Vol. 27, No.2, pp. 83-90
THE BINARY ALGORITHM FOR THE ANALYSIS OF THE SPECTRUM AMPLITUDE AND RECOVER OF HARMONIC COMPONENTS SIGNALS IN THE PRESENCE OF UNCORRELATED BACKGROUND NOISE
V. N. Yakimov, A. V. Mashkov
Samara State Technical University, Russia
The problem of detecting and distinguishing between the informative components of multicomponent signals in the presence of additive uncorrelated noises. The solution of this problem is based on harmonic analysis Fourier series using the sign-function analog-stochastic quantization of the investigated signal. This kind of quantization makes it possible to carry out extremely coarse two-level quantization without a systematic error, regardless of the statistical properties of the analyzed signals. A developed algorithm for harmonic analysis provides a high computational efficiency. The practical implementation of these algorithms basically reduces to performing the simplest arithmetic operations of summation and subtraction of the discrete values of the filtering functions calculated at the instants of time when the result of sign-function analog-stochastic quantization changes its sign. For the decision about the detection of harmonic components, we use a criterion for exceeding of a threshold value by estimating the amplitudes of the harmonic components. The results of experimental studies show that the use of sign-function analog-stochastic quantization has allowed to effectively eliminate the effect of the additive uncorrelated noise and to obtain high quality recover of continuous multicomponent signals useful.
Keywords: multicomponent signals, additive noise, harmonic analysis, sign-function analog-stochastic quantization, binary signal
REFERENСES
1. Bogdanovich V.A., Vostrecov A.G. Teoriya ustojchivogo obnaruzheniya, razlicheniya i ocenivaniya signalov [Theory of steady detection, distinction and estimation of signals]. Moscow, Fizmatlit Publ., 2004. 320 p. (In Russ.).
2. Beletckij Yu.S. Metody i algoritmy kontrastnogo obnaruzheniya signalov na fone pomekh s apriori neizvestnymi harakteristikami [Methods and algorithms of contrast detection of signals against the background of hindrances with a priori unknown characteristics]. Moscow, Radiotekhnika Publ., 2011. 436 p. (In Russ.).
3. Mirskij G.Ya. Harakteristiki stohasticheskoj vzaimosvyazi i ih izmereniya [Characteristics of stochastic interrelation and their measurement]. Moscow, Energoizdat Publ., 1982. 320 p. (In Russ.).
Contacts: Mashkov Andrey Valer'evich, [email protected]
4. Maks Zh. Metody i tekhnika obrabotki signalov pri fizi-cheskih izmereniyah. T. 1 [Methods and technology of processing of signals at physical measurements. Vol. 1]. Moscow, Mir Publ., 1983. 312 p. (In Russ.).
5. Yakimov V.N. [Mathematical representation of streams of discrete sign transformation of continuous signals]. Vestnik samarskogo gosudarstvennogo tekhnicheskogo universiteta. seriya tekhnicheskie nauki [Bulletin of the Samara State Technical University. Technical Sciences Series], 2000, vol. 8, pp. 190-192. (In Russ.).
6. Yakimov V.N. Digital harmonic analysis of multicomponent random processes. Measurement Techniques. New York, Springer, 2006, vol. 49, no. 4, pp. 341-347. Doi: 10.1007/s11018-006-0112-x.
Article received in edition: 17.04.2017
