Коррекция динамической погрешности инерционных детекторов газоаналитических приборов

Размещено на http://www.allbest.ru/

УДК 621.372.542

КОРРЕКЦИЯ ДИНАМИЧЕСКОЙ ПОГРЕШНОСТИ ИНЕРЦИОННЫХ ДЕТЕКТОРОВ ГАЗОАНАЛИТИЧЕСКИХ ПРИБОРОВ

И.А. Платонов1, В.М.Мухин1, П.К.Ланге2

1Самарский государственный аэрокосмический университет имени академика С.П. Королёва

(национальный исследовательский университет)

2Самарский государственный технический университет

Рассмотрен метод коррекции динамической погрешности инерционного детектора по теплопроводности с передаточной функцией, имеющей вид апериодического звена первого порядка. Метод основан на реализации обратной функции оператора, описывающего инерционный детектор в динамическом режиме. Метод коррекции использует параболическую сплайн - аппроксимацию дискретных значений сигнала детектора, а также его первой производной. Приведен пример коррекции сигнала детектора по теплопроводности конкретного типа. Показано, что использование математического метода коррекции на базе цифровых фильтров скользящей обработки дискретных значений сигнала позволяет в несколько раз снизить динамическую погрешность инерционного детектора. Предложено использовать описанный метод также и для коррекции сигнала детектора, оператор которого соответствует инерционному звену второго порядка. Рассмотренный метод может быть легко реализован с использованием современных микропроцессорных измерительных систем.

Ключевые слова: детектор по теплопроводности, газовый поток, передаточная функция, динамическая погрешность, сплайн - аппроксимация

инерционный детектор теплопроводность передаточный

The method of dynamic error correction for thermal conductivity inertia detector with a transfer function having the form of the first order aperiodic element. The method is based on the realization of the inverse function operator describing inertia detector in a dynamic mode. Correction method uses a parabolic spline - approximation of discrete signal detector values, as well as its first derivative. An example of correction of the signal formed by the thermal conductivity detector 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 in several times the dynamic error of inertia detector. It is proposed to use the method described above also for the correction of the signal detector, the operator of which having the form of the second order element. The considered method can be easily implemented using modern microprocessor measuring systems.

Key words: A thermal conductivity detector, gas stream, transfer function, dynamic error, spline - approximation

Введение

Известные детекторы газоаналитических приборов [1, 2] обладают повышенной инерционностью, которая затрудняет их использование в системах экспресс - анализа содержания примесей в газовой среде. Большая инерционность этих детекторов вызвана, с одной стороны, значительной величиной объема камеры детектора, с другой - длительностью процесса установления теплового равновесия в газовом потоке между термочувствительным элементом и стенкой камеры детектора. Этот фактор вызывает появление динамической погрешности определения концентрации примеси в потоке газа.

Целью работы является разработка математического метода коррекции указанной динамической погрешности.

Методика исследований

При динамических измерениях интерес представляет не выходной сигнал детектора y(t), а контролируемый входной параметр детектора x(t). Поэтому задачей коррекции сигнала y(t) является определение значений параметра x(t) с учетом оператора А, характеризующего динамические свойства детектора и его аппаратную функцию. В современных измерительных системах доступны дискретные значения сигнала y(t), формируемые аналого - цифровым преобразователем системы.

Такая задача может быть решена реализацией оператора А-1, обратного оператору А при соответствующей обработке дискретных значений сигнала y(t) [3].

В наиболее распространенном случае инерционный детектор имеет передаточную функцию апериодического звена первого порядка

, (1)

где Т - постоянная времени детектора.

При К=1 передаточная функция корректирующей цепи принимает вид

,

что соответствует реализуемому дифференциальному уравнению

.

Таким образом, это метод коррекции динамической погрешности должен реализовать функцию дифференцирования сигнала y(t) и сложения его производной с самим сигналом.

Для решения такой задачи в [3] было предложено использовать цифровой фильтр, реализующий параболическую сплайн - аппроксимацию дискретных значений сигнала y(t), а также его производной.

При использовании параболической сплайн - аппроксимации на n-м дискретном участке сигнал описывается параболической функцией

где a2 , a1 , a0 - коэффициенты аппроксимирующей функции.

Коэффициенты такой аппроксимации могут быть определены с использованием цифровых фильтров "скользящей" обработки дискретных значений сигнала y(t) в реальном темпе времени, число точек весовых функций фильтров обычно лежит в пределах 4…6.

Коэффициенты цифровых фильтров пятиточечной сплайн - аппроксимации сигнала y(t) по его дискретным значениям, определяются выражениями [4]:

,

, (2)

.

Коэффициенты пятиточечной параболической сплайн - функции, аппроксимирующей производную сигнала y(t), определяются выражениями [4]:

,

, (3)

.

Таким образом, скорректированный сигнал определяется выражением

, (4)

где

,

, (5)

.

При t=0 из (4) следует:

(6)

Как видно из (3), скорректированный сигнал z(t) формируется относительно дискретных значений сигнала y(n) детектора с задержкой в 2 дискретных интервала.

Профиль скорректированного сигнала (4) в идеальном случае должен соответствовать профилю входного параметра x(t).

Использование сплайн - аппроксимации позволяет определить значения сигнала z(t) внутри интервалов дискретизации, что является достоинством рассматриваемого метода.

Рассмотрим в качестве примера коррекцию сигнала микродетектора по теплопроводности, широко распространенного в расходомерах микропотоков газов, в газоаналитических приборах, в хроматографических анализаторах.

