Все выпуски

Рассматривается выпуклая задача оптимального управления для параболического уравнения со строго равномерно выпуклым целевым функционалом, с граничным управлением и с распределенными поточечными фазовыми ограничениями типа равенства и неравенства. Образы задающих поточечные фазовые ограничения операторов вкладываются в лебегово пространство суммируемых с $s$-й степенью функций при $s\in(1,2)$. В свою очередь, граничное управление принадлежит лебегову пространству с показателем суммируемости $r\in (2,+\infty)$. Основными результатами работы в рассматриваемой задаче оптимального управления с поточечными фазовыми ограничениями являются регуляризованные, или, другими словами, устойчивые к ошибкам исходных данных, секвенциальные принцип Лагранжа в недифференциальной форме и поточечный принцип максимума Понтрягина.

A convex optimal control problem is considered for a parabolic equation with a strictly uniformly convex cost functional, with boundary control and distributed pointwise state constraints of equality and inequality type. The images of the operators that define pointwise state constraints are embedded into the Lebesgue space of integrable with $s$-th degree functions for $s\in(1,2)$. In turn, the boundary control belongs to Lebesgue space with summability index $r\in (2,+\infty)$. The main results of this work in the considered optimal control problem with pointwise state constraints are the two stable, with respect to perturbation of input data, sequential or, in other words, regularized principles: Lagrange principle in nondifferential form and Pontryagin maximum principle.

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

This work is devoted to the construction of a characteristic polynomial of the spectral problem of a first-order differential equation on an interval with a spectral parameter in a boundary value condition with integral perturbation which is an entire analytic function of the spectral parameter. Based on the characteristic polynomial formula, conclusions about the asymptotics of the spectrum of the perturbed spectral problem are established.

Ряд задач в теории характеристических показателей Ляпунова линейных дифференциальных систем

ẋ=A(t)x,    x∈Rⁿ,    t≥0,

сводится к изучению влияния возмущений коэффициентов на характеристические показатели и другие асимптотические инварианты возмущенных систем

ẏ=A(t)y+Q(t)y,    y∈Rⁿ,    t≥0.

При этом возмущения коэффициентов предполагаются принадлежащими некоторым классам малости, то есть определенным подмножествам множества KC_n(R⁺) кусочно-непрерывных и ограниченных на положительной полуоси n×n-матриц. Обычно используемые классы возмущений, например бесконечно малые (исчезающие в бесконечности), экспоненциально убывающие либо суммируемые на полуоси, задаются конкретными аналитическими условиями, но общее определение класса малости в теории показателей отсутствует. На основе анализа свойств общепринятых классов малости нами предложено аксиоматическое определение класса малости возмущений коэффициентов линейных дифференциальных систем, которому удовлетворяет большинство таких классов, используемых в теории характеристических показателей. Это определение достаточно громоздко. Для более компактной характеристики классов малости предложено использовать следующее их свойство: множество возмущений удовлетворяет предложенному определению класса малости тогда и только тогда, когда оно является полной матричной алгеброй над произвольным нетривиальным идеалом кольца функций KC₁(R⁺) (с поточечным умножением), содержащим хотя бы одну строго положительную функцию.

A number of problems in the Lyapunov exponent theory of linear differential systems

ẋ=A(t)x,    x∈Rⁿ,    t≥0,

can be reduced to an investigation of the influence of coefficient perturbations on characteristic exponents and other asymptotic invariants of perturbed systems

ẏ=A(t)y+Q(t)y,    y∈Rⁿ,    t≥0.

Here perturbations are assumed to be in some classes of smallness, i.e. certain subsets of the space KC_n(R⁺) of piecewise continuous and bounded on the positive semiaxis n×n-matrices. Commonly used classes of perturbations, such as infinitesimal (vanishing at infinity), exponentially decaying or integrable on the positive semiaxis are defined by specific analytical conditions, but there is no general definition of the smallness class. By analyzing the desirable properties of commonly used classes, we propose an axiomatic definition for this notion, such that most of classes used in the theory of characteristic exponents satisfy this definition. Since the axioms are somewhat cumbersome, for more compact characterization we propose to use the following property of smallness classes: the set of perturbation satisfies the proposed definition if and only if it is a complete matrix algebra over an arbitrary non-trivial ideal of functional ring KC₁(R⁺) (with the pointwise multiplication) containing at least one strictly positive function.

