CC BY
10
Труды БГТУ, 2016, № 6, с. 117-121
117
УДК 681.53
И. Г. Сухорукова, Д. А. Гринюк, И. О. Оробей
Белорусский государственный технологический университет
ВЛИЯНИЕ УСЛОВИЙ ФИЛЬТРАЦИИ И СГЛАЖИВАНИЯ В ИНФОРМАЦИОННЫХ КАНАЛАХ НА КРИТЕРИЙ СЕРИЙ
В статье производится оценка возможностей использования критерия серий для адаптивной фильтрации. Выполнено моделирование работы измерительных каналов в широком диапазоне частот и в условии различных помех, которые характерны для радиоэлектроники. Тестирование производилось при помехах со следующими распределениями случайной величины: нормальным, равномерным, экспоненциальным, треугольным, трапециидальным, арктангенсуальным. Случайные величины генерировались встроенными средствами Matlab. После смешения с полезным гармоническим сигналом временной ряд сглаживался тремя типами фильтров: бегущего среднего, медианой выборки, низких частот. В процессе моделирования параметры фильтров изменялись в широком диапазоне. Длина наблюдения серий варьировалась от 8 до 64. Проводилась оценка влияния условий на основные показатели критерия: количество серий, математическое ожидание, дисперсия, полный и упрощенный коэффициенты стационарности.
На основании полученных результатов сделаны выводы о возможности использования основных параметров критерия серий для адаптивной фильтрации. Несмотря на различный тип помех серий в широком диапазоне динамики распределены по нормальному закону. Оценка коэффициентов корреляции между параметрами критерия серий и полезным сигналом позволяет рекомендовать использование количества серий и упрощенного критерия стационарности только при высокой динамике. Из-за случайного характера основных параметров критерия серий для адаптации лучше применять триггерные варианты. Предложено несколько вариантов схем адаптивной фильтрации.
Ключевые слова: адаптивная фильтрация, критерий серий, стационарность.
I. G. Suhorukova, D. A. Hryniuk, I. O. Orobei
Belarusian State Technological University
INFLUENCE OF CONDITIONS OF FILTERING AND SMOOTHING IN INFORMATION CHANNELS RUNS TESTS
The article considers the runs tests for the adaptive filtering. Simulation of the measurement channels operation in a wide range of frequencies and in different interference conditions that are in different interference conditions that are typical for electronics, was performed. Working conditions were created at various noises as in. Testing was made in case of noises from distribution of a random variable: normal, uniform, exponential, triangular, trapezoidal; arctangent. Random variables were generated by the built-in means of Matlab. After mixing with the useful harmonic signal the time series was filtered with three types of filters: the running average, median filter, the filter of low frequency. In the course of simulation parameters of filters changed in a broad range. The length of observation of series varied from 8 to 64. The impact assessment of conditions on the main indices of the criterion was carried out: quantity of series, mathematical expectation, dispersion, the complete and simplified stationarity coefficient. Influence of quantization on level on quantity of series is carried out.
Based on the results of the analysis the conclusions are made about the possibility of using of the main parameters of the criterion for adaptive filtering. Various types of noise lead to a series of normal distribution. The normal distribution is destroyed by high dynamics of deterministic time series. The correlation coefficients show that at low frequencies to adapt better to use a mathematical expectation, dispersion or complete stationarity, but at high frequencies is used quantity of series. Because of the accidental character of key parameters of the criterion of series for adaptation it is better to use trigger options. Several versions of schemes for the adaptive filtering are offered.
Key words: adaptive filtering, runs tests, stationary.
Введение. Динамика временных рядов в цифровой фильтрации может быть обнаружена критерием серий или критерием Вальда - Вольфови-ца (ВВ) [1]. Данный критерий относится к непараметрическим и широко используется в биологиче-
ских, физических, социальных и экономических исследованиях. В работе [2] были проведены исследования использования ВВ с целью оптимизации его работы, а в [3] выполнена апробация на промышленном логическом контроллере.
Одним из актуальных вопросов является работа принципа адаптации при условии не га-уссовского распределения шумов на входе фильтра. В метрологической практике при анализе работы измерительных приборов используют ограниченный ряд стандартных аппроксимирующих функций распределения [4]: нормальную, равномерную, треугольную, арксину-соидальную, экспоненциальную и Рэлея.
Математическое моделирование производилось в Ма1;1аЬ. При постановке математического эксперимента случайные данные выбранных распределений генерировались таким образом, чтобы их мощности на фиксированном периоде времени отличались не более чем на 1%.
Помехи смешивались с гармоническим сигналом (ГС) с частотами от 0 до 0,03 Гц и пропускались через фильтры. Соотношение между амплитудой полезного сигнала и величиной дисперсии помехи не ниже 10 по нормальному распределению. В качестве фильтров анализировался фильтр бегущего среднего и с медианной выборкой [1] с длиной от 3 до 7 и цифровые фильтры 1-3 порядков, фильтр Баттерворта. Производилась оценка влияния эффекта квантования сигнала по уровню на параметры критерия серий.
Количество серий N варьировалось от 8 до 64. Для всех вариантов определялись количества серий г, математическое ожидание ц, дисперсия с и оценки стационарности у1, у2:
2 N N- 1
ц=-+1;
N
= 2N N~ (2NN - N) _
N 2( N -1)
(
Y = Yi = exp
(r K a2
2 Л
Y = Y2 = -,
(1)
(2)
(3)
где N - количество положительных и отрицательных разниц до и после фильтра; К - коэффициент сжатия-растяжения вероятности статистической независимости.
Результаты и обсуждение. Для медианного (МФ), бегущего среднего (БС), фильтров первого (ФНЧ1) и третьего (ФНЧ3) порядков в отсутствии ГС и на малых частотах, при минимальных параметрах фильтров (МтТипе) имеются отклонения от нормального распределения больше, чем при максимальных (МахТипе) (рис. 1-3). Графики (рис. 1-4) построены функцией распределения плотности вероятностей (РПВ) методом непараметрического сглаживания к^еп8йу(). В целом, чем больше сдвиг фаз фильтра, тем больше форма распределения случайных данных соответствует идеальному
колоколу Гаусса (рис. 3 и 4). Наличие у помехи (как, например, при экспоненциальном распределении) ненулевого математического ожидания влияет и на ц серий.
0,2 0,15 0,1 0,05 0
-0,05 -0,1
■ N / 2 / ■ МФ ■ БС Ль. - 1 \ —2' — з - pirf^b,"..... 4 -
ФНЧ1
20
25
30
35
40
45
Рис. 1. РПВ без ГС с МтТипе при N = 64. Часть 1: 1 - нормальное; 2 - равномерное; 3 - экспоненциальное; 4 - релейное
0,25 0,2 0,15 0,1 0,05 0
-0,05 -0,1
■ N / 2 / ■ МФ ^У \ — 2 ■ \ — 3 - к™'"" 4 -
' БС
■ ФНЧ1
20 25 30 35 40 45 . РПВ без ГС с MinTune при N = 64. Часть 2 1 - нормальное; 2 - треугольное; трапециидальное; 4 - арктангенсуальное
0,15 0,1 0,05 0
-0,05 -0,1
- ' N/2^ — 2 ■ N. — 3 -
■ мф ^
■ -¿Г БС ........ 4 -
/¿Г
-rf—.
20
25
30
35
40
45
Рис. 3. РПВ без ГС с МахТипе при N = 64. Часть 1: 1 - нормальное; 2 - равномерное; 3 - экспоненциальное; 4 - релейное
0,15
0,1 • МФ
0,05 ' БС
0
-0,05 . ФНЧ1
-0,1 ■ ФНЧ3
20
25
30
35
40
45
Рис. 4. РПВ без ГС с MaxTune при N = 64. Часть 2: 1 - нормальное; 2 - треугольное; 3 - трапециидальное; 4 - арктангенсуальное
r
И. Г. Сухорукова, Д. А. Гринюк, И. О. Оробей
119
Кривая для МФ, работа которого основана на случайном механизме, наиболее приближена к кривой нормального распределения и ц = N / 2. Такой эффект происходит и при увеличении параметров фильтров.
При повышении частоты математическое ожидание серий уменьшается (рис. 5 и 6). На некоторой частоте, значение которой зависит от соотношения между дисперсией помехи и амплитудой полезного сигнала, ц приближается к нулю и колокол Гаусса «разрушается». Первоначально уход от гауссовского распределения происходит у цифровых фильтров по причине проявления фазовых сдвигов. При этом равномерное распределение, трапеции-дальное и арктангенсуальное, которые имеют большую плотность в области больших отклонений, дольше сохраняют нормальное распределение.
0,3 0,25 0,2 0,15 0,1 0,05 0
-0,05 -0,1
-10 0
10 20 30 40
Рис. 5. РПВ при средних частотах с МахТипе: 1 - нормальное; 2 - равномерное; 3 - экспоненциальное; 4 - релейное
0,2 0,15 0,1 0,05 0
-0,05 -0,1
0
5
10 15 20 25 30 35 40
Рис. 6. РПВ при больших частотах с МшТипе: 1 - нормальное; 2 - равномерное;
3 - экспоненциальное; 4 - релейное
Если при максимальной частоте полезного сигнала при длине 3 МФ и БС значение математического ожидания еще около 32, то уже при длине 4 происходит существенное отклонение, а при 7 значение математического ожидания смещается к 22-23.
При снижении г до 48 графики РПВ деформируются подобно кривым на рис. 6 для ФНЧ3. При длине г в 16 провалы опускаются до нуля. В то же время ц = N / 2. Это хорошо согласуется с положением, что при моделировании случайных процессов рекомендуется выбирать длину выборки не менее 100. Поэтому алгоритм использования критерия серий, работа которого исходит из нормального распределения серий, должен иметь длину наблюдения 64 и выше.
Изменения параметров серий для различной динамики, типов фильтров и их настроек можно видеть на рис. 7-11. Кривая 5 на графиках приведена без масштаба для визуального сравнения с остальными параметрами.
40 30 20 10 0 1
0,5
0
20
20,1 20,2 20,3 20,4 /, кс
Рис. 7. Изменение показателей критерия серий при слабой динамике 5 для БС
30 20 10
20,0 20,04 20,08 20,12 20,16 кс
Рис. 8. Изменение показателей критерия серий при короткой длине БС
30 20 10
0,5
20,0 20,04 20,08 20,12 20,16 /, кс
Рис. 9. Изменение показателей критерия серий при максимальной длине БС
30 20
20,0 20,04 20,08 20,12 20,16 кс
Рис. 10. Изменение показателей критерия серий при минимальной длине МФ
30
20
10 0
0,8
0,4 0
V ц /
''.Л \ 4 Г-' . . . ^ .
убн 1ц) "Ч.....ш /------------------------
20,0 20,04 20,08 20,12 20,16 /, кс
Рис. 11. Изменение показателей критерия серий для ФНЧ
Графики показывают, что в той или иной степени г, ц, с, у! и у2 могут выступать в качестве параметра оценки динамики тренда для адаптации. Установить же аналитическую зависимость между ними и динамикой сигнала проблематично. В то же время можно рекомендовать использование интервальной адаптации для ступенчатой перенастройки фильтров.
Влияние коэффициента сжатия-растяжения вероятности статистической независимости на у! отражено на рис. 12. Выбор значения К зависит от структуры адаптации и той гипотезы (стационарный процесс или нестационарный), которая в данный момент интересует.
К = 20 К = 2,5
К = 1 К = 0,75
К = 0,3
20 20,1 20,2 20,3 20,4 Рис. 12. Влияние К на у1
t, кс
Увеличение К склоняет у1 в сторону гипотезы стационарности, а уменьшение - нестационарности. Путем изменения К можно реализовать триг-
герную структуру с зоной нечувствительности для уменьшения зависимости от шума у1.
При одинаковом входном сигнале фильтры по-разному формируют серии и их показатели. Даже близкие по принципу работы МФ и БС дают разную картину. Цифровые фильтры имеют более узкий диапазон частот, при котором можно говорить о стабильном определении динамики ввиду большей зависимости фазы от частоты. БС и МФ имеют больший диапазон детектирования динамики, но хуже сглаживают. Возможен вариант построения наблюдателя слежения за динамикой для своевременной перенастройки фильтра основного канала. Для этого можно использовать и ФНЧ с автоподстраи-ваемой в зависимости от динамики малой постоянной времени по сравнению с основным каналом. Данный подход вместе с изменением К, как было указано выше, скорее всего окажется более выигрышным, чем поиск зависимости между уь у2, ц, г, с и 5 динамикой тренда.
Производилась оценка корреляции между Уь у2, ц, г, с и 5. Коэффициент корреляции (КК) выше между 5 и уь ц, с в среднем на четверть, чем между у2, г и 5. При этом для МФ и БС КК при частотах 10-4 начинает быстро подниматься и остается на одном уровне до 10-2, а затем резко возрастает, когда количество серий г сдвигается ближе к нулю, чем к значению N / 2. Если длина БС на среднем диапазоне частот незначительно влияет на корреляцию, то при больших частотах при некоторой длине БС коэффициент корреляции достигает максимума для ц и с, а затем начинает падать, тогда как для остальных параметров КК начинает резко расти. Для ФНЧ при частотах 10-4 в зависимости КК от постоянной времени фильтра имеется пологий максимум. По мере роста частоты максимум становится резким и смещается в сторону малых постоянных времен. На больших частотах КК для у2, г становятся выше, чем для остальных параметров.
Для трендов г на среднем значении частот для ФНЧ, где параметр колеблется в широком диапазоне, выполнено разложение в ряд Фурье. Основная гармоника больше по амплитуде остальных, однако вторая и последующие имеют сравнимое с ней значение, что говорит о нелинейности преобразования. В целом тренд г для ФНЧ при некоторых частотах напоминает модуль гармонического сигнала.
При прохождении через квантователь по уровню критерий стационарности у1 и у2 сдвигается в сторону 1 со снижением амплитуды вариации, его зависимость от динамики уменьшается. Тренд г смещается в сторону N / 2, снижает амплитуду и меньше реагирует на характеристику шума. Тренды ц и с также очищаются от шума,
И. Г. Сухорукова, Ä. А. Гринюк, И. О. Оробей
121
но увеличивают амплитуду, становятся более подобны на модуль гармонического сигнала. Случайные процессы с экспоненциальным распределением ввиду больших отклонений менее подвержены влиянию квантования.
Заключение. Критерий серий показал слабую зависимость от распространенных распределений шумов в радиоэлектронике. Ненулевое ма-
тематическое ожидание шума оказывает воздействие на параметры распределения серий. Вследствие подверженности основных показателей серий влиянию шума адаптацию лучше строить на основании триггерных схем, которые были предложены выше. В условиях работы с малым количеством разрядов АЦП на диапазоне шума лучше использовать ц или с для адаптации фильтров.
Литература
1. Айвазян С. А., Енюков И. С., Мешалкин Л. Д. Прикладная статистика: Основы моделирования и первичная обработка данных. М.: Финансы и статистика, 1983. 471 с.
2. Гринюк Д. А., Оробей И. О., Сухорукова И. Г. Эффективность работы адаптации фильтра на критерии серий // Труды БГТУ. 2011. № 6: Физ.-мат. науки и информатика. С. 107-111.
3. Гринюк Д. А., Оробей И. О., Сухорукова И. Г. Адаптация критерия серий к применению в управлении технологическими процессами // Труды БГТУ. 2014. № 6: Физ.-мат. науки и информатика. С. 92-95.
4. Басов В. Г. Информационно-измерительные системы. Основы проектирования. В 3 ч. Ч. 2. Минск: БГУИР, 2007. 78 с.
References
1. Ayvazyan S. A., Enyukov I. S., Meshalkin L. D. Prikladnaya statistika: Osnovy modelirovaniya i pervichnaya obrabotka dannykh [Applied Statistics: Basics of modeling and primary data processing]. Мoscow, Finansy i statistika Publ., 1983. 471 p.
2. Hryniuk D. A., Orobei I. O., Suhorukova I. G. Efficiency of the filter to adapt to a runs test. Trudy BGTU [Proceedings of BSTU], 2011, no. 6: Physical-mathematical sciences and informatics, pp. 107-111 (In Russian).
3. Hryniuk D. A., Orobei I. O., Suhorukova I. G. Adaptation runs test for use in process control. Trudy BGTU [Proceedings of BSTU], 2014, no. 6: Physical-mathematical sciences and informatics, pp. 92-95 (In Russian).
4. Basov V. G. Informatsionno-izmeritel'nyye sistemy. Osnovy proektirovaniya. V 3 ch. Ch. 2 [Information-measuring system. Fundamentals of design. In 3 parts. Part 2], Minsk, BGUIR Publ., 2007. 78 p.
Информация об авторах
Сухорукова Ирина Геннадьевна - ассистент кафедры информационных систем и технологий. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). Е-mail: [email protected]
Гринюк Дмитрий Анатольевич - кандидат технических наук, доцент, доцент кафедры автоматизации производственных процессов и электротехники. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). Е-mail: [email protected]
Оробей Игорь Олегович - кандидат технических наук, доцент, доцент кафедры автоматизации производственных процессов и электротехники. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). Е-mail: [email protected]
Information about the authors
Suhorukova Irina Gennadyevna - assistant lecturer, the Department of Information Systems and Technology. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Hryniuk Dmitry Anatolyevich - PhD (Engineering), Assistant Professor, Assistant Professor, the Department of Automation of Production Processes and Electrical Engineering. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Orobei Igor Olegovich - PhD (Engineering), Assistant Professor, Assistant Professor, the Department of Automation of Production Processes and Electrical Engineering. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: [email protected]
Поступила 02.03.2016
