Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
Известия Российской академии наук. Серия математическая
RUS  ENG    ЖУРНАЛЫ   ПЕРСОНАЛИИ   ОРГАНИЗАЦИИ   КОНФЕРЕНЦИИ   СЕМИНАРЫ   ВИДЕОТЕКА   ПАКЕТ AMSBIB  
Общая информация
Последний выпуск
Скоро в журнале
Архив
Импакт-фактор
Правила для авторов
Загрузить рукопись

Поиск публикаций
Поиск ссылок

RSS
Последний выпуск
Текущие выпуски
Архивные выпуски
Что такое RSS



Изв. РАН. Сер. матем.:
Год:
Том:
Выпуск:
Страница:
Найти






Персональный вход:
Логин:
Пароль:
Запомнить пароль
Войти
Забыли пароль?
Регистрация


Известия Российской академии наук. Серия математическая, 2022, том 86, выпуск 1, страницы 219–236
DOI: https://doi.org/10.4213/im9125
(Mi im9125)
 

Эта публикация цитируется в 2 научных статьях (всего в 2 статьях)

Экстремальная интерполяция с наименьшим значением нормы второй производной в пространстве Lp(R)

В. Т. Шевалдин

Институт математики и механики им. Н. Н. Красовского Уральского отделения РАН, г. Екатеринбург
Список литературы:
Аннотация: В работе в терминах разделенных разностей формулируется общая задача экстремальной функциональной интерполяции действительных функций одного переменного (для конечных разностей это задача Яненко–Стечкина–Субботина). Требуется вычислить наименьшее значение n-й производной в пространстве Lp(R), 1p, на классе функций, интерполирующих любую заданную бесконечную последовательность действительных чисел на произвольной, бесконечной в обе стороны сетке узлов на числовой оси R для класса интерполируемых последовательностей, у которых последовательность разделенных разностей n-го порядка принадлежит пространству lp(Z). В настоящей работе эта задача решается в случае n=2. Указанная величина оценивается сверху и снизу через наибольший и наименьший шаги сетки узлов.
Библиография: 12 наименований.
Ключевые слова: интерполяция, разделенная разность, сплайн, разностное уравнение.
Финансовая поддержка Номер гранта
Министерство науки и высшего образования Российской Федерации 075-02-2021-1383
Работа выполнена в рамках исследований, проводимых в Уральском математическом центре при финансовой поддержке Министерства науки и высшего образования Российской Федерации (номер соглашения 075-02-2021-1383).
Поступило в редакцию: 18.11.2020
Исправленный вариант: 06.12.2020
Англоязычная версия:
Izvestiya: Mathematics, 2022, Volume 86, Issue 1, Pages 203–219
DOI: https://doi.org/10.1070/IM9125
Реферативные базы данных:
Тип публикации: Статья
УДК: 519.65
MSC: Primary 41A05; Secondary 41A15, 41A50, 65D07

§ 1. Введение

Пусть на числовой оси R=(;+) задана бесконечная в обе стороны сетка узлов Δ={xk}kZ вида

<xk<xk+1<xk+2<,
и пусть hk=xk+1xk – шаги этой сетки, причем h_=inf, \overline{h}=\sup_k h_k. В настоящей работе всюду далее полагаем, что \underline{h}>0 и \overline{h}<+\infty. Отсюда, в частности, следует, что \lim_{k\to -\infty} x_k=-\infty, \lim_{k\to +\infty} x_k=+\infty. Для функции f, определенной на \mathbb R, положим
\begin{equation*} f(x_k)=y_k,\qquad k\in \mathbb Z, \end{equation*} \notag
где y=\{y_k\}_{k\in \mathbb Z} – произвольная последовательность действительных чисел. Разделенная разность порядка n для функции f на сетке \Delta определяется рекуррентно при помощи равенств
\begin{equation*} \begin{gathered} \, [y_k]=f[x_k]=y_k,\qquad [y_{k+1},y_k]=f[x_{k+1},x_k]=\frac{[y_{k+1}]-[y_k]}{x_{k+1}-x_k}, \\ \begin{split} [y_{k+2},y_{k+1},y_k]&=f[x_{k+2},x_{k+1},x_k] =\frac{[y_{k+2},y_{k+1}]-[y_{k+1},y_k]}{x_{k+2}-x_k},\quad\dots, \\ [y_{k+n},\dots,y_k]&=f[x_{k+n},\dots,x_k] =\frac{[y_{k+n},\dots,y_{k+1}]-[y_{k+n-1},\dots,y_k]}{x_{k+n}-x_k},\qquad k\in \mathbb Z. \end{split} \end{gathered} \end{equation*} \notag
Хорошо известно, что если исходная сетка узлов равномерная с шагом h>0 (т. е. x_m=x+mh, m\in \mathbb Z), то разделенная разность n-го порядка с точностью до постоянного множителя совпадает с обычной конечной разностью
\begin{equation*} \Delta_h^{n}f(x)=\sum_{m=0}^n (-1)^{n-m} C_n^m f(x+mh) \end{equation*} \notag
n-го порядка функции f с шагом h, а именно, имеет место равенство
\begin{equation*} f[x_{k+n},\dots,x_k]=\frac{1}{(n!)h^n}\Delta_h^{n}f(x_k). \end{equation*} \notag
Если исходная функция f дифференцируема n раз на оси \mathbb R, то имеют место простые формулы
\begin{equation*} \begin{aligned} \, f[x_{k+n},\dots,x_k]&=\frac{f^{(n)}(\xi)}{n!},\qquad \xi\in (x_k;x_{k+n}), \\ f[x_{k+n},\dots,x_k]&=\int_{D}f^{(n)}(t_0x_k+t_1x_{k+1}+\dots+t_nx_{k+n})\, dt_0\cdots dt_n, \end{aligned} \end{equation*} \notag
где D – область, удовлетворяющая следующим ограничениям:
\begin{equation*} t_0+t_1+\dots+t_n=1,\qquad t_i\geqslant 0,\qquad i=0,1,\dots,n, \end{equation*} \notag
связывающие разделенные разности и соответствующие производные (см., например, [1; гл. 1]). Более интересна другая формула (см. [2; с. 15]), в неявном виде она была впервые выписана Ж. Фаваром в 1940 г. в работе [3], а именно,
\begin{equation} f[x_{k+n},\dots,x_k]=\frac1{n!}\int_{x_k}^{x_{k+n}}M_{k,n}(t)f^{(n)}(t)\,dt, \end{equation} \tag{1.1}
где M_{k,n}(t) – нормализованный в L полиномиальный сплайн степени (n-1) с узлами x_k,\dots,x_{k+n}.

Пусть l_p=l_p(\mathbb Z) – пространство числовых последовательностей Z=\{Z_k\}_{k\in \mathbb Z} с нормой

\begin{equation*} \|Z\|_{l_p}=\|\{Z_k\}\|_{l_p}=\begin{cases} {\displaystyle \biggl( \sum_{k\in \mathbb Z}|Z_k|^p\biggr)^{1/p}}, &1\leqslant p<\infty, \\ {\displaystyle \sup_{k}|Z_k|}, &p=\infty. \end{cases} \end{equation*} \notag
Рассмотрим класс последовательностей следующего вида:
\begin{equation*} Y_{p,n}=\bigl\{ y=\{y_k\}_{k\in \mathbb Z}\colon \| \{[y_{k+n},\dots,y_k]\}\|_{l_p}\leqslant 1\bigr\}. \end{equation*} \notag
Для каждой последовательности y\in Y_{p,n} введем класс функций
\begin{equation*} F_{p,n}(y)=\{ f\colon f^{(n-1)}\in \mathrm{AC},\ f^{(n)}\in L_p(\mathbb R),\ f(x_k)=y_k,\ k\in \mathbb Z\}. \end{equation*} \notag
Здесь \mathrm{AC} – класс локально абсолютно непрерывных функций и L_p=L_p(\mathbb R), 1\leqslant p\leqslant \infty, – класс функций с обычным определением нормы
\begin{equation*} \|f\|_{L_p}=\|f\|_{L_p(\mathbb R)}=\begin{cases} {\displaystyle \biggl( \int_{\mathbb R}|f(t)|^p\,dt\biggr)^{1/p}}, & 1\leqslant p<\infty, \\ {\displaystyle \operatorname*{ess\,sup}_{t\in \mathbb R}|f(t)|}, &p=\infty. \end{cases} \end{equation*} \notag

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

\begin{equation} A_{p,n}(\Delta)=\sup_{y\in Y_{p,n}}\inf_{f\in F_{p,n}(y)}\|f^{(n)}\|_{L_p}. \end{equation} \tag{1.2}

Дифференциальный оператор взятия n-й производной и разностный оператор разделенной разности n-го порядка имеют одно и то же ядро – пространство алгебраических многочленов степени (n-1). Легко видеть, что если n-я производная некоторой действительной функции f ограничена сверху по модулю для всех x\in \mathbb R положительной константой M, то абсолютная величина разделенной разности n-го порядка на любой сетке узлов \Delta не превосходит числа M/(n!). Поэтому задачу о вычислении величин A_{p,n}(\Delta) можно считать обратной к отмеченному свойству разделенных разностей и более общей, чем известная задача Яненко–Стечкина–Субботина для конечных разностей. Она имеет богатую историю в случае равномерной сетки узлов (т. е. для конечных разностей), а для разделенных разностей оказалась гораздо труднее. Отметим, что случай n=1 является тривиальным, поскольку он сводится к локальной интерполяции ломаными.

