CC BY
80
5. Вольфовский Б.Н., Шибаев С.С., Роздобудько В.В. Алгоритмы вычисления частоты в акустооптических измерителях параметров радиосигналов // Известия высших учебных заведений России. Радиоэлектроника. - 2008. - № 4. - С. 38-46.
6. Вольфовский Б.Н., Шибаев С.С., Роздобудько В.В. Алгоритмический способ уменьшения погрешностей в акустооптических измерителях // Известия высших учебных заведений России. Радиоэлектроника. - 2011. - № 5. - С. 51-56.
7. Роздобудько В.В., Пелипенко М.И. Быстродействующий измеритель параметров СВЧ-радиосигналов // Специальная техника. - 2006. - № 1. - С. 28-36.
8. Балакший В.И., Парыгин В.Н., Чирков Л.Е. Физические основы акустооптики. - М.: Радио и связь, 1985. - 280 с.
9. Гудмен Дж. Введение в Фурье-оптику: Пер. с англ. - М.: Мир, 1970. - 364 с.
Статью рекомендовал к опубликованию д.т.н., профессор В.И. Марчук.
Шибаев Станислав Сергеевич - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: sshib75@mail.ru; 347922, г. Таганрог, ул. Шевченко, 2; тел.: 88634312482; лаборатория оптоэлектроники; с.н.с.; к.т.н.
Shibaev Stanislaw Sergeyevich - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: sshib75@mail.ru;
2, Shevchenko, Taganrog, 347922, Russia; phone: +78634312482; laboratory of optoelectronics; senior scientist; cand. of eng. sc.
УДК 621.82: 621.397
Л.К. Самойлов, И.И. Турулин, С.А. Киракосян, А.Д. Вартенков
ПОГРЕШНОСТИ ВОССТАНОВЛЕНИЯ СИГНАЛОВ В СИСТЕМАХ УПРАВЛЕНИЯ И КОНТРОЛЯ
Частота дискретизации системы управления и контроля выбирается в результате решения прямой или обратной задачи распределения погрешностей всей системы. При решении прямой задачи выделяются погрешности устройства дискретизации и погрешности устройства восстановления. Соответственно различают частоту дискретизации с точки зрения процедуры дискретизации и процедуры восстановления. Оптимальным вариантом следует считать ситуацию, когда частота дискретизации с точки зрения процедуры дискретизации будет равна частоте дискретизации с точки зрения процедуры восстановления. Такое решение может быть найдено в итерационном цикле с перераспределением погрешностей. В статье рассматриваются составляющие погрешности восстановления, что позволит организовать такой итерационный цикл. Погрешность устройства восстановления состоит из методической погрешности восстановления, погрешности цифрового восстановления (передискретизации), погрешности экстраполяции нулевого порядка на основе ЦАП и погрешности аналоговых фильтров нижних частот. На выбор частоты дискретизации также влияют инструментальные погрешности устройства восстановления. В системах управления добавляется погрешность за счет задержки информации в устройстве восстановления. В статье рассматриваются указанные погрешности трехступенчатого устройства восстановления. Показано, что значение рабочей частоты устройства цифрового восстановления не влияет на величину методической погрешности восстановления. Увеличение рабочей частоты устройства цифрового восстановления уменьшает задержку сигнала в экстраполяторе на основе ЦАП, что косвенно снижает погрешность за счет задержки информации. Даются оценки погрешностей устройств, позволяющие выбрать требуемую частоту дискретизации.
Частота дискретизации; восстановление сигналов; запаздывание сигналов; граничная частота; системы управления.
L.K. Samoilov, I.I. Turulin, S.A. Kirakosyan, A.D. Vartenkov
THE ERRORS OF SIGNALS RESTORATION IN SYSTEMS MANAGEMENT
AND CONTROL
Frequency of digitization of a control system and the control gets out as a result of the decision of a direct or return problem of distribution of errors of all system. At the decision of a direct problem errors of the device of digitization and an error of the device of restoration are allocated. Accordingly distinguish frequency of digitization from the point of view ofprocedure of digitization and restoration procedure. It is necessary to consider as an optimum variant a situation when frequency of digitization from the point of view of digitization procedure will be equal to frequency of digitization from the point of view of restoration procedure. Such decision can be found in an iterative cycle with redistribution of errors. In article making errors of restoration that will allow to organise such iterative cycle are considered. The error of the device of restoration consists of a methodical error of restoration, an error of digital restoration (redigitization), an error of extrapolation of a zero order on the basis of DAC and errors of analogue filters of the bottom frequencies. The choice of frequency of digitization is influenced also by tool errors of the device of restoration. In control systems the error at the expense of an information delay in the restoration device is added. In article the specified errors of the three-stage device of restoration are considered. It is shown, that value of working frequency of the device of digital restoration does not influence size of a methodical error of restoration. The increase in working frequency of the device of digital restoration reduces a signal delay in eksirapolyator on the basis of DAC, that indirectly reduces an error at the expense of an information delay. Estimations of errors of the devices, allowing to choose demanded frequency of digitization are given.
Sampling frequency; signal restoration; signal delay; cutoff frequency; control systems.
В системах управления и контроля частота дискретизации сигналов выбирается в процессе решения прямой и обратной задач распределения погрешностей системы [1]. При этом выделяется погрешность на процедуру дискретизации и на процедуру восстановления информации. Соответственно различают частоту дискретизации (/д) с точки зрения процесса дискретизации () и процесса восстановления
информации (/I) [1]. Из этих двух частот выбирают или наибольшую или в итерационном процессе выравнивают эти частоты [2]. Устройства временной дискретизации сигналов характеризуются тремя параметрами, определяющими / :
♦ погрешностью наложения спектров;
♦ инструментальной погрешностью устройств дискретизации;
♦ временем задержки информации.
Устройства восстановления сигналов также определяются тремя параметрами:
♦ методической погрешностью восстановления информации;
♦ инструментальной погрешностью устройств восстановления;
♦ временем задержки информации.
Современные устройства восстановления могут иметь три ступени восстановления (рис. 1):
♦ цифрового восстановления (ЦВ);
♦ экстраполяции нулевого порядка на основе ЦАП;
♦ аналоговой фильтрации с помощью ФНЧ.
Особенностью цифрового восстановления сигнала с частотой дискретизации / является тот факт, что на выходе устройства или программы восстановления будет тоже цифровой сигнал. Но отсчеты этого сигнала будут идти с частотой восстановления (/ ), которая существенно (в десятки и сотни раз) выше / . Значение частоты / должно быть кратно значению частоты / . Увеличение частоты отсчетов значительно облегчает процедуру получения аналогового сигнала.
Цифровой
сигнал
Цифровой
экстраполятор;
Цифровой интерполятор; Цифровой ФНЧ
Цифровой
экстраполятор
нулевого
порядка
ФНЧ
Аналоговый
сигнал
и А
Рис. 1. Структурная схема трехступенчатого устройства восстановления
Произвольная непрерывная функция и(/) после равномерной дискретизации по аргументу ^ представляет собой решетчатую функцию (?), которая
совпадает с исходной в точках дискретизации.
В устройстве восстановления при интерполяции по т отсчетам решетчатой функции *¥д (?) восстанавливается исходный сигнал в промежутках между отсчетами. В точках отсчетов восстановленный сигнал совпадает с решетчатой функцией. Такое восстановление называется интерполяцией (т - 1)-го порядка.
При экстраполяции по т отсчетам решетчатой функции д (?) восстанавли-
вается исходный сигнал в промежутке между т-м и т+1 отсчетами. В точках отсчетов восстановленный сигнал совпадает с решетчатой функцией. Такое восстановление называется экстраполяцией (т - 1)-го порядка.
При использовании в ЦВ интерполирующих фильтров на входе устройства ЦВ стоит квантователь, который заполняет промежутки между двумя отсчетами с частотой /д и нулевыми отсчетами с частотой /в. В интерполирующем фильтре
производится восстановление недостающих отсчетов.
Временная диаграмма работы первой и второй ступеней восстановления при использовании интерполяции в ЦВ приведена на рис. 2.
Л
Рис. 2. Временная диаграмма восстановленного сигнала на выходе двухступенчатого устройства восстановления
Методическая погрешность устройства цифрового восстановления на основе экстраполяторов и интерполяторов с порядком восстановления т > 1.
3S
(уМ1) представляет собой погрешность непрерывного устройства восстановления. Величина этой погрешности зависит от способа восстановления (интерполяция, экстраполяция), от вида и степени (т) восстанавливающего полинома, а также от величины частоты дискретизации с точки зрения погрешности восстановления /д .
На практике чаще всего устройства восстановления проектируются с использованием полинома Лагранжа [3].
Пусть имеется непрерывная функция и(1) и соответствующая ей решетчатая
функция /(/), заданная в т+1 точке /(^) , /^1), ... , /(^т).
В этих точках решетчатая функция /(/) совпадает с непрерывной функцией
Процесс получения непрерывной функции иапр (г) из решетчатой называется восстановлением или аппроксимацией.
Значения непрерывной функции в промежутках между отсчетами решетчатой функции можно аппроксимировать многочленом, который называется интерполяционной формулой Лагранжа:
Если аппроксимация решетчатой функции происходит внутри участка
[ 10, Хт ], то восстанавливающее устройство называется интерполятором.
Если решетчатая функция восстанавливается на шаг вперед на участке
[ 1п, 1т+1 ], то восстанавливающее устройство называется экстраполятором.
Абсолютное значение максимальной величины методической погрешности восстановления определяется разностью
Максимальное значение приведенной методической погрешности восстановления равно
плитудной спектральной плотности сигнала датчика при нулевой частоте.
Если в устройстве цифрового восстановления используется полином Лагранжа, то максимальное значение приведенной методической погрешности для случая экстраполяции определяется как [3]:
и(г).
т
(1)
при т > 1
1 при т = 0
при т иапр(г) = и (г).
(2)
ли (г) = и(г)-иап(г).
(3)
(4)
где А0 - максимальное значение непрерывной функции и(() или значение ам-
(5)
где Тд - интервал дискретизации; А - максимальное значение непрерывной
функции U(t) ;
M m+1
мак
- абсолютное значение максимальной величины т+1 про-
изводной восстанавливаемого сигнала датчика; т - порядок восстанавливающего полинома.
Обратим внимание, что значение частоты сигналов цифрового восстановления первой ступени / не влияет на величину методической погрешности устройства восстановления. Значение этой частоты влияет на величину методической погрешности цифроаналогового устройства восстановления на основе экстраполя-торов нулевого порядка, который должен присутствовать на второй ступени восстановления.
Для случая применения интерполяции в устройстве восстановления может быть предложен удобный и наглядный способ вычисления максимальной приведенной погрешности [4]:
(гМаU = К(т)-!|мМак, „ .
А0 (6)
fm+1\ грт+1 1 мак д
где К(т) - поправочный коэффициент.
Проведенные исследования [4] показали. что наилучшее приближение с допустимой погрешностью не более 1 % при т < 6 дает аппроксимация К(т) в виде
гг 1 л( -483,94914Е-3) - (498,07533Е-3)т + (84,64261Е-3)т2 - (6,17667Е-3)т3
К (т)= Ш .
Если необходимо использовать порядки интерполяции выше пятого, то можно воспользоваться зависимостью, которая дает аппроксимацию К(т) с максимальной погрешностью 2 % при 5 <т< 16:
К(т) = К)(-581,01970Е-3) - (378,02360Е-3)т + (46,33452Е-3)т2 -( 2,97881Е-3)п3 + (72,92118Е-6)т4 ^
Полученные выражения для интерполяционных полиномов позволяют находить аналитически значения погрешностей и частоты дискретизации для двух ступеней восстановления практически используемого ряда порядков полиномов.
Методическая погрешность интерполирующего фильтра (уЦТ ) при условии,
что частота среза фильтра (а ) больше 1 , будет определяться неравномерно' 2/д
стью амплитудно-частотной характеристики фильтра в полосе пропускания.
Величина среднеквадратической приведенной инструментальной погрешности первой ступени определяется погрешностью округления при вычислении восстановленных значений.
Приведенная инструментальная погрешность ЦАП определяется его разрядностью и равна половине младшего разряда.
Время задержки в рассматриваемом цифровом устройстве восстановления
(^) определяется временем задержки выбранного алгоритма восстановления
/ ^алг \ / ,выч \
((м ) и временем задержки при вычислениях ():
. _ .алг , .выч (9)
1М 1М ^ 1М • (9)
Если в цифровом устройстве выбран экстраполяционный алгоритм восстановления, то ^ _ 0. При выборе интерполяционного алгоритма время задержки
равно
С = (т-\)Тд. (10)
Время задержки при вычислениях зависит от числа команд при реализации программы восстановления и от скорости работы вычислительного устройства.
Время задержки второй ступени определяется параметрами ЦАП (^п ).
Третья ступень восстановления в виде ФНЧ (рис. 1) предназначена для выделения постоянной составляющей переменного ступенчатого напряжения на выходе ЦАП и уменьшения уровня пульсаций выходного напряжения с частотой / (рис. 3).
При этом происходит уменьшение методической погрешности второй ступени восстановления.
Диз = 0,5Ди2. (11)
Л
Рис. 3. Временная диаграмма сигнала на выходе трехступенчатого устройства
восстановления
Уровень пульсаций и методическая погрешность второй ступени уменьшаются пропорционально величине коэффициента передачи АЧХ фильтра третьей ступени ч(р)) на частоте р = 2ж{ъ.
Тогда можно записать, что максимальное приведенное значение методической погрешности трехступенчатого устройства восстановления при использовании экстраполяции в первой ступени (уМ )жс будет равно
(Мм Хс = \{М:+А-ТГ + 0,5[7 + Кфнч (р)] \мЦ-Те}. (12)
А0
Методическая погрешность восстановления при использовании интерполяции в первой ступени (мм)инт будет равна
(Мм )инт = -\{К(т) • •Тдт+\ + 0,5[\ + Кфнч (Р )] • \}. (13)
Методическая погрешность восстановления при использовании интерполирующего ФНЧ в первой ступени (мм)фнч будет равна
(Мм )фнч = мММт +-\{0,5[\ + Кфнч (Р )]<№Ц • Тв)}. (14)
А0
Как видно из рис. 3, постоянная составляющая ступенчатого напряжения на выходе ЦАП (жирная непрерывная линия) повторяет напряжение на выходе первой ступени, но с задержкой, равной половине такта цифрового восстановления [5]:
С=Р <15>
•Iв
п ~ ,вт
В этом случае, если рассматривать восстановленный сигнал с задержкой Хзад, то методическая погрешность восстановления может быть представлена как:
(мм =-\-{мтх\• тт+\+0,5[Кфнч \ )]| \Ц\};
А0
(Мм )нт = -^{К (т) • Мт:1\• Тт+\ + 0,5[ Кфт р)] \иЦ • Тв};
А0
(Мм )фНч =мМт +^{0,5[Кфнч (Рв )](\мЦ •Тв)}. (16)
Инструментальная погрешность третьей ступени характеризуется величиной среднеквадратической приведенной инструментальной погрешности фильтра
( Мф ). Инструментальная погрешность фильтра имеет две составляющие:
Мху - среднеквадратическая приведенная погрешность, определяемая нестабильностью коэффициента передачи фильтра как четырехполюсника;
Мдр - среднеквадратическая приведенная погрешность за счет дрейфа и смещения нуля.
Эти погрешности независимы друг от друга, что позволяет при оценке погрешности фильтра использовать геометрическое суммирование отдельных погрешностей:
(Мф)2 = (Мку)2 + (Мдр )2 . (14)
Время задержки информации в трехступенчатом устройстве восстановления определяется временем задержки в каждой ступени. Кроме того, как видно из рис. 3, при использовании третьей ступени в виде ФНЧ появляется дополнительная задержка второй ступени. Величина этой задержки равна половине такта цифрового восстановления [5]:
* ет =_\_ (15)
1 зад _ г . (15)
21 в
Тогда время задержки сигнала в трехступенчатом устройстве восстановления будет равно
> _<алг .^выч
1М = 1М ^ 1М ^ 1ЦАП ^ 1 зад ^ 1ФНЧ , (16)
где 1фнч - время задержки сигнала в ФНЧ третьей ступени.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Самойлов Л.К. Итерационные алгоритмы выбора частоты дискретизации аналоговых сигналов в цифровых системах управления и контроля // Известия ЮФУ. Технические науки. - 2012. - № 2 (127). - С. 43-53.
2. Самойлов Л.К., Жуков А.В. Выбор частоты дискретизации реальных сигналов // Приборы. Контроль. Управление. - 2009. - № 1. - С. 29-31.
3. Баранов Л.А. Квантование по уровню и временная дискретизация в цифровых системах управления. - М.: Энергоатомиздат, 1990. - 304 с.
4. Самойлов Л.К., ЧерновА.М. Аналитическое представление восстанавливающего оператора при интерполяции по Лагранжу // Известия ЮФУ. Технические науки. - 2009.
- № 2 (91). - С. 233-236.
5. Maheshwari R., Bharadia M., Gupta M. Multirate DSP and its technique to reduce the cost of the analog signal conditioning filters // International Journal of Computer Applications.
- 2010. - № 10. - P. 27-34.
Статью рекомендовал к опубликованию д.т.н., профессор Я.Е. Ромм.
Самойлов Леонтий Константинович - Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Южный федеральный университет»; e-mail: samoilov@fep.tti.sfedu.ru; 347928, г. Таганрог, пер. Некрасовский, 44; тел.: +78634311193; кафедра информационных измерительных технологий и систем; профессор.
Турулин Игорь Ильич - e-mail: iiturulin@sfedu.ru; кафедра информационных измерительных технологий и систем; профессор.
Киракосян Степан Айрапетович - e-mail: asni@fep.tti.sfedu.ru; кафедра информационных измерительных технологий и систем; студент.
Вартенков Андрей Дмитриевич - e-mail: asni@fep.tti.sfedu.ru; кафедра информационных измерительных технологий и систем; аспирант.
Samoilov Leonty Konstantinovich - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”; e-mail: samoilov@fep.tti.sfedu.ru; 44, Nekrasovskiy, Taganrog, 347928, Russia; phone: +78634311193; the department of information measuring technologies and systems; professor.
Turulin Igor Il’ich - e-mail: iiturulin@sfedu.ru; the department of information measuring technologies and systems; professor.
Kirakosyan Stepan Ajrapetovich - e-mail: asni@fep.tti.sfedu.ru; the department of information measuring technologies and systems; student.
Vartenkov Andrey Dmitrievich - e-mail: asni@fep.tti.sfedu.ru; the department of information measuring technologies and systems; postgraduate student.