Пусть $U$ — множество допустимых управлений, $T>0$ и задана шкала банаховых пространств $W[0;\tau]$, $\tau\in(0;T]$, такая, что множество сужений функций из $W=W[0;T]$ на $[0;\tau]$ совпадает с $W[0;\tau]$; $F[.;u]\colon W\to W$ — управляемый вольтерров оператор, $u\in U$. Ранее для операторного уравнения $x=F[x;u]$, $x\in W$, автором была введена система сравнения в форме функционально-интегрального уравнения в пространстве $\mathbf{C}[0;T]$. Было установлено, что для сохранения (относительно малых вариаций правой части) глобальной разрешимости операторного уравнения достаточно сохранения глобальной разрешимости указанной системы сравнения, а также установлены соответствующие достаточные условия. В данной статье рассматриваются дальнейшие примеры приложения этой теории: нелинейное волновое уравнение, сильно нелинейное волновое уравнение, нелинейное уравнение теплопроводности, сильно нелинейное параболическое уравнение.

Let $U$ be the set of admissible controls, $T>0$, and let $W[0;\tau]$, $\tau\in(0;T]$, be a scale of Banach spaces such that the set of restrictions of functions from $W=W[0;T]$ to $[0;\tau]$ coincides with $W[0;\tau]$; let $F[.;u]\colon W\to W$ be a controlled Volterra operator, $u\in U$. Earlier, for the operator equation $x=F[x;u]$, $x\in W$, the author introduced a comparison system in the form of a functional integral equation in the space $\mathbf{C}[0;T]$. It was established that to preserve (under small perturbations of the right-hand side) the global solvability of the operator equation, it is sufficient to preserve the global solvability of the specified comparison system, and the corresponding sufficient conditions were established. In this paper, further examples of application of this theory are considered: nonlinear wave equation, strongly nonlinear wave equation, nonlinear heat equation, strongly nonlinear parabolic equation.

Проводится исследование динамической эволюции шести моделей рассеянных звездных скоплений по данным о фазовых координатах звезд, полученных при численном интегрировании уравнений движения звезд. Для этой цели используются фазовые координаты звезд для 100 равноотстоящих моментов времени от начального t=0 до t_m≅5.1τ_vr (τ_vr - начальное время бурной релаксации скопления). На этом интервале времени ошибки, связанные с округлением и экспоненциальным нарастанием возмущений в исходных координатах звезд, существенно не сказываются на статистических выводах о характере движения звезд скопления. Метод исследования основан на вычислениях взаимных корреляционных функций C_1,2=C_1,2(τ,r) (τ - временная задержка, r - расстояние между точками) для флуктуаций фазовой плотности и применении Фурье-преобразования функций C_1,2 для расчета спектра частот и дисперсионных соотношений. Анализ графиков функций C_1,2, спектров частот и дисперсионных кривых подтверждает существование в моделях волн фазовой плотности, позволяет установить полный спектр радиальных колебаний фазовой плотности, отделить устойчивые колебания от неустойчивых, рассчитать периоды колебаний фазовой плотности и инкременты нарастания неустойчивых колебаний фазовой плотности. Подтверждены теоретические оценки периодов известных неустойчивых гомологических колебаний ядер моделей скоплений. Указываются некоторые астрофизические приложения полученных результатов: возникновение иррегулярных структур в рассеянных скоплениях, слабая турбулентность в движениях звезд скоплений.

The investigation of dynamical evolution of 6 open cluster models is carried out on data about phase coordinates of stars received by numerical integration of stellar motion equations. To attain the aim the phase coordinates of stars for 100 equidistant moments of time from the initial t=0 to t_m≅5.1τ_vr (τ_vr is the initial time of cluster violent relaxation), are used. Over the interval of time the rounding-off errors and errors because of exponential growth of initial coordinates perturbations do not affect statistical conclusions about motion behavior of cluster stars. The investigation method is based on calculations of mutual correlation functions C_1,2=C_1,2(τ,r) (τ  is the time delay, r is the distance between the points) for phase density fluctuations and application of Fourier transformations of functions C_1,2 in order to calculate frequency spectra and dispersion relations. The analysis of graphics C_1,2, frequency spectra and dispersion curves confirms the existence of phase density waves in cluster models, allows to get a complete spectrum of phase density radial oscillations, to separate stable and unstable oscillations, to calculate the periods of phase density oscillations and increments of unstable phase density oscillations. The theoretical estimations of periods of known unstable homological core oscillations of cluster models are confirmed. Pointed out are some astrophysical applications of results received: the origin of irregular structures in open clusters, weak turbulence of cluster star motions.

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

