ИДЕНТИФИКАЦИЯ ЛИНЕЙНО МЕНЯЮЩЕЙСЯ ЧАСТОТЫ СИНУСОИДАЛЬНОГО СИГНАЛА
ИДЕНТИФИКАЦИЯ ЛИНЕЙНО МЕНЯЮЩЕЙСЯ ЧАСТОТЫ …
3 АНАЛИЗ И СИНТЕЗ СЛОЖНЫХ СИСТЕМ
УДК 681.51.015
ИДЕНТИФИКАЦИЯ ЛИНЕЙНО МЕНЯЮЩЕЙСЯ ЧАСТОТЫ СИНУСОИДАЛЬНОГО СИГНАЛА
С.В. Арановский, А.А. Бобцов, А.А. Пыркин
Рассмотрена задача идентификации линейно меняющейся частоты синусоидального сигнала. Показана практическая значимость рассматриваемой задачи, получена динамическая модель сигнала и предложен алгоритм идентификации линейно меняющейся частоты. Приведен численный пример, иллюстрирующий работоспособность предложенного алгоритма. Ключевые слова: идентификация, линейно меняющаяся частота, синусоидальный сигнал.
Введение
В работе рассматривается задача идентификации линейно меняющейся частоты синусоидального сигнала вида y(t) Asin(t2 2 ) . Несмотря на кажущуюся математическую экзотичность задачи (тра-
диционно строится алгоритм идентификации для сигналов Asin(t ) ), необходимо отметить, что проблема идентификации линейно меняющейся частоты синусоидального сигнала имеет реальное прикладное значение и может возникать, например, в задачах компенсации возмущений. Проблема идентификации постоянной частоты гармонического сигнала была хорошо изучена за последнее десятилетие, разработано большое число функционирующих в реальном времени алгоритмов [1–8]. В то же время задача идентификации переменной частоты синусоидального сигнала в литературе практически не рассматривается. На данный момент авторам не известны методы, позволяющие в реальном времени оценивать частоту такого сигнала. Существующие методы, в частности, алгоритм [8], позволяют оценивать частоту гармонического сигнала, имеющего трапециевидную форму (рис. 1), но не обеспечивают нулевую ошибку оценивания.
Рис. 1. Оценка частоты ˆ при изменении частоты методом [8]
В работе ставится цель разработать алгоритм идентификации линейно меняющейся частоты гармонического сигнала.
Возмущения с линейно меняющейся частотой при управлении телескопом
Рассмотрим, как связана обсуждаемая в работе задача с прикладной проблемой компенсации возмущений при управлении телескопом. Широкое распространение в телескопостроении получили вентильные двигатели [9, 10]. Такие машины имеют много достоинств, но, к сожалению, не лишены недостатков. Одним из недостатков является наличие пульсирующего зубцового момента, который для трехфазных машин может быть описан как
6 M p () M a sin( pn 2 p ) ,
(1)
где – угол поворота двигателя; M p () – значение зубцового момента в угле ; M a – амплитуда
пульсаций момента; pn – число пар полюсов; p – фаза зубцового момента. Очевидно, что если ось
28 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2012, № 1 (77)
С.В. Арановский, А.А. Бобцов, А.А. Пыркин
телескопа вращается с постоянной скоростью const , то во временной области зубцовый момент опи-
сывается гармонической функцией M p (t) Ma sin(pt p ) с частотой
p
pn
6 2
const
.
Проблема компенсации зубцового момента в приводах телескопов не нова. Она решалась, напри-
мер, при построении систем управления группировкой телескопов VLT, Чили [11]. Однако при этом рас-
сматривалась задача компенсации пульсаций момента только для той оси, приводами которой они вы-
званы. Естественно, что при этом значение угла доступно прямому измерению, как и значение , и
задача компенсации сводится к определению только амплитуды и фазы сигнала при известной частоте.
Традиционно оптические телескопы включают в себя две оси – азимутальную и ось угла места,
при этом в фотоприемном устройстве используются светочувствительные матрицы. Однако в ряде ситу-
аций оправдано использование вместо матрицы линейки светочувствительных элементов, размещенной
на третьей оси телескопа – оси сканирования. В этом случае ось сканирования осуществляет периодиче-
ские колебательные движения, а поступающие со светочувствительной линейки столбцы формируют
полный кадр изображения. Недостатком такой конструкции является то, что, в отличие от азимутальной
и угломестной осей, которые всегда ортогональны, угол между осью сканирования и азимутальной мо-
жет меняться от 0º до 90º и совпадает с углом наклона оси угла места к горизонту. Соответственно мо-
менты, создаваемые приводами одной оси, оказывают возмущающее воздействие на вторую ось с коэф-
фициентом пропорциональности, равным косинусу угла между осями. Если угломестная ось ориентиро-
вана в зенит, то плоскость сканирования перпендикулярна азимутальной плоскости, и взаимное влияние
осей отсутствует. Если же ось угла места ориентирована в горизонт, то взаимное влияние осей макси-
мально, момент, создаваемый приводами азимутальной оси, полностью передается на ось сканирования,
и наоборот.
Рассмотрим задачу компенсации пульсаций зубцового момента приводов азимутальной оси, пере-
даваемых на сканирующую ось. Первый возможный подход заключается в передаче информации об угле
поворота одной оси в систему управления другой осью. В этом случае задача компенсации становится
тривиальной, но расплачиваться за это приходится увеличением информационного потока в каналах свя-
зи, появлением задержек и снижением максимальной частоты дискретизации систем управления. Второй
путь – идентификация параметров возмущающего момента с его последующей компенсацией. Это по-
требует внедрения алгоритмов идентификации, но снизит объем передаваемой по каналам связи инфор-
мации. Как было отмечено выше, при движении с постоянной скоростью зубцовый момент представляет
собой гармоническую функцию времени с постоянной частотой. В этом случае для идентификации па-
раметров возмущения (частоты, амплитуды и фазы) разработано большое число функционирующих в
реальном времени алгоритмов [1–8], позволяющих решить задачу компенсации возмущения.
Более сложная задача идентификации возникает, если предположить, что траектория движения
оси телескопа состоит из участков движения с постоянной скоростью (сопровождение медленных объек-
тов) и участков равноускоренного движения (переход между различными объектами, сопровождение
быстро движущихся объектов, проекция траекторий которых на азимутальную ось описывается парабо-
лой). В этом случае развитие зубцового момента (1) во времени будет описываться гармонической функ-
цией с линейно меняющейся частотой, и для компенсации такого возмущения потребуется решить зада-
чу идентификации линейно меняющейся частоты.
Математическая постановка задачи
Рассматривается синусоидальный сигнал вида
y(t) Asin(t2 / 2 ) ,
(2)
где A , и – неизвестные параметры. Ставится задача синтеза алгоритма идентификации неизвест-
ной частоты .
Основной результат
Продифференцируем уравнение (2) трижды: y (3ty t2 y) или p3 y (3ty t2 py) ,
(3)
где p d / dt – оператор дифференцирования; 2 – неизвестный параметр. Для оценки неизвестного параметра произведем параметризацию модели (3) следующим образом:
p( p 1)2 y (2 p2 p) y (3ty t2 py) ,
причем t 2 py pt2 y 2ty . Запишем:
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2012, № 1 (77)
29
ИДЕНТИФИКАЦИЯ ЛИНЕЙНО МЕНЯЮЩЕЙСЯ ЧАСТОТЫ …
py
2p2 p ( p 1)2
y
(
p
3 1)2
ty
(
p
p 1)2
t2
y
(
p
2 1)2
ty
.
Введем следующие обозначения
1
1 ( p 1)2
y,
2
1 ( p 1)2
ty
,
3
1 ( p 1)2
t2 y
,
при подстановке которых, выражение (4) примет вид
py (2 p2 p)1 2 3 . Тогда можно записать параметризованную модель в виде
(4)
y 21 1 2 3 21 1 (2 3 ) 1 2 ,
(5)
где 1 21 1 и 2 2 3 – измеряемые сигналы.
Теперь построим алгоритм идентификации параметра . Временно предположим, что параметр известен. Запишем алгоритм идентификации следующим образом:
ˆ k22ˆ k22 ,
(6)
где коэффициент k 0 . Для доказательства работоспособности алгоритма (6) рассмотрим ошибку
ˆ . Дифференцируя , получаем
ˆ ˆ k22 . Рассмотрим функцию Ляпунова вида
V 2 / 2 . Дифференцируя (7), получаем
(7)
V k222 k22V 0 .
(8)
Из (8) следует устойчивость положения равновесия 0 . Интегрируя (8), получаем
V (t) V (t0 )ek(t) ,
t
где
t0
– начальный момент времени;
(t
)
t0
22
d
.
Легко
видеть,
что
(t)
при
t ,
а,
следова-
тельно, lim (t) 0 .
t
Таким образом, алгоритм идентификации (6) обеспечивает сходимость настраиваемого параметра ˆ (t) к . Теперь построим реализуемый алгоритм, в котором параметр не измеряется. Из уравнения
(5) находим
y 1 2 . Подставляя (9) в (6), имеем
(9)
ˆ k22ˆ k2 y k21 . Поскольку переменная y неизвестна, то введем в рассмотрение сигнал
(10)
ˆ k2 y .
(11)
Дифференцируя (11), получаем реализацию алгоритма идентификации вида
k22ˆ k2 y k21 ,
(12)
ˆ k2 y .
(13)
Очевидно, что при дифференцировании (13) будет получено выражение (10). Алгоритм идентифи-
кации (12), (13) обеспечивает сходимость ˆ (t) к . Для поиска параметра достаточно провести сле-
дующую несложную процедуру ˆ ˆ (t) , где ˆ – оценка частоты .
Пример
Рассмотрим сигнал вида (2) с параметрами 0, 5 ; A 2 ; 1 . На рис. 2 представлена времен-
ная диаграмма сигнала y(t) . На рис. 3, а, б, представлены результаты моделирования алгоритма идентификации параметра с коэффициентами k 1 и k 0, 01 соответственно. На рис. 4, а, б, представлены графики переходных процессов оценки ˆ при коэффициентах алгоритма идентификации k 0, 001 и k 0, 2 для сигнала y(t) с параметрами 2 ; A 5 ; 2 .
30 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2012, № 1 (77)
С.В. Арановский, А.А. Бобцов, А.А. Пыркин
Рис. 2. Временная диаграмма сигнала y(t)
ˆ
ˆ
аб Рис. 3. Переходные процессы для переменной ˆ при 0, 5 : k 1(а); k 0, 01(б)
ˆ ˆ
аб
Рис. 4. Переходные процессы для переменной ˆ при 2 : k 0, 001 (а); k 0, 2 (б)
Заключение
Рассмотрена задача идентификации частоты синусоидального сигнала с линейно меняющейся частотой. Получена динамическая модель генератора сигнала вида y(t) Asin(t2 / 2 ) , предложен ал-
горитм идентификации вида (12), (13), обеспечивающий сходимость настраиваемого параметра ˆ к истинному значению . Предложенный алгоритм имеет настроечный коэффициент k, варьируя который можно изменять быстродействие алгоритма.
Работа выполнена при финансовой поддержке ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2007–2013 годы» (контракт № 11.519.11.4007).
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2012, № 1 (77)
31
УСОВЕРШЕНСТВОВАННЫЙ АЛГОРИТМ АНАЛИЗА СИГНАЛОВ …
Литература
1. Арановский С.В., Бобцов А.А., Кремлев А.С., Лукьянова Г.В., Николаев Н.А. Идентификация частоты смещенного синусоидального сигнала // Автоматика и телемеханика. – 2008. – № 9. – C. 3–9.
2. Арановский С.В., Бобцов А.А., Пыркин А.А. Адаптивный наблюдатель неизвестного синусоидального выходного возмущения для линейного объекта // Автоматика и телемеханика. – 2009. – № 11. – C. 108–116.
3. Арановский С.В., Бобцов А.А., Кремлев А.С. Компенсация конечномерного квазигармонического возмущения для нелинейного объекта // Известия РАН. Теория и системы управления. – 2006. – № 4. – C. 14–21.
4. Арановский С.В., Бобцов А.А., Никифоров В.О. Синтез наблюдателя для нелинейного объекта в условиях гармонического возмущения, приложенного к выходной переменной // Научно-технический вестник СПбГУ ИТМО. – 2010. – № 3. – C. 32–39.
5. Marino R., Tomei R. Global Estimation of Unknown Frequencies // IEEE Transactions on Automatic Control. – 2002. – № 47. – P. 1324–1328.
6. Hou M. Amplitude and frequency estimator of a sinusoid // IEEE Transactions on Automatic Control. – 2005. – V. 50. – P. 855–858.
7. Marino R., Santosuosso G.L. Output regulation for a class of nonlinear systems with unknown sinusoidal disturbances // Proceedings of the 44th IEEE Conference on Decision and Control. – 2005. – C. 3321–3326.
8. Stanislav Aranovskiy, Alexey Bobtsov, Artem Kremlev, Nikolay Nikolaev, Olga Slita. Identification of frequency of biased harmonic signal // European Journal of Control. – 2010. – № 4. – Р. 129–139.
9. Васильев В.Н., Томасов В.С., Шаргородский В.Д., Садовников М.А. Состояние и перспективы развития прецизионных электроприводов комплексов высокоточных наблюдений // Изв. вузов. Приборостроение. – 2008. – № 6. – C. 5–12.
10. Аракелян А.К., Афанасьев А.А. Вентильные электрические машины. – М.: Энергоатомиздат, 1997. – 236 с.
11. Erm T., Sandrock S. Adaptive correction of periodic errors improves telescope performance // Proceedings of the American Control Conference. – 2005. – V. 6. – Р. 3776– 3777.
Арановский Станислав Владимирович Бобцов Алексей Алексеевич
Пыркин Антон Александрович
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат технических наук, ст. научный сотрудник, s.aranovskiy@gmail.com
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, доктор технических наук,
профессор, декан, bobtsov@mail.ifmo.ru – Санкт-Петербургский национальный исследовательский университет ин-
формационных технологий, механики и оптики, кандидат технических наук, мл. научный сотрудник, a.pyrkin@gmail.com
32 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2012, № 1 (77)
3 АНАЛИЗ И СИНТЕЗ СЛОЖНЫХ СИСТЕМ
УДК 681.51.015
ИДЕНТИФИКАЦИЯ ЛИНЕЙНО МЕНЯЮЩЕЙСЯ ЧАСТОТЫ СИНУСОИДАЛЬНОГО СИГНАЛА
С.В. Арановский, А.А. Бобцов, А.А. Пыркин
Рассмотрена задача идентификации линейно меняющейся частоты синусоидального сигнала. Показана практическая значимость рассматриваемой задачи, получена динамическая модель сигнала и предложен алгоритм идентификации линейно меняющейся частоты. Приведен численный пример, иллюстрирующий работоспособность предложенного алгоритма. Ключевые слова: идентификация, линейно меняющаяся частота, синусоидальный сигнал.
Введение
В работе рассматривается задача идентификации линейно меняющейся частоты синусоидального сигнала вида y(t) Asin(t2 2 ) . Несмотря на кажущуюся математическую экзотичность задачи (тра-
диционно строится алгоритм идентификации для сигналов Asin(t ) ), необходимо отметить, что проблема идентификации линейно меняющейся частоты синусоидального сигнала имеет реальное прикладное значение и может возникать, например, в задачах компенсации возмущений. Проблема идентификации постоянной частоты гармонического сигнала была хорошо изучена за последнее десятилетие, разработано большое число функционирующих в реальном времени алгоритмов [1–8]. В то же время задача идентификации переменной частоты синусоидального сигнала в литературе практически не рассматривается. На данный момент авторам не известны методы, позволяющие в реальном времени оценивать частоту такого сигнала. Существующие методы, в частности, алгоритм [8], позволяют оценивать частоту гармонического сигнала, имеющего трапециевидную форму (рис. 1), но не обеспечивают нулевую ошибку оценивания.
Рис. 1. Оценка частоты ˆ при изменении частоты методом [8]
В работе ставится цель разработать алгоритм идентификации линейно меняющейся частоты гармонического сигнала.
Возмущения с линейно меняющейся частотой при управлении телескопом
Рассмотрим, как связана обсуждаемая в работе задача с прикладной проблемой компенсации возмущений при управлении телескопом. Широкое распространение в телескопостроении получили вентильные двигатели [9, 10]. Такие машины имеют много достоинств, но, к сожалению, не лишены недостатков. Одним из недостатков является наличие пульсирующего зубцового момента, который для трехфазных машин может быть описан как
6 M p () M a sin( pn 2 p ) ,
(1)
где – угол поворота двигателя; M p () – значение зубцового момента в угле ; M a – амплитуда
пульсаций момента; pn – число пар полюсов; p – фаза зубцового момента. Очевидно, что если ось
28 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2012, № 1 (77)
С.В. Арановский, А.А. Бобцов, А.А. Пыркин
телескопа вращается с постоянной скоростью const , то во временной области зубцовый момент опи-
сывается гармонической функцией M p (t) Ma sin(pt p ) с частотой
p
pn
6 2
const
.
Проблема компенсации зубцового момента в приводах телескопов не нова. Она решалась, напри-
мер, при построении систем управления группировкой телескопов VLT, Чили [11]. Однако при этом рас-
сматривалась задача компенсации пульсаций момента только для той оси, приводами которой они вы-
званы. Естественно, что при этом значение угла доступно прямому измерению, как и значение , и
задача компенсации сводится к определению только амплитуды и фазы сигнала при известной частоте.
Традиционно оптические телескопы включают в себя две оси – азимутальную и ось угла места,
при этом в фотоприемном устройстве используются светочувствительные матрицы. Однако в ряде ситу-
аций оправдано использование вместо матрицы линейки светочувствительных элементов, размещенной
на третьей оси телескопа – оси сканирования. В этом случае ось сканирования осуществляет периодиче-
ские колебательные движения, а поступающие со светочувствительной линейки столбцы формируют
полный кадр изображения. Недостатком такой конструкции является то, что, в отличие от азимутальной
и угломестной осей, которые всегда ортогональны, угол между осью сканирования и азимутальной мо-
жет меняться от 0º до 90º и совпадает с углом наклона оси угла места к горизонту. Соответственно мо-
менты, создаваемые приводами одной оси, оказывают возмущающее воздействие на вторую ось с коэф-
фициентом пропорциональности, равным косинусу угла между осями. Если угломестная ось ориентиро-
вана в зенит, то плоскость сканирования перпендикулярна азимутальной плоскости, и взаимное влияние
осей отсутствует. Если же ось угла места ориентирована в горизонт, то взаимное влияние осей макси-
мально, момент, создаваемый приводами азимутальной оси, полностью передается на ось сканирования,
и наоборот.
Рассмотрим задачу компенсации пульсаций зубцового момента приводов азимутальной оси, пере-
даваемых на сканирующую ось. Первый возможный подход заключается в передаче информации об угле
поворота одной оси в систему управления другой осью. В этом случае задача компенсации становится
тривиальной, но расплачиваться за это приходится увеличением информационного потока в каналах свя-
зи, появлением задержек и снижением максимальной частоты дискретизации систем управления. Второй
путь – идентификация параметров возмущающего момента с его последующей компенсацией. Это по-
требует внедрения алгоритмов идентификации, но снизит объем передаваемой по каналам связи инфор-
мации. Как было отмечено выше, при движении с постоянной скоростью зубцовый момент представляет
собой гармоническую функцию времени с постоянной частотой. В этом случае для идентификации па-
раметров возмущения (частоты, амплитуды и фазы) разработано большое число функционирующих в
реальном времени алгоритмов [1–8], позволяющих решить задачу компенсации возмущения.
Более сложная задача идентификации возникает, если предположить, что траектория движения
оси телескопа состоит из участков движения с постоянной скоростью (сопровождение медленных объек-
тов) и участков равноускоренного движения (переход между различными объектами, сопровождение
быстро движущихся объектов, проекция траекторий которых на азимутальную ось описывается парабо-
лой). В этом случае развитие зубцового момента (1) во времени будет описываться гармонической функ-
цией с линейно меняющейся частотой, и для компенсации такого возмущения потребуется решить зада-
чу идентификации линейно меняющейся частоты.
Математическая постановка задачи
Рассматривается синусоидальный сигнал вида
y(t) Asin(t2 / 2 ) ,
(2)
где A , и – неизвестные параметры. Ставится задача синтеза алгоритма идентификации неизвест-
ной частоты .
Основной результат
Продифференцируем уравнение (2) трижды: y (3ty t2 y) или p3 y (3ty t2 py) ,
(3)
где p d / dt – оператор дифференцирования; 2 – неизвестный параметр. Для оценки неизвестного параметра произведем параметризацию модели (3) следующим образом:
p( p 1)2 y (2 p2 p) y (3ty t2 py) ,
причем t 2 py pt2 y 2ty . Запишем:
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2012, № 1 (77)
29
ИДЕНТИФИКАЦИЯ ЛИНЕЙНО МЕНЯЮЩЕЙСЯ ЧАСТОТЫ …
py
2p2 p ( p 1)2
y
(
p
3 1)2
ty
(
p
p 1)2
t2
y
(
p
2 1)2
ty
.
Введем следующие обозначения
1
1 ( p 1)2
y,
2
1 ( p 1)2
ty
,
3
1 ( p 1)2
t2 y
,
при подстановке которых, выражение (4) примет вид
py (2 p2 p)1 2 3 . Тогда можно записать параметризованную модель в виде
(4)
y 21 1 2 3 21 1 (2 3 ) 1 2 ,
(5)
где 1 21 1 и 2 2 3 – измеряемые сигналы.
Теперь построим алгоритм идентификации параметра . Временно предположим, что параметр известен. Запишем алгоритм идентификации следующим образом:
ˆ k22ˆ k22 ,
(6)
где коэффициент k 0 . Для доказательства работоспособности алгоритма (6) рассмотрим ошибку
ˆ . Дифференцируя , получаем
ˆ ˆ k22 . Рассмотрим функцию Ляпунова вида
V 2 / 2 . Дифференцируя (7), получаем
(7)
V k222 k22V 0 .
(8)
Из (8) следует устойчивость положения равновесия 0 . Интегрируя (8), получаем
V (t) V (t0 )ek(t) ,
t
где
t0
– начальный момент времени;
(t
)
t0
22
d
.
Легко
видеть,
что
(t)
при
t ,
а,
следова-
тельно, lim (t) 0 .
t
Таким образом, алгоритм идентификации (6) обеспечивает сходимость настраиваемого параметра ˆ (t) к . Теперь построим реализуемый алгоритм, в котором параметр не измеряется. Из уравнения
(5) находим
y 1 2 . Подставляя (9) в (6), имеем
(9)
ˆ k22ˆ k2 y k21 . Поскольку переменная y неизвестна, то введем в рассмотрение сигнал
(10)
ˆ k2 y .
(11)
Дифференцируя (11), получаем реализацию алгоритма идентификации вида
k22ˆ k2 y k21 ,
(12)
ˆ k2 y .
(13)
Очевидно, что при дифференцировании (13) будет получено выражение (10). Алгоритм идентифи-
кации (12), (13) обеспечивает сходимость ˆ (t) к . Для поиска параметра достаточно провести сле-
дующую несложную процедуру ˆ ˆ (t) , где ˆ – оценка частоты .
Пример
Рассмотрим сигнал вида (2) с параметрами 0, 5 ; A 2 ; 1 . На рис. 2 представлена времен-
ная диаграмма сигнала y(t) . На рис. 3, а, б, представлены результаты моделирования алгоритма идентификации параметра с коэффициентами k 1 и k 0, 01 соответственно. На рис. 4, а, б, представлены графики переходных процессов оценки ˆ при коэффициентах алгоритма идентификации k 0, 001 и k 0, 2 для сигнала y(t) с параметрами 2 ; A 5 ; 2 .
30 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2012, № 1 (77)
С.В. Арановский, А.А. Бобцов, А.А. Пыркин
Рис. 2. Временная диаграмма сигнала y(t)
ˆ
ˆ
аб Рис. 3. Переходные процессы для переменной ˆ при 0, 5 : k 1(а); k 0, 01(б)
ˆ ˆ
аб
Рис. 4. Переходные процессы для переменной ˆ при 2 : k 0, 001 (а); k 0, 2 (б)
Заключение
Рассмотрена задача идентификации частоты синусоидального сигнала с линейно меняющейся частотой. Получена динамическая модель генератора сигнала вида y(t) Asin(t2 / 2 ) , предложен ал-
горитм идентификации вида (12), (13), обеспечивающий сходимость настраиваемого параметра ˆ к истинному значению . Предложенный алгоритм имеет настроечный коэффициент k, варьируя который можно изменять быстродействие алгоритма.
Работа выполнена при финансовой поддержке ФЦП «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2007–2013 годы» (контракт № 11.519.11.4007).
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2012, № 1 (77)
31
УСОВЕРШЕНСТВОВАННЫЙ АЛГОРИТМ АНАЛИЗА СИГНАЛОВ …
Литература
1. Арановский С.В., Бобцов А.А., Кремлев А.С., Лукьянова Г.В., Николаев Н.А. Идентификация частоты смещенного синусоидального сигнала // Автоматика и телемеханика. – 2008. – № 9. – C. 3–9.
2. Арановский С.В., Бобцов А.А., Пыркин А.А. Адаптивный наблюдатель неизвестного синусоидального выходного возмущения для линейного объекта // Автоматика и телемеханика. – 2009. – № 11. – C. 108–116.
3. Арановский С.В., Бобцов А.А., Кремлев А.С. Компенсация конечномерного квазигармонического возмущения для нелинейного объекта // Известия РАН. Теория и системы управления. – 2006. – № 4. – C. 14–21.
4. Арановский С.В., Бобцов А.А., Никифоров В.О. Синтез наблюдателя для нелинейного объекта в условиях гармонического возмущения, приложенного к выходной переменной // Научно-технический вестник СПбГУ ИТМО. – 2010. – № 3. – C. 32–39.
5. Marino R., Tomei R. Global Estimation of Unknown Frequencies // IEEE Transactions on Automatic Control. – 2002. – № 47. – P. 1324–1328.
6. Hou M. Amplitude and frequency estimator of a sinusoid // IEEE Transactions on Automatic Control. – 2005. – V. 50. – P. 855–858.
7. Marino R., Santosuosso G.L. Output regulation for a class of nonlinear systems with unknown sinusoidal disturbances // Proceedings of the 44th IEEE Conference on Decision and Control. – 2005. – C. 3321–3326.
8. Stanislav Aranovskiy, Alexey Bobtsov, Artem Kremlev, Nikolay Nikolaev, Olga Slita. Identification of frequency of biased harmonic signal // European Journal of Control. – 2010. – № 4. – Р. 129–139.
9. Васильев В.Н., Томасов В.С., Шаргородский В.Д., Садовников М.А. Состояние и перспективы развития прецизионных электроприводов комплексов высокоточных наблюдений // Изв. вузов. Приборостроение. – 2008. – № 6. – C. 5–12.
10. Аракелян А.К., Афанасьев А.А. Вентильные электрические машины. – М.: Энергоатомиздат, 1997. – 236 с.
11. Erm T., Sandrock S. Adaptive correction of periodic errors improves telescope performance // Proceedings of the American Control Conference. – 2005. – V. 6. – Р. 3776– 3777.
Арановский Станислав Владимирович Бобцов Алексей Алексеевич
Пыркин Антон Александрович
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, кандидат технических наук, ст. научный сотрудник, s.aranovskiy@gmail.com
– Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, доктор технических наук,
профессор, декан, bobtsov@mail.ifmo.ru – Санкт-Петербургский национальный исследовательский университет ин-
формационных технологий, механики и оптики, кандидат технических наук, мл. научный сотрудник, a.pyrkin@gmail.com
32 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2012, № 1 (77)