Для конечных разностей вида \Delta_{h}^{n}f(x) эта задача для любой индивидуальной последовательности y=\{y_k\}_{k\in \mathbb Z} возникла в исследованиях академика Н. Н. Яненко при построении разностных методов решения дифференциальных уравнений. После бесед с Н. Н. Яненко профессор С. Б. Стечкин в начале 60-х гг. предложил своему ученику (ныне члену-корреспонденту РАН) Ю. Н. Субботину рассмотреть данную задачу в экстремальной постановке (1.2) для равномерной сетки узлов (обозначим ее через \overline{\Delta}) и всего класса последовательностей Y_{p,n}. К тому времени оценки сверху (неточные) для величины A_{p,n}(\overline{\Delta}) были получены В. С. Рябеньким (p=\infty) [4], [5] и С. Л. Соболевым [6]. Ю. Н. Субботин в нескольких своих работах (см. [7]–[9]) успешно справился с этой задачей, вычислив точно величину A_{n,p}(\overline{\Delta}) для равномерной сетки \overline{\Delta} при всех n\in \mathbb N и 1\leqslant p\leqslant \infty. Важным моментом решения было то, что экстремальными функциями в данной задаче оказались полиномиальные сплайны (с правильными узлами “склейки”) и их обобщения. Работы Ю. Н. Субботина дали мощный толчок к исследованию аппроксимативных и экстремальных свойств этих кусочно-полиномиальных функций. Многочисленные обобщения и применения результатов Ю. Н. Субботина изложены в обзорной статье [10].

Близкой по постановке к задаче (1.2) является интерполяционная задача Ж. Фавара [3] (полную библиографию и некоторые обобщения см., например, в [11]). Ж. Фавар рассматривал только случай p=\infty и равенства (1.1) для конечного набора чисел k, а именно, для k=0,1,\dots,m-n (m – любое фиксированное натуральное число, большее n) (см. [3]) и любых (m+1) точек x_0<x_1<\dots<x_m, в которых заданы значения y_0=f(x_0), \dots, y_m=f(x_m) некоторой функции f. Пусть имеют место неравенства

\begin{equation*} |f[x_{k+n},\dots,x_k]|\leqslant M,\qquad k=0,1,\dots,m-n, \end{equation*} \notag
где M – фиксированное положительное число. При таких ограничениях Ж. Фавар рассмотрел задачу нахождения функции f на интервале (x_0;x_m), у которой n-я производная в пространстве L_{\infty} является наименьшей.

В свою очередь, при формулировке задачи (1.2) мы предполагаем, что равенства (1.1) выполнены для всех целых значений параметра k, и при p=\infty все эти разделенные разности ограничены (см. определение класса Y_{\infty,n}), и минимизируем n-ю производную не локально, а на всей оси \mathbb R.

В недавней работе С. И. Новикова и автора [12] нам удалось продвинуться в решении задачи (1.2) при n=2 и p=\infty для произвольной сетки узлов \Delta (не требуя ограничений вида \underline{h}>0 и \overline{h}<+\infty). А именно, из результатов работы [12], в частности, следует неравенство

\begin{equation} 2\leqslant A_{\infty,2}(\Delta)\leqslant 18. \end{equation} \tag{1.3}
Видимо, задачу (1.2) в общей постановке (т. е. для разделенных разностей на произвольной сетке узлов) больше никто не рассматривал.

Основной результат настоящей работы состоит в следующем.

Теорема. Пусть 1\leqslant p\leqslant \infty и \underline{h}=\inf_k h_k>0, \overline{h}=\sup_k h_k<+\infty. Имеет место двойное неравенство

\begin{equation*} C_1(p)\underline{h}^{1/p}\leqslant A_{p,2}(\Delta)\leqslant C_2(p)\overline{h}^{1/p}, \end{equation*} \notag
где
\begin{equation*} \begin{aligned} \, C_1(p)&=\begin{cases} \biggl(\dfrac{2p-1}{p-1}\biggr)^{(p-1)/p}2^{1/p}, &1<p<\infty, \\ 2, &p=1, \\ 2, &p=\infty, \end{cases} \\ C_2(p)&=\begin{cases} \dfrac{18p}{p-1} \biggl(\dfrac{2p-1}{p-1}\biggr)^{(p-1)/p}2^{1/p}, &1<p<\infty, \\ 2, &p=1, \\ 4, &p=\infty. \end{cases} \end{aligned} \end{equation*} \notag

Отметим, что Ю. Н. Субботин [7]–[9], решая задачу (1.2) для равномерной сетки узлов \overline{\Delta}=\{kh\}_{k\in \mathbb Z}, рассматривал только случай h=1. Но если в его работах сделать соответствующую замену переменных и рассмотреть равномерную сетку с произвольным шагом h, то в случае n=2 основной результат Ю. Н. Субботина (с учетом связи между конечными и разделенными разностями) может быть переформулирован следующим образом:

\begin{equation} A_{p,2}(\overline{\Delta})=\begin{cases} 2\biggl(\dfrac{2p-1}{p-1}\biggr)^{(p-1)/p}h^{1/p}, &1<p<\infty, \\ 2h, &p=1, \\ 4, &p=\infty. \end{cases} \end{equation} \tag{1.4}

Структура настоящей работы следующая. В § 2 и § 3 будут получены соответственно оценки снизу и сверху для величины A_{p,2}(\Delta) при 1<p<\infty. В § 4 рассмотрен случай p=1, и, наконец, в § 5 уточнена оценка сверху в неравенстве (1.3).

§ 2. Оценка снизу величины A_{p,2}(\Delta)

Функцию f\in F_{p,2}(y), 1\leqslant p\leqslant \infty, с помощью формулы Тейлора запишем в виде

\begin{equation*} f(x)=f(x_k)+f'(x_k)(x-x_k)+\int_{x_k}^{x}(x-t)f''(t)\,dt. \end{equation*} \notag
Отсюда легко получить (см. [12]) следующее представление разделенной разности второго порядка:
\begin{equation} [y_{k+2},y_{k+1},y_k]=\frac{1}{h_k+h_{k+1}}\biggl[ \int_{x_{k+1}}^{x_{k+2}}\frac{x_{k+2}-t}{h_{k+1}}f''(t)\,dt+ \int_{x_{k}}^{x_{k+1}}\frac{t-x_{k}}{h_{k}}f''(t)\,dt\biggr]. \end{equation} \tag{2.1}
Заметим, что равенство (2.1) является частным случаем равенства (1.1) при n=2.

Пусть числа p и q удовлетворяют равенству

