Модель для оптимизации настроечных коэффициентов программно-аппаратного измерительного канала Текст научной статьи по специальности «Компьютерные и информационные науки»

Модель для оптимизации настроечных коэффициентов программно-аппаратного измерительного канала

В.П. Шевчук, Е.Ю. Силаева

Волжский политехнический институт (филиал) Волгоградского государственного

технического университета

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

Введение. При разработке крупных информационно-измерительных и управляющих систем особое внимание необходимо уделить взаимодействию различных блоков системы друг с другом. Для этого применяется аппаратно-программное моделирование, при этом блоки реализуются в виде моделей [1, 2].

Структура программно-аппаратного измерительного канала газоанализатора, показанная на рисунке 1, представляет собой последовательное соединение четырех преобразователей:

1) динамическое преобразование аналоговой измерительной информации (первичный преобразователь);

2) преобразователь формы сигнала, осуществляющий процесс измерения (аналого-цифровой преобразователь);

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

4) цифро-аналоговое преобразование [3].

На рисунке 1 приняты следующие обозначения: Х(1) - измеряемая физическая величина; У(1) - выходная величина программно-аппаратного

измерительного канала; Ир(¿) - эквивалентная весовая функция аналоговой части измерительного канала газоанализатора; МВА - модуль ввода аналогового сигнала; к,(1Т§) - эквивалентная весовая функция программной части измерительного канала газоанализатора; кг(1) - весовая функция восстанавливающего элемента; ЕС1(1-Т^) - ошибка квантования по уровню в аналого-цифровом преобразователе; Т$ - период опроса газоанализатора, с; г - дискретный момент времени.

Рис. 1. - Структура программно-аппаратного измерительного канала Постановка задачи. Анализ рис. 1 показывает, что динамическая характеристика программно-аппаратного измерительного канала газоанализатора, кп(1), может быть представлена сверткой четырех динамических характеристик вышеупомянутых элементов измерительного канала:

кп (г ■ Тс ) = Ир (г ■ Тс ■ Тс - г ■ Т5 )• к, (г ■ Т5 ) • кг (г ■ Т5 ), (1)

где Тс - время измерения, 5(1Тс - 1-Т$)- динамическая характеристика аналого-цифрового преобразования (функция окно - модель процесса измерения).

Весовая функция большинства первичных преобразователей имеет вид: к —

кр (() = ^ ■ еТР, (2)

р

где Тр - постоянная времени; кр - коэффициент передачи.

:

Математическую модель периодического отбора целесообразно представить в виде стационарной случайной периодизированной последовательности [4, 5]:

Ы-Х ■ ), (3)

где - знак последовательности; N - количество измерений на к-

той циклограмме.

Имитационное моделирование программно-аппаратного измерительного канала.

Если предположить, что ошибка от квантования по уровню является аддитивно наложенной коррелированной помехой с уровнем шума £=0.5: Еч (■ Т ) = гпё(кр (■ Т5 )-0.01)~ Ир (г ■ Т5 )■ 0.005), и автокорреляционной функцией

экспоненциального вида: ЯП(¡■Т5)= 0.241 ■е~2'5'гТ , то сигнал на выходе модуля ввода аналоговой информации можно описать математической моделью: Хя (гТ8 ) = кр (Т ) + Еч (Т ).

Компенсацию динамических свойств первичного преобразователя целесообразно осуществлять методом экспоненциального сглаживания. Алгоритм экспоненциального сглаживания представляет собой рекуррентное соотношение, определяющее текущее значение программного динамического преобразователя У(■Т8) через текущее значение измеряемой физической величины и предыдущего значения программного преобразователя:

У (Т ) = У ((г - 1).Т8) + ^ [[ (г-Т8 )-У ( - 1)Т )],

Т 1/

где Т^- - постоянная времени интегрирования программного преобразователя. При этом динамический компенсатор имеет весовую функцию вида [4]:

К ) =

Т

Т 1/