Time-periodic perturbations of an asymmetric pendulum-type equation close to an integrable standard equation of a mathematical pendulum are considered. For an autonomous equation, the problem of limit cycles, which reduces to the study of the Poincaré-Pontryagin generating functions, is solved. A partition of the parameter plane into domains with different behavior of the phase curves is constructed. Basic phase portraits for each domain of the obtained partition are given. For a nonautonomous equation, the question of the structure of the resonance zones, to which the solution of the problem of synchronization of oscillations leads, is studied. Averaged equations of the pendulum type, describing the behavior of solutions of the original equation in individual resonance zones, are calculated and analyzed. The global behavior of solutions in cells that do not contain small neighborhoods of unperturbed separatrices is ascertained. Using the analytical Melnikov method and numerical modeling, the basic bifurcations of the nonautonomous equation associated with the appearance of nonrough homoclinic curves are studied. On the plane of the main parameters, a bifurcation diagram for the Poincaré map generated by the original equation, describing different types of homoclinic tangencies of the separatrices of the saddle fixed point, is constructed. Homoclinic zones (those domains of parameters for which homoclinic trajectories to the saddle fixed point exist) with nonsmooth bifurcation boundaries are found.

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

This paper investigates the rolling motion of a spherical top with an axisymmetric mass distribution on a smooth horizontal plane performing periodic vertical oscillations. For the system under consideration, equations of motion and conservation laws are obtained. It is shown that the system admits two equilibrium points corresponding to uniform rotations of the top about the vertical symmetry axis. The equilibrium point is stable when the center of mass is located below the geometric center, and is unstable when the center of mass is located above it. The equations of motion are reduced to a system with one and a half degrees of freedom. The reduced system is represented as a small perturbation of the problem of the Lagrange top motion. Using Melnikov’s method, it is shown that the stable and unstable branches of the separatrix intersect transversally with each other. This suggests that the problem is nonintegrable. Results of computer simulation of the top dynamics near the unstable equilibrium point are presented.

Исследуется эволюция угла наклона оси вращения планеты в поле притяжения звезды и внешних планет, входящих в планетную систему. Считаем, что исследуемая планета является динамически-симметричным твердым телом $(A = B)$. Полагаем также, что сама планета и внешние планеты движутся по кеплеровским эллипсам вокруг звезды со средними движениями $\omega$ и $\omega_2,\ldots ,\omega_N$, где $N$ - число небесных тел, воздействующих на планету. В переменных Депри-Андуайе получена функция Гамильтона задачи в рамках спутникова приближения. Проведено осреднение функции Гамильтона по быстрым переменным вращательного и орбитального движений при условии отсутствия резонансов между быстрыми частотами указанных движений. Показано, что осредненная функция Гамильтона содержит, помимо классических параметров, параметры $D_i$, являющиеся функционалами на семействе орбит исследуемой планеты и внешних планет. Показано, что осредненная функция Гамильтона допускает разделение переменных и, как следствие, существует три первых интеграла в инволюции. При рассмотрении гравитационных моментов от внешних планет как малых возмущений, получены, с помощью интеграла энергии осредненных уравнений, явные приближенные формулы для угла нутации исследуемой планеты. Получены также приближенные формулы для возмущенного периода прецессии планеты. Проведены расчеты размаха колебаний по углу нутации планеты, возмущенного периода ее прецессии для частного случая планетной системы, состоящей из звезды, самой планеты и массивной внешней планеты (подобной Юпитеру) с симметрично расположенными орбитами, плоскости которых пересекаются под углом $\gamma$.