\begin{equation*} \frac{1}{p}+\frac{1}{q}=1 \end{equation*} \notag
(при p=1 полагаем q=\infty, а при p=\infty считаем, что q=1). При получении оценки снизу величины A_{p,2}, 1\leqslant p<\infty, будем действовать методом Ю. Н. Субботина [7]–[9]. Рассмотрим произвольную последовательность y^*=\{ y_k^*\}_{k\in \mathbb Z}\in Y_{p,2}, удовлетворяющую условию
\begin{equation} [y_{k+2}^*,y_{k+1}^*,y_k^*]=\begin{cases} (-1)^k(2N+1)^{-1/p}, &|k|\leqslant N, \\ 0, &|k|> N, \end{cases} \end{equation} \tag{2.2}
где N – произвольное натуральное число, N\geqslant 2. Принимая во внимание (2.1) и (2.2), имеем
\begin{equation*} \begin{aligned} \, &(2N+1)^{1/q}=(2N+1)(2N+1)^{-1/p}=\sum_{k}(-1)^k[y_{k+2}^*,y_{k+1}^*,y_k^*] \\ &\quad =\sum_{k=-N}^N\frac{(-1)^k}{h_k+h_{k+1}}\biggl[ \int_{x_{k+1}}^{x_{k+2}}\frac{x_{k+2}-t}{h_{k+1}}f''(t)\,dt+ \int_{x_{k}}^{x_{k+1}}\frac{t-x_{k}}{h_{k}}f''(t)\,dt\biggr]=A+B, \end{aligned} \end{equation*} \notag
где
\begin{equation} \begin{aligned} \, A &=\frac{(-1)^N}{h_N+h_{N+1}}\int_{x_{N+1}}^{x_{N+2}}\frac{x_{N+2}-t}{h_{N+1}}f''(t)\,dt \\ &\qquad + \frac{(-1)^{-N}}{h_{-N}+h_{-N+1}} \int_{x_{-N}}^{x_{-N+1}}\frac{t-x_{-N}}{h_{-N}}f''(t)\,dt, \\ B&=\sum_{k=-N+1}^{N}(-1)^{k-1}\int_{x_k}^{x_{k+1}}\varphi_k(t)f''(t)\,dt. \end{aligned} \end{equation} \tag{2.3}
Функции
\begin{equation} \varphi_k(t)=\frac{x_{k+1}-t}{h_k(h_{k-1}+h_k)}-\frac{t-x_k}{h_k(h_k+h_{k+1})},\qquad t\in [x_k;x_{k+1}],\quad k=-N+1,\dots, N, \end{equation} \tag{2.4}
играют важную роль в дальнейших рассуждениях. При любом k функция \varphi_k(t) является линейной на отрезке [x_k;x_{k+1}], причем имеют место равенства
\begin{equation} \varphi_k(x_k)=\frac{1}{h_{k-1}+h_k},\qquad \varphi_k(x_{k+1})=-\frac{1}{h_{k}+h_{k+1}}. \end{equation} \tag{2.5}
Кроме того, \varphi_k(t_k)=0 при
\begin{equation} t_k=\frac{x_k(h_{k-1}+h_k)+x_{k+1}(h_k+h_{k+1})}{h_{k-1}+2h_k+h_{k+1}}\in (x_k;x_{k+1}). \end{equation} \tag{2.6}
При 1<p<\infty оценим абсолютные величины чисел A и B (см. (2.3)). В силу неравенства Гёльдера получим
\begin{equation*} \begin{aligned} \, |A| &\leqslant \|f''\|_{L_p(\mathbb R)}\biggl\{ \frac{1}{h_{N+1}(h_N+h_{N+1})}\|x_{N+2}-t\|_{L_q[x_{N+1};x_{N+2}]} \\ &\qquad+\frac{1}{h_{-N}(h_{-N}+h_{-N+1})}\|t-x_{-N}\|_{L_q[x_{-N};x_{-N+1}]}\biggr\} \\ &= \frac{\|f''\|_{L_p(\mathbb R)}}{(q+1)^{1/q}}\biggl\{ \frac{h_{N+1}^{1/q}}{h_N+h_{N+1}}+\frac{h_{-N}^{1/q}}{h_{-N}+h_{-N+1}}\biggr\}. \end{aligned} \end{equation*} \notag
Поскольку 1<p<\infty, то q>1 и из выписанной оценки следует неравенство
\begin{equation} |A|\leqslant \dfrac{h^{1/q-1}_{N+1}+h^{1/q-1}_{-N}}{(q+1)^{1/q}}\|f''\|_{L_p(\mathbb R)} \leqslant 2\max\{1,\underline{h}^{-1/p}\}\|f''\|_{L_p(\mathbb R)},\qquad 1<p<\infty. \end{equation} \tag{2.7}
Из (2.3) имеем
\begin{equation} |B|\leqslant \sum_{k=-N+1}^N\int_{x_k}^{x_{k+1}}|\varphi_k(t)|\cdot |f''(t)|\,dt=\overline{B}+\overline{\overline{B}}, \end{equation} \tag{2.8}
где
\begin{equation} \begin{aligned} \, \overline{B}&=\int_{t_N}^{x_{N+1}}|\varphi_N(t)|\cdot |f''(t)|\,dt+ \int_{x_{-N+1}}^{t_{-N+1}}|\varphi_{-N+1}(t)|\cdot |f''(t)|\,dt, \\ \overline{\overline{B}}&=\int_{t_{-N+1}}^{t_{N}}\psi(t) |f''(t)|\,dt. \end{aligned} \end{equation} \tag{2.9}
Здесь
\begin{equation} \psi(t)=\begin{cases} -\varphi_k(t), & t\in [t_k;x_{k+1}], \\ \varphi_{k+1}(t), &t\in [x_{k+1};t_{k+1}], \end{cases} \qquad k=-N+1,\dots,N-1. \end{equation} \tag{2.10}
Из (2.5) следует, что \varphi_{k+1}(x_{k+1})=-\varphi_{k}(x_{k+1}), и поэтому функция \psi(t)\geqslant 0 является непрерывной функцией на отрезке [t_{-N+1};t_N]. График этой функции похож на “пилу” с вершинами в точках x_{-N+2},\dots, x_N и нулями в точках t_{-N+1},\dots, t_N (см. (2.6)). Оценим сверху величины \overline{B} и \overline{\overline{B}} при 1<p<\infty с помощью неравенства Гёльдера. Получим
\begin{equation} \overline{B}\leqslant 2\max\{1,\underline{h}^{-1/p}\}\|f''\|_{L_p(\mathbb R)},\qquad \overline{\overline{B}}\leqslant \|f''\|_{L_p(\mathbb R)}\|\psi\|_{L_q[t_{-N+1};t_N]}. \end{equation} \tag{2.11}
Значит, из (2.7), (2.8) и (2.11) имеем
\begin{equation*} (2N+1)^{1/q}\leqslant |A|+|B|\leqslant \|f''\|_{L_p(\mathbb R)}\bigl(\|\psi\|_{L_q[t_{-N+1};t_N]}+O(1)\bigr). \end{equation*} \notag
Поэтому для любого натурального числа N\geqslant 2 имеет место оценка
\begin{equation*} \|f''\|_{L_p(\mathbb R)}\geqslant \frac{(2N+1)^{1/q}}{\|\psi\|_{L_q[t_{-N+1};t_N]}+O(1)},\qquad 1<p<\infty. \end{equation*} \notag
Переходя в этом неравенстве к пределу при N\to \infty, получаем следующую оценку величины A_{p,2}(\Delta), 1<p<\infty:
\begin{equation} A_{p,2}(\Delta)\geqslant \varlimsup_{N\to \infty}\biggl(\frac{1}{2N}\int_{t_{-N+1}}^{t_N}\psi^q(t)\,dt \biggr)^{-1/q}. \end{equation} \tag{2.12}
В этом неравенстве присутствует функция \psi(t) (см. (2.4) и (2.10)). Пересчитаем интеграл в правой части неравенства (2.12) в терминах h_k – длин шагов сетки \Delta, используя то, что функция \psi кусочно линейна. Из (2.4) с учетом (2.6) имеем
\begin{equation} \begin{aligned} \, Q_k&=\int_{t_k}^{t_{k+1}}\psi^q(t)\,dt=\int_{t_k}^{x_{k+1}}\biggl[ \frac{t-x_k}{h_k(h_k+h_{k+1})}-\frac{x_{k+1}-t}{h_k(h_{k-1}+h_k)}\biggr]^q\,dt \\ &\qquad+\int_{x_{k+1}}^{t_{k+1}}\biggl[ \frac{x_{k+2}-t}{h_{k+1}(h_{k}+h_{k+1})}-\frac{t-x_{k+1}}{h_{k+1}(h_{k+1}+h_{k+2})}\biggr]^q\,dt. \end{aligned} \end{equation} \tag{2.13}
Вычисляя интегралы в (2.13) с помощью линейных замен переменных, получаем
\begin{equation} Q_k=\frac{1}{(q+1)(h_k+h_{k+1})^q}\biggl[ \frac{h_k(h_{k-1}+h_{k})}{h_{k-1}+2h_k+h_{k+1}} +\frac{h_{k+1}(h_{k+1}+h_{k+2})}{h_k+2h_{k+1}+h_{k+2}}\biggr]. \end{equation} \tag{2.14}
Таким образом, с учетом (2.13), (2.14) оценка (2.12) принимает вид
\begin{equation} A_{p,2}(\Delta)\geqslant \varlimsup_{N\to \infty}\biggl( \frac{1}{2N} \sum_{k=-N+1}^{N-1} Q_k\biggr)^{-1/q}, \qquad 1<p<\infty. \end{equation} \tag{2.15}
В дальнейшем мы загрубим полученную оценку с учетом условий теоремы (мы рассматриваем только сетки \Delta, удовлетворяющие условиям \underline{h}>0 и \overline{h}<+\infty). Но интересно отметить, что оценка (2.15) величины A_{p,2}(\Delta) является точной для равномерной сетки узлов \overline{\Delta}, т. е. в случае h_k=h. В этом случае для любого числа k имеет место равенство
\begin{equation*} Q^{-1/q}_k=2(q+1)^{1/q}h^{(q-1)/q}=2\biggl(\frac{2p-1}{p-1}\biggr)^{(p-1)/p}h^{1/p}, \end{equation*} \notag
и поэтому
\begin{equation*} A_{p,2}(\overline{\Delta})\geqslant 2\biggl(\frac{2p-1}{p-1}\biggr)^{(p-1)/p}h^{1/p} \end{equation*} \notag
(см. равенство (1.4)).

Из (2.14) следует, что

