ПОВЫШЕНИЕ ТОЧНОСТИ ОЦЕНКИ ЧАСТОТЫ СИНУСОИДАЛЬНОГО СИГНАЛА С ИСПОЛЬЗОВАНИЕМ НЕЛИНЕЙНОГО ФИЛЬТРА
СИСТЕМЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ
УДК 681.5.03
А. А. БОБЦОВ, А. В. КРЫЛОВ, А. А. ПЫРКИН
ПОВЫШЕНИЕ ТОЧНОСТИ ОЦЕНКИ ЧАСТОТЫ СИНУСОИДАЛЬНОГО СИГНАЛА С ИСПОЛЬЗОВАНИЕМ НЕЛИНЕЙНОГО ФИЛЬТРА
Рассматривается новая схема идентификации частоты измеряемого синусоидального сигнала с использованием нелинейного фильтра. В отличие от большинства известных аналогов, предлагаемый алгоритм позволяет получить более точные оценки за меньшее время.
Ключевые слова: идентификация, гармонический сигнал, фильтрация.
Введение. Актуальность проблемы идентификации частоты измеряемого синусоидаль-
ного сигнала в условиях влияния внешних возмущений и шумов подтверждается наличием
большого числа публикаций, посвященных идентификации частоты синусоидального или
частот мультигармонического сигналов [1—17]. Однако в большинстве указанных работ не
обсуждается задача качества идентификации, а именно ускорения вычислительной процеду-
ры и повышения точности оценок при неучтенных возмущениях, присутствующих в канале
измерения полезного сигнала. Анализ и рекомендации по увеличению скорости сходимости
оценки частоты синусоидального сигнала к истинному значению были опубликованы в рабо-
тах [11, 14—16], а их робастность к неучтенным возмущениям и помехам исследована в ра-
ботах [13—16]. В частности, в работе [15] была предложена новая схема идентификации час-
тоты измеряемого синусоидального сигнала в условиях неучтенных возмущающих воздейст-
вий и шумов, присутствующих в канале измерения полезного сигнала. Одно из отличий схе-
мы, предложенной в работе [15], от известных аналогов заключается в возможности компен-
сации неучтенных возмущений путем подбора коэффициентов алгоритма идентификации.
Однако такой подход характеризуется потерей скорости параметрической сходимости при
устранении влияния возмущений и шумов. По мнению авторов настоящей статьи, назрела не-
обходимость совершенствования имеющихся методов с учетом сохранения одновременно
двух показателей качества, т.е. скорости параметрической сходимости и точности оценок.
В данной статье на основе ранних работ Первозванского [18] предлагается новый под-
ход к идентификации частоты синусоидального сигнала с использованием нелинейного
фильтра, что позволит решить задачу сохранения качества.
Постановка задачи. Рассмотрим, как и в работе [15], измеряемый сигнал
y(t) = y(t) + δ(t) ;
(1)
y(t) = σ sin(ωt + ϕ) ,
(2)
где амплитуда σ , частота ω и фаза ϕ — неизвестные постоянные величины, а ограниченный
гладкий сигнал δ(t) характеризует неучтенное возмущающее воздействие или шум в канале
измерения полезного сигнала y(t) .
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
8 А. А. Бобцов, А. В. Крылов, А. А. Пыркин
Ставится задача синтеза схемы идентификации частоты ω, обеспечивающей сходимость оценки частоты ω к ее истинному значению ω с минимальной ошибкой. Иными словами, необходимо обеспечить следующее целевое условие:
lim ω − ω(t) ≤ δ при δ(t) ≠ 0 ,
t→∞
(3)
где число δ = δ(δ(t)) .
Схема идентификации частоты синусоидального сигнала с использованием нелинейного фильтра. Для синтеза схемы идентификации частоты ω сначала рассмотрим сигнал (1) при отсутствии возмущения δ(t) , т.е.
y = y = σsin(ωt + ϕ) .
(4)
Кратко повторим основные положения работы [15]. Известно [13], что для моделирования сигнала (4) можно использовать следующие дифференциальные уравнения:
x1 = x2 ;
(5)
x2 = −ω2 x1; y = k1x1 + k2 x2 = k1x1 + k2 x1, где k1 и k2 — строго положительные постоянные коэффициенты.
(6) (7)
Пусть искомый параметр θ = −ω2 , тогда идеальный алгоритм идентификации параметра θ может быть представлен следующим образом:
θ = kx12 (θ − θ) ,
(8)
где θ — текущая оценка параметра θ . Алгоритм вида (8) основан на классическом методе настройки неизвестных параметров
и при выполнении условий предельной интегральной невырожденности гарантирует асимп-
тотическую сходимость θ к θ (см., например, [7, 13]). Однако схема идентификации (8) тех-
нически нереализуема, так как содержит неизмеряемый сигнал x1 и неизвестный параметр θ .
Для того чтобы получить реализуемый алгоритм, проведем следующие преобразования. Из уравнения (7) имеем
x1 = k2−1(−k1x1 + y) = −ax1 + by ,
(9)
где a = k1 / k2 , b = k2−1 . Тогда, пренебрегая экспоненциально затухающими членами, вызванными ненулевыми
начальными условиями x1(0) , для восстановления переменной x1 будем использовать урав-
нение (9) при x1(0) = 0 . Для того чтобы компенсировать неопределенность x12θ в выражении (8), воспользуемся уравнением (6), при этом алгоритм (8) принимает вид
θ = −kx12θ + kx1x2 . Из выражений (5) и (9) получаем уравнение
(10)
подставляя которое в (10), имеем
x2 = −ax1 + by ,
(11)
θ = −kx12θ + kx1(−ax1 + by) ,
(12)
где функции x1 и x1 определяются из уравнения (9). Очевидно, что алгоритм (12) по-прежнему нереализуем, так как содержит неизвестную
функцию y . Для компенсации этой неизвестной составляющей введем в рассмотрение новую
переменную
ς = θ − kbx1 y .
(13)
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
Повышение точности оценки частоты синусоидального сигнала
9
Дифференцируя (13), получаем реализуемый алгоритм идентификации параметра θ :
ς = θ − kbx1 y − kbx1 y = −kx12θ − kax1 x1 − kbx1 y ,
(14)
θ = ς + kbx1 y , ω = θ .
(15)
Итак, уравнения (9), (14) и (15) представляют собой схему идентификации частоты синусоидального сигнала вида (4), где в случае ненулевого неучтенного возмущения δ(t) в
канале измерения полезного сигнала вместо y(t) используется y(t) = y(t) + δ(t) . Заметим, что
предложенная схема идентификации является робастной относительно неучтенных возмущений. При наличии внешнего возмущения δ(t) в канале измерений алгоритм оценивания час-
тоты позволяет обеспечить выполнение целевого условия (3), где величина δ зависит от δ(t) .
Следует отметить, что уравнение (9) представляет собой низкочастотный фильтр первого порядка: это, в свою очередь, позволяет подавлять возмущение δ(t) за счет выбора коэффици-
ентов a и b . С другой стороны, для повышения точности оценивания частоты ω необходимо уменьшать параметры a и b , что влечет за собой увеличение времени оценивания частоты.
В целях повышения точности и быстродействия оценивания рассмотрим нелинейный
фильтр, структурная схема которого представлена на рис. 1.
Выходная переменная фильтра ω f определяется соотношением
t
ω f = ∫ ϑ(τ)d τ , 0
где сигнал ϑ(τ) является выходом нелинейного звена „насыщение“:
(16)
( )ωf
⎧c
=
⎪ ⎨
d
ω−ωf
, ω−ωf
≤ d;
⎪⎩c sgn (ω − ω f ), ω − ω f > d,
(17)
где c и d — настроечные параметры.
На вход фильтра поступает сигнал оценки частоты ω, рассчитанной по закону (15).
Суть фильтра заключается в том, чтобы удалить из сигнала ω высокочастотную составляю-
щую, связанную с возмущением δ(t) . В работе [18] показано, что пропускная способность
нелинейного звена „насыщение“ понижается с ростом дисперсии входного сигнала: действи-
тельно, чем больше амплитуда на входе, тем меньше коэффициент передачи звена.
ω
c –d
1 ωf
ds
–c
Рис. 1
Как видно из рис. 1, нелинейное звено (17) в фильтре ограничивает скорость роста переменной ω f благодаря константам c и d . Настройка фильтра сводится к определению констант c и d , которые, в свою очередь, определяют максимальную скорость изменения градиента полезного сигнала. Исходя из этого следует, что при настройке фильтра (16), (17) должны учитываться динамические свойства входного процесса.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
10 А. А. Бобцов, А. В. Крылов, А. А. Пыркин
Если входной сигнал содержит компонент с заведомо большей скоростью, чем полез-
ный сигнал, то этот компонент будет подавлен фильтром (16), (17). Если в возмущении δ(t)
присутствует большой по амплитуде выброс, вызванный, например, сбоем датчика, то нели-
нейный фильтр (16), (17) в силу свой структуры не пропустит его, в отличие от любого ли-
нейного фильтра.
Результат работы алгоритма идентификации частоты (формулы (9), (14) и (15)) без ис-
пользования фильтра (16), (17) зависит только от выбора параметров k , a и b [15]. С ростом
значений этих параметров увеличится скорость сходимости оценки частоты в установившую-
ся область, но размер этой области, как и дисперсия сигнала, будет увеличиваться. Наличие
фильтра (16), (17) обеспечивает существенное повышение точности оценивания при сохране-
нии быстродействия алгоритма идентификации частоты.
Для иллюстрации работоспособности предложенной схемы идентификации рассмотрим
пример.
Пример. Сравним работу схемы идентификации, предложенной в настоящей статье, со
схемой, рассмотренной в работе [15]. На рис. 2 и 3 приведены результаты моделирования
схем идентификации (9), (14)—(17).
а)
y (t) ,
о.е. 10
Измеряемый сигнал y (t)
б)
ω, рад/с 1,3
Оценка частоты ω=1,3 рад/с
51
0ω ωf
–5 0,5
–10 –15
0
а) y (t) , о.е.
5
5 10 15 Измеряемый сигнал y (t)
t, с Рис. 2
б)
0
ω, рад/с 1,2
1
5 10 15 t, с Оценка частоты ω=1 рад/с
0,8 0
0,6
–5 0,4 ω1 0,2 ω2
–10 ω f
0 5 10 15 t, с
0
5 10 15 t, с
Рис. 3
Возмущение δ(t) типа белый шум моделировалось как последовательность случайных
чисел с нормальным распределением, сменяющихся на каждом интервале времени t0 . На рис. 2, а представлена временная диаграмма для измеряемого сигнала y(t) , где
y(t) = 8sin(1,3t) . При моделировании были выбраны следующие параметры сигнала δ(t) :
мощность N = 0, 001, интервал квантования t0 = 0, 001 c . На рис. 2, б представлены результаты
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
Повышение точности оценки частоты синусоидального сигнала
11
оценивания частоты при следующих параметрах алгоритма идентификации: a = 0,5 , b = 0,3 ,
k = 0, 2 , c = 0, 3 , d = 0, 3 . Оценка частоты ω, полученная по алгоритму (9), (14), (15), соответ-
ствует результату, приведенному в работе [15]. При использовании дополнительного фильтрующего устройства (16) получена оценка ω f . Очевидно, что оценка ω f точнее ω при том
же времени переходного процесса. В статье [15] рассматривается возможность компенсации возмущения δ(t) за счет вы-
бора коэффициентов a и b без использования дополнительных фильтрующих устройств. Рис. 3 иллюстрирует, что нелинейный фильтр (16), (17) позволяет получить более точную оценку за меньшее время. На рис. 3, а представлена временная диаграмма для измеряемого сигнала y(t) , где y(t) = 8sin(1t + 1) . В этом случае при моделировании были выбраны сле-
дующие параметры сигнала δ(t) : мощность N = 0, 05 , интервал квантования t0 = 0, 05 c . На рис. 3, б представлены результаты оценивания частоты при различных параметрах схемы идентификации: оценки ω1(t) и ω2 (t) получены на основе алгоритма (9), (14), (15) при
a1 = 0, 2 , b1 = 0, 2 , k1 = 0, 7 и a2 = 0,1, b2 = 0,1, k2 = 0, 2 соответственно; оценка ω f (t) полу-
чена на основе сигнала ω1(t) при c = 1 , d = 1 . Здесь, как и на рис. 2, видно, что использование нелинейного фильтра значительно повышает точность оценивания частоты. Для обеспечения такой же точности оценивания без использования фильтра (16), (17), следуя работе [15], были выбраны другие параметры схемы идентификации. Сравнение времени переходного процесса для ω2 (t) и ω f (t) показывает очевидное преимущество метода оценивания час-
тоты с использованием нелинейного фильтра (16), (17). Заключение. Предложен новый метод оценивания частоты сигнала, содержащего по-
лезную составляющую в форме синусоидальной функции времени, и возмущающего сигнала, вызванного помехами в канале измерения. Достоинством представленной схемы является обеспечение желаемой точности идентификации частоты за меньшее время в сравнении с известными мировыми аналогами.
Работа выполнена при поддержке Российского фонда фундаментальных исследований, грант № 09-08-00139-а.
СПИСОК ЛИТЕРАТУРЫ
1. Bodson M., Douglas S. C. Adaptive algorithms for the rejection of periodic disturbances with unknown frequencies // Automatica. 1997. Vol. 33. P. 2213—2221.
2. Hsu L., Ortega R., Damm G. A globally convergent frequency estimator // IEEE Transact. on Automatic Control. 1999. Vol. 46. P. 967—972.
3. Mojiri M., Bakhshai A. R. An adaptive notch filter for frequency estimation of a periodic signal // IEEE Transact. on Automatic Control. 2004. Vol. 49. P. 314—318.
4. Marino R., Tomei R. Global estimation of unknown frequencies // IEEE Transact. on Automatic Control. 2002. Vol. 47. P. 1324—1328.
5. Xia X. Global frequency estimation using adaptive identifiers // IEEE Transact. on Automatic Control. 2002. Vol. 47. P. 1188—1193.
6. Obregon-Pulido G., Castillo-Toledo B., Loukianov A. A. Globally convergent estimator for n-frequencies // IEEE Transact. on Automatic Control. 2002. Vol. 47. P. 857—863.
7. Bobtsov A., Lyamin A., Romasheva D. Algorithm of parameter’s identification of polyharmonic function // 15th IFAC World Congress on Automatic Control. Barcelona, Spain, 2002.
8. Бобцов А. А., Кремлев А. С. Адаптивная идентификация частоты смещенного синусоидального сигнала // Изв. вузов. Приборостроение. 2005. Т. 48, № 4. С. 22—26.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
12 А. А. Бобцов, А. В. Крылов, А. А. Пыркин
9. Hou M. Amplitude and frequency estimator of a sinusoid // IEEE Transact. on Automatic Control. 2005. Vol. 50. P. 855—858.
10. Арановский С. В., Бобцов А. А., Кремлев А. С., Лукьянова Г. В. Робастный алгоритм идентификации частоты синусоидального сигнала // Изв. РАН. Сер. Теория и системы управления. 2007. № 3. С. 1—6.
11. Арановский С. В., Бобцов А. А., Кремлев А. С. и др. Идентификация частоты смещенного синусоидального сигнала // Автоматика и телемеханика. 2008. № 9. С. 3—9.
12. Aranovskiy S., Bobtsov A., Kremlev A. et al. Identification of frequency of biased harmonic signal // IFAC Workshop on Adaptation and Learning in Control and Signal Processing (ALCOSP 07). St. Petersburg, 2007.
13. Bobtsov A. New approach to the problem of globally convergent frequency estimator // Intern. Journal of Adaptive Control and Signal Processing. 2008. N 3. P. 306—317.
14. Aranovskiy S., Bobtsov A., Kremlev A. et al. Identification of frequency of biased harmonic signal // European Journal of Control. 2010. N 2.
15. Бобцов А. А., Николаев Н. А., Слита О. В. Новая схема идентификации частоты синусоидального сигнала // Мехатроника, автоматизация, управление. 2010. № 11. С. 2—4.
16. Пыркин А. А. Адаптивный алгоритм компенсации параметрически неопределенного смещенного гармонического возмущения для линейного объекта с запаздыванием в канале управления // Автоматика и телемеханика. 2010. № 8. С. 62—78.
17. Бобцов А. А., Колюбин С. А., Пыркин А. А. Компенсация неизвестного мультигармонического возмущения для нелинейного объекта с запаздыванием по управлению // Там же. 2010. № 11. С. 136—148.
18. Первозванский А. А. Случайные процессы в нелинейных автоматических системах. М.: Физматгиз, 1962. 352 с.
Алексей Алексеевич Бобцов
Андрей Валентинович Крылов Антон Александрович Пыркин
Сведения об авторах — д-р техн. наук, профессор; Санкт-Петербургский государственный
университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: bobtsov@mail.ifmo.ru — канд. техн. наук; ЗАО „НАВИС“, Санкт-Петербург; начальник отдела; E-mail: a.krylov@navisincontrol.com — канд. техн. наук; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: a.pyrkin@gmail.com
Рекомендована кафедрой систем управления и информатики СПбГУ ИТМО
Поступила в редакцию 18.01.11 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
УДК 681.5.03
А. А. БОБЦОВ, А. В. КРЫЛОВ, А. А. ПЫРКИН
ПОВЫШЕНИЕ ТОЧНОСТИ ОЦЕНКИ ЧАСТОТЫ СИНУСОИДАЛЬНОГО СИГНАЛА С ИСПОЛЬЗОВАНИЕМ НЕЛИНЕЙНОГО ФИЛЬТРА
Рассматривается новая схема идентификации частоты измеряемого синусоидального сигнала с использованием нелинейного фильтра. В отличие от большинства известных аналогов, предлагаемый алгоритм позволяет получить более точные оценки за меньшее время.
Ключевые слова: идентификация, гармонический сигнал, фильтрация.
Введение. Актуальность проблемы идентификации частоты измеряемого синусоидаль-
ного сигнала в условиях влияния внешних возмущений и шумов подтверждается наличием
большого числа публикаций, посвященных идентификации частоты синусоидального или
частот мультигармонического сигналов [1—17]. Однако в большинстве указанных работ не
обсуждается задача качества идентификации, а именно ускорения вычислительной процеду-
ры и повышения точности оценок при неучтенных возмущениях, присутствующих в канале
измерения полезного сигнала. Анализ и рекомендации по увеличению скорости сходимости
оценки частоты синусоидального сигнала к истинному значению были опубликованы в рабо-
тах [11, 14—16], а их робастность к неучтенным возмущениям и помехам исследована в ра-
ботах [13—16]. В частности, в работе [15] была предложена новая схема идентификации час-
тоты измеряемого синусоидального сигнала в условиях неучтенных возмущающих воздейст-
вий и шумов, присутствующих в канале измерения полезного сигнала. Одно из отличий схе-
мы, предложенной в работе [15], от известных аналогов заключается в возможности компен-
сации неучтенных возмущений путем подбора коэффициентов алгоритма идентификации.
Однако такой подход характеризуется потерей скорости параметрической сходимости при
устранении влияния возмущений и шумов. По мнению авторов настоящей статьи, назрела не-
обходимость совершенствования имеющихся методов с учетом сохранения одновременно
двух показателей качества, т.е. скорости параметрической сходимости и точности оценок.
В данной статье на основе ранних работ Первозванского [18] предлагается новый под-
ход к идентификации частоты синусоидального сигнала с использованием нелинейного
фильтра, что позволит решить задачу сохранения качества.
Постановка задачи. Рассмотрим, как и в работе [15], измеряемый сигнал
y(t) = y(t) + δ(t) ;
(1)
y(t) = σ sin(ωt + ϕ) ,
(2)
где амплитуда σ , частота ω и фаза ϕ — неизвестные постоянные величины, а ограниченный
гладкий сигнал δ(t) характеризует неучтенное возмущающее воздействие или шум в канале
измерения полезного сигнала y(t) .
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
8 А. А. Бобцов, А. В. Крылов, А. А. Пыркин
Ставится задача синтеза схемы идентификации частоты ω, обеспечивающей сходимость оценки частоты ω к ее истинному значению ω с минимальной ошибкой. Иными словами, необходимо обеспечить следующее целевое условие:
lim ω − ω(t) ≤ δ при δ(t) ≠ 0 ,
t→∞
(3)
где число δ = δ(δ(t)) .
Схема идентификации частоты синусоидального сигнала с использованием нелинейного фильтра. Для синтеза схемы идентификации частоты ω сначала рассмотрим сигнал (1) при отсутствии возмущения δ(t) , т.е.
y = y = σsin(ωt + ϕ) .
(4)
Кратко повторим основные положения работы [15]. Известно [13], что для моделирования сигнала (4) можно использовать следующие дифференциальные уравнения:
x1 = x2 ;
(5)
x2 = −ω2 x1; y = k1x1 + k2 x2 = k1x1 + k2 x1, где k1 и k2 — строго положительные постоянные коэффициенты.
(6) (7)
Пусть искомый параметр θ = −ω2 , тогда идеальный алгоритм идентификации параметра θ может быть представлен следующим образом:
θ = kx12 (θ − θ) ,
(8)
где θ — текущая оценка параметра θ . Алгоритм вида (8) основан на классическом методе настройки неизвестных параметров
и при выполнении условий предельной интегральной невырожденности гарантирует асимп-
тотическую сходимость θ к θ (см., например, [7, 13]). Однако схема идентификации (8) тех-
нически нереализуема, так как содержит неизмеряемый сигнал x1 и неизвестный параметр θ .
Для того чтобы получить реализуемый алгоритм, проведем следующие преобразования. Из уравнения (7) имеем
x1 = k2−1(−k1x1 + y) = −ax1 + by ,
(9)
где a = k1 / k2 , b = k2−1 . Тогда, пренебрегая экспоненциально затухающими членами, вызванными ненулевыми
начальными условиями x1(0) , для восстановления переменной x1 будем использовать урав-
нение (9) при x1(0) = 0 . Для того чтобы компенсировать неопределенность x12θ в выражении (8), воспользуемся уравнением (6), при этом алгоритм (8) принимает вид
θ = −kx12θ + kx1x2 . Из выражений (5) и (9) получаем уравнение
(10)
подставляя которое в (10), имеем
x2 = −ax1 + by ,
(11)
θ = −kx12θ + kx1(−ax1 + by) ,
(12)
где функции x1 и x1 определяются из уравнения (9). Очевидно, что алгоритм (12) по-прежнему нереализуем, так как содержит неизвестную
функцию y . Для компенсации этой неизвестной составляющей введем в рассмотрение новую
переменную
ς = θ − kbx1 y .
(13)
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
Повышение точности оценки частоты синусоидального сигнала
9
Дифференцируя (13), получаем реализуемый алгоритм идентификации параметра θ :
ς = θ − kbx1 y − kbx1 y = −kx12θ − kax1 x1 − kbx1 y ,
(14)
θ = ς + kbx1 y , ω = θ .
(15)
Итак, уравнения (9), (14) и (15) представляют собой схему идентификации частоты синусоидального сигнала вида (4), где в случае ненулевого неучтенного возмущения δ(t) в
канале измерения полезного сигнала вместо y(t) используется y(t) = y(t) + δ(t) . Заметим, что
предложенная схема идентификации является робастной относительно неучтенных возмущений. При наличии внешнего возмущения δ(t) в канале измерений алгоритм оценивания час-
тоты позволяет обеспечить выполнение целевого условия (3), где величина δ зависит от δ(t) .
Следует отметить, что уравнение (9) представляет собой низкочастотный фильтр первого порядка: это, в свою очередь, позволяет подавлять возмущение δ(t) за счет выбора коэффици-
ентов a и b . С другой стороны, для повышения точности оценивания частоты ω необходимо уменьшать параметры a и b , что влечет за собой увеличение времени оценивания частоты.
В целях повышения точности и быстродействия оценивания рассмотрим нелинейный
фильтр, структурная схема которого представлена на рис. 1.
Выходная переменная фильтра ω f определяется соотношением
t
ω f = ∫ ϑ(τ)d τ , 0
где сигнал ϑ(τ) является выходом нелинейного звена „насыщение“:
(16)
( )ωf
⎧c
=
⎪ ⎨
d
ω−ωf
, ω−ωf
≤ d;
⎪⎩c sgn (ω − ω f ), ω − ω f > d,
(17)
где c и d — настроечные параметры.
На вход фильтра поступает сигнал оценки частоты ω, рассчитанной по закону (15).
Суть фильтра заключается в том, чтобы удалить из сигнала ω высокочастотную составляю-
щую, связанную с возмущением δ(t) . В работе [18] показано, что пропускная способность
нелинейного звена „насыщение“ понижается с ростом дисперсии входного сигнала: действи-
тельно, чем больше амплитуда на входе, тем меньше коэффициент передачи звена.
ω
c –d
1 ωf
ds
–c
Рис. 1
Как видно из рис. 1, нелинейное звено (17) в фильтре ограничивает скорость роста переменной ω f благодаря константам c и d . Настройка фильтра сводится к определению констант c и d , которые, в свою очередь, определяют максимальную скорость изменения градиента полезного сигнала. Исходя из этого следует, что при настройке фильтра (16), (17) должны учитываться динамические свойства входного процесса.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
10 А. А. Бобцов, А. В. Крылов, А. А. Пыркин
Если входной сигнал содержит компонент с заведомо большей скоростью, чем полез-
ный сигнал, то этот компонент будет подавлен фильтром (16), (17). Если в возмущении δ(t)
присутствует большой по амплитуде выброс, вызванный, например, сбоем датчика, то нели-
нейный фильтр (16), (17) в силу свой структуры не пропустит его, в отличие от любого ли-
нейного фильтра.
Результат работы алгоритма идентификации частоты (формулы (9), (14) и (15)) без ис-
пользования фильтра (16), (17) зависит только от выбора параметров k , a и b [15]. С ростом
значений этих параметров увеличится скорость сходимости оценки частоты в установившую-
ся область, но размер этой области, как и дисперсия сигнала, будет увеличиваться. Наличие
фильтра (16), (17) обеспечивает существенное повышение точности оценивания при сохране-
нии быстродействия алгоритма идентификации частоты.
Для иллюстрации работоспособности предложенной схемы идентификации рассмотрим
пример.
Пример. Сравним работу схемы идентификации, предложенной в настоящей статье, со
схемой, рассмотренной в работе [15]. На рис. 2 и 3 приведены результаты моделирования
схем идентификации (9), (14)—(17).
а)
y (t) ,
о.е. 10
Измеряемый сигнал y (t)
б)
ω, рад/с 1,3
Оценка частоты ω=1,3 рад/с
51
0ω ωf
–5 0,5
–10 –15
0
а) y (t) , о.е.
5
5 10 15 Измеряемый сигнал y (t)
t, с Рис. 2
б)
0
ω, рад/с 1,2
1
5 10 15 t, с Оценка частоты ω=1 рад/с
0,8 0
0,6
–5 0,4 ω1 0,2 ω2
–10 ω f
0 5 10 15 t, с
0
5 10 15 t, с
Рис. 3
Возмущение δ(t) типа белый шум моделировалось как последовательность случайных
чисел с нормальным распределением, сменяющихся на каждом интервале времени t0 . На рис. 2, а представлена временная диаграмма для измеряемого сигнала y(t) , где
y(t) = 8sin(1,3t) . При моделировании были выбраны следующие параметры сигнала δ(t) :
мощность N = 0, 001, интервал квантования t0 = 0, 001 c . На рис. 2, б представлены результаты
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
Повышение точности оценки частоты синусоидального сигнала
11
оценивания частоты при следующих параметрах алгоритма идентификации: a = 0,5 , b = 0,3 ,
k = 0, 2 , c = 0, 3 , d = 0, 3 . Оценка частоты ω, полученная по алгоритму (9), (14), (15), соответ-
ствует результату, приведенному в работе [15]. При использовании дополнительного фильтрующего устройства (16) получена оценка ω f . Очевидно, что оценка ω f точнее ω при том
же времени переходного процесса. В статье [15] рассматривается возможность компенсации возмущения δ(t) за счет вы-
бора коэффициентов a и b без использования дополнительных фильтрующих устройств. Рис. 3 иллюстрирует, что нелинейный фильтр (16), (17) позволяет получить более точную оценку за меньшее время. На рис. 3, а представлена временная диаграмма для измеряемого сигнала y(t) , где y(t) = 8sin(1t + 1) . В этом случае при моделировании были выбраны сле-
дующие параметры сигнала δ(t) : мощность N = 0, 05 , интервал квантования t0 = 0, 05 c . На рис. 3, б представлены результаты оценивания частоты при различных параметрах схемы идентификации: оценки ω1(t) и ω2 (t) получены на основе алгоритма (9), (14), (15) при
a1 = 0, 2 , b1 = 0, 2 , k1 = 0, 7 и a2 = 0,1, b2 = 0,1, k2 = 0, 2 соответственно; оценка ω f (t) полу-
чена на основе сигнала ω1(t) при c = 1 , d = 1 . Здесь, как и на рис. 2, видно, что использование нелинейного фильтра значительно повышает точность оценивания частоты. Для обеспечения такой же точности оценивания без использования фильтра (16), (17), следуя работе [15], были выбраны другие параметры схемы идентификации. Сравнение времени переходного процесса для ω2 (t) и ω f (t) показывает очевидное преимущество метода оценивания час-
тоты с использованием нелинейного фильтра (16), (17). Заключение. Предложен новый метод оценивания частоты сигнала, содержащего по-
лезную составляющую в форме синусоидальной функции времени, и возмущающего сигнала, вызванного помехами в канале измерения. Достоинством представленной схемы является обеспечение желаемой точности идентификации частоты за меньшее время в сравнении с известными мировыми аналогами.
Работа выполнена при поддержке Российского фонда фундаментальных исследований, грант № 09-08-00139-а.
СПИСОК ЛИТЕРАТУРЫ
1. Bodson M., Douglas S. C. Adaptive algorithms for the rejection of periodic disturbances with unknown frequencies // Automatica. 1997. Vol. 33. P. 2213—2221.
2. Hsu L., Ortega R., Damm G. A globally convergent frequency estimator // IEEE Transact. on Automatic Control. 1999. Vol. 46. P. 967—972.
3. Mojiri M., Bakhshai A. R. An adaptive notch filter for frequency estimation of a periodic signal // IEEE Transact. on Automatic Control. 2004. Vol. 49. P. 314—318.
4. Marino R., Tomei R. Global estimation of unknown frequencies // IEEE Transact. on Automatic Control. 2002. Vol. 47. P. 1324—1328.
5. Xia X. Global frequency estimation using adaptive identifiers // IEEE Transact. on Automatic Control. 2002. Vol. 47. P. 1188—1193.
6. Obregon-Pulido G., Castillo-Toledo B., Loukianov A. A. Globally convergent estimator for n-frequencies // IEEE Transact. on Automatic Control. 2002. Vol. 47. P. 857—863.
7. Bobtsov A., Lyamin A., Romasheva D. Algorithm of parameter’s identification of polyharmonic function // 15th IFAC World Congress on Automatic Control. Barcelona, Spain, 2002.
8. Бобцов А. А., Кремлев А. С. Адаптивная идентификация частоты смещенного синусоидального сигнала // Изв. вузов. Приборостроение. 2005. Т. 48, № 4. С. 22—26.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6
12 А. А. Бобцов, А. В. Крылов, А. А. Пыркин
9. Hou M. Amplitude and frequency estimator of a sinusoid // IEEE Transact. on Automatic Control. 2005. Vol. 50. P. 855—858.
10. Арановский С. В., Бобцов А. А., Кремлев А. С., Лукьянова Г. В. Робастный алгоритм идентификации частоты синусоидального сигнала // Изв. РАН. Сер. Теория и системы управления. 2007. № 3. С. 1—6.
11. Арановский С. В., Бобцов А. А., Кремлев А. С. и др. Идентификация частоты смещенного синусоидального сигнала // Автоматика и телемеханика. 2008. № 9. С. 3—9.
12. Aranovskiy S., Bobtsov A., Kremlev A. et al. Identification of frequency of biased harmonic signal // IFAC Workshop on Adaptation and Learning in Control and Signal Processing (ALCOSP 07). St. Petersburg, 2007.
13. Bobtsov A. New approach to the problem of globally convergent frequency estimator // Intern. Journal of Adaptive Control and Signal Processing. 2008. N 3. P. 306—317.
14. Aranovskiy S., Bobtsov A., Kremlev A. et al. Identification of frequency of biased harmonic signal // European Journal of Control. 2010. N 2.
15. Бобцов А. А., Николаев Н. А., Слита О. В. Новая схема идентификации частоты синусоидального сигнала // Мехатроника, автоматизация, управление. 2010. № 11. С. 2—4.
16. Пыркин А. А. Адаптивный алгоритм компенсации параметрически неопределенного смещенного гармонического возмущения для линейного объекта с запаздыванием в канале управления // Автоматика и телемеханика. 2010. № 8. С. 62—78.
17. Бобцов А. А., Колюбин С. А., Пыркин А. А. Компенсация неизвестного мультигармонического возмущения для нелинейного объекта с запаздыванием по управлению // Там же. 2010. № 11. С. 136—148.
18. Первозванский А. А. Случайные процессы в нелинейных автоматических системах. М.: Физматгиз, 1962. 352 с.
Алексей Алексеевич Бобцов
Андрей Валентинович Крылов Антон Александрович Пыркин
Сведения об авторах — д-р техн. наук, профессор; Санкт-Петербургский государственный
университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: bobtsov@mail.ifmo.ru — канд. техн. наук; ЗАО „НАВИС“, Санкт-Петербург; начальник отдела; E-mail: a.krylov@navisincontrol.com — канд. техн. наук; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра систем управления и информатики; E-mail: a.pyrkin@gmail.com
Рекомендована кафедрой систем управления и информатики СПбГУ ИТМО
Поступила в редакцию 18.01.11 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 6