We investigate the evolution of the obliquity of a planet in the gravitational field of a star and other planets comprising a planetary system. The planet is assumed to be an axially symmetric rigid body ($A=B$). This planet and other planets move around the star along Keplerian ellipses with frequencies $\omega$ and $\omega_2,\ldots,\omega_N$, respectively, where $N$ is the number of celestial bodies (material points) affecting the planet. We derive Hamiltonian for the problem in the Depri-Andoyer variables in the satellite approximation. The Hamiltonian is averaged over the fast variables of the rotational and orbital motions, assuming that the motions are not resonant. The averaged Hamiltonian involves, in addition to the classic parameters, parameters $D_i$, that can be considered as functionals on the family of orbits of celestial bodies comprising the planetary system. The averaged Hamiltonian admits separation of variables, which implies the existence of three first integrals in involution. Regarding the gravitational torques of the other planets as small perturbations, we obtain from the energy integral of the averaged equations explicit approximate expressions for obliquity of the planet and its perturbed period of precession. We investigate numerically the amplitude of oscillations of the planet's obliquity and it's perturbed period of precession for a planetary system involving a star, the planet itself and another massive planet (similar to Jupiter), whose orbits satisfy certain symmetry conditions and orbital planes intersect at angle $\gamma$.

Краевые задачи теории функции комплексных переменных эффективно используются при исследовании равновесия однородных упругих сред. Наиболее сложные системы краевых задач соответствуют случаю, когда упругое тело обладает анизотропными свойствами. Анизотропия среды приводит к появлению в краевых условиях функции сдвига, которая в общем случае нарушает аналитичность искомых функций. В работе проводится исследование систем краевых задач со сдвигом для аналитических векторов, соответствующих трем основным задачам теории упругости (первая, вторая и смешанная задачи). Системы аналитических векторов со сдвигом сводятся к равносильным системам из краевых задач Гильберта для аналитических функций, содержащих интегральные члены со слабой особенностью. Полученное общее решение основных краевых задач анизотропной теории упругости позволяет проверить указанные задачи на устойчивость относительно возмущений краевых условий и формы контура. Такое исследование актуально в связи с необходимостью применения приближенных численных методов к решению краевых задач со сдвигом. Основным результатом работы следует считать доказательство устойчивости систем векторных краевых задач со сдвигом для аналитических функций на пространстве Гёльдера, соответствующих основным задачам теории упругости для анизотропных тел относительно изменения краевых условий и формы контура.

The boundary problems of the complex-variable function theory are effectively used while investigating equilibrium of homogeneous elastic mediums. The most complicated systems of the boundary value problems correspond to the case when an elastic body exhibits anisotropic properties. Anisotropy of the medium results in the drift of boundary conditions of the function that in general disrupts analyticity of the functions of interest. The paper studies systems of the boundary value problems with drift for analytic vectors corresponding to the primal elastic problems (first, second and mixed problems). Systems of analytic vectors with drift are reduced to equivalent systems of Hilbert boundary value problems for analytic functions with weak singularity integrators. The obtained general solution of the primal boundary value problems for the anisotropic theory of elasticity allows us to check the above problems for stability with respect to perturbations of boundary value conditions and contour shape. The research is relevant as there is necessity to apply approximate numerical methods to the boundary value problems with drift. The main research result comes to be a proof of stability of the systems of the vector boundary value problems with drift for analytic functions on the H\"older space corresponding to the primal problems of the elastic theory for anisotropic bodies in the case of change in the boundary value conditions and contour shape.

Конечно-разностные модели дифференциальных уравнений в частных производных, такие как уравнения Лапласа или Пуассона, приводят к конечной сети. Дискретизированное уравнение на неограниченном множестве на плоскости или в пространстве приводит к бесконечной сети. В бесконечной сети оператор Шрёдингера (возмущенный оператор Лапласа, $q$-оператор Лапласа) определяется для развития теории дискретного потенциала, которая имеет модель в уравнении Шрёдингера в евклидовых пространствах. Исследуется связь между $\Delta$-теорией оператора Лапласа и $\Delta_q$-теорией. В $\Delta_q$-теории уравнение Пуассона решается, если сеть является деревом, и в общем случае получается каноническое представление для неотрицательных $q$-супергармонических функций.

Finite-difference models of partial differential equations such as Laplace or Poisson equations lead to a finite network. A discretized equation on an unbounded plane or space results in an infinite network. In an infinite network, Schrödinger operator (perturbed Laplace operator, $q$-Laplace) is defined to develop a discrete potential theory which has a model in the Schrödinger equation in the Euclidean spaces. The relation between Laplace operator $\Delta$-theory and the $\Delta_q$-theory is investigated. In the $\Delta_q$-theory the Poisson equation is solved if the network is a tree and a canonical representation for non-negative $q$-superharmonic functions is obtained in general case.