Информационные технологии
УДК 621.372.542
КОРРЕКЦИЯ ДИНАМИЧЕСКОЙ ПОГРЕШНОСТИ ИНЕРЦИОННЫХ ИЗМЕРИТЕЛЬНЫХ ПРЕОБРАЗОВАТЕЛЕЙ С ПЕРЕДАТОЧНОЙ ФУНКЦИЕЙ ВТОРОГО ПОРЯДКА
К.Л. Куликовский, П.К. Ланге
Самарский государственный технический университет 443100, г. Самара, ул. Молодогвардейская, 244
Рассмотрен метод коррекции динамической погрешности инерционных измерительных преобразователей с передаточной функцией 2-го порядка на основе коррекции аппаратной функции преобразователей с использованием сплайн - аппроксимации дискретных значений сигнала, а также его первой и второй производных.
Ключевые слова: измерительный преобразователь, аппаратная функция, сплайн-аппроксимация, динамическая погрешность.
Большинство измерительных преобразователей (ИП) характеризуется динамической погрешностью, которая, собственно, и определяет верхнюю границу их частотной характеристики. Динамическая погрешность существенно искажает получаемую с первичных измерительных преобразователей информацию. К таким преобразователям относятся различные термопреобразователи, акселерометры, электрохимические датчики гидрохимических параметров и т. д. Динамические измерительные преобразователи характеризуются передаточной функцией обычно первого или второго порядков.
В приборах для измерения анализа состава и свойств веществ (хроматографах, спектрометрах) наличие аппаратной функции вызывает уширение спектральных линий, в результате чего ухудшается разрешение таких пиков, близко расположенных друг к другу.
Аналогичная проблема возникает и при реализации многих других аналитических методов.
В хроматографических и электрохимических анализаторах, например, используются различного рода детекторы (первичные ИП), формирующие собственно аналитический сигнал в функции от времени и характеризующиеся инерционностью, что также вызывает расширение аналитических пиков и смещение их положений на спектрограмме.
Динамическая погрешность зависит от постоянных времени ИП, а также скорости изменения измеряемой величины и может быть скорректирована с использованием аппаратных и программных методов [1].
Константин Лонгинович Куликовский (д.т.н., проф.), профессор, каф. информационноизмерительной техники.
Петр Константинович Ланге (д.т.н., проф.), профессор, каф. информационно-измерительной техники.
62
Применение аппаратных методов сопряжено с решением ряда проблем. Снижение динамической погрешности путем уменьшения постоянной времени ИП, как правило, снижает его надежность и долговременную стабильность. Кроме того, уменьшение инерционности ИП связано со значительными технологическими трудностями.
В связи с этим все большее распространение получает коррекция динамических погрешностей с использованием программных методов.
В ряде случаев коррекция динамической погрешности представляет собой обратную задачу [1], то есть задачу восстановления входного сигнала по известной информации об операторе А (аппаратной функции) ИП.
В частности, коррекция аппаратной функции ИП первого порядка на основе сплайн-аппроксимации дискретных значений измерительного сигнала и его производной рассмотрена в [2].
Рассмотрим задачу измерения мгновенных значений параметра х{(), который преобразуется ИП в сигналу(0, формируемый на его выходе (рис. 1).
x(t) . А y(t) А-1 2^)
ИП КЗ
Р и с. 1. Структурная схема тракта прохождения сигнала
При динамических измерениях, как правило, интерес представляет не выходной сигнал ИП у(0, а параметр х(0 . Поэтому задачей обработки результатов измерений является определение значений параметра х(0 по выходному сигналу у(0 и оператору А, описывающему динамические свойства ИП, что представляет собой решение задачи коррекции его аппаратной функции. Наиболее просто такая задача решается реализацией оператора А-;, обратного оператору А, с использованием корректирующего звена КЗ (см. рис. 1) в аппаратном или программном виде, который обрабатывает сигнал у(0. Однако по сущности своей такая задача является некорректной, поскольку обратный оператор должен реализовывать функцию предсказания сигнала, что физически реализовать невозможно [1].
В связи с этим корректное решение обратной задачи при измерениях динамических параметров может быть выполнено, если предусмотреть определенное запаздывание в формировании значений сигнала z(t) на выходе корректирующего звена, что не требует реализации функции предсказания.
Рассмотрим коррекцию аппаратной функции преобразователей, имеющих передаточную функцию второго порядка вида
К(Р) = (/, р + ф2р +1) ’ где К0 - статический коэффициент преобразования ИП, Т\, Т2 - постоянные времени ИП.
При К0=1 передаточная функция корректирующего звена принимает вид
W ■1( р) = ^4 = (т1 р + 1)(г2 р +1)=Т1Т1 р2 +Т + Т2 )р + 1, у( р)
что соответствует реализуемому им дифференциальному уравнению
*2 у (ґ) + (т, + т2)*у(')
(1)
л2 х 1 27 *
Таким образом, это звено должно реализовать функцию одно- и двукратного дифференцирования сигнала у(ґ) и сложения его производных с самим сигналом (рис. 2).
Р и с. 2. Структурная схема корректирующей цепи
Для решения такой задачи предлагается использовать цифровой фильтр, реализующий сплайн-аппроксимацию дискретных значений сигнала у^), а также его производной.
При использовании параболической сплайн-аппроксимации на п-м дискретном участке сигнал описывается параболической функцией (рис. 3)
У,
г (ґ) = а 2 [п ]2 + а1 Я ] + а о [я]
(2)
где а2 , а\ , а0 - коэффициенты аппроксимации.
Коэффициенты пятиточечной параболической сплайн-функции, аппроксимирующей сигнал у^) по его дискретным значениям, определяются следующими выражениями [2]:
а0п ] = Х6 (- у [п - 2]+4 у[п -1]+10 у[п ]+4 у[п +1] - у[п + 2]), ^
а1[п] = Х/8 (у П - 2]-6 у [п -1]+ 6 у [п +1]-у [п + 2]), >- (3)
а2[п ] = Ц6 (- у[п - 2] + 7 у[п -1] - 6 у[п ] - 6 у[п +1] + 7 у[п + 2] - у[п + 3]). „
Коэффициенты пятиточечной параболической сплайн-функции, аппроксимирующей производную сигнала у^), определяются выражениями [2]:
(4)
ь0[п] = Х2 (у[п - 2] - 8у[п -1] + 8у[п +1] - у[п + 2]),
Ь1 [п] = % (- у[п - 2] + 10у[п -1] -18у[п] + 10у[п + 1] - у[п + 2]),
ь2 [п] = 112 (- у[п - 2] - 9у[п -1] + 22у[п] - 22у[п +1] + 9у[п + 2]) - у[п + 3].
Сплайн-аппроксимация второй производной сигнала у(^ на п-м участке его дискретизации определяется выражением
у"п 0) = Л 2 М 2 + Л1 №) + Л 0 И •
При использовании пятиточечного параболического сплайн-фильтра коэффициенты Л2[п], Л\[п], Л0[п] записываются [3] в виде
Л0 [п] = 14 (- у[п - 2] + 2у[п]- у[п + 2]), >,
^[«1 = 29(- у[п - 2] + 2у[п -1]-2у\п +1] + у[п + 2], „ (5)
Л2 [п] = 19 (у[п - 2] - 3у[п - 1] + 2у[п] + 2у[п +1] - 3у[п + 2] + у[п + 3]). ^
Уn(t)
Г[п-2] Г [п-1] Г [п] Г [п+1] Г [п+2] t [п+3] Г
Р и с. 3. Аппроксимация выходного сигнала измерительного преобразователя параболической сплайн-функцией
Таким образом, сигнал zK(Г) на выходе корректирующего фильтра на п-м интервале дискретизации описывается соотношениями
zk (Г) = а2 [п] Г2 + а1 [п] Г + {Т1 + Т2 )(ь2 [п] Г2 + Ь1 [п] Г + Ь0 [п]) +
+ Т1Т2 (Л2 [п] Г2 + Л1 [п] Г + Л0 [п])=(а0 [п]+(Т1 + Т2 )Ь2 [п] + Т1Т2Л2 [п]) Г2 +
+ (а1 [п]+(Т1 + Т2 )Ь1 [п] + Т1Т2Л2 И Г + (Т1 + Т2 К [п] + Т1Т2Л0 [п].
Следовательно, выходной сигнал zK(Г) также определяется параболической сплайн - функцией
гк (г) = с2[п] Г2 + С1[п] Г + с0[п1 (6)
где С2 [п] = а2 [п]+(Т1 + Т2 )ь2 [п] + Т1Т2Л2 [п],
С1 [п] = а1 [п] + (Т1 + Т2 )^ [п] + ТТ ^1 [п],
с0 [п] = а0 [п] + (Т1 + Т2 К [п] + Т1Т2Л0 [п] .
Как видно из (6), дискретные значения этого сигнала определяются выражением
^0 [п] = а0 [п] + (Т1 + Т2 К [п] + Т1Т2Л0 [п] . (7)
Из (5) следует, что сигнал zK(Г) формируется с запаздыванием в три дискретных интервала.
Аналогичным образом могут быть определены выражения (7) и при использовании параболических сплайн-функций, определенных на четном числе точек дискретизации, например четырехточечных сплайн-функций.
Свойства разработанных корректирующих сплайн-фильтров удобно рассмотреть на основе анализа их частотных характеристик.
Частотные характеристики сплайн-фильтров, аппроксимирующих значения сигнала, а также его первой и второй производных, определяются из выражений для (3), (4) и (5) для коэффициентов а0, Ь0, с0 соответственно при использовании их 2-
преобразований и при подстановке гп = в]тп :
5 - 4cos та- соэ2та
ао (Іта) = “ ьо (Іта) = І dо (Іта) =
8
8sin та - sin2та
6
1 - cos 2та
2
Амплитудно-частотная характеристика корректирующего фильтра в этом случае определяется выражением
Кк (Іта) = ^(та)+ І1 (та),
где
/ ч 5 + 4ТТ cos т 1 + 4Т(Т^ Т( \ Т + Т2 / . ч
Л(т) =----------------------------cos2m , I(т) =------------^тт - sln2т).
8 2 8 6
Фазочастотная характеристика (ФЧХ) корректирующего фильтра определяется выражением
/ ч I(т)
Фк (т) = ■
Качество коррекции гармонического сигнала удобно оценивать частотными характеристиками всей цепи преобразования, состоящей из детектора аналитического сигнала и корректирующего фильтра. АЧХ такой цепи
К£(т) = KD (т)кк (т) , где KD (т) - АЧХ инерционного измерительного преобразователя:
KD (т)= I----------Г 1 I =Т.
V1 + (Т1т)2 V1 + (Т2т)2
В качестве примера можно привести АЧХ всей цепи преобразования сигнала для случая, когда детектор представляет собой апериодическое звено второго порядка (1) для постоянных времени детектора Т]=10 , Т2 = 20 и использования пятиточечных и четырехточечных параболических сплайн-фильтров. Графики таких характеристик в зависимости от относительной частоты / (/ = 1N, где N - число интервалов дискретизации на периоде гармонического сигнала) приведены на рис. 4.
Из рассмотрения этого рисунка видно, что использование корректирующего преобразования на основе параболических сплайн-фильтров позволяет в несколько раз расширить частотный диапазон преобразования сигнала в случае, когда детектор аналитического прибора представляет собой инерционное апериодическое звено второго порядка.
Качество восстановления формы реальных измерительных сигналов необходимо рассмотреть на примере сигналов определенной формы. Известно, что реальные аналитические параметры х(Г), преобразуемые инерционным ИП в сигнал у(Г), во
многих случаях представляют собой симметричные пики колоколообразной формы. В связи с этим представляет интерес задача восстановления формы таких пиков с помощью описанных выше алгоритмов.
Рассмотрим в качестве примера прохождение через инерционный детектор аналитического прибора единичного сигнала гауссовой формы
х(ґ) = ехр
ґ - 5 1.5
Весовая функция ИП, имеющего передаточную функцию (1) второго порядка, определяется выражением
к(ґ) =
1
Ті - т2
ехр
ґ
ехр
11 У
ґ
12 У
• 1(ґ).
(8)
/
Р и с. 4. Амплитудно-частотные характеристики всей цепи преобразования для детектора сигнала, представляющего собой апериодическое звено второго порядка:
Аі(і) - апериодического звена второго порядка; А04(і) - всей цепи преобразования сигнала при использовании в корректирующем фильтре четырехточечной сплайн-функции; А05(і) - всей цепи преобразования сигнала при использовании в корректирующем фильтре пятититочечной сплайн-функции
При входном сигнале х^) гауссовой формы (8) сигнал у(^ на выходе детектора анализатора определяется функцией свертки
Т - Т
11 2
ехр
-(ґ- г - 5)2 1.5
ехр( - Т) - ехР( - Т~} Т1 Т2
1(ґ) сіг.
Оценка возможностей работы корректирующего фильтра второго порядка при использовании вышеописанного алгоритма сплайн-аппроксимации сигнала у(^, его первой и второй производных осуществляется аналогичным образом.
При постоянных времени измерительного преобразователя ИП Т1 = 10 и Т2=1 графики входного сигнала х(£), его выходного сигнала у(^ и сигнала на выходе корректирующего фильтра z(t) для пятиточечного алгоритма сплайн-аппроксимации сигнала у(^ и его производной приведены на рис. 5.
Из рис. 5 видно, что разработанный алгоритм корректирующей цифровой фильтрации достаточно эффективно восстанавливает по форме сигнал х(£), действующий на входе инерционного детектора. Запаздывание скорректированного сиг-
1
нала z(t) относительно нескорректированного х^) составляет 3 дискретных интервала. Для удобства восприятия это запаздывание не отражено на рис. 5. Корректирующий сплайн-фильтр восстанавливает форму пика сигнала с погрешностью порядка 20 % .
1*0)
^A0
I " \ *~v I \, K0
-
о 5 10 15 20 25 30 35 40
Р и с. 5. Коррекция сигнала x(t) гауссовой формы, искаженного инерционным преобразователем, имеющим передаточную функцию второго порядка:
y(t) - сигнал на выходе детектора; z(t) - скорректированный сигнал
Таким образом, использование сплайн-аппроксимации измерительного сигнала, формируемого инерционными измерительными преобразователями, позволяет снизить их динамические погрешности и тем самым восстановить аналитические сигналы, действующие на их входах.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Солопченко Г.Н. Обратные задачи в измерительных процедурах // Измерение, контроль, автоматизация. - М: ЦНИИТЭИприборостроения. - 1983. - Вып. 2(46). - С. 32-46.
2. Ланге П.К., Карпова Н.Е. Аппроксимативные методы коррекции динамической погрешности // Датчики и системы. - Ин-т проблем управления РАН. - М: 2001. - №12. - С. 10-12.
3. Ланге П. К. Дифференцирующие цифровые фильтры дискретизированных сигналов со сплайн-аппроксимацией их второй производной // Тематич. сборник науч. трудов «Информационноизмерительные и управляющие системы». Самара: Поволж. регион. науч.-техн. центр Метрологической академии РФ. 2003. - Вып. 12. - С. 43-51.
Статья поступила в редакцию 2 ноября 2011 г.
CORRECTION OF THE DYNAMIC ERROR OF THE INERTIAL SENSORS HAVING THE SECOND ORDER TRANSFER FUNCTION
K.L. Kulikovskij, P.K. Lange
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100
In article is considered a method for the dynamic error correction of the inertial sensors having the second order transfer function on the basis of the correction of the sensor instrumental function using a spline-approximation of the discrete values of the signal and its first and second derivatives.
Keywords: sensor, instrumental function, spline approximation, dynamic error.
Konstantin L. Kulikovskij (Dr. Sci. (Techn.)), Professor. Petr K. Lange (Dr. Sci. (Techn.)), Professor.