Переходная характеристика одного из таких детекторов [5] изображена на рис. 1. В соответствии со стандартной методикой [6] определения инерционности детектора по его переходной характеристике постоянная времени Т детектора в выражении (1) оказалась равной 0,1с.

Рис. 1. Диаграмма выходного сигнала UY микродетектора по теплопроводности при скачкообразном входном параметре x(t) = 1(t).

Известно [6,7], что концентрация введенной в газовый поток примеси в виде импульса распределяется в соответствии с функцией, близкой к Гауссовой. Это объясняется протеканием ряда физических процессов, основным из которых является процесс диффузии. При протекании газового потока через камеру детектора по теплопроводности детектор формирует сигнал y(t) в соответствии с величиной мгновенной концентрации примеси в потоке, однако инерционность детектора вызывает появление динамической погрешности.

Рассмотрим коррекцию выходного сигнала такого детектора при его входном параметре, изменяющемся в соответствии с Гауссовой функцией единичной высоты (рис. 2)

. (7)

Ширина такой функции составляет примерно 0,05с.

При передаточной функции (1) детектора

,

его переходная характеристика определяется выражением .

Выражение для сигнала y(t) на выходе ИП может быть получено с использованием интеграла Дюамеля:

Для конкретного значения постоянной времени Т=0,1с передаточной функции (1) детектора график этой функции изображен на рис. 2. Как видно из рассмотрения графика, такой детектор характеризуется очень большой динамической погрешностью.

Рис. 2. Сигналы на входе и выходе инерционного детектора.

х - график изменения параметра на входе детектора,

y - сигнал на выходе детектора.

Функция скорректированного сигнала z(t) определяется при аппроксимации дискретных значений сигнала y(t) параболической сплайн - функцией с использованием выражений (2) - (3) для коэффициентов аппроксимации, а также выражений (4) - (5) для метода коррекции.

Эти функции изображены на рис. 3, при рассмотрении которого видно, что рассмотренный метод достаточно эффективно восстанавливает по форме сигнал x(t), действующий на входе инерционного детектора, с запаздыванием в 2 интервала дискретизации.

Рис. 3. Диаграммы сигналов.

1 - параметр x(t) на входе инерционного детектора,

2 - сигнал y(t) на выходе детектора,

3 - скорректированный сигнал z(t-2), сдвинутый на 2 дискретных интервала влево,

n - номера дискретных отсчетов сигнала y(t).

Результаты

При наличии в потоке газа нескольких примесей, размещенных в потоке на некотором расстоянии друг от друга, детектор регистрирует последовательно два пика, форма которых близка к Гауссовой функции. Такая ситуация характерна для различных хроматографических и ряда других анализаторов. В этом случае инерционность детектора может вызвать существенную динамическую погрешность.

На рис. 5 представлены диаграммы сигналов: входного параметра x(t) микродетектора по теплопроводности, сигнала y(t), формируемого на выходе детектора, скорректированного сигнала z(t) с использованием описанного алгоритма коррекции. Входной параметр x(t) в данном случае представляет собой сумму двух Гауссовых функций вида (7), смещенных друг относительно друга на интервал 0,02с. Ширина Гауссовых функций составляет примерно 0,06с. Алгоритм коррекции использует параболическую сплайн - аппроксимацию дискретных значений сигнала детектора, интервал дискретизации составляет 0,01 с. Из рассмотрения диаграмм видно, что динамическая погрешность, в данном случае составляющая десятки процентов, с использованием описанного алгоритма коррекции снижается до нескольких процентов.

Рис. 4. Коррекция сигнала детектора, представляющего собой две смещенные Гауссовые функции.

Заключение

Метод коррекции динамической погрешности инерционного детектора по теплопроводности показывает хорошие результаты в случае, когда достаточно точно известна передаточная функция детектора. В том случае, когда передаточная функция представляет собой апериодическое звено первого порядка, динамическая погрешность снижается в значительной степени. Описанный подход может быть применен и в том случае, когда передаточная функция детектора соответствует апериодическому или колебательному звену второго порядка. Метод коррекции, описываемый выражениями (2) - (6), достаточно прост и может быть легко реализован современными микропроцессорными системами.

Библиографический перечень

1. Высокоэффективная газовая хроматография / Под ред. К.Хайвер, пер. с англ. - М.: Мир. - 1993. - 288 с.

2. Berezkin V.G. Chromatographic Adsorption Analysis: Selected Works of M.S. Tswett, Ellis Horwood, New York, 1990, see: L.S. Ettre, Milestones in Chromatography, LC GS North America. - 2003. - р.p. 458-468.

3. Коррекция динамической погрешности измерительных преобразователей на основе сплайн - аппроксимации сигнала / П.К.Ланге // Известия Самарского научного центра РАН. 2003. №2. Т.5. С.162-168.

4. Сплайн - аппроксимация дискретных значений сигналов с применением методов цифровой фильтрации / П.К. Ланге // Сб. труд. Самарского гос. тех. ун-та. Серия "Физ.-матем. науки". 2003. Вып.19. С. 134-138.

5. Динамические характеристики микродетектора теплопроводности для газоаналитических приборов / И.А. Платонов, П.К.Ланге, И.Н.Колесниченко, В.И.Платонов // Измерительная техника. 2015. №6. С. 71-73.

6. Бражников В.В. Детекторы для хроматографии. М.: Машиностроение. 1992. 320 с.

7. Е. Лейбнитц, Х. Штруппе; [Под ред. А.А. Жуховицкого]. Руководство по газовой хроматографии. М.: Мир, 1969. 503с.

Размещено на Allbest.ru