\begin{equation*} \begin{aligned} \, Q_k&=\frac{1}{(q+1)(h_k+h_{k+1})^{q-1}}\biggl[ \frac{1+{h_{k-1}}/{h_k}}{(1+{h_{k+1}}/{h_k})({h_{k-1}}/{h_k}+2+{h_{k+1}}/{h_k})} \\ &\qquad+\frac{1+{h_{k+2}}/{h_{k+1}}}{(1+{h_{k}}/{h_{k+1}})({h_{k}}/{h_{k+1}}+2+ {h_{k+2}}/{h_{k+1}})}\biggr] \\ &<\frac{1}{(q+1)(h_k+h_{k+1})^{q-1}}\biggl[\frac{1}{1+{h_{k+1}}/{h_k}} +\frac{1}{1+{h_k}/{h_{k+1}}} \biggr] \\ &=\frac{1}{(q+1)(h_k+h_{k+1})^{q-1}}\leqslant \frac{1}{(q+1)2^{q-1}\underline{h}^{q-1}}, \end{aligned} \end{equation*} \notag
поскольку \underline{h}=\inf_k h_k. Поэтому из (2.15) получим
\begin{equation} A_{p,2}(\Delta)\geqslant \varlimsup_{N\to \infty}\biggl( \frac{2N}{2N-1} (q+1)2^{q-1}\underline{h}^{q-1} \biggr)^{1/q}=(q+1)^{1/q}2^{1/p}\underline{h}^{1/p},\qquad 1<p<\infty. \end{equation} \tag{2.16}
Напомним, что числа p и q связаны соотношением 1/p\,{+}\,1/q=1. Таким образом, при 1<p<\infty в теореме, сформулированной во введении, доказана оценка снизу с константой
\begin{equation*} C_1(p)=\biggl(\frac{2p-1}{p-1}\biggr)^{(p-1)/p}2^{1/p}. \end{equation*} \notag

§ 3. Оценка сверху величины A_{p,2}(\Delta)

Пусть 1<p<\infty. Покажем, что для любой последовательности y\in Y_{p,2} существует функция f\in F_{p,2}(y) – обобщенный параболический сплайн с узлами \{ t_k\}_{k\in \mathbb Z} (см. (2.6)), причем f(t_k)=0, k\in \mathbb Z, который интерполирует значения последовательности y в точках \{x_k\}_{k\in \mathbb Z}, т. е. f(x_k)=y_k, k\in \mathbb Z, и для него справедлива оценка

\begin{equation*} \|f''\|_{L_p(\mathbb R)}\leqslant C_2(p)\overline{h}^{1/p} \end{equation*} \notag
с некоторой положительной константой C_2(p), где \overline{h}=\sup_k h_k. Построим этот сплайн с помощью введенной в § 2 функции \psi(t) (см. (2.4) и (2.10)). Для любого числа k\in \mathbb Z полагаем
\begin{equation} f''(t)=\begin{cases} Z_k\psi^{q-1}(t), &x_k\leqslant t<t_k, \\ Z_{k+1}\psi^{q-1}(t), &t_k\leqslant t<t_{k+1}, \\ Z_{k+2}\psi^{q-1}(t), &t_{k+1}\leqslant t<x_{k+2}, \end{cases} \qquad k\in \mathbb Z. \end{equation} \tag{3.1}
Числа \{Z_k\}_{k\in \mathbb Z} подлежат дальнейшему определению. Поскольку f(x_k)=y_k, то
\begin{equation*} f[x_{k+2},x_{k+1},x_k]=[y_{k+2},y_{k+1},y_k], \end{equation*} \notag
и можно воспользоваться представлением (2.1). Подставляя (3.1) в равенство (2.1), для определения чисел \{Z_k\}_{k\in \mathbb Z} получаем разностное уравнение
\begin{equation} A_k Z_{k+2}+B_k Z_{k+1}+C_kZ_k=[y_{k+2},y_{k+1},y_k],\qquad k\in \mathbb Z, \end{equation} \tag{3.2}
где
\begin{equation} \begin{aligned} \, A_k&=\frac{1}{h_k+h_{k+1}}\int_{t_{k+1}}^{x_{k+2}}\frac{x_{k+2}-t}{h_{k+1}}\psi^{q-1}(t)\,dt, \\ B_k&=\frac{1}{h_k+h_{k+1}}\biggl[ \int_{x_{k+1}}^{t_{k+1}}\frac{x_{k+2}-t}{h_{k+1}}\psi^{q-1}(t)\,dt+ \int_{t_{k}}^{x_{k+1}}\frac{t-x_{k}}{h_{k}}\psi^{q-1}(t)\,dt\biggr], \\ C_k&=\frac{1}{h_k+h_{k+1}}\int_{x_{k}}^{t_{k}}\frac{t-x_{k}}{h_{k}}\psi^{q-1}(t)\,dt. \end{aligned} \end{equation} \tag{3.3}
Вычисляя интегралы (3.3) с помощью замен переменных, после элементарных преобразований получаем
\begin{equation*} \begin{aligned} \, A_k &=\frac{h_{k+1}(h_k+h_{k+1})}{q(q+1)(h_{k+1}+h_{k+2})^{q-1}(h_k+2h_{k+1}+h_{k+2})^2}, \\ B_k &=\frac{h_{k}(h_{k-1}+h_{k})^2}{(h_{k}+h_{k+1})^{q-1}(h_{k-1}+2h_{k}+h_{k+1})^2} \biggl[\frac{1}{(q+1)(h_k+h_{k+1})}+\frac{1}{q(h_{k-1}+h_k)} \biggr] \\ &\ +\frac{h_{k+1}(h_{k+1}\,{+}\,h_{k+2})^2}{(h_{k}\,{+}\,h_{k+1})^{q-1}(h_{k}\,{+}\,2h_{k+1}\,{+}\,h_{k+2})^2} \biggl[\frac{1}{(q\,{+}\,1)(h_k\,{+}\,h_{k+1})}\,{+}\,\frac{1}{q(h_{k+1}\,{+}\,h_{k+2})} \biggr], \\ C_k&=\frac{h_{k}(h_{k}+h_{k+1})}{q(q+1)(h_{k-1}+h_{k})^{q-1}(h_{k-1}+2h_{k}+h_{k+1})^2}. \end{aligned} \end{equation*} \notag
С помощью замены переменных
\begin{equation} \widetilde{Z}_k=\frac{Z_k}{q(q+1)(h_{k-1}+h_k)^{q-1}},\qquad k\in \mathbb Z, \end{equation} \tag{3.4}
упростим разностное уравнение (3.2):
\begin{equation} a_k \widetilde{Z}_{k+2}+b_k\widetilde{Z}_{k+1}+c_k\widetilde{Z}_{k}=[y_{k+2},y_{k+1},y_k], \end{equation} \tag{3.5}
где
\begin{equation} \begin{gathered} \, a_k=\frac{h_{k+1}(h_k+h_{k+1})}{(h_k+2h_{k+1}+h_{k+2})^2},\qquad c_k=\frac{h_{k}(h_k+h_{k+1})}{(h_{k-1}+2h_{k}+h_{k+1})^2}, \\ \begin{aligned} \, b_k=b_k(q)&=\frac{h_{k}(h_{k-1}+h_{k})}{(h_{k-1}+2h_{k}+h_{k+1})^2} \biggl[\frac{q}{h_k+h_{k+1}}+\frac{q+1}{h_{k-1}+h_k} \biggr] \\ &\qquad+\frac{h_{k+1}(h_{k+1}+h_{k+2})}{(h_k+2h_{k+1}+h_{k+2})^2} \biggl[\frac{q}{h_k+h_{k+1}}+\frac{q+1}{h_{k+1}+h_{k+2}} \biggr]. \end{aligned} \end{gathered} \end{equation} \tag{3.6}
С помощью элементарных преобразований (3.6) нетрудно проверить равенство
\begin{equation} a_k+b_k(1)+c_k=1. \end{equation} \tag{3.7}
Свойства чисел a_k, b_k(1) и c_k были исследованы в совместной работе С. И. Новикова и автора [12; лемма 1 и следствие к этой лемме]. С учетом (3.7) справедливо следующее утверждение.

Лемма 3.1 (см. [12]). При любом k\in \mathbb Z имеют место неравенства

\begin{equation*} \begin{gathered} \, 0<a_k<\frac14,\qquad 0<c_k<\frac14,\qquad a_k+c_k<\frac49, \\ \frac59< b_k(1)<1,\qquad \frac{1}{9}<b_k(1)-a_k-c_k<1. \end{gathered} \end{equation*} \notag

Из определения чисел b_k(q) легко видеть, что b_k(q)\geqslant b_k(1) при q\geqslant 1, и поэтому из (3.7) следует неравенство

\begin{equation} a_k+b_k(q)+c_k\geqslant 1,\qquad q\geqslant 1,\quad k\in \mathbb Z. \end{equation} \tag{3.8}

Лемма 3.2. Для любой последовательности y\in Y_{p,2}, 1<p<\infty, разностное уравнение (3.5) имеет решение \widetilde{Z}=\{\widetilde{Z}_k\}_{k\in \mathbb Z}\in l_p, и это решение единственно.

