Все выпуски

Рассмотрено применение барицентрического метода для численного решения задач Дирихле и Неймана для уравнения Гельмгольца в ограниченной односвязной области $\Omega\subset\mathbb{R}^2$. Основное допущение в решении заключается в задании границы $\Omega$ в кусочно-линейном представлении. Отличительная особенность барицентрического метода состоит в порядке формирования глобальной системы векторных базисных функций для $\Omega$ через барицентрические координаты. Установлены существование и единственность решения задач Дирихле и Неймана для уравнения Гельмгольца барицентрическим методом и определена оценка скорости сходимости. Уточнены особенности алгоритмической реализации метода.

The application of the barycentric method for the numerical solution of Dirichlet and Neumann problems for the Helmholtz equation in the bounded simply connected domain $\Omega\subset\mathbb{R}^2$ is considered. The main assumption in the solution is to set the $\Omega$ boundary in a piecewise linear representation. A distinctive feature of the barycentric method is the order of formation of a global system of vector basis functions for $\Omega$ via barycentric coordinates. The existence and uniqueness of the solution of Dirichlet and Neumann problems for the Helmholtz equation by the barycentric method are established and the convergence rate estimate is determined. The features of the algorithmic implementation of the method are clarified.

В первой части определено и исследовано нелинейное метрическое пространство $\langle\overline{\rm G}^\infty[a,b],d\rangle$, состоящее из функций, действующих из отрезка $[a,b]$ в расширенную числовую ось $\overline{\mathbb R}$. По определению предполагается, что для любых $x\in\overline{\rm G}^\infty[a,b]$ и $t\in(a,b)$ существуют предельные числа $x(t-0),x(t+0)\in\overline{\mathbb R}$ (и числа $x(a+0),x(b-0)\in\overline{\mathbb R}$). Доказана полнота пространства. Оно является замыканием пространства ступенчатых функций в метрике $d$. Во второй части работы определено и исследовано нелинейное пространство ${\rm RL}[a,b]$. Всякая кусочно-гладкая функция, определенная на $[a,b]$, содержится в ${\rm RL}[a,b]$. Всякая функция $x\in{\rm RL}[a,b]$ имеет ограниченное изменение. Для нее определены все односторонние производные (со значениями в метрическом пространстве $\langle\overline{\mathbb R},\varrho\rangle$). Функция левосторонних производных непрерывна слева, а функция правосторонних производных непрерывна справа. Обе функции, доопределенные на весь отрезок $[a,b]$, принадлежат пространству $\overline{\rm G}^\infty[a,b]$. В заключительной части работы определены и исследованы два подпространства пространства ${\rm RL}[a,b]$. В подпространствах сформулированы и обсуждены перспективные постановки для простейших вариационных задач.

In the first part of the paper, the nonlinear metric space $\langle\overline{\rm G}^\infty[a,b],d\rangle$ is defined and studied. It consists of functions defined on the interval $[a,b]$ and taking the values in the extended numeric axis $\overline{\mathbb R}$. For any $x\in\overline{\rm G}^\infty[a,b]$ and $t\in(a,b)$ there are limit numbers $x(t-0),x(t+0) \in\overline{\mathbb R}$ (and numbers $x(a+0),x(b-0)\in\overline{\mathbb R}$). The completeness of the space is proved. It is the closure of the space of step functions in the metric $d$. In the second part of the work, the nonlinear space ${\rm RL}[a,b]$ is defined and studied. Every piecewise smooth function defined on $[a,b]$ is contained in ${\rm RL}[a,b]$. Every function $x\in{\rm RL}[a,b]$ has bounded variation. All one-sided derivatives (with values in the metric space $\langle\overline{\mathbb R},\varrho\rangle$) are defined for it. The function of left-hand derivatives is continuous on the left, and the function of right-hand derivatives is continuous on the right. Both functions extended to the entire interval $[a,b]$ belong to the space $\overline{\rm G}^\infty[a,b]$. In the final part of the paper, two subspaces of the space ${\rm RL}[a,b]$ are defined and studied. In subspaces, promising formulations for the simplest variational problems are stated and discussed.

