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

Раздел 6. Инженерная экология и смежные вопросы. Аппроксимация нелинейных статических характеристик оптико-абсорбционных газоанализаторов

д.т.н. проф. Рылов В.А.

Университет машиностроения л>а гу1ол>(с1)таИ. г и

Аннотация. Исследована возможность аппроксимации статических характеристик оптико-абсорбционных газоанализаторов дробно-степенной функцией.

Ключевые слова: линеаризация статических характеристик, оптико-абсорбционные газоанализаторы, аппроксимация дробно-степенной функцией. Одним из основных методов диагностики и регулировки двигателей автомобильного транспорта является контроль состава выхлопных газов. Для этих целей используются многоканальные автоматические газоанализаторы для определения основных компонентов: оксидов углерода (СО, СОг), суммы углеводородов (СпНт), оксидов азота (Ж)х), кислорода (Ог). Нормы содержания вредных примесей для автотранспорта устанавливаются международными и национальными стандартами и постоянно ужесточаются с целью улучшения экологической обстановки в мегаполисах. Для определения СО, СОг и СпНт в анализируемой смеси традиционно используется избирательный оптико-абсорбционный метод, являющийся одним из направлений молекулярной спектроскопии.

Метод основан на измерении поглощения излучения в области среднего ИК диапазона (1,5 - 10 мкм), где большинство газов имеет характерные полосы поглощения. Для обеспечения необходимой избирательности приметаются интерференционные оптические фильтры, которые пропускают ИК излучение только в узкой области, совпадающей со спектральной полосой поглощения выбранного компонента. Оптический блок газоанализатора содержит малоинерционный модулируемый источник ИК излучения, кювету с анализируемой пробой и детекторы излучения с интерференционными фильтрами.

К особенностям такого прибора следует отнести то, что для измерения трех различных компонентов применяется общая кювета заданной длины, что создает определенные трудности при выборе оптимальной оптической плотности для каждого компонента. Это связано с тем, что каждый из них отличается от других как диапазоном измерения (СпНт - до 1%; СО -до 5%; СОг - до 20%), так и интенсивностью полосы поглощения. Следствием этого являются нелинейные статические характеристики для отдельных компонентов. Пример нелинейной статической характеристики для измерительного преобразователя на СОг представлен на рисунке 1

Рисунок 1 - Статическая характеристика измерительного оптико-абсорбционного преобразователя на СО2 (верхняя кривая) и критерий нелинейности (нижняя кривая) в зависимости от оптической толщины слоя анализируемой пробы

По горизонтальной оси на графике отложена оптическая толщина газового слоя - произведение концентрации определяемого компонента (%) на длину оптической кюветы (мм); статические характеристики будут одинаковы для одного преобразователя СОг с диапазоном измеряемых концентраций 20 % и длиной кюветы 10 мм и для другого - с диапазоном 100 % и длиной кюветы 2 мм.

По вертикальной оси отложена величина относительного поглощения у, равная отношению уменьшения выходного сигнала II к начальному сигналу при С = 0:

у

Щ 0) 1 ;

Для количественной оценки нелинейности статической характеристики используется критерий нелинейности, который равен отношению значений производной в начале и конце характеристики:

У 1с=о

Ne = —^.

У 1с=с,

шах

Обычно не применяют нелинейные статические характеристики, для которых критерий превышает величину 10-15; для этого используют более короткие газовые кюветы; при более высоких значениях критерия существенно возрастают погрешности в конце диапазона.

С применением микропроцессоров в аналитических приборах одной из функций микропроцессора стала линеаризация статической характеристики. Общую схему преобразования можно представить в виде нескольких последовательных стадий (блоков):

• получение электрического аналогового сигнала на выходе первичного измерительного преобразователя С/(С);

• усиление аналогового сигнала и преобразование его в цифровой сигнал у(С);

• нелинейное преобразование цифрового сигнала (линеаризация); ЛХу);

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

Экспериментальные данные представляют значения выходного сигнала при различных концентрациях диоксида углерода в смеси

Таблица 1

Концентрация, % Выходной сигнал Отн. изм. сигн.

СО 0,00 U0 2618 У 0,0000

С1 1,82 U1 2429 у! 0,0722

С2 4,92 U2 2231 у2 0,1478

сз 10,00 из 2030 yi 0,2246

С4 15,00 U4 1878 у4 0,2826

С5 20,10 U5 1763 у5 0,3266

С6 47,7 (?) U6 1425 уб 0,4557

С7 100 U7 1142 у! 0,5638

На рисунке 2а представлены результаты аппроксимации функции линеаризации ('(>') по стандартному алгоритму полиномиальной аппроксимации Excel. Функция линеаризации С(у) описывается полиномом 4 степени и содержит 5 параметров (коэффициентов):