Доказательство. Разностное уравнение (3.5) перепишем в следующем виде:
\begin{equation*} \widetilde{Z}_{k+1}=\frac{1}{a_k+b_k(q)+c_k}\bigl( [y_{k+2},y_{k+1},y_k]+a_k\bigl(\widetilde{Z}_{k+1}-\widetilde{Z}_{k+2}\bigr) +c_k\bigl(\widetilde{Z}_{k+1} -\widetilde{Z}_{k}\bigr)\bigr). \end{equation*} \notag
Рассмотрим нелинейный оператор T, который ставит в соответствие произвольной последовательности \widetilde{Z}=\{\widetilde{Z}_{k+1}\}_{k\in \mathbb Z}\in l_p последовательность
\begin{equation*} \biggl\{ \frac{1}{a_k+b_k(q)+c_k}\bigl( [y_{k+2},y_{k+1},y_k]+a_k\bigl(\widetilde{Z}_{k+1}-\widetilde{Z}_{k+2}\bigr)+ c_k\bigl(\widetilde{Z}_{k+1}-\widetilde{Z}_{k}\bigr) \bigr)\biggr\}_{k\in \mathbb Z}\in l_p. \end{equation*} \notag
В силу леммы 3.1 и неравенства (3.8) имеем
\begin{equation*} \begin{aligned} \, &\bigl\| T\widetilde{Z}^{(1)}-T\widetilde{Z}^{(2)}\bigr\|_{l_p}= \biggl\| \biggl\{ \frac{a_k+c_k}{a_k+b_k(q)+c_k}\bigl(\widetilde{Z}_{k+1}^{(1)} -\widetilde{Z}_{k+1}^{(2)}\bigr) \\ &\ \quad +\frac{a_k}{a_k+b_k(q)+c_k}\bigl(\widetilde{Z}_{k+2}^{(2)} -\widetilde{Z}_{k+2}^{(1)}\bigr)+\frac{c_k}{a_k+b_k(q)+c_k}\bigl(\widetilde{Z}_{k}^{(2)} -\widetilde{Z}_{k}^{(1)}\bigr)\biggr\}\biggr\|_{l_p} \\ &\ \leqslant \biggl\| \biggl\{ \frac{a_k+c_k}{a_k+b_k(q)+c_k} \bigl(\widetilde{Z}_{k+1}^{(1)}-\widetilde{Z}_{k+1}^{(2)}\bigr)\biggr\}\biggr\|_{l_p}{+}\, \biggl\|\biggl\{\frac{a_k}{a_k+b_k(q)+c_k} \bigl(\widetilde{Z}_{k+2}^{(2)}-\widetilde{Z}_{k+2}^{(1)}\bigr)\biggr\}\biggr\|_{l_p} \\ &\ \quad +\biggl\|\biggl\{\frac{c_k}{a_k+b_k(q)+c_k} \bigl(\widetilde{Z}_{k}^{(2)}-\widetilde{Z}_{k}^{(1)}\bigr)\biggr\}\biggr\|_{l_p} \\ &\ \leqslant\biggl( \frac49+\frac14+\frac14\biggr) \bigl\|\widetilde{Z}^{(1)}-\widetilde{Z}^{(2)}\bigr\|_{l_p} =\frac{17}{18}\bigl\|\widetilde{Z}^{(1)}-\widetilde{Z}^{(2)}\bigr\|_{l_p}. \end{aligned} \end{equation*} \notag
Значит, оператор T является сжимающим оператором в полном метрическом пространстве l_p=l_p(\mathbb Z) с константой сжатия 17/18<1. Здесь \widetilde{Z}^{(1)}=\{\widetilde{Z}_{k+1}^{(1)}\}_{k\in \mathbb Z} и \widetilde{Z}^{(2)}=\{\widetilde{Z}_{k+1}^{(2)}\}_{k\in \mathbb Z}. Поэтому согласно теореме о сжимающем операторе уравнение T\widetilde{Z}=\widetilde{Z} (т. е. разностное уравнение (3.5)) имеет решение \widetilde{Z}\in l_p, и это решение единственно. Лемма 3.2 доказана.

Лемма 3.3. Для решения \widetilde{Z}=\{\widetilde{Z}_k\}_{k\in \mathbb Z} разностного уравнения (3.5) справедлива оценка

\begin{equation*} \bigl\|\widetilde{Z}\bigr\|_{l_p}\leqslant 18,\qquad 1<p<\infty. \end{equation*} \notag

Доказательство. В лемме 3.2 было доказано, что разностное уравнение (3.5) имеет единственное решение \widetilde{Z}=\{\widetilde{Z}_{k+1}\}\in l_p, но при этом не был найден явный вид этого решения. Для доказательства леммы 3.3 будем рассуждать от противного. Пусть \|\widetilde{Z}\|_{l_p}>18. Поскольку y\in Y_{p,2}, то
\begin{equation} \begin{aligned} \, 1&\geqslant \|\{ [y_{k+2},y_{k+1},y_k]\}\|_{l_p}=\bigl\| \bigl\{ a_k\widetilde{Z}_{k+2}+b_k(q)\widetilde{Z}_{k+1}+c_k\widetilde{Z}_k \bigr\}\bigr\|_{l_p} \\ &\geqslant \bigl|\, \bigl\| \bigl\{ b_k(q)\widetilde{Z}_{k+1}\bigr\}\bigr\|_{l_p}- \bigl\| \bigl\{ a_k\widetilde{Z}_{k+2}+c_k\widetilde{Z}_{k}\bigr\}\bigr\|_{l_p} \bigr|. \end{aligned} \end{equation} \tag{3.9}
Покажем, что выражение под знаком модуля в последнем неравенстве больше 1. В самом деле, в силу леммы 3.1 имеем
\begin{equation*} \begin{aligned} \, &\bigl\| \bigl\{ b_k(q)\widetilde{Z}_{k+1}\bigr\} \bigr\|_{l_p}- \bigl\| \bigl\{ a_k\widetilde{Z}_{k+2}+c_k\widetilde{Z}_{k}\bigr\}\bigr\|_{l_p}\geqslant \frac59 \bigl\| \bigl\{ \widetilde{Z}_{k+1}\bigr\} \bigr\|_{l_p} \\ &\qquad-\bigl\| \bigl\{ a_k\widetilde{Z}_{k+2}\bigr\} \bigr\|_{l_p}- \bigl\| \bigl\{ c_k\widetilde{Z}_{k}\bigr\} \bigr\|_{l_p}\geqslant \biggl(\frac59-\frac14-\frac14\biggr) \bigl\|\widetilde{Z}\bigr\|_{l_p}=\frac{1}{18}\bigl\|\widetilde{Z}\bigr\|_{l_p}>1, \end{aligned} \end{equation*} \notag
что противоречит неравенству (3.9). Лемма 3.3 доказана.

Оценим теперь норму второй производной в пространстве L_p(\mathbb R) построенной функции f\in F_{p,2}(y). Из (3.1) имеем

\begin{equation} \begin{aligned} \, \|f''\|_{L_p(\mathbb R)}&=\biggl( \sum_{k\in \mathbb Z}\int_{t_k}^{t_{k+1}}|Z_{k+1}|^p \psi^{(q-1)p}(t)\,dt\biggr)^{1/p} \\ &=\biggl( \sum_{k\in \mathbb Z}|Z_{k+1}|^p\int_{t_k}^{t_{k+1}} \psi^{q}(t)\,dt\biggr)^{1/p}. \end{aligned} \end{equation} \tag{3.10}
Принимая во внимание (2.13), (2.14) и (3.4), из (3.10) получаем
\begin{equation*} \begin{aligned} \, \|f''\|_{L_p(\mathbb R)} &=\biggl( \sum_{k\in \mathbb Z}\bigl|\widetilde{Z}_{k+1}\bigr|^p (h_k+h_{k+1})^{(q-1)p}q^p(q+1)^p \frac{M_k}{(q+1)(h_k+h_{k+1})^q}\biggr)^{1/p} \\ &=\frac{q(q+1)}{(q+1)^{1/p}} \biggl( \sum_{k\in \mathbb Z}\bigl|\widetilde{Z}_{k+1}\bigr|^p M_k\biggr)^{1/p}, \end{aligned} \end{equation*} \notag
где
\begin{equation*} M_k=\frac{h_k(h_{k-1}+h_k)}{h_{k-1}+2h_k+h_{k+1}}+ \frac{h_{k+1}(h_{k+1}+h_{k+2})}{h_{k}+2h_{k+1}+h_{k+2}}<h_k+h_{k+1}\leqslant 2\overline{h}. \end{equation*} \notag
Отсюда с учетом леммы 3.3 получаем неравенство
\begin{equation*} \|f''\|_{L_p(\mathbb R)}\leqslant 18\cdot 2^{1/p}q(q+1)^{(p-1)/p}\,\overline{h}^{1/p}. \end{equation*} \notag
Переходя от переменной q к переменной p, из доказанного неравенства выводим, что
\begin{equation*} A_{p,2}(\Delta)\leqslant C_2(p)\,\overline{h}^{1/p},\qquad 1<p<\infty, \end{equation*} \notag
где
\begin{equation*} C_2(p)=2^{1/p}\frac{18p}{p-1}\biggl(\frac{2p-1}{p-1}\biggr)^{(p-1)/p}. \end{equation*} \notag
Теорема, сформулированная во введении, доказана при 1<p<\infty.

§ 4. Случай p=1

Для оценки снизу величины A_{1,2}(\Delta) рассмотрим последовательность y^*=\{y_k^*\}_{k\in \mathbb Z}, удовлетворяющую условию

