Все выпуски

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

The influence of additive and parametrical noise on attractors of the one-dimensional system governed by the stochastic differential Ito equation is investigated. It is shown that unlike additive, parametrical disturbances lead to the shift of extrema of probability density function. For the value of this shift, a decomposition on small parameter of noise intensity is obtained. It is shown that the influence of the parametrical noise can change not only the arrangement, but also the quantity of extrema of probability density function. The corresponding noise-induced phenomena are studied for three dynamical models in detail. An analysis of the error for the different order estimations of the shift of extrema for the probability density function is presented by the example of a linear model. Two scenarios of the transition between unimodal and bimodal forms of the stochastic attractor are investigated for systems with different types of cubic nonlinearity.

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

We study the stability of linear autonomous scalar difference equations with complex coefficients. For an equation with an arbitrary number of delays, we propose a simple proof of the linear connectivity of the stability region in the space of coefficients. This result allows us to assert that the stability region of the equation in the space of coefficients is the region of the $D$-decomposition of this space containing the origin of coordinates. Further, we consider some equations with two delays and complex coefficients, for which we give detailed analytic and geometric descriptions of the regions of uniform and exponential stability.

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

We study the variational approach to setting and solving of the function approximating problem by quasipolynomials which are the solutions of the homogeneous autonomous linear difference or differential equations.

Для автономного линейного разностного уравнения четвертого порядка исследуется область устойчивости $S_4$ в четырехмерном пространстве, параметрами которого являются коэффициенты уравнения. Исследованы и построены в трехмерном пространстве сечения области $S_4$ гиперплоскостями четырехмерного пространства при фиксированном значении четвертого параметра. Показано, что эти сечения диффеоморфны области устойчивости $S_3$ линейного разностного уравнения третьего порядка.

For an autonomous fourth-order linear difference equation, the stability domain $S_4$ in four-dimensional space is investigated, the parameters of which are the equation coefficients. Sections of the domain $S_4$ by hyperplanes of four-dimensional space are studied and constructed in three-dimensional space for a fixed value of the fourth parameter. It is shown that these sections are diffeomorphic to the stability domain $S_3$ of the third-order linear difference equation.

Изучаются свойства дискретной вариационной задачи динамической аппроксимации в комплексном евклидовом (L + 1)-мерном пространстве E. Она обобщает известные задачи среднеквадратической полиномиальной аппроксимации функций, заданных своими отсчетами в конечном интервале. В рассматриваемой задаче аппроксимация последовательности y = {y_i}_L⁰ отсчетов функции y(t) ∈ L²[0, T], T = Lh на сетке I_h осуществляется решениями однородных линейных дифференциальных или разностных уравнений заданного порядка n с постоянными, но, возможно, неизвестными коэффициентами. Тем самым показано, что в последнем случае задача аппроксимации включает в себя и задачу идентификации. Анализ ее особенностей - основная тема статьи. Ставится задача нахождения вектора коэффициентов  разностного уравнения Σ_n⁰ ŷ_i+k α_i = 0, где k = 0,L − n. Оптимизируются коэффициенты  и начальные условия переходного процесса y этого уравнения. Цель оптимизации - наилучшая аппроксимация исследуемого динамического процесса y ∈ E. Критерий аппроксимации  минимум величины ||y − ŷ||²_E. Показано, что изучаемая вариационная задача сводится к задачам проектирования в E вектора y на ядра разностных операторов с неизвестными коэффициентами α ∈ ω ⊂ S ⊂ Eⁿ⁺¹. Здесь α - направление, S - сфера или гиперплоскость. Показана связь изучаемой задачи с задачами дискретизации и идентифицируемости. Тогда координаты вектора y ∈ E есть точное решение дифференциального уравнения на сетке I_h и y = ŷ. Дано сравнение изучаемой задачи вариационной идентификации с алгебраическими методами идентификации. Показано, что ортогональные дополнения к ядрам разностных операторов всегда имеют теплицев базис. Это приводит к быстрым проекционным алгоритмам вычислений. Показано, что задача нахождения оптимального вектора α сводится к задаче безусловной минимизации функционала идентификации, зависящего от направления  в Eⁿ⁺¹. Предложена итерационная процедура его минимизации на сфере с широкой областью и высокой скоростью сходимости. Изучаемую вариационную задачу можно применять при математическом моделировании в управлении и научных исследованиях. При этом на конечных интервалах может использоваться, в частности, возможность кусочно-линейной динамической аппроксимации сложных динамических процессов разностными и дифференциальными уравнениями указанного типа.