Понятие равномерной полной управляемости линейной системы, введенное Р. Калманом, играет ключевую роль в задачах управления асимптотическими характеристиками линейных систем управления, замкнутых по принципу линейной обратной связи. Е.Л. Тонков установил необходимое и достаточное условие равномерной полной управляемости для систем с кусочно-непрерывными и ограниченными коэффициентами. Критерий Тонкова можно положить в основу определения равномерной полной управляемости. Если условия на коэффициенты системы ослабить, то определения Калмана и Тонкова перестают совпадать. Здесь установлены необходимые условия и достаточные условия равномерной полной управляемости по Калману и по Тонкову для систем с измеримыми, локально суммируемыми коэффициентами. Введено определение равномерной полной управляемости, которое обобщает определение Тонкова и совпадает с определением Калмана, если матрица $B(\cdot)$ ограничена. Доказаны некоторые известные результаты об управляемости линейных систем, в которых можно ослабить требования на коэффициенты. Доказано, что если линейная управляемая система $\dot x=A(t)x+B(t)u$, $x\in\mathbb{R}^n$, $u\in\mathbb{R}^m$, с измеримой ограниченной матрицей $B(\cdot)$ равномерно вполне управляема в смысле Калмана, то для любой измеримой и интегрально ограниченной $m\times n$-матричной функции $Q(\cdot)$ система $\dot x=(A(t)+B(t)Q(t))x+B(t)u$ равномерно вполне управляема по Калману.

The notion of uniform complete controllability of linear system introduced by R. Kalman plays a key role in problems of control of asymptotic properties for linear systems closed by linear feedback control. E.L. Tonkov has found a necessary and sufficient condition of uniform complete controllability for systems with piecewise continuous and bounded coefficients. The Tonkov criterion can be considered as the definition of uniform complete controllability. If the coefficients of the system satisfy weak conditions then the definitions of Kalman and Tonkov are not coincide. We obtain necessary conditions and sufficient conditions for uniform complete controllability in the sense of Kalman and Tonkov for systems with measurable and locally integrable coefficients. We introduce a new definition of uniform complete controllability that extends the definition of Tonkov and coincides with the definition of Kalman providing the matrix $B(\cdot)$ is bounded. We prove some known results on the controllability of linear systems that allow the weakening of the requirements on the coefficients. We prove that if a linear control system $\dot x=A(t)x+B(t)u$, $x\in\mathbb{R}^n$, $u\in\mathbb{R}^m$, with measurable and bounded matrix $B(\cdot)$ is uniformly completely controllable in the sense of Kalman then for any measurable and integrally bounded $m\times n$-matrix function $Q(\cdot)$ the system $\dot x=(A(t)+B(t)Q(t))x+B(t)u$ is also uniformly completely controllable in the sense of Kalman.