\begin{equation} [y_{k+2}^*,y_{k+1}^*,y_k^*]=\begin{cases} 1, &k=0, \\ 0, &k\ne 0, \end{cases} \end{equation} \tag{4.1}
и покажем, что для любой функции f\in F_{1,2}(y^*) имеет место неравенство
\begin{equation} \|f''\|_{L_1(\mathbb R)}\geqslant 2\underline{h}. \end{equation} \tag{4.2}
Пусть N – произвольное натуральное число, N\geqslant 2. Из (2.1) и (4.1) имеем
\begin{equation} \begin{aligned} \, 1&=\sum_{k=-N}^N (-1)^k [y_{k+2}^*,y_{k+1}^*,y_k^*] \\ &=\sum_{k=-N}^{N}\frac{(-1)^k}{h_k+h_{k+1}}\biggl[ \int_{x_{k+1}}^{x_{k+2}}\frac{x_{k+2}-t}{h_{k+1}}f''(t)\,dt+ \int_{x_{k}}^{x_{k+1}}\frac{t-x_{k}}{h_{k}}f''(t)\,dt\biggr]. \end{aligned} \end{equation} \tag{4.3}
Для любой функции f\in F_{1,2}(y^*), как и в случае 1<p<\infty, правую часть равенства (4.3) представим в виде A+B, где A и B определены равенствами (2.3). Заметим, что в силу условий теоремы \underline{h}>0 и \overline{h}<+\infty, и поэтому если f''\in L_1(\mathbb R), то при N\to \infty
\begin{equation*} \int_{x_{N+1}}^{x_{N+2}}|f''(t)|\,dt\to 0,\qquad \int_{x_{-N}}^{x_{-N+1}}|f''(t)|\,dt\to 0. \end{equation*} \notag
В силу этого обстоятельства легко видеть, что A=o(1) и аналогично \overline{B}=o(1) при N\to \infty (см. (2.8), (2.9)). Из равенств (2.4), (2.5) и (2.10) имеем
\begin{equation*} \psi(t)\leqslant \sup_{k}(h_k+h_{k+1})^{-1}\leqslant (2\underline{h})^{-1}. \end{equation*} \notag
Поэтому величину
\begin{equation*} \overline{\overline{B}}=\int_{t_{-N+1}}^{t_{N}}\psi(t) |f''(t)|\,dt \end{equation*} \notag
(см. (2.9)) можно оценить следующим образом:
\begin{equation*} \overline{\overline{B}}\leqslant (2\underline{h})^{-1}\|f''\|_{L_1(\mathbb R)}. \end{equation*} \notag
Тогда из равенства (4.3) получим
\begin{equation*} 1\leqslant |A|+|B|=|A|+\overline{B}+\overline{\overline{B}}\leqslant \biggl( \frac{1}{2\underline{h}}+o(1)\biggr)\|f''\|_{L_1(\mathbb R)}. \end{equation*} \notag
Переходя к пределу при N\to \infty в этом неравенстве, получаем неравенство (4.2), и, значит,
\begin{equation} A_{1,2}(\Delta)\geqslant 2\underline{h}. \end{equation} \tag{4.4}

Для доказательства оценки сверху величины A_{1,2}(\Delta) для любого числа \varepsilon>0 и любой последовательности y\in Y_{1,2} построим функцию f(t)=f_{\delta}(t)\in F_{1,2}(y) такую, что

\begin{equation} \|f''\|_{L_1(\mathbb R)}\leqslant (2+\varepsilon)\overline{h}. \end{equation} \tag{4.5}
Положим
\begin{equation} f''(t)=\begin{cases} \dfrac{Z_k}{2\delta}, &t\in [x_k-\delta;x_k+\delta], \\ 0, &t\notin [x_k-\delta;x_k+\delta], \end{cases} \qquad k\in \mathbb Z, \end{equation} \tag{4.6}
где \delta – произвольное положительное число, \delta< (1/2)\underline{h}. Числа \{Z_k\}_{k\in \mathbb Z} подлежат дальнейшему определению. Подставляя (4.6) в равенство (2.1), получаем разностное уравнение относительно чисел \{Z_k\}_{k\in \mathbb Z}:
\begin{equation} \overline{A}_k Z_{k+2}+\overline{B}_k Z_{k+1}+\overline{C}_kZ_k=[y_{k+2},y_{k+1},y_k],\qquad k\in \mathbb Z, \end{equation} \tag{4.7}
где
\begin{equation} \begin{gathered} \, \overline{A}_k=\frac{\delta}{4h_{k+1}(h_k+h_{k+1})},\qquad \overline{B}_k=\frac{1}{4(h_{k}+h_{k+1})}\biggl(4-\frac{\delta}{h_{k+1}}-\frac{\delta}{h_k} \biggr), \\ \overline{C}_k=\frac{\delta}{4h_{k}(h_k+h_{k+1})}. \end{gathered} \end{equation} \tag{4.8}
Обоснуем, что для любой последовательности y\in Y_{1,2} уравнение (4.7) имеет решение и затем оценим норму этого решения в пространстве l_1. Уравнение (4.7) с учетом (4.8) можно переписать в виде
\begin{equation} Z_{k+1} =(h_k+h_{k+1})\bigl( [y_{k+2},y_{k+1},y_k] +\overline{A}_k(Z_{k+1}-Z_{k+2})+\overline{C}_k(Z_{k+1}-Z_k)\bigr), \qquad k\in \mathbb Z. \end{equation} \tag{4.9}
Рассмотрим нелинейный оператор T, который ставит в соответствие каждой последовательности Z=\{Z_{k+1}\}_{k\in \mathbb Z} последовательность, стоящую в правой части равенства (4.9). Поскольку \delta<(1/2)\underline{h}, то из неравенства
\begin{equation*} \begin{aligned} \, &\| TZ^{(1)}-TZ^{(2)}\|_{l_1}=\bigl\| \bigl\{ \bigl( \overline{A}_k+\overline{C}_k\bigr)(h_k+h_{k+1})\bigl(Z_{k+1}^{(1)}-Z_{k+1}^{(2)}\bigr) \\ &\ \quad +\overline{A}_k(h_k+h_{k+1})\bigl(Z_{k+2}^{(2)}-Z_{k+2}^{(1)}\bigr)+\overline{C}_k(h_k+h_{k+1}) \bigl(Z_{k}^{(2)}-Z_{k}^{(1)}\bigr)\bigr\}\bigr\|_{l_1} \\ &\ \leqslant \biggl\| \biggl\{ \frac{\delta(h_k+h_{k+1})}{4h_kh_{k+1}} \bigl(Z_{k+1}^{(1)}-Z_{k+1}^{(2)}\bigr)\biggr\}\biggr\|_{l_1}+ \biggl\| \biggl\{ \frac{\delta}{4h_{k+1}} \bigl(Z_{k+2}^{(2)}-Z_{k+2}^{(1)}\bigr)\biggr\}\biggr\|_{l_1} \\ &\ \quad +\biggl\| \biggl\{ \frac{\delta}{4h_{k}}\bigl(Z_{k}^{(2)}-Z_{k}^{(1)}\bigr)\biggr\}\biggr\|_{l_1} \leqslant \biggl( \frac14+\frac18+\frac18\biggr)\bigl\| Z^{(1)}\,{-}\,Z^{(2)}\bigr\|_{l_1}\,{=}\,\frac12\| Z^{(1)}\,{-}\,Z^{(2)}\|_{l_1} \end{aligned} \end{equation*} \notag
следует, что T является сжимающим оператором в полном метрическом пространстве l_1=l_1(\mathbb Z). Следовательно, по теореме о сжимающем операторе уравнение (4.9) (а потому и уравнение (4.7)) имеет решение Z\in l_1, и это решение единственно. Оценим норму этого решения в пространстве l_1.

Перепишем уравнение (4.7) в виде

