Все выпуски

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

In parametrical family of subspaces of space of regulated functions the concept of the adjoint integral (in everyone subspace own integral is applied) is defined. In subspace, representing their crossing, the concept of the adjoint integral also is defined. This subspace includes the space of functions of the bounded variation. In any subspace on the basis of the adjoint integral the concept of the generalized regulated function and its adjoint generalized derivative is defined. Solvability of linear impulse systems in terms of adjoint generalized functions is proved.

В статье рассматривается задача устойчивой реконструкции неизвестного входа системы по результатам неточных измерений ее решения. Суть задачи состоит в следующем. Имеется система, описываемая распределенным уравнением второго порядка, решение которой зависит от входа, меняющегося со временем. Как вход, так и решение заранее не известны. В дискретные моменты времени измеряется решение уравнения. Результаты измерения неточны. Требуется построить алгоритм приближенного восстановления входа, обладающий свойствами динамичности и устойчивости. Свойство динамичности означает, что текущие значения приближений входа вычисляются в реальном времени (он-лайн). Свойство устойчивости — что приближения являются достаточно точными, при хорошей точности измерений. Задача относится к классу обратных задач. Представленный в статье алгоритм основан на конструкциях теории устойчивого динамического обращения в комбинации с методами некорректных задач и позиционного управления.

In this paper, we consider the stable reconstruction problem of the unknown input of a distributed system of second order by results of inaccurate measurements of its solution. The content of the problem considered is as follows. We consider a distributed equation of second order. The solution of the equation depends on the input varying in the time. The input, as well as the solution, is not given in advance. At discrete times the solution of the equation is measured. These measurements are not accurate in general. It is required to design an algorithm for approximate reconstruction of the input that has dynamical and stability properties. The dynamical property means that the current values of approximations of the input are produced on-line, and the stability property means that the approximations are arbitrarily accurate for a sufficient accuracy of measurements. The problem refers to the class of inverse problems. The algorithm presented in the paper is based on the constructions of a stable dynamical inversion and on the combination of the methods of ill-posed problems and positional control theory.

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

A mathematical model of competition under conditions of biological invasion, written in the form of a system of nonlinear parabolic equations, is considered. The competition of two closely related species — resident and invader — is studied. The dynamics of populations in a heterogeneous area is determined by local interaction and diffusion. For the invader population, interspecific taxis and directed migration caused by heterogeneity of living conditions are taken into account. In computational experiments, sets of migration parameters corresponding to various invasion scenarios are determined. An analysis of the influence of initial distributions on competitive exclusion and coexistence of species is given.

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

We consider a control system defined by a linear time-invariant system of differential equations with lumped and distributed delays in the state variable. We construct a controller for the system as linear static output feedback with lumped and distributed delays in the same nodes. We study a finite spectrum assignment problem for the closed-loop system. One needs to construct gain coefficients such that the characteristic function of the closed-loop system becomes a polynomial with arbitrary preassigned coefficients. We obtain conditions on coefficients of the system under which the criterion was found for solvability of the finite spectrum assignment problem. Corollaries on stabilization by linear static output feedback with several delays are obtained for the closed-loop system.

Статья посвящена исследованию эффективности применения технологии параллельных вычислений на многопроцессорных системах с общей памятью для задач приближенного расчета множеств достижимости нелинейных управляемых систем в конечномерном евклидовом пространстве. В рамках исследования предложен параллельный алгоритм приближенного построения множеств достижимости, основанный на пошаговой вычислительной схеме с использованием узлов «кубических» сеток для аппроксимации множеств. Предложенный алгоритм предназначен для проведения расчетов на ЭВМ архитектуры SMP и решает вопросы разделения задачи на отдельные подзадачи, синхронизации работы параллельных частей алгоритма и равномерного распределения нагрузки между процессорами. Численное моделирование примеров на ЭВМ с двумя 4-ядерными процессорами с использованием предложенного в статье параллельного алгоритма показало высокую эффективность применения технологии параллельных вычислений для расчета множеств достижимости сеточными методами.

