CC BY
34
КРАТКИЕ СООБЩЕНИЯ
УДК 517.9 DOI: 10.14529/mmp140111
ЗАДАЧА ОПТИМАЛЬНОГО ИЗМЕРЕНИЯ ДЛЯ МОДЕЛИ ИЗМЕРИТЕЛЬНОГО УСТРОЙСТВА С ДЕТЕРМИНИРОВАННЫМ МУЛЬТИПЛИКАТИВНЫМ ВОЗДЕЙСТВИЕМ И ИНЕРЦИОННОСТЬЮ
А.В. Келлер, М.А. Сагадеева
В последнее время результаты теории уравнений соболевского типа активно применяются для измерения динамически искаженных сигналов. В данной работе рассматривается задача оптимального измерения для системы, на которую произведено известное мультипликативное воздействие, которое имеет вид скалярной функции переменной t. Построены точное и приближенное решения задачи оптимального измерения для указанной системы.
Статья состоит из двух частей. В первой части формулируется постановка задачи оптимального измерения для системы с детерминированным мультипликативным воздействием, а во второй приводятся формулы точных и приближенных решений рассматриваемой задачи.
Ключевые слова: оптимальное измерение; система леонтьевского типа; модель Шестакова - Свиридюка.
1. Задача оптимального измерения
При моделировании динамических процессов возникают математические и инженерные задачи в различных областях технических наук [1, 2]. Среди них особо выделяют проблематику динамических измерений, это обусловлено изменением требований к результатам измерений, которое является следствием выдвигаемых требований к качеству испытаний и эффективности производства. Так, одной их наиболее значимых в теории динамических измерений является проблема восстановления измеряемого сигнала. Отметим, что эта задача математически некорректна, поэтому для ее решения А.Л. Шестаковым и его учениками были предложены технически обоснованные гипотезы [2], решения которых были воплощены «в металл>. Затем А.Л. Шестаковым и Г.А. Свиридюком для решения задачи восстановления динамически искаженного сигнала было предложено использовать методы теории оптимального управления [3], а полученную задачу - называть задачей оптимального измерения. Построенная математическая модель - модель Шестакова - Свиридюка -позволила начать численные исследования задачи оптимального измерения [4], опираясь на результаты о численных решениях задачи Коши для систем леонтьевского типа [5].
Для постановки задачи оптимального измерения для модели измерительного устройства с учетом детерминированного мультипликативного воздействия и инерционности введем в рассмотрение пространство состояний ^ = {x € L2((0,t); К”) : X € L2((0,t); К”)}, пространство измерений U = {u € L2((0,т); К”) : up+1 € L2((0,t); R”)} и пространство наблюдений Y = N[^] при некотором фиксированном т € К+. Выделим в U замкнутое и
выпуклое подмножество Яд - множество допустимых измерений. Требуется найти оптимальное измерение V € Яд почти всюду на (0,т), удовлетворяющее системе леонтьевского типа
Ьх(£) = а(£)Мх(£) + Ои(£), (1)
у(1) = Мх(г), (2)
при начальных условиях Шоуолтера - Сидорова [6]
[Кьа(Ы)]р+1(х(0) - хо)=0, (3)
минимизирущее значение функционала
Здесь х = (х1,х2,---,хга) и X = (X 1,Х2,...,хп) - вектор-функции состояния и скорости изменения состояния измерительного устройства соответственно; и = (П1,П2, ... ,ип) и У = (У1,У2,... ,Ут) - вектор-функции измерений и наблюдений измерительного устройства соответственно; у^ = (у^1,у^2,...,У<1т) - наблюдение, полученное в ходе натурного эксперимента; п - число параметров состояния системы; Ь и М - квадратные матрицы порядка п (ёе1 Ь = 0), предстваляющие собой взаимовлияние скоростей изменения состояния и собственно состояния измерительного устройства, при этом скалярная функция а : (0, т) ^ М+. описывает изменение во времени параметров системы этих взаимовлияний; квадратная матрица О порядка п и прямоугольная матрица N размера тхп, характеризуют взаимовлияние параметров измерения и связь между состоянием системы и наблюдением соответственно, функция а : (0, т) ^ Мп будет описана ниже.
Отметим, что представленная модель датчика, описываемая системой леонтьевского типа (1), содержит мультипликативное воздействие а(Ь) на измерительное устройство, то есть система (1) нестационарна (разрешимость таких систем в более общем случае рассматривается, например, в [7]). Кроме того, системы леонтьевского типа являются конечномерным аналогом уравнений соболевского типа [8], и при изучении задачи оптимального измерения воспользуемся результатами о существовании решения задачи оптимального управления для (1), полученными в [9].
2. Точное и приближенное решения задачи оптимального измерения
Пусть Ь и М — квадратные матрицы порядка п, ёе1 Ь = 0. Следуя [8], будем называть множества рь(М) = [ц € С : (1еЛ(цЬ — М) =0} и аь(М) = С \ рь(М) соответственно Ь-резольвентным множеством, и Ь-спектром матрицы М. Нетрудно показать, что либо рь(М) = 0, либо Ь-спектр матрицы М состоит из конечного множества точек [8]. Кроме того, множества рь(М) и аь(М) не изменяются при переходе к другим базисам.
Для комплексной переменной ц € С определим матрично-значные функции (цЬ — М) — Е^(М) = (цЬ — М)-1Ь, Ь£(М) = Ь(цЬ — М)-1 с областью определения рь(М) и будем их называть соответственно Ь-резольвентой, правой и левой Ь-резольвентами матрицы М.
Определение 1. Матрица М называется Ь-регулярной, если рь(М) = 0 и (Ь,р)-регу-лярной, при р равном порядку полюса в ж для функции ёе1(цЬ — М)-1.
Замечание 1. Если бесконечность является устранимой особой точкой Ь-резольвенты матрицы М, то р = 0. Для квадратных матриц параметр р не может превосходить размерности пространства п.
В силу результатов [9] справедлива следующая
Теорема 1. Пусть матрица М (Ь,р)-регулярна и ёе1 М = 0. Тогда для любых Хо Є К”, и Є и и а Є Ср+1([0, т]; К+), отделенной от нуля, существует единственное решение х Є ^ задачи Шоуолтера - Сидорова (3) для (1), имеющее вид
x(u,t) = lim xk(u,t) = lim
k^-ж k^-ж
t ( / * N -1
X0 +
11 і } Y' ^
L — kM j a(()d( j L + j ¡L — 1 M J a(Z)d( j L ^L — 1 M^ (kLk(M))pDu(s)ds+
0
p+i\ { 1 d\ q Du(t)
+ g (м-'((к4(М)Г — l.)L)'M-'(U — (^(МГ1) ^1(І) dtj a(t)
(5)
причем выражение (? в последнем слагаемом означает последовательное приме-\a(t) dt j
нение q раз данного оператора.
Выражением (5) определены точное x(u,t) и приближенное Xk(u,t) решения задачи (1), (3). Решением задачи оптимального измерения является функция v € Uq такая, что
J (v) = / Л™, J (u) (6)
(u,x(u))£Ug хН
где функционал J(и) имеет вид (4), и пара (u,x(u)) € Uq х ^ удовлетворяет системе (1), (2) с начальным условием (3).
В силу результатов, полученных в [10], справедлива следующая
Теорема 2. Пусть матрица M (L,p)-регулярна, p € {0}UN и det M = 0. Тогда при любых Xo € Rn, и € U и a € Cp+1([0,rj;R+), отделенной от нуля, существует единственная пара (v, x(v)) такая, что x(v) — сильное решение задачи Шоуолтера - Сидорова (3) для системы (1), (2), а функция v минимизирует (6) с функционалом (4), причем они связаны формулой
x(v) = lim xk(v, t).
k^-ж
Следуя [10], опишем приближенные решения задачи оптимального измерения. Заменим пространство управлений U на конечномерное пространство U£ = Hj+1(Rra) вектор-
(£ £ £ \
tj ,j2c2itj ,...Y.] Cnj tj I . j=o j=o j=0 )
Учитывая вид (5), необходимо чтобы £ > p. Подставляя и£ вместо и в (5) и (4), будем
рассматривать задачу минимизации функционала J(v£) = min J(и£), в результате получим
решение (v£,x£), причем x£ = x(v£,t) = lim Xk(v£,t).
k^-ж
Приближенным решением задачи оптимального измерения (1) - (4), (6) является функ-у£к (¿), минимизирующая функционал
1
Jk(u£)=¿ /|lNxkq)(u£,t)—ydq)(t)ll2dt
q=0 о
Сходимость таких приближенных решений к точному доказана в [10].
Литература
1. Куропатенко, В.Ф. О моделировании динамических процессов в сферических и цилиндрических оболочках / В.Ф. Куропатенко, Ю.Н. Андреев // Вычислительная механика сплошных сред. - 2010. - Т. 3, № 4. - С. 53-67.
2. Шестаков, А.Л. Коррекция динамической погрешности измерительного преобразователя линейным фильтром на основе модели датчика / А.Л. Шестаков // Известия высших учебных заведений. Приборостроение. - 1991. - Т. 34, № 4. - С. 8-13.
3. Шестаков, А.Л. Новый подход к измерению динамически искаженных сигналов / А.Л. Шестаков, Г.А. Свиридюк // Вестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2010. - № 16 (192), вып. 5. - C. 116-120.
4. Келлер, А.В. Задача оптимального измерения: численное решение, алгоритм программы / А.В. Келлер, Е.И. Назарова // Известия Иркутского государственного университета. Серия: Математика. - 2011. - Т. 4, № 3. - C. 74-82.
5. Свиридюк, Г.А. Численное решение систем уравнений леонтьевского типа / Г.А. Свиридюк, С.В. Брычев // Известия вузов. Математика. - 2003. - № 8. - C. 46-52.
6. Свиридюк, Г.А. Задача Шоуолтера - Сидорова как феномен уравнений соболевского типа / Г.А. Свиридюк, С.А. Загребина // Известия Иркутского государственного университета. Серия: Математика. - 2010. - Т. 3, № 1. - С. 104-125.
7. Сагадеева, М.А. Исследование устойчивости решений линейных уравнений соболевского типа: дис. ... канд. физ.-мат. наук / М.А. Сагадеева. - Челябинск, 2006. - 119 с.
8. Свиридюк, Г.А. К общей теории полугрупп операторов // Успехи математичеких наук.
- 1994. - Т. 49, № 4. - С. 47-74.
9. Сагадеева, М.А. Оптимальное управление решениями нестационарных уравнений соболевского типа специального вида в относительно секториальном случае / М.А. Сагадеева, А.Д. Бадоян // Вестник МаГУ. Математика. - 2013. - Вып. 15. - C. 68-80.
10. Келлер, А.В. Численное решение задач оптимального и жесткого управления для одной нестационарной системы леонтьевского типа / А.В. Келлер, М.А. Сагадеева // Научные ведомости Белгородского государственного университета. Серия: Математика. Физика.
- 2013. - Т. 32, № 19. - C. 57-66.
Алевтина Викторовна Келлер, доктор физико-математических наук, доцент, заведующий кафедрой «Математическое моделирование^ Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].
Минзиля Алмасовна Сагадеева, кандидат физико-математических наук, доцент, кафедра «Информационно-измерительная техника>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), [email protected].____________________________
Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming & Computer Software",
2014, vol. 7, no. 1, pp. 134-138.
MSC 47D06, 49J15, 93A30 DOI: 10.14529/mmp140111
The Optimal Measurement Problem for the Measurement Transducer Model with a Deterministic Multiplicative Effect and Inertia
A.V. Keller, South Ural State University, Chelyabinsk, Russian Federation. [email protected],
M.A. Sagadeeva, South Ural State University, Chelyabinsk, Russian Federation, [email protected]
The results of the theory of Sobolev-type equations are extensively used to measure of dynamically distorted signals recently. In this paper the authors consider the optimal measurement for the system where the well-known multiplicative effect was produced which in its turn has the form of a scalar function of the variable t. The authors develop the exact and approximate solutions of the optimal measurement problem for the specified system.
The paper consists of two parts. The statement of the problem is formulated in the first part as an optimal measurement for the system with a deterministic multiplicative effect, and the second part presents the formulas of exact and approximate solutions of the problem.
Keywords: optimal measurement; Leontiev type system; Shestakov-Sviridyuk model.
References
1. Kuropatenko V.F., Andreev Yu.N. Simulation of Dynamic Processes in Spherical and Cylindrical Shells [O modelirovanii dinamicheskih processov v sfericheskih i cilindricheskih obolochkah]. Vychislitel’naya mehanika sploshnyh sred [Computational Continuum Mechanics], 2010, vol. 3, no. 4, pp. 53-67.
2. Shestakov A.L. Dynamic Error Correction Transducer Linear Filter-based Sensor Model [Korrektsiya dinamicheskoy pogreshnosti izmeritel’nogo preobrazovatelya lineynym fil’trom na osnove modeli datchika]. Izvestiya VUZ. Priborostroenie, 1991, vol. 34, no. 4, pp. 8-13.
3. Shestakov A.L., Sviridyuk G.A. A New Approach to Measuring Dynamically Distorted Signals [Novyy podkhod k izmereniyu dinamicheski iskazhennykh signalov]. Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming & Computer Software", 2010, no. 16 (192), issue 5, pp. 116-120.
4. Keller A.V., Nazarova E.I. Optimal Measuring Problem: the Computation Solution,
the Program Algorithm [Zadacha optimal’nogo izmereniya: chislennoe reshenie, algoritm programmy]. Izvestiya Irkutskogo gosudarstvennogo universiteta. Seriya "Matematika" [The Bulletin of Irkutsk State University. Series "Mathematics"], 2011, vol. 4, no. 3, pp. 74-82.
5. Sviridyuk G.A., Brychev S.V. Numerical solution of Leontiev type systems [Chislennoe reshenie sistem uravneniy leont’evskogo tipa]. Izvestiya VUZ. Matematika [Russian Mathematics], 2003, no. 8, pp. 46-52.
6. Sviridyuk G.A., Zagrebina S.A. The Showalter-Sidorov Problem as a Phenomena of the Sobolev Type Equations [Zadacha Shouoltera - Sidorova kak fenomen uravneniy sobolevskogo tipa]. Izvestiya Irkutskogo gosudarstvennogo universiteta. Seriya "Matematika" [The Bulletin of Irkutsk State University. Series "Mathematics"], 2010, vol. 3, no. 1, pp. 133-137.
7. Sagadeeva M.A. Investigation of Solutions Stability for Linear Sobolev Type Equations [Issledovanie ustoychivosti resheniy lineynykh uravneniy sobolevskogo tipa: dis...kand. fiz.-mat. nauk]. Chelyabinsk, 2006.
8. Sviridyuk G.A. On the General Theory of Operator Semigroups. Russian Mathematical Surveys, 1994, vol. 49, no. 4, pp. 45-74. DOI: 10.1070/RM1994v049n04ABEH002390
9. Sagadeeva M.A., Badoyan A.D. The Optimal Control over Solutions of Special Form of Nonstacionary Sobolev Type Equations in Relatively Spectral Case [Optimal’noe upravlenie resheniyami nestatsionarnykh uravneniy sobolevskogo tipa spetsial’nogo vida v otnositel’no sektorial’nom sluchae]. Vestnik Magnitogorskogo gosudarstvennogo universiteta. Matematika [Bulletin of Magnitogorsk State University. Mathematics], 2013, no. 15, pp. 68-80.
10. Keller A.V., Sagadeeva M.A. The Numerical Solution of Optimal and Hard Control for Nonstationary System of Leontiev Type [Chislennoe reshenie zadach optimal’nogo i zhestkogo upravleniya dlya odnoy nestacionarnoy sistemy leont’evskogo tipa]. Nauchnye vedomosti BelGU. Seriya: Matematika i fizika [Belgorod State University Scientific Bulletin. Mathematics & Physics], 2013, vol. 32, no. 19, pp. 57-66.
Поступила в редакцию 15 ноября 2013 г.