Some properties of the discrete variational problem of the dynamic approximation in the complex Euclidean (L + 1)-dimensional space are studied here. It generalizes familiar problems of the mean square polynomial approximation of the functions given on the finite interval in accordance with their references. In the problem under consideration sequence approximation y = {y_i}_L⁰ of the references of the function y(t) ∈ L²[0, T], T = Lh on the lattice I_h is achieved by solving homogeneous linear differential equations or difference equations of the given order n with constant but possibly unknown coefficients. Thus, it is shown that in the latter case the approximation problem also includes the identification problem. The analysis of its properties is the main subject of the article. The problem is set to find vector of coefficients  of difference equation Σ_n⁰ ŷ_i+k α_i = 0, where k = 0,L − n. Coefficients  and initial conditions of the transient process by of this equation are optimized. The optimization purpose is to achieve the best approximation of the dynamic process y ∈ E being considered here. The approximation criterion is a minimum of the quantity ||y − ŷ||²_E. The variational problem under study is shown to be reduced to the problem of projecting vector y in E on the kernels of the difference operators with unknown coefficients  α ∈ ω ⊂ S ⊂ Eⁿ⁺¹, where  is a direction, S is a sphere or a hyperplane. The problem under study is shown to be related to the problems of the discretization and identifiability. In this case vector coordinates y ∈ E is an exact solution of differential equation on the lattice I_h and y = ŷ. The problem of the variational identification is compared with algebraic methods of identification. The orthogonal complement to the kernels of the difference operators are shown to always have Toeplitz basis. This results in fast projecting algorithms of computation. The problem of finding optimal vector α is shown to be reduced to the problem of the absolute minimization of the identification functional depending on the direction  in Eⁿ⁺¹. The iterative procedure of its minimization on a sphere with wide domain and high speed of convergence is presented here. The variational problem considered here can be applied in mathematical modeling for control problem and research purposes. On the finite intervals, for example, it is possible to use piecewise-linear dynamic approximations of the complex dynamic processes with difference and differential equations of the specified type.

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

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.

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

Two problems of nonlinear guaranteed estimation for states of dynamical systems are considered. It is supposed that unknown measurable in $t$ disturbances are linearly included in the equation of motion and are additive in the measurement equations. These disturbances are constrained by nonlinear integral functionals, one of which is analog of functional of the generalized work. The studied problem consists in creation of the information sets according to measurement data containing the true position of the trajectory. The dynamic programming approach is used. If the first functional requires solving a nonlinear equation in partial derivatives of the first order which is not always possible, then for functional of the generalized work it is enough to find a solution of the linear Lyapunov equation of the first order that significantly simplifies the problem. Nevertheless, even in this case it is necessary to impose additional conditions on the system parameters in order for the system trajectory of the observed signal to exist. If the motion equation is linear in state variable, then many assumptions are carried out automatically. For this case the issue of mutual approximation of information sets via inclusion for different functionals is discussed. In conclusion, the most transparent linear quadratic case is considered. The statement is illustrated by examples.

Выведена формула, связывающая фундаментальное решение и матрицу Коши линейного автономного скалярного уравнения нейтрального типа.

The paper presents simple formulas relating the fundamental solution of a linear differential-difference equation of neutral type to its Cauchy function.

В данной работе исследуются различные разновидности показателей колеблемости (верхние или нижние, сильные или слабые) нулей, корней, гиперкорней, строгих и нестрогих знаков ненулевых решений линейных однородных дифференциальных систем на положительной полуоси. На множестве ненулевых решений систем установлены соотношения между этими показателями колеблемости. Доказано, что все сильные показатели колеблемости (в отличие от частот Сергеева смен знаков, нулей и корней, а также всех слабых показателей колеблемости), рассматриваемые как функции на множестве решений линейных однородных дифференциальных систем с непрерывными на полуоси коэффициентами, не являются остаточными (т.е. могут меняться при изменении решения на конечном отрезке). Кроме того, при любом наперед заданном натуральном $n\ge2$ приводится пример $n$-мерной дифференциальной системы, у которой все сильные показатели колеблемости некоторого решения не совпадают с соответствующими слабыми показателями. При этом все слабые и все сильные показатели на выбранном решении совпадают соответственно между собой. При доказательстве результатов настоящей работы отдельно рассмотрены случаи четности и нечетности $n$.

In this paper, we study various types of exponents of oscillation (upper or lower, strong or weak) of zeros, roots, hyperroots, strict and non-strict signs of non-zero solutions of linear homogeneous differential systems on the positive semi-axis. On the set of non-zero solutions of autonomous systems the relations between these exponents of oscillation are established. It is proved that all strong exponents of oscillations (unlike Sergeev's frequencies of sign changes, zeros and roots, as well as all the weak exponents of oscillations) considered as functions on the set of solutions to linear homogeneous $n$-dimensional differential systems with continuous coefficients on the semi-line are not residual (i.e. can be changed when changing solution on a finite interval). Besides, at any beforehand given natural $n\ge2$ we give the example of $n$-dimensional differential system, for some solution of which all strong oscillation exponents differ from corresponding weak exponents. In this case, all weak and all strong exponents on the chosen solution coincide with each other, respectively. When proving the results of this work, the case of parity and odd $n$ are considered separately.

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

In this paper, we study the sections of the generating series for solutions to a linear multidimensional difference equation with constant coefficients and find recurrent relations for these sections. As a consequence, a multidimensional analogue of Moivre's theorem on the rationality of sections of the generating series depending on the form of the initial data of the Cauchy problem for a multidimensional difference equation is proved. For problems on the number of paths on an integer lattice, it is shown that the sections of their generating series represent the well-known sequences of polynomials (Fibonacci, Pell, etc.) with a suitable choice of steps.