The paper investigates the effectiveness of shared memory parallel programming approach for constructing approximate attainable sets of nonlinear control systems in a finite-dimensional Euclidean space. In this study, we propose a parallel iterative algorithm for constructing approximate attainable sets employing a regular Cartesian grid for spatial discretization. The proposed algorithm has been designed for implementation on SMP systems and handles such issues as data decomposition, threads synchronization and distribution of work between multiple threads. Numerical experiments on a system with two quad-core processors confirmed a high efficiency of shared memory parallel programming approach for applying grid-based methods to construct approximate attainable sets.

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

The stable sequential Pontryagin maximum principle or, in other words, the regularized Pontryagin maximum principle in iterative form is formulated for the optimal control problem of a linear parabolic equation with distributed, initial and boundary controls and operator semiphase equality constraint. The main difference between it and the classical Pontryagin maximum principle is that, firstly, it is formulated in terms of minimizing sequences, secondly, the iterative process occurs in dual space, and thirdly, it is resistant to error of raw data and gives a minimizing approximate solution in the sense of J. Warga. So it is a regularizing algorithm. The proof of the regularized Pontryagin maximum principle in iterative form is based on the dual regularization methods and iterative dual regularization. The results of model calculations of the concrete optimal control problem illustrating the work of the algorithm based on the regularized iterative Pontryagin maximum principle are presented. The problem of finding a control triple with minimal norm under a given equality constraint at the final instant of time or, in other words, the inverse final observation problem of finding a normal solution is used as a concrete model optimal control problem.

В настоящей работе проведено исследование модели деформаций системы из $n$ стилтьесовских струн, расположенных вдоль геометрического графа-звезды, с нелинейным условием в узле. Соответствующая граничная задача имеет вид $$ \left\{\begin{array}{lll} -\left(p_iu_i^\prime\right)(x)+\displaystyle{\int_{0}^{x}}u_i\,dQ_i=F_i(x)-F_i(+0)-(p_iu_i')(+0),\quad i=1,2, \ldots, n,\\ \sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0),\\ u_1(0)=u_2(0)=\ldots=u_n(0)=u(0),\\ (p_iu_i')(l_i-0)+u_i(l_i)\Delta Q_i(l_i)=\Delta F_i(l_i),\quad i=1,2,\ldots, n. \end{array} \right. $$ Здесь функции $u_i(x)$ определяют деформации каждой из струн; $F_i(x)$ описывают распределение внешней нагрузки; $p_i(x)$ характеризуют упругость струн; $Q_i(x)$ описывают упругую реакцию внешней среды. Скачок $\Delta F_i(l_i)$ равняется сосредоточенной в точке $l_i$ внешней силе; скачок $\Delta Q_i(l_i)$ совпадает с жесткостью упругой опоры (пружины), прикрепленной к точке $l_i$. Условие $\sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0)$ возникает за счет наличия в узле ограничителя, представленного отрезком $[-m,m]$, на перемещение струн под воздействием внешней нагрузки, то есть предполагается, что $|u(0)|\leq m$. Здесь через $N_{[-m,m]}u(0)$ обозначен нормальный конус к $[-m,m]$ в точке $u(0)$. В работе проведен вариационный вывод модели; доказаны теоремы существования и единственности решения; проанализированы критические нагрузки, при которых происходит соприкосновение струн с ограничителем; приведена явная формула представления решения.