\begin{equation} \begin{aligned} \, &\frac{\delta}{h_{k+1}}Z_{k+2}+\biggl(4-\frac{\delta}{h_{k+1}}-\frac{\delta}{h_k}\biggr)Z_{k+1}+ \frac{\delta}{h_k}Z_k \\ &\qquad=4(h_k+h_{k+1})[y_{k+2},y_{k+1},y_k],\qquad k\in \mathbb Z, \end{aligned} \end{equation} \tag{4.10}
и докажем, что для любого числа \varepsilon >0 имеет место неравенство
\begin{equation} \|Z\|_{l_1}\leqslant (2+\varepsilon)\overline{h}. \end{equation} \tag{4.11}
Будем рассуждать от противного. Пусть \|Z\|_{l_1}>(2+\varepsilon)\overline{h} для некоторого числа \varepsilon>0. Поскольку y\in Y_{1,2}, то
\begin{equation} \| \{4(h_k+h_{k+1})[y_{k+2},y_{k+1},y_k]\} \|_{l_1}\leqslant 8\overline{h}, \end{equation} \tag{4.12}
т. е. норма в пространстве l_1 последовательности, стоящей в правой части разностного уравнения (4.10), не превосходит числа 8\overline{h}. Покажем, что та же норма последовательности, стоящей в левой части этого уравнения, будет больше 8\overline{h} при подходящем выборе числа \delta (это число фигурирует в равенстве (4.6)). Имеем
\begin{equation*} \begin{aligned} \, &\biggl\| \biggl\{ \frac{\delta}{h_{k+1}}Z_{k+2}+\biggl(4-\frac{\delta}{h_{k+1}}-\frac{\delta}{h_{k}}\biggr)Z_{k+1}+ \frac{\delta}{h_{k}}Z_k\biggr\}\biggr\|_{l_1} \\ &\qquad\geqslant \biggl|\, \biggl\| \biggl\{ \biggl(4-\frac{\delta}{h_{k+1}}-\frac{\delta}{h_{k}}\biggr)Z_{k+1} \biggr\}\biggr\|_{l_1}- \biggl\| \biggl\{ \frac{\delta}{h_{k+1}}Z_{k+2} +\frac{\delta}{h_{k}}Z_{k}\biggr\}\biggr\|_{l_1} \biggr|. \end{aligned} \end{equation*} \notag
Оценим снизу выражение (обозначим его через Q), стоящее в правой части этого неравенства под знаком модуля с учетом того, что \delta<(1/2)\underline{h}. Получим
\begin{equation} \begin{aligned} \, Q &\geqslant \biggl\| \biggl\{ \biggl( 4-\frac{\delta}{h_{k+1}}-\frac{\delta}{h_{k}}\biggr)Z_{k+2} \biggr\}\biggr\|_{l_1}- \biggl\| \biggl\{ \frac{\delta}{h_{k+1}}Z_{k+2}\biggr\}\biggr\|_{l_1}- \biggl\| \biggl\{ \frac{\delta}{h_{k}}Z_{k}\biggr\}\biggr\|_{l_1} \\ &\geqslant \biggl( 4-\frac{2\delta}{\underline{h}} -\frac{2\delta}{\underline{h}}\biggr)\|Z\|_{l_1}= \frac{4(\underline{h}-\delta)}{\underline{h}}\|Z\|_{l_1}> \frac{4(2+\varepsilon)(\underline{h}-\delta)\overline{h}}{\underline{h}}. \end{aligned} \end{equation} \tag{4.13}
Выберем число \delta\,{<}\min\{ (1/2)\underline{h},\varepsilon\underline{h}/(2+\varepsilon)\}. Тогда 4(2+\varepsilon)(\underline{h}-\delta)\,{>}\,8\underline{h}, и поэтому из (4.13) следует, что Q>8\overline{h}. Получили противоречие с неравенством (4.12).

Оценим теперь норму функции f'' в пространстве L_1=L_1(\mathbb R). Из (4.6) имеем

\begin{equation*} \|f''\|_{L_1(\mathbb R)}=\sum_{k\in \mathbb Z}\int_{x_k-\delta}^{x_k+\delta}\biggl|\frac{Z_k}{2\delta} \biggr|\,dt=\sum_{k\in \mathbb Z}|Z_k|=\|Z\|_{l_1}\leqslant (2+\varepsilon)\overline{h} \end{equation*} \notag
для любого числа \varepsilon>0, и неравенство (4.5) доказано. Устремляя число \varepsilon к нулю (тогда \delta также стремится к нулю), отсюда получим, что
\begin{equation} A_{1,2}(\Delta)\leqslant 2\overline{h}. \end{equation} \tag{4.14}
Таким образом, из (4.4) и (4.14) следует, что
\begin{equation*} 2\underline{h}\leqslant A_{1,2}(\Delta)\leqslant 2\overline{h}. \end{equation*} \notag
Теорема, сформулированная во введении, в случае p=1 также доказана.

§ 5. Уточнение оценки сверху величины A_{\infty,2}(\Delta)

В настоящем параграфе в случае p=\infty не будем накладывать на сетку \Delta ограничения вида \underline{h}>0 и \overline{h}<+\infty. Пусть

\begin{equation*} x_{k+1/2}=\frac{x_k+x_{k+1}}{2},\qquad k\in \mathbb Z. \end{equation*} \notag
Для любой последовательности y\in Y_{\infty,2} построим параболический сплайн f\in F_{\infty,2}(y) с узлами \{ x_{k+1/2}\}_{k\in \mathbb Z}, полагая
\begin{equation} f''(t)=\begin{cases} Z_k, & x_k\leqslant t<x_{k+1/2}, \\ Z_{k+1}, & x_{k+1/2}\leqslant t<x_{k+3/2}, \\ Z_{k+2}, & x_{k+3/2}\leqslant t<x_{k+2}, \end{cases} \qquad k\in \mathbb Z, \end{equation} \tag{5.1}
где числа \{Z_k\}_{k\in \mathbb Z} подлежат дальнейшему определению. Из (2.1) и (5.1) имеем
\begin{equation} [y_{k+2},y_{k+1},y_k]=\widetilde{A}_kZ_{k+2}+\widetilde{B}_kZ_{k+1}+\widetilde{C}_kZ_{k}, \qquad k\in \mathbb Z, \end{equation} \tag{5.2}
где
\begin{equation} \widetilde{A}_k=\frac{h_{k+1}}{8(h_k+h_{k+1})},\qquad \widetilde{B}_k=\frac38,\qquad \widetilde{C}_k=\frac{h_{k}}{8(h_k+h_{k+1})}. \end{equation} \tag{5.3}
По схеме, изложенной в предыдущих параграфах, докажем существование решения \{Z_k\}_{k\in \mathbb Z} разностного уравнения (5.2). Это уравнение с учетом (5.3) перепишем в виде
\begin{equation} \begin{aligned} \, Z_{k+1}&=2[y_{k+2},y_{k+1},y_k]+\frac{h_{k+1}}{4(h_k+h_{k+1})}(Z_{k+1}-Z_{k+2}) \\ &\qquad+\frac{h_{k}}{4(h_k+h_{k+1})}(Z_{k+1}-Z_{k}). \end{aligned} \end{equation} \tag{5.4}
Нелинейный оператор T, ставящий в соответствие каждой последовательности Z=\{Z_{k+1}\}_{k\in \mathbb Z} последовательность, стоящую в правой части равенства (5.4), является сжимающим, поскольку из (5.4) легко следует, что
\begin{equation*} \|TZ^{(1)}-TZ^{(2)}\|_{l_{\infty}}<\frac12\|Z^{(1)}-Z^{(2)}\|_{l_{\infty}}. \end{equation*} \notag
Значит, по известной теореме о сжимающем операторе в полном метрическом пространстве l_{\infty}=l_{\infty}(\mathbb Z) получаем, что разностное уравнение (5.2) имеет единственное решение Z=\{Z_{k+1}\}_{k\in \mathbb Z}\in l_{\infty}.

Лемма 5.1. Для решения разностного уравнения (5.2) справедлива оценка

\begin{equation*} \|Z\|_{l_{\infty}}=\sup_k |Z_k|\leqslant 4. \end{equation*} \notag

Доказательство. Доказательство леммы 5.1, приведенное ниже, практически повторяет доказательство аналогичной леммы 3 [12] и представлено здесь в основном для полноты изложения.

Пусть Z=\{Z_{k+1}\}_{k\in \mathbb Z} – единственное ограниченное решение уравнения (5.2). Обозначим \alpha=\sup_k|Z_k| и будем рассуждать от противного. Пусть \alpha>4. Тогда существует положительное число \beta такое, что \alpha>\beta+4 (например, можно взять \beta=(1/2)(\alpha-4)). Поскольку \alpha=\sup_k|Z_k|, то для любого целого числа k имеет место неравенство |Z_k|\leqslant \alpha, и для любого числа \varepsilon>0 найдется такое целое число m, что

\begin{equation*} \alpha-\varepsilon<|Z_{m+1}|\leqslant \alpha. \end{equation*} \notag
Не ограничивая общности, будем далее считать, что Z_{m+1}\,{>}\,0 и поэтому \alpha-\varepsilon<Z_{m+1}\leqslant \alpha. Оценим величину
\begin{equation*} D_m=\bigl| \widetilde{A}_mZ_{m+2}+\widetilde{B}_mZ_{m+1}+\widetilde{C}_mZ_{m} \bigr| \end{equation*} \notag
снизу и сверху и придем тем самым к противоречию. Поскольку y\in Y_{\infty,2}, то
\begin{equation} D_m=|[y_{m+2},y_{m+1},y_m]|\leqslant 1. \end{equation} \tag{5.5}
С другой стороны,
\begin{equation} D_m\geqslant\bigl| \widetilde{B}_mZ_{m+1}-\bigl|\widetilde{A}_mZ_{m+2}+\widetilde{C}_mZ_{m}\bigr| \bigr|. \end{equation} \tag{5.6}
Рассмотрим выражение, стоящее под знаком внешнего модуля в правой части неравенства (5.6). Имеем
\begin{equation} \begin{aligned} \, &\widetilde{B}_mZ_{m+1}-\bigl|\widetilde{A}_mZ_{m+2}+\widetilde{C}_mZ_{m} \bigr|> \widetilde{B}_m(\alpha-\varepsilon)-\bigl(\widetilde{A}_m+\widetilde{C}_m\bigr)\alpha \\ &\qquad=\alpha\bigl(\widetilde{B}_m-\widetilde{A}_m-\widetilde{C}_m\bigr) -\widetilde{B}_m\varepsilon> (4+\beta)\bigl(\widetilde{B}_m-\widetilde{A}_m-\widetilde{C}_m \bigr)-\widetilde{B}_m\varepsilon \\ &\qquad=4\bigl( \widetilde{B}_m-\widetilde{A}_m-\widetilde{C}_m\bigr)+\beta\bigl( \widetilde{B}_m-\widetilde{A}_m-\widetilde{C}_m\bigr)-\widetilde{B}_m\varepsilon. \end{aligned} \end{equation} \tag{5.7}
Поскольку \widetilde{B}_m=3/8 и \widetilde{B}_m-\widetilde{A}_m-\widetilde{C}_m=1/4, то из неравенств (5.6) и (5.7) следует, что
\begin{equation*} D_m>1+\frac{\beta}{4}-\frac{3\,\varepsilon}{8}. \end{equation*} \notag
Устремляя положительное число \varepsilon к нулю, получаем, что для некоторого числа m=m(\varepsilon)\in \mathbb Z имеет место неравенство D_m>1, что противоречит (5.5). Лемма 5.1 доказана.

