Аннотация:
Построены семейства оптимальных методов восстановления аналитических в полосе функций и их производных по неточно заданному следу преобразования Фурье этих функций на вещественной оси. При этом от методов дополнительно требуется, чтобы они были точны на подпространствах целых функций.
Библиография: 12 названий.
Одна из распространенных идей при построении численных методов состоит в том, что ищутся методы, точные на некотором подпространстве функций. При этом исходят из естественных соображений, основанных на том, что если исходная функция приближается с достаточной точностью элементами этого подпространства, то и соответствующий метод (как правило, являющийся некоторым линейным функционалом или оператором от этой функции) будет иметь приемлемую погрешность. Характерным примером здесь являются квадратурные формулы, построенные из условия их точности на алгебраических многочленах фиксированной степени, а наиболее ярким примером являются квадратурные формулы Гаусса (см., например, [1]).
Другой подход к построению численных методов, а в более широком смысле – к аппроксимации в целом, связан с идеями А. Н. Колмогорова. В этом случае вначале фиксируется некоторая априорная информация о функциях – некоторое множество функций (класс), для которых затем строится оптимальный (наилучший) метод из условия его минимальной погрешности на этом классе функций. Здесь также одним из характерных примеров являются квадратурные формулы, впервые построенные в такой постановке С. М. Никольским (см. [2]).
В работе [3] было предложено совместить эти два подхода: один, идущий от Гаусса и основанный на построении методов, точных на подпространствах, и другой, идущий от А. Н. Колмогорова, который основан на построении методов, оптимальных на данном классе. Иначе говоря, предлагалось искать методы, которые были бы оптимальны на классе и в то же время точны на некотором фиксированном подпространстве. В рамках такого подхода в работах [4] и [5] были решены некоторые задачи о восстановлении решений уравнений математической физики.
В настоящей работе рассматриваются задачи построения оптимальных методов восстановления аналитических в полосе функций и их производных по неточно заданному следу преобразования Фурье этих функций на вещественной оси. Причем от оптимальных методов дополнительно требуется, чтобы они были точны на подпространствах целых функций.
§ 2. Постановка задачи
Пусть $X$ – линейное пространство, $Y,Z$ – линейные нормированные пространства и $A\colon X\to Z$, $I\colon X\to Y$ – линейные операторы. Рассматривается задача об оптимальном восстановлении значений оператора $A$ на множестве $W\subset X$ по неточно заданным значениям оператора $I$ на элементах этого множества. Считается, что для каждого элемента $x\in W$ известно значение $y\in Y$ такое, что $\|Ix-y\|_Y\leqslant\delta$, где $\delta>0$ – число, характеризующее погрешность исходной информации об элементах множества $W$. Требуется по значению $y$ восстановить наиболее точным образом значение $Ax$. Любой метод восстановления представляет из себя отображение $m\colon Y\to Z$, которое элементу $y\in Y$ ставит в соответствие значение $m(y)\in Z$, принимаемое за приближенное значение $Ax$.
а методы, на которых достигается нижняя грань, называются оптимальными на множестве $W$. Сформулированная задача относится к теории оптимального восстановления. Более подробные сведения об этой теории и задачах, рассматриваемых в рамках этой теории, можно найти в обзорной статье [6] и монографиях [7]–[10].
Пусть $L\subset X$ – линейное подпространство $X$. Будем говорить, что метод $m\colon Y\to Z$ точен на $L$, если $Ax=m(Ix)$ для всех $x\in L$. Рассмотрим множество $\mathcal E_L$, состоящее из линейных операторов $m\colon Y\to Z$, точных на $L$. Положим
Предложение (см. [4]). Пусть $L\subset X$ – линейное подпространство $X$ и $m^*\colon Y\to Z$ – линейный оператор, являющийся оптимальным методом восстановления оператора $A$ на множестве $W+L$. Тогда
функций и их производных при условии, что методы восстановления точны на подпространстве целых функций $\mathcal B_{\sigma,2}(\mathbb R)$, являющемся подпространством в $L_2(\mathbb R)$, образованном сужениями на $\mathbb R$ целых функций экспоненциального типа $\sigma$.
Перейдем к точной постановке задачи. Пространством Харди $\mathcal H_2^\beta$ называется множество функций $f$, аналитических в полосе $S_\beta$, для которых
Через $\mathcal H_2^{r,\beta}$ (пространство Харди–Соболева) будем обозначать множество аналитических в полосе $S_\beta$ функций, для которых $f^{(r)}\in\mathcal H_2^\beta$.
Обозначим через $H_2^{r,\beta}$ множество функций $f\in\mathcal H_2^{r,\beta}\cap L_2(\mathbb R)$, для которых $\|f^{(r)}\|_{\mathcal H_2^\beta}\leqslant1$. Если $\sigma>0$, то через $\mathcal B_{\sigma,2}(\mathbb R)$ обозначим подпространство в $L_2(\mathbb R)$, образованное сужениями на $\mathbb R$ целых функций экспоненциального типа $\sigma$. Как хорошо известно, $f\in\mathcal B_{\sigma,2}(\mathbb R)$ тогда и только тогда, когда носитель преобразования Фурье $Ff$ принадлежит отрезку $\Delta_\sigma=[-\sigma,\sigma]$. По определению $\mathcal B_{0,2}(\mathbb R)=\{0\}$.
Рассмотрим задачу оптимального восстановления $k$-й производной функции $f\in H_2^{r,\beta}+\mathcal B_{\sigma,2}(\mathbb R)$, $k\leqslant r$, по следу на $\Delta_{\sigma_1}$, $\sigma_1>0$, ее преобразования Фурье, заданному с погрешностью в метрике $L_2(\Delta_{\sigma_1})$, т.е. считается, что вместо следа на $\Delta_{\sigma_1}$ функции $Ff$ известна функция $y\in L_2(\Delta_{\sigma_1})$ такая, что
Иными словами, мы хотим найти оптимальные методы восстановления $k$-й производной на классе $H_2^{r,\beta}$ среди точных на подпространстве целых функций $\mathcal B_{\sigma,2}(\mathbb R)$. Без требования точности методов на $\mathcal B_{\sigma,2}(\mathbb R)$ эта задача рассматривалась в работе [11].
§ 3. Основные результаты
Рассмотрим функцию $y=s(x)$, $x\geqslant0$, которая задается параметрически
В силу вогнутости функции $s$ найдется точка $t_0$ такая, что касательная к $s$ в точке $(x(t_0),y(t_0))$ будет параллельна прямой $y=\lambda_2x$. Тем самым точка $t_0$ находится из уравнения
В силу того, что функция в правой части (3.2) при $t_0\in[0,+\infty)$ монотонно возрастает от нуля до $+\infty$, для любого числа $\lambda_1>0$ найдется точка $t_0>0$ такая, что касательная к $s$ в точке $(x(t_0),y(t_0))$ будет проходить через точку $(0,\lambda_1)$. Обозначим такую точку $t_0$ через $h_1(\lambda_1)$.
Функция $t^r\sqrt{\operatorname{ch}2\beta t}$ монотонно возрастает при $t\in\mathbb R_+$ от $0$ до $+\infty$. Поэтому для любого $x\in\mathbb R_+$ существует единственное решение уравнения
принадлежащее интервалу $[0,+\infty)$. Обозначим его через $\mu_{r\beta}(x)$.
Через $\widehat\sigma_1$ обозначим значение параметра $t$, при котором $t_0=\widehat t_0=\mu_{r\beta}(\sqrt{2\pi}/\delta)$, т.е. $x(\widehat t_0)=2\pi/\delta^2$. Положим $\widehat\sigma=h(\widehat\sigma_1)$. Уравнение касательной, проходящей через точку $(x(\widehat t_0),y(\widehat t_0))$, будет иметь вид $y=\widehat\lambda_1+\widehat\lambda_2x$, где
Через $\Theta(\sigma,\sigma_1)$ обозначим множество измеримых функций $\theta$ на $[-\sigma_1,-\sigma)\cup(\sigma,-\sigma_1]$, для которых $|\theta(t)|\leqslant1$ для п.в. $\sigma<|t|\leqslant\sigma_1$.
Доказательство. По основной теореме о представлении аналитических функций над трубчатыми областями (см. [12]) следует, что $f\in\mathcal H_2^\beta$ в том и только том случае, если она имеет вид
Докажем, что $f\in\mathcal H_2^{r,\beta}\cap L_2(\mathbb R)$ принадлежит классу $H_2^{r,\beta}+\mathcal B_{\sigma,2}(\mathbb R)$ тогда и только тогда, когда
Действительно, если $f\in H_2^{r,\beta}+\mathcal B_{\sigma,2}(\mathbb R)$, то $f=f_1+f_2$, где $f_1\in H_2^{r,\beta}$, а $f_2\in\mathcal B_{\sigma,2}(\mathbb R)$. Тогда, учитывая, что $Ff_2$ сосредоточено на отрезке $\Delta_{\sigma}$, имеем
Обратно, пусть $f\in\mathcal H_2^{r,\beta}\cap L_2(\mathbb R)$ такова, что выполнено неравенство (3.9). Обозначим через $f_2\in L_2(\mathbb R)$ функцию, для которой $Ff_2=\chi_{\sigma}Ff$, где $\chi_{\sigma}$ – характеристическая функция отрезка $\Delta_\sigma$. Тогда ясно, что $f_2\in\mathcal B_{\sigma,2}(\mathbb R)$. Положим $f_1=f-f_2$. Очевидно, что $f_1\in\mathcal H_2^{r,\beta}\cap L_2(\mathbb R)$, и в силу (3.8) (учитывая, что $Ff_1=0$ на $\Delta_{\sigma}$) будем иметь
т.е. $f=f_1+f_2\in H_2^{r,\beta}+\mathcal B_{\sigma,2}(\mathbb R)$.
Пусть $f\in\mathcal H_2^{r,\beta}\cap L_2(\mathbb R)$ такова, что $\|Ff\|_{L_2(\Delta_{\sigma_1})}\leqslant\delta$, и выполнено неравенство (3.9). Для любого метода $m\colon L_2(\Delta_{\sigma_1})\to L_2(\mathbb R)$ имеем
Как было отмечено, в силу вогнутости функции $s$ найдется точка $t_0$ такая, что касательная к $s$ в точке $(x(t_0),y(t_0))$ будет параллельна прямой $y=\lambda_2x$. Сама касательная, проходящая через точку $(x(t_0),y(t_0))$, будет иметь вид $y= \lambda_1+ \lambda_2x$, где
Из того, что $\sigma_1\leqslant\widehat\sigma_1$, вытекает, что $t_0\leqslant\widehat t_0$. Следовательно, $t_0^{2r}\operatorname{ch}2\beta t_0\leqslant2\pi/\delta^2$. Для каждого $n\in\mathbb N$ такого, что $h(\sigma_1)<t_0-1/n$, рассмотрим функцию $f_n$, для которой
Пусть $(\sigma_1,\sigma)\in\Sigma_4$, а $t_0$ – точка, которая определялась при рассмотрении случая, когда $(\sigma_1,\sigma)\in\Sigma_2$. Положим $\xi=h_1(\sigma^{2k})$. В силу того, что $\sigma\leqslant h(\sigma_1)$, получаем, что $\xi\leqslant t_0<\sigma_1$. Кроме того, так как $\sigma\geqslant\widehat\sigma$, то $\xi\geqslant\widehat t_0$. Следовательно, $\xi^{2r}\operatorname{ch}2\beta\xi\geqslant2\pi/\delta^2$. Для каждого $n\in\mathbb N$ такого, что $1/n<\sigma$ и $\xi-1/n>\sigma$, рассмотрим функцию $f_n$, для которой
Будем искать оптимальные методы восстановления $m_a\colon L_2(\Delta_{\sigma_1})\to L_2(\mathbb R)$ среди отображений, которые в образах Фурье представляются в виде
Отсюда следует, что если функция $a$ не равна п.в. единице на $\Delta_\sigma$, то, поскольку $\mathcal B_{\sigma,2}(\mathbb R)$ – линейное пространство, значение задачи (3.16) (и тем самым задачи (3.15)) равно бесконечности, т.е. погрешность метода с таким $a$ бесконечна.
Пусть $a\equiv1$ на $\Delta_\sigma$. Оценим сверху максимизируемый функционал в (3.16), представив его для этого в виде суммы трех слагаемых:
Пусть $(\sigma_1,\sigma)\in\Sigma_4$. Тогда $\lambda_2$ – угловой коэффициент касательной к $s$ в точке $(x(\xi),y(\xi))$, а величина $\sigma_1^{2(k-r)}\operatorname{ch}^{-1}2\beta\sigma_1$ – угловой коэффициент касательной к $s$ в точке $(x(t_0),y(t_0))$ (точка $t_0$ определялась при рассмотрении оценки снизу для случая $(\sigma_1,\sigma)\in\Sigma_2$). Поскольку $\xi\leqslant t_0$, а функция $s$ вогнутая, то
Таким образом, оценка (3.20) справедлива во всех областях.
Если предположить, что функция $a$ такова, что $S_a\leqslant1$, то, складывая неравенства (3.17), (3.19) и (3.20), получаем следующую оценку для максимизируемого функционала в задаче (3.16):
Покажем, что функции $a$, для которых $S_a\leqslant1$, существуют. Заметим (“выделяя полный квадрат”), что условие $S_a\leqslant1$ равносильно тому, что для п.в. $\sigma<|t|\leqslant\sigma_1$ выполнено неравенство
при $\sigma<|t|\leqslant\sigma_1$, то такие $a$, очевидно, существуют и описываются равенством (3.5).
Если $(\sigma_1,\sigma)\in\Sigma_1$, то прямая $y=\lambda_1+\lambda_2x$ параллельна касательной, проведенной к $s$ в точке $(x(t_0),y(t_0))$, где $t_0$ определяется равенством
(см. (3.1)), и в силу условия $\sigma\geqslant h(\sigma_1)$ проходит не ниже касательной. Следовательно, в силу вогнутости $s$ для всех $x\geqslant0$ выполняется неравенство $\lambda_1+\lambda_2x\geqslant s(x)$. Отсюда вытекает условие (3.21). В остальных трех случаях прямые $y=\lambda_1+\lambda_2x$ являются касательными к $s$, и по тем же соображениям условие (3.21) выполнено.
3) Пусть $k=0$, $\sigma_1>0$, $\sigma\geqslant0$ и $\sigma\leqslant\sigma_1$. Как было показано выше, функции $f_n$, определенные равенством (3.13), являются допустимыми в задаче (3.11). Следовательно,
Будем искать оптимальные методы восстановления $m_a\colon L_2(\Delta_{\sigma_1})\to L_2(\mathbb R)$ среди отображений, которые в образах Фурье представляются в виде (3.14) с $k=0$. Рассуждая по той же схеме, которая была приведена выше, считаем, что $a\equiv1$ на $\Delta_\sigma$, и будем оценивать сверху максимизируемый функционал в (3.16) (при $k=0$), представив его для этого в виде суммы трех слагаемых:
Если предположить, что функция $a$ такова, что $\widetilde S_a\leqslant1$, то, учитывая оценки, полученные для $I_2$ и $I_3$, получаем следующую оценку для максимизируемого функционала в задаче (3.16) (при $k=0$):
Такие $a$, очевидно, существуют и описываются равенством (3.6).
Теорема доказана.
§ 4. Обсуждение оптимальных методов
При восстановлении $f^{(k)}$ на классе $H_2^{r,\beta}+\mathcal B_{\sigma,2}(\mathbb R)$ по неточно заданному преобразованию Фурье функции $f$ на отрезке $[-\sigma_1,\sigma_1]$ возникают следующие два вопроса:
Иными словами, нас интересует, не является ли часть информации, которая нам дается о функции $f$, лишней и нельзя ли среди получаемого семейства оптимальных методов найти тот, который был бы точен на более широком подпространстве и не увеличивал бы погрешность оптимального восстановления?
Будем рассматривать случай, когда $k\geqslant1$. Ответы на эти вопросы зависят от того, в какую из областей $\Sigma_j$, $j= 1,2,3,4$, попадает точка $(\sigma_1,\sigma)$.
В случае, когда $(\sigma_1,\sigma)\in\Sigma_1$, из (3.4) видно, что при уменьшении $\sigma_1$ или увеличении $\sigma$ погрешность оптимального восстановления увеличивается. Тем самым ответы на эти вопросы в этом случае отрицательные.
Если $(\sigma_1,\sigma)\in\Sigma_2$, то погрешность оптимального восстановления не изменится для точки $(\sigma_1,h(\sigma_1)$. Это означает, что можно расширить исходное подпространство $\mathcal B_{\sigma,2}(\mathbb R)$ до подпространства $\mathcal B_{h(\sigma_1),2}(\mathbb R)$, не увеличивая при этом погрешность оптимального восстановления.
При $(\sigma_1,\sigma)\in\Sigma_4$ можно уменьшить отрезок, на котором задается информация о функции $f$, сузив его до отрезка $[-\sigma_1',\sigma_1']$, где $\sigma_1'$ таково, что $h(\sigma_1')=\sigma$.
Наконец, если $(\sigma_1,\sigma)\in\Sigma_3$, то можно и сузить отрезок, на котором задается информация о функции $f$, до отрезка $[-\widehat\sigma_1,\widehat\sigma_1]$, и расширить подпространство до $\mathcal B_{\widehat\sigma,2}(\mathbb R)$.
С. М. Никольский, Квадратурные формулы, 4-е изд., Наука, М., 1988, 256 с. ; исп. пер. 3-го изд.: S. Nikolski, Fórmulas de cuadratura, Editorial Mir, Moscow, 1990, 293 pp.
2.
С. М. Никольский, “К вопросу об оценках приближений квадратурными формулами”, УМН, 5:2(36) (1950), 165–177
3.
Г. Г. Магарил-Ильяев, К. Ю. Осипенко, “Точность и оптимальность методов восстановления функций по их спектру”, Функциональные пространства, теория приближений, смежные разделы математического анализа, Сборник статей. К 110-летию со дня рождения академика Сергея Михайловича Никольского, Труды МИАН, 293, МАИК «Наука/Интерпериодика», М., 2016, 201–216; англ. пер.: G. G. Magaril-Il'yaev, K. Yu. Osipenko, “Exactness and optimality of methods for recovering functions from their spectrum”, Proc. Steklov Inst. Math., 293 (2016), 194–208
4.
Е. А. Балова, К. Ю. Осипенко, “Оптимальные методы восстановления решений задачи Дирихле, точные на подпространствах сферических гармоник”, Матем. заметки, 104:6 (2018), 803–811; англ. пер.: E. A. Balova, K. Yu. Osipenko, “Optimal recovery methods for solutions of the Dirichlet problem that are exact on subspaces of spherical harmonics”, Math. Notes, 104:6 (2018), 781–788
5.
С. А. Унучек, “Оптимальные методы восстановления решения уравнения теплопроводности, точные на тригонометрических полиномах”, Матем. заметки, 113:1 (2023), 118–131; англ. пер.: S. A. Unuchek, “Optimal recovery methods exact on trigonometric polynomials for the solution of the heat equation”, Math. Notes, 113:1 (2023), 116–128
6.
C. A. Micchelli, T. J. Rivlin, “A survey of optimal recovery”, Optimal estimation in approximation theory (Freudenstadt, 1976), The IBM Research Symposia Series, Plenum, New York, 1977, 1–54
7.
Дж. Трауб, X. Вожьняковский, Общая теория оптимальных алгоритмов, Мир, М., 1983, 384 с. ; пер. с англ.: J. F. Traub, H. Woźniakowski, A general theory of optimal algorithms, ACM Monogr. Ser., Academic Press, Inc., New York–London, 1980, xiv+341 с.
8.
L. Plaskota, Noisy information and computational complexity, Cambridge Univ. Press, Cambridge, 1996, xii+308 pp.
9.
K. Yu. Osipenko, Optimal recovery of analytic functions, Nova Science Publ., Inc., Huntington, NY, 2000, 220 pp.
10.
К. Ю. Осипенко, Введение в теорию оптимального восстановления, Лань, СПб., 2022, 388 с.
11.
К. Ю. Осипенко, “Неравенство Харди–Литтлвуда–Полиа для аналитических функций из пространств Харди–Соболева”, Матем. сб., 197:3 (2006), 15–34; англ. пер.: K. Yu. Osipenko, “The Hardy–Littlewood–Pólya inequality for analytic functions in Hardy–Sobolev spaces”, Sb. Math., 197:3 (2006), 315–334
12.
И. Стейн, Г. Вейс, Введение в гармонический анализ на евклидовых пространствах, Мир, М., 1974, 336 с. ; пер. с англ.: E. M. Stein, G. Weiss, Introduction to Fourier analysis on Euclidean spaces, Princeton Math. Ser., 32, Princeton Univ. Press, Princeton, NJ, 1971, x+297 с.
Образец цитирования:
К. Ю. Осипенко, “О восстановлении аналитических функций, точном на подпространствах целых функций”, Матем. сб., 215:3 (2024), 100–118; K. Yu. Osipenko, “Recovery of analytic functions that is exact on subspaces of entire functions”, Sb. Math., 215:3 (2024), 383–400