( Т \

1 - ^ Т

,при г = 1,2,...,N (4)

0,при г = 0

Далее сигнал в виде решётчатой функции У(гТ$ ) подаётся на вторичный прибор для визуализации (монитор, планшет, дисплей компьютера и т.п.) весовая функция которого имеет вид [4]:

Применение восстанавливающего элемента нулевого порядка (5) (элемент запоминания) позволяет преобразовать решетчатую функцию У(]'Т§) в ступенчатую функцию выходного сигнала измерительного канала. [6].

При этом методическая составляющая погрешности обработки информации возникающая в программно-аппаратном измерительном канале газоанализатора вычисляется как: Л(г • Тс) = Ир (г • Тс) - У (г • Тс). Статистические характеристики погрешности обработки информации будет иметь следующие значения: оценка математического ожидания ошибки МЛ = -0.16, оценка дисперсии ошибки БЛ « 0 .23 (см. рис.2).

Рис 2. Погрешность программно-аппаратного измерительного канала газоанализатора (Тс = 1с, Кр = 100, Тр = 20 с, Тв = 3 с, Тр = 100, а=1/ТР=0.05 Гц).

Как видно, из анализа математических моделей элементов измерительного канала (2-5), программное обеспечение измерительного

(5)

100 I сел

канала имеет семь настроек: время измерения физической величины - Тс, диапазон дросселирования сигнала газоанализатора (коэффициент передачи умноженный на 100%) - КР, постоянная времени интегрирования первичного преобразователя - ТР, период опроса первичного преобразователя - Т8, постоянная времени интегрирования программного преобразователя - Т, скорость изменения физической величины - а=1/ТР, Гц, и длина реализации -N=50'Т8 [3]. Из них, четыре отражают свойства измеряемой физической величины и аппаратной части измерительного канала (а, ТР, КР, и Тс), а три другие (Т8, Т, и N являются настройками программного обеспечения системы измерительного канала.

Таким образом, используя модель формирования динамической погрешности процессорного измерительного канала, подробно описанную в работах [3, 4, 7, 8], можно получить модель для оптимизации настроек программного обеспечения, которая будет иметь вид:

Тг + е~г'/Тр

Т

2 ■ Т1 - Т8

8

а ■ ТР

у2■ Тр ■ Вх 1-(а■ Тр)

2

Р/ у

■ Т - Т8)

Т, - е

Т /Т

+

Т/ + е

-аг.

(Тг - Т,)

■ (Тг - Т5) 1 - (а ■ Тр)) Тг - еа ■ Т - Т8)

-шт\к02( 0),К'Т2( 0)]+1

вх

при ограничениях вида:

Л Т ■В

Кк(0) = - Т В

1 - а ■Тр

-аТЛ

2 ■ а ■Тр

-ТЛ/Т,,

-аТ,

Т /Т

к / - '(Т/ -Г8) 1 + а■Тр г/ -е ■ (Т, -Т5);