Из леммы 5.1 и равенства (5.1) выводим неравенство

\begin{equation*} \|f''\|_{L_{\infty}(\mathbb R)}=\sup_k|Z_k|\leqslant 4 \end{equation*} \notag
и поэтому A_{\infty,2}(\Delta)\leqslant 4.

Из результатов работы [12] следует, что для любой сетки узлов \Delta имеет место неравенство

\begin{equation*} 2\leqslant A_{\infty,2}(\Delta)\leqslant 18. \end{equation*} \notag
Предыдущие рассуждения позволяют в этом неравенстве уменьшить оценку сверху. Итак, для любой сетки узлов \Delta справедливо неравенство
\begin{equation*} 2\leqslant A_{\infty,2}(\Delta)\leqslant 4. \end{equation*} \notag
При этом можно утверждать, что
\begin{equation*} \inf_{\Delta}A_{\infty,2}(\Delta)=2,\qquad \sup_{\Delta}A_{\infty,2}(\Delta)=\max_{\Delta}A_{\infty,2}=4, \end{equation*} \notag
поскольку в [12] было еще доказано, что для геометрической сетки узлов вида \Delta_r=\{ r^kh\}_{k\in \mathbb Z}, r\geqslant 1, h>0, имеет место равенство
\begin{equation*} A_{\infty,2}(\Delta_r)=\frac{2(r+1)^2}{r^2+1},\qquad r\geqslant 1, \end{equation*} \notag
и при этом A_{\infty,2}(\Delta_1)=A_{\infty,2}(\overline{\Delta})=4, \lim_{r\to \infty}A_{\infty,2}(\Delta_r)=2.

Список литературы

1. А. О. Гельфонд, Исчисление конечных разностей, 3-е изд., Наука, М., 1967, 375 с.  mathscinet  zmath; англ. пер.: A. O. Gelfond, Calculus of finite differences, Int. Monogr. Adv. Math. Phys., Hindustan Publishing Corp., Delhi, 1971, vi+451 с.  mathscinet  zmath
2. С. Б. Стечкин, Ю. Н. Субботин, Сплайны в вычислительной математике, Наука, М., 1976, 248 с.  mathscinet  zmath
3. J. Favard, “Sur l'interpolation”, J. Math. Pures Appl. (9), 19:9 (1940), 281–306  mathscinet  zmath
4. В. С. Рябенький, “Необходимые и достаточные условия хорошей обусловленности краевых задач для систем обыкновенных разностных уравнений”, Ж. вычисл. матем. и матем. физ., 4:2 (1964), 242–255  mathnet  mathscinet  zmath; англ. пер.: V. S. Ryaben'kii, “Necessary and sufficient conditions for good definition of boundary value problems for systems of ordinary difference equations”, U.S.S.R. Comput. Math. Math. Phys., 4:2 (1964), 43–61  crossref
5. В. С. Рябенький, А. Ф. Филиппов, Об устойчивости разностных уравнений, Гостехиздат, М., 1956, 171 с.  mathscinet  zmath; нем. пер.: V. S. Rjabenki, A. F. Filippow, Über die Stabilität von Differenzengleichungen, Mathematik für Naturwissenschaft und Technik, 3, VEB Deutscher Verlag der Wissenschaften, Berlin, 1960, viii+136 pp.  mathscinet  zmath
6. С. Л. Соболев, Лекции по теории кубатурных формул, Ч. 2, Изд-во Новосибирского ун-та, Новосибирск, 1965, 293 с.  mathscinet
7. Ю. Н. Субботин, “О связи между конечными разностями и соответствующими производными”, Экстремальные свойства полиномов, Сборник работ, Тр. МИАН СССР, 78, Наука, М., 1965, 24–42  mathnet  mathscinet  zmath
8. Ю. Н. Субботин, “Функциональная интерполяция в среднем с наименьшей n-й производной”, Приближение функций в среднем, Сборник работ, Тр. МИАН СССР, 88, 1967, 30–60  mathnet  mathscinet  zmath; англ. пер.: Yu. N. Subbotin, “Functional interpolation in the mean with smallest n derivative”, Proc. Steklov Inst. Math., 88 (1967), 31–63
9. Ю. Н. Субботин, “Экстремальные задачи функциональной интерполяции и интерполяционные в среднем сплайны”, Тр. МИАН СССР, 138 (1975), 118–173  mathnet  mathscinet  zmath; англ. пер.: Yu. N. Subbotin, “Extremal problems of functional interpolation, and mean interpolation splines”, Proc. Steklov Inst. Math., 138 (1977), 127–185
10. Ю. Н. Субботин, С. И. Новиков, В. Т. Шевалдин, “Экстремальная функциональная интерполяция и сплайны”, Тр. ИММ УрО РАН, 24, № 3, 2018, 200–225  mathnet  crossref  mathscinet
11. Th. Kunkle, “Favard's interpolation problem in one or more variables”, Constr. Approx., 18:4 (2002), 467–478  crossref  mathscinet  zmath
12. С. И. Новиков, В. Т. Шевалдин, “О связи между второй разделенной разностью и второй производной”, Тр. ИММ УрО РАН, 26, № 2, 2020, 216–224  mathnet  crossref  mathscinet

Образец цитирования: В. Т. Шевалдин, “Экстремальная интерполяция с наименьшим значением нормы второй производной в пространстве L_p(\mathbb R)”, Изв. РАН. Сер. матем., 86:1 (2022), 219–236; Izv. Math., 86:1 (2022), 203–219
Цитирование в формате AMSBIB
\RBibitem{She22}
\by В.~Т.~Шевалдин
\paper Экстремальная интерполяция с~наименьшим значением нормы второй производной в~пространстве $L_p(\mathbb R)$
\jour Изв. РАН. Сер. матем.
\yr 2022
\vol 86
\issue 1
\pages 219--236
\mathnet{http://mi.mathnet.ru/im9125}
\crossref{https://doi.org/10.4213/im9125}
\mathscinet{http://mathscinet.ams.org/mathscinet-getitem?mr=4461231}
\zmath{https://zbmath.org/?q=an:1489.41001}
\adsnasa{https://adsabs.harvard.edu/cgi-bin/bib_query?2022IzMat..86..203S}
\transl
\jour Izv. Math.
\yr 2022
\vol 86
\issue 1
\pages 203--219
\crossref{https://doi.org/10.1070/IM9125}
\isi{https://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcApp=Publons&SrcAuth=Publons_CEL&DestLinkType=FullRecord&DestApp=WOS_CPL&KeyUT=000772176900001}
\scopus{https://www.scopus.com/record/display.url?origin=inward&eid=2-s2.0-85128167519}
Образцы ссылок на эту страницу:
  • https://www.mathnet.ru/rus/im9125
  • https://doi.org/10.4213/im9125
  • https://www.mathnet.ru/rus/im/v86/i1/p219
  • Эта публикация цитируется в следующих 2 статьяx:
    1. C. Э. Нохрин, В. Т. Шевалдин, “О достаточных условиях существования решения бесконечно-разностного уравнения с переменными коэффициентами”, Чебышевский сб., 25:2 (2024), 243–250  mathnet  crossref
    2. Sergey I. Novikov, “Interpolation with minimum value of L_{2}-norm of differential operator”, Ural Math. J., 10:2 (2024), 107–120  mathnet  crossref
    Citing articles in Google Scholar: Russian citations, English citations
    Related articles in Google Scholar: Russian articles, English articles
    Известия Российской академии наук. Серия математическая Izvestiya: Mathematics
    Статистика просмотров:
    Страница аннотации:413
    PDF русской версии:57
    PDF английской версии:32
    HTML русской версии:160
    Список литературы:72
    Первая страница:15
     
      Обратная связь:
    math-net2025_03@mi-ras.ru
     Пользовательское соглашение  Регистрация посетителей портала  Логотипы © Математический институт им. В. А. Стеклова РАН, 2025