С

Стандартный алгоритм аппроксимации Excel компенсирует погрешности (и ошибки) экспериментальных данных дополнительными членами полинома более высоких степеней; в результате чего получается очень высокий коэффициент корреляции (5 девяток). Однако эта функция аппроксимации не является «гладкой», о чем свидетельствует наличие членов полинома с отрицательными коэффициентами.

Другой алгоритм основан на аппроксимации экспериментальных данных «гладкими» функциями, в качестве которых могут быть использованы гиперболические степенные

функции типа:

С(У) = К2

1

<1-к3-у)г

-1

1 и

; у = I--; г - «степень» гиперболы

Un

Одно из основных свойств «гладких» функций заключается в том, что во всем диапазоне изменений аргумента (0, 1/Кз) их производные (любого порядка) не обращаются в нуль. Примерами таких «гладких» функций может служить гипербола:

1

С{у)=К2

\-К3-у

-1

1-К3-у

или экспонента

С

Интуитивно понятно, что нелинейная статическая характеристика реального измерительного преобразователя является «гладкой» функцией (как согнутая упругая линейка), а отклонения от номинальной гладкой функции связано с погрешностями эксперимента. Поэтому такая аппроксимация выявляет погрешности эксперимента: увеличение погрешности в одной экспериментальной точке в несколько раз по сравнению с остальными свидетельствует об ошибочном экспериментальном результате, требующем корректировки.

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

На рисунке 26 представлены результаты аппроксимации приведенных экспериментальных данных степенной гиперболической функцией (г = 6):

120 100 80 60 40 20

у = 0,166637х4  0,906340х3 + 2,926234х2 + 0,598157х +.

0,035732

R = 0,999992

С02 0 - 100 V.

11= Л /.

t г = зг у.

II з = -.53 у. Л -1 = - .21 X у 15= .13 у/ В 6 - ,09 у. !7 = - .32 у.

0.1 0.2 0.3 0.1 0.5

ГрадуироБочные данные и характеристика линеаризатора ССц)

а б

Рисунок 2 - Определение параметров функции линеаризации при аппроксимации полиномом 4 степени

Особенность данного эксперимента заключалась в том, что для всех экспериментальных точек, кроме С6, использовались ПГС с известной погрешностью аттестации; смесь С6 получалось на неаттестованой установке динамического смешения, поэтому значение концентрации возможно было оценить только с большой погрешностью. В результате обработки данных необходимо было определить 4 параметра: показатель степени г, масштабные коэффициенты К2, К3 и значение концентрации Св.

С учетом этого была принята следующая методика нахождения параметров линеаризующей функции:

1) задавалось значение показателя степени гиперболы г;

2) по методу наименьших квадратов определялись значения параметров Кз и значение СКО для данной аппроксимирующей кривой;

3) затем проводилась корректировка (подгонка) значения концентрации С'б, при котором значение СКО достигало минимального значения; это значение Сб принималось как оптимальное для данного значения показателя степени г. В результате «подгонки» при значении С'б = 47,7 % получилось минимальное значение СКО = 0,165 %, при этом значения параметров составили К2 = 1,132; Кз = 0,6845

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

В сводной таблице 2 приведены значения СКО и С 6 для функций линеаризации, лежащих в диапазоне от гиперболы до экспоненты. Минимальное значение СКО, определяющее выбор оптимальной функции (параметр г), находится на уровне 0,17 % и соответствует «гиперболам 5-7 степени»

г С 6 СКО % Къ

со 49,6 0,285 6,198

8 48,2 0,177 0,5611

7 48,0 0,168 0,6166

6 47,7 0,165 0,6845

5 47,4 0,170 0,7682

4 46,9 0,198 0,8751

3 46,2 0,259 1,0147

2 45,0 0,384 1,2035

1 42,7 0,630 1,462

Таблица 2

Экспонента

Гипербола

На основании проведенной обработки результатов для конкретного измерительного

преобразователя можно сделать следующие общие выводы:

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

• функции содержат два масштабных коэффициента и показатель степени;

• наличие одного параметра (К2) в виде сомножителя существенно упрощает алгоритм определения оптимальных параметров функции аппроксимации по методу наименьших квадратов;

• метод не является критичным по отношению к выбору показателя степени г , как следует из таблицы, при изменении параметра г в диапазоне от 4 до 8 СКО погрешности аппроксимации не превышает 0,2 %, что существенно меньше основной приведенной погрешности аналитических измерительных преобразователей (обычно 2-5 %).

Методы интенсификации процессов кристаллизации при криогранулировании жидкофазных суспензий в жидком азоте. Исследование окрестности фронта смены режимов кипения

к.т.н. проф. Белуков C.B., Малышев P.E.

Университет машиностроения 8(499)267-07-14

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

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