ш1п[к°Т2(0),К1Т2(0)\-

\К0Ш(0), если К°Т2(0) < К\2(0), К\2(0), если К°Г2(0) > К'Т2(0),

Ту > Г > 1; Тя < Тг < 10ТЯ;

,(6)

где Вх - оценка дисперсии измеряемой физической величины; а - скорость изменения физической величины; 8 - уровень помехи; Кух(1Т8) - оценка

1

1

2

е

е

взаимокорреляционной функции истинных значений измеряемой величины Y[iTS] и показаний программно-аппаратного Z(iTS) измерительного канала.

Предлагаемая математическая модель позволяет найти настроечные коэффициенты программного обеспечения измерительного канала, которые обеспечивают минимальные потери информации в программно-аппаратном измерительном канале, у таких вторичных приборов, как газоанализаторы. Как видно, из модели (6), относительная погрешность измерения и визуализации зависит как от параметров, которые отражают свойства измеряемой физической величины и аппаратной части измерительного канала (а, S, TP, КР, и TC), так и параметров, которые являются настройками программного обеспечения измерительного канала (Ts, Tf, и N ).

Заключение. Модель для оптимизации настроечных коэффициентов программно-аппаратного измерительного канала (6) имеет минимум и может быть использована в системах автоматического управления для текущей идентификации настроек регуляторов. Подобные решения по идентификации параметров приведены в работах [9, 10].

Литература:

1. Шайкин А. С., Шайкина Е. В. Применение комплекса полунатурного моделирования в процессе проектирования информационно-измерительных и управляющих систем // Инженерный вестник Дона, 2014, №1 URL: ivdon.ru/ru/magazine/archive/n1y2014/2248.

2. Шевчук В.П. Классификация информационно-измерительных систем по типу уравнения измерения для определения критерия управления // Метрология. 2008. №12. С. 3 - 16.

3. Данилов С.И. Параметрический синтез измерительных каналов в автоматизированной системе управления технологическим процессом: автореф. дис. ... канд. тех. наук: 05.11.16. Волгоград, 2000. 16 с.

4. Шевчук В.П. Моделирование метрологических характеристик интеллектуальных измерительных приборов и систем. М.: ФИЗМАТЛИТ, 2011. 320 с.

5. Шевчук В.П. Расчет динамических погрешностей интеллектуальных измерительных систем. М.: ФИЗМАТЛИТ, 2008. 288 с.

6. Шевчук В.П. Информативность и эффективность интеллектуальной измерительной техники // Метрология. 2012. № 1. С. 12-21.

7. Шевчук В.П., Авдеюк О.А. и др. Проблемы определения погрешности визуализации обобщенных критериев управления // Инженерный вестник Дона, 2017, №4 URL: ivdon.ru/ru/magazine/archive/n4y2017/4390.

8. Stroet, P., Measuring Phase and Delay Errors Accurately in I/Q Modulators. Linear Technology, 2005. 102. Date Views 18.11.2017 URL: cds.linear.com/docs/en/application-note/an102f.pdf.

9. Браганец С.А., Гольцов А.С., Савчиц А.В. Идентификация математической модели главного золотника для системы диагностики и адаптивного управления открытием направляющего аппарата // Инженерный вестник Дона, 2013, №4 URL: ivdon.ru/ru/magazine/archive/n4y2013/1906.

10. Astrom K.J. Advanced PID control. -ISA. Triangle Park, 2006. 446 p.

References:

1. Shaykin A. S., Shaykina E. V. Inzenernyj vestnik Dona (Rus), 2014, №1. URL: ivdon.ru/ru/magazine/archive/n1y2014/2248.

2. Shevchuk V.P. Metrologija. 2008. №12. pp. 3 - 16.

3. Danilov S.I. Parametricheskiy sintez izmeritel'nykh kanalov v avtomatizirovannoy sisteme upravleniya tekhnologicheskim protsessom [Parametric synthesis of measuring channels in an automated process control system]: avtoref. dis. ... kand. tekh. nauk: 05.11.16.. Volgograd, 2000. 16 p.

4. Shevchuk V.P. Modelirovanie metrologicheskikh kharakteristik intellektual'nykh izmeritel'nykh priborov i sistem [Modeling of metrological characteristics of intelligent measuring instruments and systems]. M.: FIZMATLIT, 2011. 320 p.

5. Shevchuk V.P. Raschet dinamicheskikh pogreshnostey intellektual'nykh izmeritel'nykh sistem [Calculation of dynamic errors in intelligent measuring systems]. M.: FIZMATLIT, 2008. 288 p.

6. Shevchuk V.P. Metrologija. 2012. №1. pp. 12-21.

7. Shevchuk V.P., Avdeyuk O.A. i dr. Inzenernyj vestnik Dona (Rus), 2017, №4. URL: ivdon.ru/ru/magazine/archive/n4y2017/4390.

8. Stroet, P., Measuring Phase and Delay Errors Accurately in I/Q Modulators. Linear Technology, 2005. 102. Date Views 18.11.2017 URL: cds.linear.com/docs/en/application-note/an102f.pdf.

9. Braganets S.A., Gol'tsov A.S., Savchits A.V. Inzenernyj vestnik Dona (Rus), 2013, №4. URL: ivdon.ru/ru/magazine/archive/n4y2013/1906.

10. Astrom K.J. Advanced PID control. ISA. Triangle Park, 2006. 446 p.