способ идентификации линейного объекта

Формула изобретения

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

Описание изобретения к патенту

Изобретение относится к технической кибернетике и предназначено для идентификации линейных динамических объектов с переменным запаздыванием. Способ может быть применен для определения (аппроксимации) времени переменного запаздывания на основе информации о сигналах объекта на входе и выходе и реализован с использованием ЭВМ в автоматическом режиме, в реальном масштабе времени.

Известна система идентификации параметров нестационарного объекта с переменным запаздыванием (патент РФ 2027214, МПК G 05 В 23/00, опубликованный 20.01.95).

Изобретение относится к системам идентификации переменных параметров динамических объектов, в том числе и переменного запаздывания, зависящего от времени. Способ определения переменного запаздывания основан на подстройке параметров системы по сигналу означающему наличие отклонений параметров объекта и модели. При параметры объекта и модели равны. Лежащее в основе изобретения представление x(t-) в ряд Тейлора в точке t обладает очень низкой точностью и существенно зависит от числа членов разложения, к тому же автор ограничивается двумя членами разложения. Следует отметить, ведение дифференциаторов еще более усугубляет положение, так как вносит в систему импульсные воздействия, что вряд ли может гарантировать высокое качество идентификации.

Наиболее близким к предлагаемому способу является способ идентификации линейного объекта (патент РФ 2146063, МПК G 05 В 17/02, опубликованный 27.02.2000) путем определения входного и выходного сигналов объекта, подачи их на идентификатор непрерывной дроби с последующим восстановлением дискретной математической модели объекта.

Способ заключается в следующем: по дискретным измерениям вход-выходных сигналов в равноотстоящие моменты времени с шагом дискретизации t строят дискретную передаточную функцию (ДПФ) как отношение Z-преобразования выходного сигнала к Z-преобразованию входного сигнала. Для получения ДПФ идентифицируемого объекта пользуются представлением в виде соответствующей правильной С-дроби, применив модифицированный алгоритм В.Висковатова, для которого строится идентифицирующая матрица. Если в (0)-строке матрицы конечное число k первых элементов равны нулю, а последующие отличны от нуля, то необходимо осуществить сдвиг влево на k элементов, при этом восстановленная ДПФ умножается на z^-k.

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

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

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

Реализация способа поясняется структурной схемой, которая содержит:

блок 1 изменения амплитуды входного сигнала;

блок 2 объекта идентификации;

блок 3 идентификатора непрерывной дроби;

блок 4 восстановления дискретной математической модели;

блок 5 перехода к непрерывной математической модели;

блок 6 оценки временного запаздывания;

блок 7 накопления оценок временного запаздывания;

блок 8 идентификации временного запаздывания.

Фиксированный входной сигнал x(k) из блока 1 поступает на вход блока 2 объекта идентификации и на вход блока 3 идентификатора непрерывной дроби. На блок 3 поступает также выходной сигнал y(k). Блок 3 преобразовывает значения входного и выходного сигналов объекта в идентифицирующую матрицу и формирует непрерывную дробь. Затем коэффициенты непрерывной дроби поступают на вход блока 4 восстановления дискретной математической модели объекта. Далее параметры полученной модели поступают на блок 5 перехода к непрерывной математической модели реакции объекта, в котором определяют модельные значения выходного сигнала объекта идентификации. Затем в блоке 6 из множества дискретных значений, выходного сигнала y(k) определяют время t_k, при котором данный сигнал отличен от нуля. Задают шаг дискретизации меньше наперед заданной точности и путем обратного табулирования от момента времени t_k определяют величину временного запаздывания объекта. Далее величина временного запаздывания поступает в блок 7 накопления оценок транспортного запаздывания. Потом изменяют величину амплитуды входного сигнала и повторяют процесс определения величины временного запаздывания. Накопленные значения в блоке 7 поступают на вход блока 8 идентификации переменного запаздывания объекта, в котором его определяют с помощью аппроксимации дробями Тиле с требуемой погрешностью вычислений.

Предлагаемый способ осуществляется следующим образом.

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



где Y(z) Z - преобразование выходного сигнала;

X(z) Z - преобразование входного сигнала.

Для получения ДПФ идентифицируемого объекта пользуются представлением (1) в виде соответствующей правильной С-дроби, применив модифицированный алгоритм В. Висковатова. На основании (1) расчетным путем определяется следующая матрица:



в которой (-1)-строка и (0)-строка содержат дискетные значения входного x(kt) и дискретные значения выходного сигналов y(kt) в моменты времени {nt}^N₀, и они являются начальными условиями при построении матрицы, а элементы _m(nt) последовательно определяются с помощью соотношений:



где _-1(nt) = x(nt); ₀(nt) = y(nt), m=l,2,3,...; n=0,1,2,3,... .

Тогда элементы первого столбца матрицы (2) порождают частные числители правильной С-дроби, и ДПФ будет иметь вид:



При аппроксимации дробно-рациональной функции в матрице (2) наблюдается появление нулевой строки, номер которой позволяет идентифицировать порядок функции. Если в (0)-строке матрицы (2) конечные числа k первых элементов равны нулю, а последующие отличны от нуля, то необходимо осуществить сдвиг влево на k элементов и далее определить элементы матрицы (2) по правилу (3), при этом, восстанавливая правильную С-дробь (4), элементы ₀(0) умножаются на z^-k.

Существенным для данного правила является то, чтобы первый элемент (-1)-строки был отличен от нуля, то есть входное воздействие в начальный момент времени х(0)0. Получив ДПФ в виде (4) и свернув ее в дробно-рациональную функцию

,

переходят от данного выражения к непрерывной математической модели объекта с помощью взаимно однозначного соответствия



По ней однозначно определяется непрерывная реакция объекта в виде аналитической зависимости. Начиная с первого ненулевого значения реакции объекта, осуществляется обратное табулирование полученной зависимости, при котором определяется значение аргумента, соответствующего нулевому значению реакции. Зная момент начала действия входного сигнала и полученное значение времени, определяется время запаздывания с точностью, определяемой точностью табулирования.

Приведем пример, иллюстрирующий решение задачи SP-идентификации предлагаемым методом.

Пример 1

Известно, что модель, связывающая расход топлива с давлением тяги одновального турбинно-реактивного двигателя (ТРД), имеет вид



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



а ДПФ имеет вид



в которой нули и полюса

z₁ ⁿ=0.88207, z₂ ⁿ=0.77919, z₁ ^н=0.68757, z₂ ^н=-0.656668

В соответствии с аналогом временное запаздывание равно =kt=4.375 с, и для повышения точности оценки временного запаздывания необходимо уменьшить шаг дискретизации.

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

s₁ ⁿ=1.0038, s₂ ⁿ=1.996, s₁ ^н=2.9967,

и, следовательно, непрерывная передаточная функция (НПФ)



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



где



Для нахождения транспортного запаздывания зафиксируем первое ненулевое измерение переходной характеристики y(t)=0.07486 и момент времени t=4.375 с. Составим уравнение относительно неизвестной величины временного запаздывания , которое принимает вид



Полагая ⁰ = 4.375 c, итерационно уменьшаем значение = - до тех пор, пока не получим численное решение уравнения с заданной точностью .

При =0.001 транспортное запаздывание =4.229 с, а при =0.0001 =4.3 с.

Тем самым мы полностью восстановили непрерывную модель объекта.

Данный подход можно использовать для идентификации линейных динамических объектов с переменным запаздыванием. (Например, на практике время запаздывания зависит от величины управляющего воздействия).

Пример 2

Математическая модель изменения площади проходного сечения сопла ТРД с изменением давления тяги описывается данной передаточной функцией:



зависит от величины поперечного сечения сопла, ₀ = 3, k=2.

На вход подается ступенчатое воздействие x(t)=х1(t).

I. Пусть x₁=2, тогда (х)=3+1=4 и



t=0.5 c

Тогда идентифицирующая матрица принимает вид



Восстанавливая ДПФ по первому столбцу идентифицирующей матрицы, получим



Используя принцип вариации шага дискретизации при t=0.125 с, получим ДПФ:



Определим момент времени, при котором значение выходного сигнала отлично от нуля

⁰ = 4.125 c,

тогда получим НПФ



определим непрерывную реакцию объекта в виде



где













Согласно вышеприведенной методике определим величину транспортного запаздывания с заданной точностью =0.0001: ^м = 4c.

II. Пусть х₂=4, тогда период дискретизации оставляем без изменения t=0.125 с. Выполняя преобразования, аналогичные (I), получим

⁰ = 3.625 c,

^м = 3.5 c при =0,0001

Выполняя аналогичные выкладки, на следующих этапах получаем

III. x₃=6, (х)=3.333, t=0.125c,

⁰ = 3.375 c,

^м = 3.333 c при =0,0001

IV. x₃=8, (х)=3,25, t=0.125 с,

⁰ = 3.375 c,

^м = 3.250 c при =0,0001.

V. x₃=10, (х)=3,2, t=0.125 с,

⁰ = 3.25 c,

^м = 3.20 c, при =0,0001.

На следующем этапе, используя обратные разности (дробь Тиле), построим зависимость для определения переменного запаздывания.

Получаем следующие значения в интерполяционной таблице:



Для приведенной таблицы получим интерполяционную непрерывную дробь вида



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