In the present paper we study a model of deformations for a system of $n$ Stieltjes strings located along a geometric graph-star with a nonlinear condition at the node. The corresponding boundary value problem has the form $$ \left\{\begin{array}{lll} -\left(p_iu_i^\prime\right)(x)+\displaystyle{\int_0^x}u_idQ_i=F_i(x)-F_i(+0)-(p_iu_i')(+0), \quad i=1,2, \ldots, n,\\ \sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0),\\ u_1(0)=u_2(0)=\ldots=u_n(0)=u(0),\\ (p_iu_i')(l_i-0)+u_i(l_i)\Delta Q_i(l_i)=\Delta F_i(l_i), \quad i=1,2,\ldots, n. \end{array} \right. $$ Here the functions $u_i(x)$ determine the deformations of each of the strings; $F_i(x)$ describe the distribution of the external load; $p_i(x)$ characterize the elasticity of strings; $Q_i(x)$ describe the elastic response of the environment. The jump $\Delta F_i(l_i)$ is equal to the external force concentrated at the point $l_i$; the jump $\Delta Q_i(l_i)$ coincides with the stiffness of the elastic support (spring) attached to the point $l_i$. The condition $\sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0)$ arises due to the presence of a limiter in the node represented by the segment $ [-m,m]$, on the movement of strings under the influence of an external load, thus it is assumed that $|u(0)|\leq m$. Here $N_{[-m,m]}u(0)$ denotes the normal cone to $[-m,m]$ at the point $u(0)$. In the present paper a variational derivation of the model is carried out; existence and uniqueness theorems for solutions are proved; the critical loads at which the strings come into contact with the limiter are analyzed; an explicit formula for the representation of the solution is presented.

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

We study a multidimensional case of a nonlinear reaction-diffusion system modeled by a system of two parabolic equations with power nonlinearities. Such systems can be used to simulate the process of propagation in space of interacting distributed formations of robots of two types. Such equations also describe the processes of nonlinear diffusion in reacting two-component continuous media. An original version of the reduction method is proposed, which reduces the construction of the dependence of the exact solution on spatial variables to the solution of the Helmholtz equation, and the dependence on time to the solution of a linear system of ordinary differential equations. A number of examples of multiparameter families of exact solutions given by elementary functions are constructed.

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

We study the stochastic dynamics of the two-dimensional Hindmarsh-Rose model. In the deterministic Hindmarsh-Rose model the parametric zones of coexistence of different stable attractors (equilibria and limit cycles) are possible. The emergence of high amplitude oscillations under the influence of random disturbances on the system in these zones is due to the presence of a limit cycle. However, the stochastic generation of high amplitude oscillations is possible in a parametric zone where the deterministic system has the only stable equilibrium. This article discusses this case. For a sufficiently low noise intensity values, random states concentrate near the stable equilibrium. With the increasing of the noise intensity, trajectories go far from the equilibrium making high amplitude oscillations in the neighborhood of the unstable equilibrium. This phenomenon is confirmed by changing of the probability distribution of random trajectories. This effect is analyzed using the stochastic sensitivity function technique. A method of estimation of critical values for noise intensity is proposed.

Рассматривается линейная система управления, заданная стационарным дифференциальным уравнением с одним сосредоточенным и одним распределенным запаздыванием. В системе на вход подается линейная комбинация из $m$ сигналов и их производных до порядка $n-p$ включительно, а выход представляет собой $k$-мерный вектор линейных комбинаций состояния и его производных до порядка не более $p-1$. Для этой системы исследуется задача управления спектром с помощью линейной статической обратной связи по выходу с сосредоточенным и распределенным запаздываниями. Получены необходимые и достаточные условия разрешимости задачи произвольного размещения спектра посредством статической обратной связи по выходу, имеющей тот же вид, что и система. Получены следствия о стабилизации системы.

A linear control system defined by a stationary differential equation with one lumped and one distributed delay is considered. In the system, the input is a linear combination of $m$ variables and their derivatives of order not more than $n-p$ and the output is a $k$-dimensional vector of linear combinations of the state and its derivatives of order not more than $p-1$. For this system, a spectrum assignment problem by linear static output feedback with delays is studied. Necessary and sufficient conditions are obtained for solvability of the arbitrary spectrum assignment problem by static output feedback controller of the same type as the system. Corollaries on stabilization of the system are obtained.