ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2017. № 2 (54)
Приборостроение, метрология и информационно-измерительные приборы и системы
УДК 621.372.542
КОРРЕКЦИЯ ДИНАМИЧЕСКОЙ ПОГРЕШНОСТИ ИНЕРЦИОННОГО ИЗМЕРИТЕЛЬНОГО ПРЕОБРАЗОВАТЕЛЯ*
П.К. Ланге
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
Рассмотрен метод коррекции динамической погрешности инерционного измерительного преобразователя с передаточной функцией, имеющей вид апериодического звена первого порядка. Метод основан на реализации обратной функции оператора, описывающего инерционный измерительный преобразователь в динамическом режиме. Метод коррекции использует параболическую сплайн-аппроксимацию дискретных значений измерительного сигнала, формируемого измерительным преобразователем, а также его первой производной. Приведен пример коррекции сигнала измерительного преобразователя конкретного типа. Показано, что использование математического метода коррекции на базе цифровых фильтров скользящей обработки дискретных значений сигнала позволяет в несколько раз снизить динамическую погрешность инерционного измерительного преобразователя. Предложено использовать описанный метод также и для коррекции сигнала измерительного преобразователя, передаточная функция которого соответствует инерционному звену второго порядка. Рассмотренный метод может быть легко реализован с использованием современных микропроцессорных измерительных систем.
Ключевые слова: измерительный преобразователь, инерционность, частотная характеристика, передаточная функция, измеряемый параметр, динамическая погрешность, сплайн-аппроксимация.
Измерительные преобразователи (ИП), как известно, обладают большей или меньшей инерционностью, которая затрудняет их использование в системах экспресс-анализа контролируемого параметра. Такие преобразователи, как термопары, терморезисторы, измерители электрохимических параметров жидкостей, измерители параметров индуктивных промышленных установок, могут характеризоваться длительностью переходного процесса от нескольких секунд до нескольких минут. Большая инерционность таких преобразователей вызвана, с одной стороны, длительностью процессов установления равновесия по окончании физических и химических процессов внутри конструкции измерительного преобразователя, с другой - инерционностью его чувствительных элементов. Этот фак-
*Работа выполнена при поддержке гранта РФФИ № 16-08-00252 А.
Петр Константинович Ланге (д.т.н., проф.), профессор кафедры «Информационно-измерительная техника». 58
тор вызывает появление динамической погрешности при регистрации быстро изменяющихся измерительных параметров.
Целью данной работы является разработка математического метода коррекции указанной динамической погрешности.
При динамических измерениях интерес представляет не выходной сигнал преобразователя у((), а контролируемый входной параметр х(1). Поэтому задачей коррекции сигнала у (О является определение значений параметра х(0 с учетом оператора К, характеризующего динамические свойства преобразователя и его аппаратную функцию. В современных измерительных системах аналоговый выходной сигнал у(() ИП преобразуется аналого-цифровым преобразователем в его дискретные значения у(п).
Задача коррекции может быть решена реализацией оператора К1, обратного оператору К, при соответствующей обработке дискретных значений сигнала
у(0 [1].
В наиболее распространенном случае инерционный ИП имеет передаточную функцию апериодического звена первого порядка
(1)
где Т - постоянная времени ИП.
При К0= 1 передаточная функция корректирующей цепи принимает вид
К-\р) ■ ■ Тр+Т, У ( р)
что соответствует реализуемому дифференциальному уравнению
2(1) ■ Тш + у(*). (2)
ш
Таким образом, такой метод коррекции динамической погрешности должен реализовать функцию дифференцирования сигнала у(() и сложения его производной с самим сигналом.
Фактически при этом решается обратная задача, однако для возможности ее физической реализации выходной сигнал корректирующей цепи должен иметь определенное запаздывание по отношению к параметру, контролируемому инерционным ИП.
Для решения задачи коррекции в [2] было предложено использовать цифровой фильтр, реализующий параболическую сплайн-аппроксимацию дискретных значений сигнала у(0, а также его производной.
При использовании параболической сплайн-аппроксимации на п-м дискретном участке сигнал описывается параболической функцией
Уи (* )■ Ь2 Н *2 + ЬТ Н * + Ь0 И ,
где а2 , а1, а0 - коэффициенты аппроксимирующей параболической функции.
Коэффициенты такой аппроксимации могут быть определены с использованием цифровых фильтров «скользящей» обработки дискретных значений сигнала у (О в реальном темпе времени. Число точек весовых функций таких фильтров обычно лежит в пределах 4.. .6.
В частности, коэффициенты цифровых фильтров пятиточечной сплайн-аппроксимации сигнала у(0 по его дискретным значениям у(п) определяются выражениями [2]:
Ь И = Хб (" У [п " 2] + 4 У [п " Ч +10 У Н + 4 У [п + !]" У [п + 2]);
Ь [п] = ^ (У [п - 2]- 6 у [п -1] + 6 у [п +1]- у [и + 2]) ; Ь [п]= У16(-У[п -2] + 7у[п -1]-6у[п]-6у[п +1] + 7у[п + 2]-у[п + 3] ).
(3)
'16*
Коэффициенты пятиточечной параболической сплайн-функции, аппроксимирующей производную сигнала у(0, определяются выражениями [3]:
С [п] = Ц2 (У [п - 2]- 8у [п -1] + 8 у [п +1]- у [п + 2]); С [п] = %(-х[п - 2] + 10х[п -1]- 18х[п] + 10х[п +1]-х[п + 2]); с2 [п] = Уп (У [п - 2]-9 у [п -1] + 22 у [п]- 22 у [п +1] + 9 у [п + 2]-у [п + 3]).
>
(4)
Таким образом, скорректированный сигнал определяется выражением
) = ^ [п] t2 + [п^ + ^ [п], (5)
где
¿о[ п] = ьо [п]+т ■ со [п];
[п] = Ь [п] + Т ■ с [п];
[п] = Ь2 [п] + Т ■ с [п].
(6)
При I = 0 из (4) следует
х(п) = а0[п] = Ьо [п] + Т■ Со [п]. (7)
Как видно из (3), скорректированный сигнал г^) формируется относительно дискретных значений выходного сигнала у(п) измерительного преобразователя с задержкой в 2 дискретных интервала.
Форма скорректированного сигнала (4) в идеальном случае должна соответствовать форме входного параметра х(0.
Использование сплайн-аппроксимации позволяет определить значения сигнала г(0 внутри интервалов дискретизации, что является достоинством рассматриваемого метода.
Рассмотрим в качестве примера коррекцию сигнала детектора по теплопроводности, широко распространенного в расходомерах микропотоков газов, в газоаналитических приборах, в хроматографических анализаторах и других приборах для анализа состава веществ.
Рассмотрим возможность коррекции динамической погрешности такого детектора с постоянной времени Т в выражении (1), равной 0,1 с.
Во многих практических случаях [5] контролируемый с помощью такого ИП параметр имеет вид функции, близкой к гауссовому закону.
Рассмотрим коррекцию выходного сигнала такого ИП при измеряемом параметре, изменяющемся в соответствии с гауссовой функцией единичной высоты (рис. 1) с шириной примерно 0,05 с:
г ^ - о.о5)2
х^) = ехр[---]. (8)
о.о15
При передаточной функции (1) ИП
W (p ) =
Tp +1
его переходная характеристика определяется выражением
_ t
h(t) = (1 - ) .
Выражение для сигнала y(t) на выходе ИП может быть получено с использованием интеграла Дюамеля:
1 t
y(t) =1 i
, (t-г-0.05)2
exp(--)
0.015
exp(-—
Для конкретного значения постоянной времени Т = 0,1 с передаточной функции (1) ИП график этой функции изображен на рис. 1. Как видно из рассмотрения графика, такой ИП характеризуется очень большой динамической погрешностью.
Рис. 1. Сигналы на входе и выходе инерционного ИП:
х - график изменения параметра на входе ИП; у - сигнал на выходе ИП
Функция скорректированного сигнала определяется при аппроксимации дискретных значений сигнала у(0 параболической сплайн-функцией с использованием выражений (3) - (4) для коэффициентов аппроксимации, а также выражений (5) - (6) для метода коррекции.
Графики этих функций представлены на рис. 2, при рассмотрении которого видно, что предложенный алгоритм достаточно эффективно восстанавливает по форме сигнал х(0, действующий на входе инерционного детектора, однако скорректированный сигнал г^) запаздывает на 2 интервала дискретизации по отношению к контролируемому параметру х(0.
При наличии в потоке газа нескольких примесей, размещенных в потоке на некотором расстоянии друг от друга, детектор регистрирует последовательно
соответствующие им пики. Такая ситуация характерна для различных хромато-графических анализаторов и ряда газоаналитических приборов.
X
0.5
0
0 4 8 12 16 20 24 28 32 36 п 0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 t с
Рис. 2. Диаграммы сигналов при коррекции динамической погрешности.
1 - параметр х() на входе инерционного детектора,
2 - сигнал у^) на выходе детектора,
3 - скорректированный сигнал г^-2), сдвинутый на 2 дискретных интервала влево,
п - номера дискретных отсчетов сигнала у(7).
На рис. 3 представлены диаграммы сигналов: входного параметра х(0 инерционного ИП, сигнала у(0, формируемого на его выходе, скорректированного сигнала г^) с использованием описанного алгоритма коррекции. Входной параметр х(0 в данном случае представляет собой сумму двух гауссовых функций вида (7), смещенных друг относительно друга на интервал 0,02 с.
Ширина гауссовых функций составляет примерно 0,06 с. Алгоритм коррекции использует параболическую сплайн-аппроксимацию дискретных значений сигнала ИП, интервал дискретизации составляет 0,01 с. Скорректированный сигнал г{(), так же как и в вышеописанном примере, запаздывает на два дискретных интервала относительно входного контролируемого параметра х(0, его график для удобства сравнения на рис. 3 сдвинут влево на величину этого запаздывания. Из рассмотрения диаграмм видно, что динамическая погрешность инерционного ИП, в данном случае составляющая десятки процентов, с использованием описанного алгоритма коррекции снижается до нескольких процентов.
Необходимо иметь в виду, что при реализации описанного алгоритма в связи с использованием в нем процедуры дифференцирования при наличии периодической помехи на выходе схемы коррекции могут возникать шумы, интенсивность которых зависит от частотной полосы шумов. Однако в том случае, когда период дискретизации сигнала кратен периоду помехи, шумы не возникают, что характерно для большинства алгоритмов цифровой фильтрации.
1 3 /
2
их,
0.
0.
0
0 4 8 12 16 20 24 28 32 36 п
Рис. 3. Коррекция сигнала инерционного ИП, представляющего собой сумму двух гауссовых функций
Метод коррекции динамической погрешности инерционного ИП показывает хорошие результаты в случае, когда достаточно точно известна передаточная функция детектора, и если передаточная функция представляет собой апериодическое звено первого порядка, то динамическая погрешность снижается в значительной степени. В принципе описанный подход коррекции динамической погрешности может быть применен и в том случае, когда передаточная функция детектора соответствует апериодическому или колебательному звену второго порядка. Метод коррекции, описываемый выражениями (3) - (7), достаточно прост и может быть легко реализован современными микропроцессорными системами.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Солопченко Г.Н. Обратные задачи в измерительных процедурах. Измерение, контроль, автоматизация / ЦНИИТЭИприборостроения. - 1983. - Вып. 2(46). - С. 32-46.
2. Ланге П.К. Сплайн-аппроксимация дискретных значений сигналов с применением методов цифровой фильтрации // Сб. труд. Самарского гос. тех. ун-та. Сер. Физ.-матем. науки. -2003. - Вып. 19. - С. 134-138.
3. Ланге П.К. Дифференцирующие цифровые фильтры со сплайн-аппроксимацией производной дискретизированных сигналов // Сб. труд. Самарского гос. тех. ун-та. Сер. Технические науки. - 2004. - Вып. 24. - С. 42-49.
4. Платонов И.А., Ланге П.К., Колесниченко И.Н., Платонов В.И. Динамические характеристики микродетектора теплопроводности для газоаналитических приборов // Измерительная техника. - 2015. - № 6. - С. 71-73.
Статья поступила в редакцию 1 февраля 2017 г.
CORRECTION OF THE DYNAMIC ERROR OF INERTIAL SENSORS
P.K. Lange
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
The method of dynamic error correction for inertia sensor having a transfer function type aperiodic operator of the first order is considered. The method is based on the realization of the inverse sensors operator in a dynamic mode. Correction method uses a parabolic spline - approximation of discrete values of sensors signal, as well as its first derivative. An example of correction of the signal formed by the inertia sensor of a particular type is described. It is shown that the use of a mathematical correction method based on digital filter in-line processing of discrete signal values allows to reduce the dynamic error of inertia sensor in several times. It is proposed to use the method described above also for the correction of the sensors signal, the operator function of which has the second order. The considered method can be easily implemented using modern microprocessor measuring systems.
Keywords: sensor, inertia, frequency response, transfer function, measured parameter, dynamic error, spline-approximation.
Petr K. Lange (Dr. Sci. (Techn.)), Professor.