CC BY
114
УДК 681.311
DOI 10.21685/2072-3040-2018-3-6
Н. П. Кривулин
ВОССТАНОВЛЕНИЕ ВХОДНЫХ СИГНАЛОВ НЕСТАЦИОНАРНЫХ ДИНАМИЧЕСКИХ СИСТЕМ1
Аннотация.
Актуальность и цели. Разработка аппаратных и программных методов восстановления входных сигналов нестационарных динамических систем и повышения их точности является актуальной проблемой, которая возникает во многих отраслях техники. Данная работа посвящена разработке аппаратных и программных методов восстановления входных сигналов нестационарных динамических систем, моделируемых дифференциальными и интегральными уравнениями. Предлагаемые в работе методы не зависят от физической природы измеряемой величины, что позволяет применять их во многих областях физики и техники и строить вычислительные алгоритмы восстановления входных сигналов измерительных преобразователей, выполненных на различных элементных базах.
Материалы и методы. Методы восстановления входных сигналов основаны на численных алгоритмах и результатах идентификации и параметрической идентификации. Предлагаемые методы позволяют осуществить коррекцию выходных сигналов динамических систем. Коррекция основана на построении обратного оператора к оператору, моделирующему измерительный преобразователь. Рассмотрена программная и аппаратная реализация обратного оператора.
Результаты. Разработанные методы восстановления входных сигналов осуществляют редукцию к идеальному прибору.
Выводы. Построены вычислительные алгоритмы, позволяющие восстанавливать с высокой точностью входные сигналы непрерывных динамических систем.
Ключевые слова: восстановление входных сигналов, редукция к идеальному прибору, идентификация динамических систем, импульсная переходная функция, вычислительные алгоритмы.
N. P. Krivulin
RECOVERING INPUT SIGNALS OF NON STATIONARY DYNAMICAL SYSTEMS
Abstract.
Background. The development of hardware and software methods for restoring input signals and increasing their accuracy is an actual problem that occurs in many branches of technology. This work is devoted to the development of methods for reconstructing input signals of non stationary dynamical systems. The methods proposed in the work do not depend on the physical nature of the measured value, which makes it possible to apply in many measurement areas and to build computa-
1 Работа поддержана РФФИ. Грант 16-01-00594.
© 2018 Кривулин Н. П. Данная статья доступна по условиям всемирной лицензии Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), которая дает разрешение на неограниченное использование, копирование на любые носители при условии указания авторства, источника и ссылки на лицензию Creative Commons, а также изменений, если таковые имеют место.
tional algorithms for restoring the input signals of measuring converters made on a different element base.
Materials and methods. Methods for restoring input signals are based on numerical algorithms and parametric identification. The proposed ones allow to correct the output signal of dynamic systems, consisting in the hardware or software implementation of the inverse operator to the operator simulating the measuring transducer.
Results. The methods of restoring input signals developed in this work allow us to construct algorithms that perform reduction to an ideal instrument.
Conclusions. The methods allow the correction of distortions by a computing device, which can be a specialized computing device made in the form of a microcontroller or a PC. Computational algorithms that allow to reconstruct the input signals of continuous systems with high accuracy are considered.
Key words: restoration of input signals, reduction to an ideal device, identification of dynamic systems, impulse response function, computational algorithms.
Введение
Проблема восстановления входных сигналов измерительных систем занимает одну из ключевых позиций во многих областях техники. Она возникает при динамических измерениях, в радиолокации, теории распознавания изображений, теории связи и т.д. Обзоры методов решения данной проблемы можно найти в работах [1-3].
Как отмечается в монографии [2], существует два подхода к решению данной проблемы. Первый, априорный, заключается в том, что путем совершенствования конструкции измерительных систем добиваются минимума искажений в системе приема сигнала; в основу второго, апостериорного, положены алгоритмы восстановления входного сигнала по выходному сигналу. Построение таких алгоритмов осуществляет редукцию к идеальному прибору.
В данной работе используется второй подход, основанный на алгоритмических методах идентификации динамических систем.
Задаче восстановления входных сигналов посвящено большое число работ (библиографию см. в [2, 4]), в которых в основном рассматриваются стационарные динамические системы, описываемые уравнениями вида
J g (t -т) x(T)d т) = f (t),
0
где x(t) - входной сигнал, подлежащий восстановлению; f (t) - известный выходной сигнал; g (t) - известная ипульсная переходная функция.
В данной работе рассматриваются методы восстановления входных сигналов нестационарных динамических систем, описываемых уравнениями вида
J g (t, т) x(t) d т = f (t),
0
где x(t) - входной сигнал, подлежащий восстановлению; f (t) - известный выходной сигнал; g(t, т) - импульсная переходная функция.
Предложено несколько подходов к восстановлению входных сигналов:
1) численный метод восстановления входного сигнала при известной импульсной переходной функции;
2) методы восстановления входного сигнала и динамических характеристик систем по результатам идентификации и параметрической идентификации.
Постановка задачи. Требуется восстановить исследуемые параметры входного сигнала в том виде, в котором они поступают на вход измерительной системы, т.е. скорректировать искажения, которые вносятся собственно измерительной системой.
1. Восстановление входных сигналов динамических систем по известным динамическим характеристикам
При известных динамических характеристиках систем одной из основных задач является восстановление входных сигналов. В данном разделе предложены методы восстановления входных сигналов по известной импульсной переходной функции.
В случае если a-priori импульсная переходная функция неизвестна, то она восстанавливается на основе методов, изложенных в работе [5].
Рассмотрим непрерывную динамическую систему, функционирование которой описывается уравнением
¡gЦ, т)х(т)dт = f ^), 0 < t <~ , (1)
0
где g ^, т) - известная импульсная переходная функция; х^) - входной сигнал, подлежащий восстановлению; / ^) - выходной сигнал.
Постановка задачи. Пусть g(^т) - известная импульсная переходная функция и известен выходной сигнал /^) системы (1). Требуется определить входной сигнал х^).
1.1. Численные методы восстановления входного сигнала
При решении прикладных задач вместо уравнения (1) приходится рассматривать уравнение вида
T
^ т) х(т^ т = f ^) + пЦ ), 0 < t < Г, (2)
0
где п(1 ) - шум; Г - момент времени такой, что при t > Г выполняется неравенство \/^) + п(})| <е, е>0 - заданная погрешность.
Замечание. Предполагается, что уровень заданной погрешности не меньше уровня внешнего шума : )| < е .
Обозначим через Е множество всевозможных выходных сигналов исследуемого измерительного преобразователя (ИП), определенных на сегменте [0,Г]. Обозначим через П множество случайных функций, характеризующих внешний шум.
Так как шум п(V) носит случайный характер, то при реализации вычислительных алгоритмов его целесообразно отнести к вычислительной погрешности задания правой части уравнения (2).
Обозначим через {фу(V)}, i = 1,2,..., множество линейно независимых функций, определенных и плотных в метрике С[0 у ]. Отметим, что это может
быть множество ортогональных полиномов, семейство сплайнов и т.д. Выделим конечную систему функций {фу(V)}, У = 1,2,...,Ы, такую, что любая функция / е ^ может быть аппроксимирована с точностью £ на сегменте ^ Т ].
Рассмотрим следующую систему интегральных уравнений:
Т
|g(V,т)х(тМт = фу(V), 0 < V <Т, У = 1,2,...,N, (3)
0
состоящую из N интегральных уравнений Фредгольма первого рода.
Известно [6], что решение уравнений Фредгольма первого рода является некорректной задачей. Однако в данном случае известны точные значения ядра уравнения и правой части и классические численные методы могут дать решение с необходимой точностью. Для решения системы уравнений (3), У = 1,2,..., N, можно использовать различные численные методы.
В частности, можно использовать метод механических квадратур. Вначале преобразуем систему уравнений (3):
, 2 -1. Проведем замену переменных: - =—1 +— V, тогда —1 <t < 1;
T
■ -/Л5*iТе ■•T
2
т = —1 +—т, тогда —1 <т < 1.
Т
2. Введем функции: g(-, т) 1 — Р g ^ + 1),~(т +1)
0 при - = ±1,
XI Т-(- +1)1 , фу(-) = фу ^ +1) I, У = 1,2,...,N .
2 11 при - *±1, ^2
x(t) =
л-р 2
Тогда система уравнений (3) примет вид 1
| т)х(ТМт = фу(-), —1 < -< 1, У = 1,2,...,N. (4)
—1
Приближенное решение системы уравнений (4) будем искать в виде функции хп (-) = 41 — -2 ип (-), где ип (-) - полином Чебышева второго рода
степени п на сегменте
[-1,1].
Подставив функции хп (¡) в уравнение (4), приравняв левые и правые
части полученного выражения в узлах ^ , к = 1,2,...,п, являющихся нулями
полинома Чебышева первого рода на сегменте [-1,1], и воспользовавшись
квадратурными формулами Гаусса наивысшей алгебраической степени точности, получаем систему метода механических квадратур:
п
£ ак§(([,Vk)ип ^к) = ф,-(¡1), I = 1,2,...,N, (5)
к=1
где {¡к}, {^к} , к = 1,2,...,N - узлы полиномов Чебышева второго и первого рода на сегменте [-1,1]; ак , к = 1,2,...,N, - коэффициенты квадратурной
формулы Гаусса с весом V1 - ¡2 .
Таблицы узлов и коэффициентов квадратурной формулы Гаусса можно найти в [7].
В случае если система уравнений (5) оказывается неустойчивой, то необходимо применить метод регуляризации. Запишем систему уравнений (5) в форме
СУ = 2 , (6)
где С = {с1к }, к, I = 1,2,...., п, сш = §(?г, V к (V к), к, I = 1,2,...., N
у ^аь..^ап}, 2 = {ф,(fl),...,ф,(1п)}.
Будем искать решение уравнения (6) в виде У = С*В, где С* -матрица, сопряженная с С, В = {¿1 ,...,Ьп } - искомый вектор.
Регуляризацию системы уравнений (6) будем проводить следующим образом.
Найдем собственные значения (,...,Хп) матрицы СС* . Пусть Х1 = 0,
X,- > 0, , = 2,3,...,п. Пусть у,- - собственный вектор матрицы СС*, отвечающий собственному значению X,-, , = 1,2,..., п.
Систему уравнений (6) разложим по собственным векторам у,, г = 1,2,...,п, и регуляризацию проведем сдвигом собственного значения Х1, заменяя Х1 на Х1 +Р1, где Р1 - достаточно малое положительное число.
Располагая решениями системы (4), решения системы (3) будут иметь
вид
х(г) = ^1 - [Тг -^Х(Тг -^.
Получив решения системы (3), вернемся к восстановлению входных сигналов ИП (2). Прежде всего проведем фильтрацию правой части уравнения (2) фильтром низких частот (настроенным таким образом, чтобы функции ф,■ (г), ■ = 1,2,...,п, не претерпели изменений). В результате получаем
функцию /(V) = /(V) + п^) . Затем функцию /(V) разлагаем в ряд по базисным функциям { ф^)}, у = 1,2,..., п :
п
/(V) = 2 У к Фк (V).
к=1
Зная решение ху (V) системы уравнений (3) для каждой функции фу (V), находим входной сигнал
п
X* ^) = 2 Укхк К).
к=1
Замечание. Предложенный метод позволяет:
1) избежать влияния шума при вычислениях;
2) создать обширную базу данных для ИП различной природы.
1.2. Восстановление входного сигнала по результатам параметрической идентификации
Во многих случаях для восстановления входного сигнала целесообразно представить модель в виде дифференциального уравнения. Это позволит дополнительно исследовать устойчивость и чувствительность модели. Аналогичный подход был реализован в работе [8] для дискретных систем, где математическая модель определялась в виде разностного уравнения с переменными коэффициентами.
Постановка задачи. Пусть известна импульсная переходная функция (ИПФ) g (V, т) динамической системы, описываемой уравнением (1).
Требуется построить математическую модель динамической системы (1) в виде дифференциального уравнения с переменными коэффициентами ап (V), Оп—), ..., а0^):
(V) /чЛп—1 / (V) .... . . .
Оп (V + Оп—1 (V)-п—Т- +... + О0^) / (V) = х^). (7)
Лп Лп 1
Здесь х^) - входной сигнал, подлежащий восстановлению; /(V) -выходной сигнал.
Отметим, что функция /(V), определяемая выражением (1) как выходной сигнал, является решением дифференциального уравнения (7).
Восстановление входного сигнала будем проводить по формуле (7). Для этого нужно определить коэффициенты ап (V), ап—1^),..., а0^) дифференциального уравнения (7) и найти производные выходного сигнала /(V):
Лп/(V) Лп—У(V) #(0
Следуя работе [9], восстановление коэффициентов ап(V), ап—1^), ..., О0(:) дифференциального уравнения (7) будем проводить, используя ИПФ g (V, т) динамической системы (1), которая предполагается известной.
Используя ИПФ § (г, т) системы (1), определим виртуальные входные сигналы %1(г), Х2(г), ...,хп+1(г) такие, чтобы им в качестве выходных сигналов соответствовала система линейно независимых функций
ф1(г), ф2(г),..., Фп+1(г). (8)
Условие линейной независимости системы функций (8) заключается в отличии от нуля вронскиана
W ((), Ф2 (t),..., )) =
) ф(и-1) (t) 92n)(t) 92n-1)(t)
Ф1 (t) 9i(t) 92(t) 92(t)
фП+^О фП+ЛО - фП+iCt) Фи+i (t)
n(«-i)
* 0. (9)
Определим искомые виртуальные входные сигналы %1(г), Х2(г), ..., хп+1 (г) из решения интегральных уравнений
jg(t,т)Xi(T)dт = ф,-(t), i = 1,2,...,n
+1.
(10)
Решая систему (10), находим для каждого выходного виртуального сигнала ф,(г), , = 1,2,..,п +1, системы (10) соответствующий виртуальный входной сигнал х, (г), , = 1,2,..., п +1.
Подставляя ф,(г) и х, (г), , = 1,2,...,п +1, в (7), приходим к системе уравнений относительно ^(г),а^),...,ап(г):
ап (г )ф(п) (г) + ап-1 (г )ф(п-1 )(г) +... + а0 ф, (г) = х, (г),, = 1,2,..., п +1. (11)
Так как определитель Вронского (9) системы линейных алгебраических уравнений (11) отличен от нуля, то решение системы (11) будет однозначным.
В частности, если в качестве линейно независимых функций (8) выступают функции ) = 1(г), У2(г) = г, ...,уп+1(г) = гп , то, как показано в работе [9], при г >0 коэффициенты дифференциального уравнения (7) определяются формулами
ац(г ) = х^г); а^г ) = х2(г) - а0(г )г;
... 5
п(г) =1 (п+1(0 -п(п -1)... • 2ап-1 (г)г -... -а0(г)гп). п! \ '
Таким образом, располагая функцией § (г, т), вычисляем коэффициенты уравнения (7).
Отметим, что для определения ап(V), ап—^), ..., а0^) не требуется подавать входные сигналы и фиксировать отклик системы на эти воздействия. Достаточно знать ИПФ системы.
Располагая функциями ак(V), к = 0,1,...,п , входной сигнал х^) восстанавливаем по формуле (7). Для этого необходимо вычислить производные от функции /(V), что, как известно, является некорректной задачей.
Для вычисления производных, следуя работам [10, 11], в ряде случаев может оказаться полезной вычислительная схема:
значения /(V).
Вычислительная схема (12) позволяет вычислять производные любого конечного порядка и имеет достаточно высокую точность и устойчивость. Параметром регуляризации в схеме (12) является к. Зависимость точности и устойчивости от параметра к исследована в [11], где предложены и другие численные методы вычисления производных, часть которых имеют большую точность, но меньшую устойчивость.
При практическом вычислении производных от функции /(V), Vе [0,Г], может быть использован следующий достаточно устойчивый метод. Пропустив функцию / (V) через фильтр низких частот, полученную функцию - (V) интерполируем по узлам Чебышева первого рода, трансформированным с сегмента [—1,1] на сегмент [0,Г ], и вычисляем производные от построенного полинома. В этом случае производные вычисляются аналитически.
Примеры, реализующие способ восстановления коэффициентов дифференциального уравнения с переменными коэффициентами, рассмотрены в работе [9].
Метод восстановления входного сигнала по аппаратной коррекции основан на обработке выходного сигнала динамических систем корректирующим устройством (КУ), на выходе которого будет наблюдаться входной сигнал измерительной системы. Для обработки выходного сигнала предлагается использовать корректирующее устройство, реализация которого может быть выполнена как аппаратно, так и программно.
Как и в случае определения ИПФ дискретных динамической системы [8], восстановление входного сигнала по коррекции выходного основано на следующем утверждении.
которая использует все измеренные в узлах
tk = Tk / N, k = 0,1,...,N,
2. Восстановление входного сигнала аппаратной коррекцией выходного сигнала
Рассмотрим интегральное уравнение
t
Jh(t,т)ф(т^т = y(t), 0< t <~ . (13)
0
Обозначим через Ф(p), Y(p) преобразования Лапласа функций 9(t) и y (t) соответственно:
Ф(р) = Z[9(t)] = J 9(t)е-ptdt, Y(p) = L[y(t)] = J y(t)e-ptdt. 00
Пусть преобразование Лапласа функции h(t, т) по переменной t удовлетворяет условию
H (p, т) = L[h(t, т)] = J h(t, т)е-ptdt = Н( p)e-Tq( p), (14)
0
где q(p) и HT(p) - аналитические функции в полуплоскости Re p >с (а = const) .
Тогда преобразование Лапласа уравнения (13) будет иметь вид
Н( ^Ф (q (p)) = Y (p). (15)
Доказательство данного утверждения и примеры восстановления ИПФ приведены в работе [12].
Постановка задачи. Пусть функционирование первичного измерительного преобразователя (ПИП) описывается уравнением
t
Jg (t, т) x(T)d т = f (t), 0 < t <~, (16)
0
где x(t) - входной сигнал ПИП, подлежащий восстановлению; f (t) - выходной сигнал ПИП; g (t, т) - ИПФ ПИП.
Требуется построить измерительную систему (ИС), состоящую из ПИП и КУ, причем функционирование последнего описывается уравнением
t
JgK (t, т) f (т)d т = x(t), 0 < t <~, (17)
0
где f (t) - входной сигнал КУ (выходной сигнал ПИП (16)); x(t) - выходной сигнал КУ (входной сигнал ПИП (16)); gK (t,т) - ИПФ КУ.
Структурная схема включения КУ в ИС приведена на рис. 1. Представим два метода построения импульсной переходной функции КУ. Первый основан на проведении серии тестовых испытаний при неизвестной импульсной переходной функции ПИП. При этом определяются как импульсная переходная функция КУ, так и импульсная переходная функция ПИП.
Рис. 1. Структурная схема ИС с аппаратной коррекцией: х(7) - входной сигнал ИС ПИП; х(7) - выходной сигнал КУ и ИС; g(7, т) - ИПФ ПИП; g (7, т) - ИПФ КУ
При втором методе восстановления импульсной переходной функции КУ предполагается, что известна импульсная переходная функция ПИП. При этом не требуется проведения дополнительных испытаний.
2.1. Восстановление импульсной переходной функции корректирующего устройства по серии тестовых испытаний
Данный метод позволяет определить одновременно: g(7,т) - ИПФ ПИП и импульсную переходную функцию gк (7,т) КУ.
Восстановление функций g (7, т), gк (7, т) будем проводить по серии из двух эталонных входных сигналов хЭ (7), 7 = 1,2, и соответсвующих
Э
выходных эталонных сигналов (7), 7 = 1,2, наблюдаемых на ПИП (16). Тогда для ИПФ g (7, т) ПИП на основании (16) получим соотношения
7
¡g (7, т) хЭ (т^т = (7), 7 = 1,2, 0 < 7 <~. (18)
0
Для ИПФ gк (7, т) КУ на основании (17) получим соотношения
г
^К (7, т)У7Э (т^т = хЭ (7), 7 = 1,2, 0 < 7 <~. (19)
0
Э
Обозначим преобразования Лапласа входных сигналов х^ (7), 7 = 1,2, и
Э
выходных сигналов (7), 7 = 1,2, в виде
ХЭ(р) = ДхЭ(7)] = \хЭ(7)в~ptd7; К?(р) = ЩЭ(7)] = \/Э(7)е-р^7, 7 = 1,2.
00
Пусть преобразования Лапласа ИПФ g ( 7, т) - ПИП (16) и gк ( 7, т) - КУ (17) по переменной 7 удовлетворяют условию (14):
G (p, т) = L[g (t, т)] = J g (t, ^e-ptdt = G( p)e-тq( p), 0
Gk (p, т) = L[ gK (t, т)] = J gK (t, т)e-ptdt = Gk (p)e-%qK (p),
0
где q(p), qK (p) и G(p), Gk (p) - подлежащие определению аналитические функции в полуплоскости Rep >а (а = const) .
Тогда преобразование Лапласа системы уравнений (18) для ПИП (16) на основании (13)-(15) примет вид
G(p)X? (q(p)) = F? (p), i = 1,2, (20)
где q(p) и G(p) - искомые аналитические функции в полуплоскости Rep > а(а = const).
Аналогично преобразование Лапласа системы уравнений (19) для КУ (17) будет иметь вид
Gk (p)F3 (qK (p)) = X? (p), i = 1,2, (21)
где qK (p) и Gk (p) - искомые аналитические функции в полуплоскости Rep >а (а = const).
Рещая систему уравнений (20) относительно функций q(p) и G(p), определим преобразование Лапласа ИПФ g (t, т) ПИП в виде
G (p, т) = G (p )e—%q(p).
Значения ИПФ g (t, т) найдем как обратное преобразование Лапласа по формуле обращения
а0
g(t,т) = i J ept-4(p)G(p)dp.
а0 —
Рещая систему уравнений (21) относительно функций qK (p) и Gk (p), определим изображение Лапласа ИПФ gK (t, т) КУ в виде
Gk (p, т) = Gk (p)e—тqк (p).
Значения ИПФ gK (t, т) найдем как обратное преобразование Лапласа по формуле обращения
а0 +i~
gK (t,т) = 2П7 J e^(p)Gk (p)dp.
а0 —
Тогда восстановление входного сигнала x(t) в режиме работы ИП определяется по формуле (17) в виде
г
х(г) = ¡£к (г, т)/(т^т.
о
Примеры восстановления импульсных переходных функций по серии тестовых испытаний рассмотрены в работах [1, 5, 12]
2.2. Определение ИПФ КУ по ИПФ ПИП
Рассмотрим ИС, состоящую из ПИП с ИПФ £ (г, т) и КУ, структурная схема которой приведена на рис. 1.
Требуется, располагая ИПФ £ (г, т) ПИП, определить ИПФ £ к (г, т) КУ. Пусть преобразование Лапласа ИПФ £ (г, т) ПИП по переменной г удовлетворяет условию (14):
О (р, т) = Ь[£ (г, т)] = ¡ £ (г, т)в-РгЛ = О (р)е-тд( р), (22)
о
где О(р), д(р) - некоторые известные (или восстановленные описанным выше методом) аналитические функции в области.
Пусть ИПФ £к (г, т) КУ удовлетворяет условию (14):
^
О (Р, т) = Ь[ £к (г, т)] = ¡ £к (г, т)е-ргЛ = Ок (р )е-тдк (р), (23)
о
где Ок (р), дк (Р) - некоторые подлежащие определению аналитические функции в области Яе р >а.
Обозначим преобразование Лапласа входного сигнала ПИП через X (р), преобразование Лапласа выходного сигнала ПИП - через Е (р).
Запишем уравнения соответствия входных и выходных сигналов для ПИП и КУ аналогично уравнениям (18), (19):
г г
¡£ (г, т)х(т)^т = /(г), ¡£к (г, т)/(т^т = х(г). (24)
оо
Применяя преобразование Лапласа к уравнениям (24), при условиях (22), (23), на основании (13)-(15) получим систему нелинейных алгебраическмх уравнений относительно Ок (р), дк (р):
О (р) X (д( р)) = Е (р), Ок (р) Е (дк (р)) = X (р). Подставляя Е (р) из первого уравнения во второе, получим
Ок (р)О (дк (р)) X (д(дк (р)) = X (р). (25)
Среди решений уравнения (25) имеется следующее: д(дк (р)) = р,
Ок (р)О(дк (р)) = 1.
Если функция д(р) имеет обратную, то
ЧК (Р) = Ч Р), °К (р) =-1- при 0(д р)) Ф 0.
о (Ч -Р))
ИПФ КУ определяется по формуле обращения:
°о ер1 -Тч-1(Р)
8К т) = 2Л7 ] О( -1( )) ¿Р . Оо-гте О(ч (р))
Вычисляя функцию gк ((, т) по сетке узлов (, X/),
tk = —В, к = 0,1,...,N X/ = /В, / = 0,1,...,N где В - достаточно большое
положительное число, методом теории аппроксимации можно построить функцию gк ^, т), аппроксимирующую с любой наперед заданной точностью £ (£>0) функцию gк^,т) .
Заключение
Предложенный в работе подход может быть распространен и на другие динамические системы: нелинейные системы, динамические системы с распределенными параметрами, эредитарные системы. При этом можно получить обощение ряда результатов работ [13, 14] для нелинейных систем, для динамических систем с распределенными параметрами [15, 16] и для эредитарных систем [17, 18].
Библиографический список
1. Бойков, И. В. Аналитические и численные методы идентификации динамических систем : монография / И. В. Бойков, Н. П. Кривулин. - Пенза : Изд-во ПГУ, 2016. - 398 с.
2. Василенко, Г. И. Теория восстановления сигналов / Г. И. Василенко. - М. : Сов. радио, 1979. - 272 с.
3. Грановский, В. А. Динамические измерения: теория и метрологическое обеспечение - вчера и сегодня / В. А. Грановский // Датчики и системы. - 2016. -№ 3 (201). - С. 57-72.
4. Сизиков, В. С. Устойчивые методы обработки результатов измерений : учеб. пособие / В. С. Сизиков. - СПб. : Специальная литература, 1999. - 240 с.
5. Бойков, И. В. Определение динамических характеристик измерительных преобразователей с распределенными параметрами / И. В. Бойков, Н. П. Кривулин // Измерительная техника. - 2000. - № 9. - С. 29.
6. Тихонов, А. Н. Методы решения некорректных задач / А. Н. Тихонов,
B. Я. Арсенин. - 2-е изд. - М. : Наука ; Главная редакция физико-математической литературы, 1979. - 386 с.
7. Кронрод, А. С. Узлы и веса квадратурных формул. Шестнадцатизначные таблицы / А. С. Кронрод. - М. : Наука, 1964. - 144 с.
8. Бойков, И. В. Восстановление входных сигналов дискретных динамических систем / И. В. Бойков, Н. П. Кривулин // Измерительная техника. - 2017. - № 11. -
C. 3-7.
9. Бойков, И. В. Восстановление параметров линейных систем, описываемых дифференциальными уравнениями с переменными коэффициентами / И. В. Бойков, Н. П. Кривулин // Измерительная техника. - 2013. - № 4. - С. 6-11.
10. Бойков, И. В. Приближенные методы вычисления сингулярных и гиперсингулярных интегралов. Ч. 2. Гиперсингулярные интегралы / И. В. Бойков. - Пенза : Изд-во ПГУ, 2009. - 252 с.
11. Бойков, И. В. Приближенные методы вычисления интегралов Адамара и решение гиперсингулярных интегральных уравнений / И. В. Бойков, H. Ф. Добрынина, Л. H. Домнин. - Пенза : Изд-во ПГТУ, 1996. - 188 с.
12. Бойков, И. В. Определение временных характеристик линейных систем с распределенными параметрами / И. В. Бойков, Н. П. Кривулин // Метрология. - 2012. -№ 8. - С. 3-14.
13. Бойков, И. В. Идентификация параметров нелинейных динамических систем / И. В. Бойков, Н. П. Кривулин // Сибирский журнал индустриальной математики. -2018. - Т. XXI, № 2 (74). - С. 17-31.
14. Щербаков, М. А. Итерационный метод оптимальной нелинейной фильтрации изображений / М. А. Щербаков // Известия высших учебных заведений. Поволжский регион. Технические науки. - 2011. - № 4 (20). - С. 43-56.
15. Кривулин, Н. П. Определение параметров физических процессов, описываемых дифференциальными уравнениями в частных производных с переменными коэффициентами / Н. П. Кривулин // Математическое и компьютерное моделирование естественнонаучных и социальных проблем : сб. ст. VIII Междунар. науч.-техн. конф. молодых специалистов, аспирантов и студентов / под ред. И. В. Бой-кова. - Пенза : Изд-во ПГУ, 2014. - С. 172-178.
16. Бойков, И. В. Параметрическая идентификация линейных динамических систем с распределенными параметрами / И. В. Бойков, Н. П. Кривулин // Метрология. - 2014. - № 7. - С. 13-23.
17. Бойков, И. В. Параметрическая идентификация систем, математические модели которых описываются дифференциальными уравнениями с производными дробных порядков / И. В. Бойков, Н. П. Кривулин // Метрология. - 2013. - № 9. -С. 3-16.
18. Бойков, И. В. Параметрическая идентификация эредитарных систем с распределенными параметрами / И. В. Бойков, Н. П. Кривулин // Известия высших учебных заведений. Поволжский регион. Технические науки. - 2013. - № 2 (26). -С. 120-129.
References
1. Boykov I. V., Krivulin N. P. Analiticheskie i chislennye metody identifikatsii dinamich-eskikh sistem: monografiya [Analytical and numerical methods of synamic system identification: monograph]. Penza: Izd-vo PGU, 2016, 398 р.
2. Vasilenko G. I. Teoriya vosstanovleniya signalov [The theory of signal restoration]. Moscow: Sov. radio, 1979, 272 р.
3. Granovskiy V. A. Datchiki i sistemy [Sensors and systems]. 2016, no. 3 (201), рр. 57-72.
4. Sizikov V. S. Ustoychivye metody obrabotki rezul'tatov izmereniy: ucheb. posobie [Stable methods of measurement results processing: tutorial]. Saint-Petersburg: Spetsial'na-ya literatura, 1999, 240 р.
5. Boykov I. V., Krivulin N. P. Izmeritel'naya tekhnika [Measurement technology]. 2000, no. 9, р. 29.
6. Tikhonov A. H., Arsenin V. Ya. Metody resheniya nekorrektnykh zadach [Methods of incorrect problem solution]. 2nd ed. Moscow: Nauka; Glavnaya redaktsiya fiziko-matematicheskoy literatury, 1979, 386 p.
7. Kronrod A. S. Uzly i vesa kvadraturnykh formul. Shestnadtsatiznachnye tablitsy [Nodes and weights of quadrature formulas. Sixteen-digit tables]. Moscow: Nauka, 1964, 144 p.
8. Boykov I. V., Krivulin N. P. Izmeritel'naya tekhnika [Measurement technology]. 2017, no. 11, pp. 3-7.
9. Boykov I. V., Krivulin N. P. Izmeritel'naya tekhnika [Measurement technology]. 2013, no. 4, pp. 6-11.
10. Boykov I. V. Priblizhennye metody vychisleniya singulyarnykh i gipersingulyarnykh in-tegralov. Ch. 2. Gipersingulyarnye integraly [Approximate methods of computing singular and hypersingular integrals. Part 2. Hypersingular integrals]. Penza: Izd-vo PGU, 2009, 252 p.
11. Boykov I. V., Dobrynina H. F., Domnin L. H. Ppiblizhennye metody vychisleniya in-tegpalov Adamapa i peshenie gipepsingulyapnykh integpal'nykh uravneniy [Approximate methods of computing Hadamard integrals and solving hypersingular integral equations]. Penza: Izd-vo PGTU, 1996, 188 p.
12. Boykov I. V., Krivulin N. P. Metrologiya [Metrology]. 2012, no. 8, pp. 3-14.
13. Boykov I. V., Krivulin N. P. Sibirskiy zhurnal industrial'noy matematiki [Siberian journal of industrial mathematics]. 2018, vol. XXI, no. 2 (74), pp. 17-31.
14. Shcherbakov M. A. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Tekhnicheskie nauki [University proceedings. Volga region. Engineering sciences]. 2011, no. 4 (20), pp. 43-56.
15. Krivulin N. P. Matematicheskoe i komp'yuternoe modelirovanie estestvennonauchnykh i sotsial'nykh problem: sb. st. VIII Mezhdunar. nauch.-tekhn. konf. molodykh spetsialis-tov, aspirantov i studentov [Mathematical and computer modeling of natural scientific and social problems: proceedings of VIII International scientific and technical conference of young specialists, undergraduate and postgraduate students]. Penza: Izd-vo PGU, 2014, pp. 172-178.
16. Boykov I. V., Krivulin N. P. Metrologiya [Metrology]. 2014, no. 7, pp. 13-23.
17. Boykov I. V., Krivulin N. P. Metrologiya [Metrology]. 2013, no. 9, pp. 3-16.
18. Boykov I. V., Krivulin N. P. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Tekhnicheskie nauki [University proceedings. Volga region. Engineering sciences]. 2013, no. 2 (26), pp. 120-129.
Кривулин Николай Петрович
кандидат технических наук, доцент, кафедра высшей и прикладной математики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
E-mail: krivulin@bk.ru
Krivulin Nikolay Petrovich Candidate of engineering sciences, associate professor, sub-department of higher and applied mathematics, Penza State University (40 Krasnaya street, Penza, Russia)
УДК 681.311 Кривулин, Н. П.
Восстановление входных сигналов нестационарных динамических систем / Н. П. Кривулин // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2018. - № 3 (47). - С. 64-78. - Б01 10.21685/2072-3040-2018-3-6.