Изучаются свойства дискретной вариационной задачи динамической аппроксимации в комплексном евклидовом (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.

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

ẋ=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.

Рассматривается линейная нестационарная управляемая система $$\dot x =A(t)x+ B(t)u, \quad x\in\mathbb{R}^n,\quad u\in\mathbb{R}^m,\quad t\in \mathbb{R}, \qquad \qquad (1)$$ с кусочно-непрерывными и ограниченными $\omega$-периодическими матрицами коэффициентов $A(\cdot)$ и $B(\cdot)$. Управление в системе (1) строится по принципу линейной обратной связи $u=U(t)x$ с кусочно-непрерывной и ограниченной матричной функцией $U(t)$, $t\in \mathbb{R}$. Для замкнутой системы $$\dot x =(A(t)+B(t)U(t))x, \quad x\in\mathbb{R}^n, \quad t\in \mathbb{R}, \qquad \qquad (2)$$ исследуется вопрос об условиях ее равномерной глобальной достижимости. Наличие последнего свойства у системы (2) означает существование такой матричной функции $U(t)$, $t\in \mathbb{R}$, которая обеспечивает для матрицы Коши $X_U(t,s)$ этой системы выполнение равенств $X_U((k+1)T,kT)=H_k$ при фиксированном $T>0$ и произвольных $k\in\mathbb{Z}$, $\det H_k>0$. Представленная задача решается в предположении равномерной полной управляемости (в смысле Калмана) системы (1), соответствующей замкнутой системе (2), т.е. при условии существования для системы (1) таких чисел $\sigma>0$ и $\alpha_i>0$, $i=\overline{1,4}$, что при всяких числе $t_0\in\mathbb{R}$ и векторе $\xi\in \mathbb{R}^n$ справедливы неравенства $$\alpha_1\|\xi\|^2\leqslant\xi^*\int\nolimits_{t_0}^{t_0+\sigma}X(t_0,s)B(s)B^*(s)X^*(t_0,s)\,ds\,\xi\leqslant\alpha_2\|\xi\|^2,$$ $$\alpha_3\|\xi\|^2\leqslant\xi^*\int\nolimits_{t_0}^{t_0+\sigma}X(t_0+\sigma,s)B(s)B^*(s)X^*(t_0+\sigma,s)\,ds\,\xi\leqslant\alpha_4 \|\xi\|^2,$$ в которых $X(t,s)$ - матрица Коши линейной системы (1) при $u(t)\equiv0.$ Доказано, что свойство равномерной полной управляемости (в смысле Калмана) периодической системы (1) является необходимым и достаточным условием равномерной глобальной достижимости соответствующей системы (2).

We consider a linear time-varying control system $$\dot x =A(t)x+ B(t)u, \quad x\in\mathbb{R}^n,\quad u\in\mathbb{R}^m,\quad t\in \mathbb{R}, \qquad \qquad (1)$$ with piecewise continuous and bounded $\omega$-periodic coefficient matrices $A(\cdot)$ and $B(\cdot).$ We construct control of the system (1) as a linear feedback $u=U(t)x$ with piecewise continuous and bounded matrix function $U(t)$, $t\in \mathbb{R}$. For the closed-loop system $$\dot x =(A(t)+B(t)U(t))x, \quad x\in\mathbb{R}^n, \quad t\in \mathbb{R}, \qquad \qquad (2)$$ the conditions of its uniform global attainability are studied. The latest property of the system (2) means existence of matrix $U(t)$, $t\in \mathbb{R}$, ensuring equalities $X_U((k+1)T,kT)=H_k$ for the state-transition matrix $X_U(t,s)$ of the system (2) with fixed $T>0$ and arbitrary $k\in\mathbb{Z}$, $\det H_k>0$. The problem is solved under the assumption of uniform complete controllability (by Kalman) of the system (1), corresponding to the closed-loop system (2), i.e. assuming the existence of such numbers $\sigma>0$ and $\alpha_i>0$, $i=\overline{1,4}$, that for any number $t_0\in\mathbb{R}$ and vector $\xi\in \mathbb{R}^n$ the following inequalities hold: $$\alpha_1\|\xi\|^2\leqslant\xi^*\int\nolimits_{t_0}^{t_0+\sigma}X(t_0,s)B(s)B^*(s)X^*(t_0,s)\,ds\,\xi\leqslant\alpha_2\|\xi\|^2,$$ $$\alpha_3\|\xi\|^2\leqslant\xi^*\int\nolimits_{t_0}^{t_0+\sigma}X(t_0+\sigma,s)B(s)B^*(s)X^*(t_0+\sigma,s)\,ds\,\xi\leqslant\alpha_4 \|\xi\|^2,$$ where $X(t,s)$ is the state-transition matrix of linear system (1) with $u(t)\equiv0.$ It is proved that the property of uniform complete controllability (by Kalman) of the periodic system (1) is a necessary and sufficient condition of uniform global attainability of the corresponding system (2).

Рассматривается уравнение в частных производных первого порядка с эффектом наследственности:

$$ \frac{\partial u(x,t)}{\partial t} + a \frac{\partial u(x,t)}{\partial x} = f ( x, t, u(x,t), u_t(x,\cdot)),$$ $$u_t(x,\cdot) = \{u(x,t+s), -\tau\leqslant s <0\}.$$

Для такого уравнения, с позиций принципа разделения конечномерной и бесконечномерной составляющих состояния, строятся сеточные методы: аналог семейства схем бегущего счета, аналог схемы Кранка-Николсон, метод аппроксимации на середину квадрата. Для учета эффекта наследственности применяются одномерная и двойная кусочно-линейная интерполяции и экстраполяция продолжением. Доказывается, что рассмотренные методы имеют порядки локальной погрешности: соответственно $O(h+\Delta)$, $O(h+\Delta^2)$ и $O(h^2+\Delta^2)$, где $h$ - шаг дискретизации по пространственной переменной, $\Delta$ - шаг дискретизации по временной переменной. Исследуются свойства двойной кусочно-линейной интерполяции. Используя результаты общей теории разностных схем, установлены условия устойчивости предложенных методов. С помощью вложения в общую схему численных методов для функционально-дифференциальных уравнений получены теоремы о порядках сходимости сконструированных алгоритмов. Приведены тестовые примеры по сравнению погрешностей методов.

We consider a first-order partial differential equation with heredity effect

$$ \frac{\partial u(x,t)}{\partial t} + a \frac{\partial u(x,t)}{\partial x} = f ( x, t, u(x,t), u_t(x,\cdot)),$$ $$u_t(x,\cdot) = \{u(x,t+s), -\tau\leqslant s <0\}.$$

For such an equation we construct grid methods using the principle of separation of finite-dimensional and infinite-dimensional state components. These grid methods are: analog of running schemes family, analog of Crank-Nicolson scheme, an approximation method to the middle of the square. The one-dimensional and double piecewise linear interpolation and the extrapolation by continuation are applied in order to account the effect of heredity. It is shown that the considered methods have orders of a local error: $O (h +\Delta) $, $O (h +\Delta^2) $ and $O (h^2 +\Delta^2)$ respectively, where $h$ is the spatial discretization interval, $\Delta$ is the time discretization interval. Properties of double piecewise linear interpolation are investigated. Using the results of the general theory of differential schemes, stability conditions of the proposed methods are established. Including them in the general scheme of numerical methods for the functional-differential equations, theorems of orders of proposed algorithms convergence are received. Test examples comparing errors of methods are given.

Пусть $X_0\subseteq\mathbb R^n$ — непустое открытое множество и $X_0\subseteq X\subseteq\overline X_0$. Допускается, что множество $X_0$ не ограничено и/или имеет счетное число компонент связности. В работе исследуются некоторые пространства функций $f\colon X\to\mathbb R$, наделенные специальной нормой $\|\cdot\|$. В определении нормы фигурирует $n$-мерный вектор $(\Delta x)^{-1}\Delta f$, являющийся аналогом отношения $\frac{\Delta f}{\Delta x}$, порождающего понятие производной функции одной переменной. Вектор $(\Delta x)^{-1}\Delta f$ можно ассоциировать с вектором $\mathrm{grad}\,f(\cdot)$. Обратимая матрица $\Delta x$ порядка $n$ состоит из специальных приращений аргумента ${x\in \mathbb R^n}$, а вектор $\Delta f$ состоит из специальных приращений функции $f$. Доказан ряд свойств вектора $(\Delta x)^{-1}\Delta f$, получена точная формула для его евклидовой нормы. Доказана полнота по специальной норме $\|\cdot\|$ пространства $\mathcal G(X)$, состоящего из непрерывных ограниченных функций $f\colon X\to\mathbb R$ и имеющих дополнительные ограничения типа ограничений Липшица–Гёльдера. Подобные функции играют важную роль при решении задач математической физики. Исследован ряд актуальных подпространств пространства $\mathcal G(X)$, доказано, что два из них банаховы, одно из них при $n=1$ и при определенных условиях является замыканием пространства кусочно-линейных функций $f\colon X\to\mathbb R$.

Let $X_0\subseteq\mathbb R^n$ be a nonempty open set and $X_0\subseteq X\subseteq\overline X_0$. We admit that the set $X_0$ is unbounded and/or has a countable number of connected components. In this paper, we study some spaces of functions $f\colon X\to\mathbb R$ endowed with a special norm $\|\cdot\|$. The definition of the norm involves an $n$-dimensional vector $(\Delta x)^{-1}\Delta f$, which is an analogue of the relation $\frac{\Delta f}{\Delta x}$ generating the concept of the derivative of a function of one variable. The vector $(\Delta x)^{-1}\Delta f$ can be associated with the vector $\mathrm{grad}\,f(\cdot)$. The invertible matrix $\Delta x$ of order $n$ consists of special increments of the argument $x\in \mathbb R^n$, and the vector $\Delta f$ consists of special increments of the function $f$. A number of properties of the vector $(\Delta x)^{-1}\Delta f$ is proved, and an exact formula for its Euclidean norm is obtained. We prove the completeness with respect to a special norm $\|\cdot\|$ of the space $\mathcal G(X)$ consisting of continuous bounded functions $f\colon X\to\mathbb R$ and having additional restrictions of the Lipschitz–Hölder type. Such functions play an important role in solving mathematical physics problems. A number of important subspaces of the space $\mathcal G(X)$ is investigated. It is proved that two of them are Banach, and one of them, for $n=1$ and under certain conditions, is the closure of the space of piecewise linear functions $f\colon X\to\mathbb R$.

Определяется понятие регулярно гладкой функции. Кусочно-гладкие функции являются регулярно гладкими, а всякая регулярно гладкая функция является липшицевой. Регулярно гладкие функции имеют конечные односторонние производные: левосторонняя производная непрерывна слева, а правосторонняя непрерывна справа. Односторонние производные порождают понятие регулярной производной. Пространство регулярно гладких функций является замыканием пространства кусочно-линейных функций по норме пространства липшицевых функций. Пространство кусочно-гладких функций всюду плотно в пространстве регулярно гладких функций. Получен аналог уравнения Эйлера для простейшей вариационной задачи в пространстве регулярно гладких функций.

The concept of regular smooth function is defined. Any piecewise smooth function is regular smooth function, and any regular smooth function is Lipschitzian. Any regular smooth function has finite one-sided derivatives: the left-side derivative is continuous at the left and the right-side derivative is continuous on the right. One-sided derivatives generate concept of a regular derivative. The space of regular smooth functions is the closure of the space of piecewise linear functions on norm of space Lipschitzian functions. The space of piecewise smooth functions is everywhere dense in space of regular smooth functions. The analogue of the equation of Euler for the elementary variational problem in space of regular smooth functions is proved.

Изучаются возможности аппроксимации произвольной кусочно-непрерывной функции на конечном отрезке линейной комбинацией $\mu$ функций Гаусса с целью дальнейшего их использования для аппроксимации управлений в сосредоточенных задачах оптимального управления. Напомним, что функция Гаусса (квадратичная экспонента) - это функция вида $\varphi(x)=\dfrac{1}{\sigma\sqrt{2\pi}} \exp\left[ -\dfrac{(x-m)^2}{2\sigma^2} \right]$. В отличие от исследований, проводившихся ранее другими авторами, рассматривается случай, когда параметры функций Гаусса (так же как и коэффициенты линейной комбинации) являются варьируемыми и подбираются, в частности, путем минимизации отклонения аппроксимации от аппроксимируемой функции либо (в том случае, когда речь идет об аппроксимации задачи оптимального управления) путем минимизации целевого функционала. Этот подход позволяет аппроксимировать задачи оптимального управления сосредоточенными системами конечномерными задачами математического программирования сравнительно небольшой размерности (в отличие от кусочно-постоянной или кусочно-линейной аппроксимации на фиксированной сетке с малым шагом, как это обычно делается). Приводятся результаты численных экспериментов, подтверждающие эффективность изучаемого подхода.

We study the opportunities of approximation of a piecewise continuous function on a finite segment by a linear combination of $\mu$ Gaussian functions, with the object of their usage for control approximation in lumped problems of optimal control. Recall that a Gaussian function (quadratic exponent) is one defined as follows $\varphi(x)=\dfrac{1}{\sigma\sqrt{2\pi}} \exp\left[ -\dfrac{(x-m)^2}{2\sigma^2} \right]$. Unlike investigations carried out by another authors, we consider the case where the parameters of a Gaussian function (with the coefficients of a linear combination) are varied and selected, in particular, by minimization of the difference between a function being approximated and its approximation, or (in the case of an optimal control problem) by minimization of the cost functional. Such an approach gives the opportunity to approximate optimal control problems with lumped parameters by finite dimensional problems of mathematical programming of comparatively small dimension (as opposed to piecewise constant or piecewise linear approximation on a fixed mesh with small width which is usually used). We present also some results of numerical experiments which substantiate efficiency of the approach under study.