Processing math: 5%

数学物理学报, 2024, 44(3): 525-538

加权 Korobov 空间中多元 L-逼近问题的指数收敛易处理性

张杰,1,*, 孙艺铭,1, 刘永平,2

1.山东建筑大学理学院 济南 250101

2.北京师范大学数学科学学院 北京 100875

EC-tractability of Multivariate L-approximation in Weighted Korobov Spaces

Zhang Jie,1,*, Sun Yiming,1, Liu Yongping,2

1. College of Science, Shandong Jianzhu University, Jinan 250101

2. Department of Mathematics, Beijing Normal University, Beijing 100875

通讯作者: Email:zhangjie1991@sdjzu.edu.com

收稿日期: 2023-07-31   修回日期: 2023-10-17  

基金资助: 国家自然科学基金(12101369)
国家自然科学基金(11871006)
山东省高等学校青创团队计划(2022KJ209)
山东建筑大学博士基金项目(X19090Z)

Received: 2023-07-31   Revised: 2023-10-17  

Fund supported: NSFC(12101369)
NSFC(11871006)
Development Plan of Youth Innovation Team of Shandong Provincial Colleges and Universities(2022KJ209)
Doctoral Foundation Project of Shandong Jianzhu University(X19090Z)

作者简介 About authors

孙艺铭:Email:sunyimingzz123@163.com;

刘永平:Email:ypliu@bnu.edu.cn

摘要

该文主要研究最坏框架下加权 Korobov 空间中多元 L-逼近问题的指数易处理性.多元逼近问题中的算法使用的信息取自由线性泛函组成的线性信息类 Λall和函数值组成的标准信息类 Λstd.该问题的指数收敛-拟多项式易处理性和指数收敛-一致弱易处理性之前并没有被研究,该文最终通过两个权参数序列给出使得多元 L-逼近问题具有这两种指数收敛易处理性的充分必要条件.

关键词: 指数收敛易处理性; Korobov 空间; 最坏框架; 多元 L-逼近问题

Abstract

In this paper we study exponential tractability of multivariate L-approximation for weighted Korobov spaces in the worst case setting. We consider all algorithms that use the class Λall of all linear functionals and the class Λstd of only function evaluations as information. We give matching necessary and sufficient conditions for notions of EC-quasi-polynomial tractability and EC-uniform weak tractability which have not been discussed before in terms of two weight parameters of the problem.

Keywords: Exponential convergence tractability; Korobov spaces; Worst case setting; Multivariate L -approximation

PDF (797KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张杰, 孙艺铭, 刘永平. 加权 Korobov 空间中多元 L-逼近问题的指数收敛易处理性[J]. 数学物理学报, 2024, 44(3): 525-538

Zhang Jie, Sun Yiming, Liu Yongping. EC-tractability of Multivariate L-approximation in Weighted Korobov Spaces[J]. Acta Mathematica Scientia, 2024, 44(3): 525-538

1 引言

多元连续问题易处理性的研究始于 1994 年[20,21],历经近三十年的发展, 越来越多的学者对该问题产生兴趣,现在已经变成非常热门的研究领域[11-13].对于给定的正整数 d, 多元连续问题通常指的是一系列解算子 S={Sd:FdGd}dN,其中 Fd 是由 d 元函数组成的函数类, Gd 是一个 Banach 空间. 事实上,多元逼近和多元积分[3,6,16,25]是其中两个最重要且被广泛研究的多元连续问题. 通常使用包含所有函数值的标准信息类 Λstd 或者包含所有连续线性泛函的线性信息类 Λall去逼近这些算子.

本文使用的算法其质量通过最坏框架下的逼近误差进行测度, 即函数空间单位球上的最大误差.由于维数 d 通常很大(成百上千), 研究算法的误差时不仅要考虑如何依赖于信息算子的数量,而且要考虑怎样依赖于维数 d. 在逼近 Sd 时,误差不超过一个常数的 ε 倍的算法中所需要信息算子的最小数量称为信息复杂度,并记作 nΛ(ε,d), 其中 Λ{Λall,Λstd}.如果该常数等于 1, 此时称为绝对误差标准. 如果该常数等于初始误差(没有使用任何信息的误差), 此时称为归一化误差标准.关于信息复杂度的基础知识可以参见文献[18].

根据多项式收敛和指数收敛, 易处理性被分成两种类型.对于代数易处理性[4,14,15],它刻画的是信息复杂度如何依赖于 dε1. 而对于指数收敛易处理性,它刻画的是信息复杂度怎样依赖于 d(1+lnε1).近些年, 指数收敛易处理性已变成多元问题易处理性中最热门的研究领域之一,参见文献[1-3,5-8,10,17,19,22,24,26].

然而, 大多数关于指数收敛易处理性的文献过去主要研究定义在 Hilbert 空间上的多元连续问题,参见文献[1-3,5,10,17,19,24,26].本文我们讨论定义在加权 Korobov 空间中多元 L-逼近问题的指数收敛易处理性.Kritzer等[7]首次研究了该问题, 得到了 Λ{Λall,Λstd}以及绝对和归一化误差标准下的一些指数收敛易处理性结果,包括指数收敛-(强)多项式易处理性, 指数收敛-弱易处理性和指数收敛-(s,t)-弱易处理性 (s=1, t>1).2020 年, 许贵桥[22]讨论了多元 Lp-逼近问题 (1p),并得到了一些指数收敛易处理性概念的充分必要条件.

本文继续对上述问题进行讨论,研究之前文献中没有解决的多元 L-逼近问题的指数收敛-拟多项式易处理性和指数收敛-一致弱易处理性.更确切地说, 我们得到了 Λ{Λall,Λstd}及绝对误差标准下指数收敛-拟多项式易处理性和 Λstd及绝对和归一化误差标准下指数收敛-一致弱易处理性的充分必要条件.

2 加权 Korobov 空间中多元 L-逼近问题

本文主要讨论 L 范数下具有指数权的 Korobov 空间中多元函数逼近问题,简称为多元 L-逼近.下面首先给出已被大量文献研究过的加权 Korobov 空间的定义,参见文献[2,5,7,22,24].

a={ak}k1b={bk}k1 为两个正实数权序列.通常假定它们满足

0<a1a2akb:=inf
(2.1)

对一个固定正实数 \omega\in(0,1), {\rm i}=\sqrt{-1}, a, b>0, 假设解析的 Korobov 核K_{d,\textbf{a},\textbf{b}} 具有如下的乘积形式

\begin{equation}\label{Kernel} K_{d,\textbf{a},\textbf{b}}(\textbf{x},\textbf{y}) =\prod_{k=1}^{d}K_{1,a_{k},b_{k}}(x_{k},y_{k}),\quad \textbf{x},\textbf{y}\in[0,1]^{d}, \end{equation}
(2.2)

其中 K_{1,a,b} 为单变量解析的 Korobov 核,

\begin{equation*} K_{1,a,b}(x,y)=\sum\limits_{h\in\mathbb{Z}} \omega^{a|h|^{b}}\exp(2\pi {\rm i}h(x-y)),\quad x,y\in[0,1]. \end{equation*}

因此

\begin{equation}\label{Korobov kernel} K_{d,\textbf{a},\textbf{b}}(\textbf{x},\textbf{y})=\sum\limits_{\textbf{h}\in \mathbb{Z} ^{d}} \omega_{\textbf{h}}\exp(2\pi {\rm i}\textbf{h}\cdot(\textbf{x}-\textbf{y})), \quad \textbf{x},\textbf{y}\in[0,1]^{d}, \end{equation}
(2.3)

这里\omega_{\textbf{h}}=\omega^{\sum_{k=1}^{d}a_{k}|h_{k}|^{b_{k}}},\textbf{h}=(h_{1},h_{2},\cdots,h_{d}), \textbf{x}=(x_{1},x_{2},\cdots,x_{d}), \textbf{y}=(y_{1},y_{2},\cdots,y_{d})

\textbf{h}\cdot (\textbf{x}-\textbf{y})=\sum\limits_{k=1}^{d}h_{k}(x_{k}-y_{k}).

经常用 H(K_{d,\textbf{a},\textbf{b}}) 来表示解析的 Korobov 空间,它是一个具有核(2.3)式的再生核 Hilbert 空间[2,5].

现在考虑多元 \mathbb{L}_{\infty}-逼近问题, 它通过如下的嵌入算子定义

\begin{equation*} \text{APP}_{d,\infty}:H(K_{d,\textbf{a},\textbf{b}})\rightarrow L_{\infty}([0,1]^{d}), \quad \text{APP}_{d,\infty}(f)=f, \end{equation*}

即解算子是 H(K_{d,\textbf{a},\textbf{b}})L_{\infty}([0,1]^{d}) 的嵌入映射. 由文献[7]可知这个嵌入是连续的, 而且

\|\text{APP}_{d,\infty}\|=\prod_{j=1}^{d}\left( 1+2\sum\limits_{h=1}^{\infty}\omega^{a_{j}h^{b_{j}}}\right)^{1/2}.

通常利用线性非自适应算法[11,18]来逼近 \text{APP}_{d,\infty},该算法使用的信息算子为线性信息或标准信息. 即函数 f\in H(K_{d,\textbf{a},\textbf{b}}) 通过如下的算法

\begin{equation} A_{n,d}(f)=\sum_{k=1}^{n}\alpha_{k}L_{k}(f) \end{equation}
(2.4)

进行逼近, 其中 L_{k}\in\Lambda^{\text{all}}=H^{*}(K_{d,\textbf{a},\textbf{b}}) (H(K_{d,\textbf{a},\textbf{b}}) 的对偶空间) 或L_{k}\in\Lambda^{\text{std}} =\{L|L(g)=g(\textbf{x}),\ \forall\ g\in H(K_{d,\textbf{a},\textbf{b}}),\textbf{x}\in[0,1]^{d}\}, 1\leq k \leq n, \alpha_{k}L_{\infty}([0,1]^{d}) 中的元素. 因为 H(K_{d,\textbf{a},\textbf{b}}) 是一个再生核 Hilbert 空间, 所以函数值也是连续线性泛函, 因此 \Lambda^{\text{std}}\subseteq \Lambda^{\text{all}}.

由 (2.4) 式的算法 A_{n,d} 进行逼近的最坏框架下的误差定义为

\begin{equation*} e^{L_\infty}(A_{n,d},H(K_{d,\textbf{a},\textbf{b}})):=\mathop{ \sup\limits_{f\in H(K_{d,\textbf{a},\textbf{b}})}}_{\|f\|_{H(K_{d,\textbf{a},\textbf{b}})}\leq1} \| f-A_{n,d}(f) \|_{L_{\infty}([0,1]^{d})}, \end{equation*}

其中 \|\cdot\|_{H(K_{d,\textbf{a},\textbf{b}})} 表示 H(K_{d,\textbf{a},\textbf{b}}) 的范数, \|\cdot\|_{L_{\infty}([0,1]^{d})} 表示本性上确界. n 阶最小最坏框架下的误差定义为

e^{L_\infty,\Lambda}(n,d)=\inf_{A_{n,d}}e^{L_\infty}(A_{n,d},H(K_{d,\textbf{a},\textbf{b}})),

这里的下确界选自所有 (2.4) 式的算法. 当 n=0 时, 用 0 逼近函数 f, 此时的初始误差为

e^{L_\infty,\Lambda}(0,d)=\|\text{APP}_{d,\infty}\|=\prod_{j=1}^{d} \left( 1+2\sum\limits_{h=1}^{\infty}\omega^{a_{j}h^{b_{j}}}\right)^{1/2}.

显然 d 非常大时初始误差可能任意大, 这意味着多元\mathbb{L}_{\infty}-逼近问题可能不是规范化的.

对于 \varepsilon\in(0,1), d\in \mathbb{N} , \Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\}, 用 n^{L_\infty,\Lambda}_{\text{abs}}(\varepsilon,d) 表示绝对误差标准下的信息复杂度, 它被定义成使得误差不大于预先给定的精度 \varepsilon 的算法需要的信息算子的最小数量, 即

n^{L_\infty,\Lambda}_{\text{abs}}(\varepsilon,d):= \min\{n\in \mathbb{N} :e^{L_\infty,\Lambda}(n,d)\leq \varepsilon\}.

类似地, 归一化误差标准下的信息复杂度定义为

n^{L_\infty,\Lambda}_{\text{nor}}(\varepsilon,d):= \min\{n\in \mathbb{N} :e^{L_\infty,\Lambda}(n,d)\leq \varepsilon e^{L_\infty,\Lambda}(0,d)\}.

由于 e^{L_\infty,\Lambda^{\text{all}}}(n,d)\leq e^{L_\infty,\Lambda^{\text{std}}}(n,d), 从而

\begin{equation} n^{L_\infty,\Lambda^{\text{all}}}_{X}(\varepsilon,d) \leq n^{L_\infty,\Lambda^{\text{std}}}_{X}(\varepsilon,d), \end{equation}
(2.5)

其中 X\in\{\text{abs},\text{nor}\}. 又因为 e^{L_\infty,\Lambda}(0,d)>1, 故对任意\Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\}

\begin{equation} n^{L_\infty,\Lambda}_{\text{nor}}(\varepsilon,d) \leq n^{L_\infty,\Lambda}_{\text{abs}}(\varepsilon,d). \end{equation}
(2.6)

下面给出指数收敛易处理性的概念. 对于 \Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\}X\in\{\text{abs},\text{nor} \}, 称多元 \mathbb{L}_{\infty}-逼近问题\text{APP}_{d,\infty}

·指数收敛-强多项式易处理性 (EC-SPT)当且仅当存在非负实数 Cp 使得

\begin{equation} n^{L_\infty,\Lambda}_{X}(\varepsilon,d)\leq C(1+\ln \varepsilon^{-1})^{p} \end{equation}
(2.7)

对所有 d\in\mathbb{N},\varepsilon\in(0,1) 成立.将满足 (2.7) 式中 p 的下确界称为 EC-SPT 的指数, 记作 p^{*}.

·指数收敛-多项式易处理性 (EC-PT)当且仅当存在非负实数 C, pq 使得

\begin{equation*} n^{L_\infty,\Lambda}_{X}(\varepsilon,d)\leq Cd^{q}(1+\ln \varepsilon^{-1})^{p} \end{equation*}

对所有 d\in\mathbb{N},\varepsilon\in(0,1) 成立.

·指数收敛-拟多项式易处理性 (EC-QPT)当且仅当存在非负实数 Ct 使得

\begin{equation} n^{L_\infty,\Lambda}_{X}(\varepsilon,d)\leq C\exp\left(t(1+\ln d)(1+\ln(1+\ln\varepsilon^{-1}))\right) \end{equation}
(2.8)

对所有 d\in\mathbb{N},\varepsilon\in(0,1) 成立.将满足(2.8)式中 t 的下确界称为 EC-QPT 的指数, 记作 t^{*}.

·指数收敛-一致弱易处理性 (EC-UWT) 当且仅当

\begin{equation*} \lim\limits_{d+\varepsilon^{-1}\rightarrow+\infty } \frac{\ln n^{L_\infty,\Lambda}_{X}(\varepsilon,d)}{d^{s}+(\ln\varepsilon^{-1})^{t}}=0 \end{equation*}

对所有 s, t>0 成立.

·指数收敛-弱易处理性 (EC-WT) 当且仅当

\begin{equation*} \lim\limits_{d+\varepsilon^{-1}\rightarrow \infty} \frac{\ln n^{L_\infty,\Lambda}_{X}(\varepsilon,d)}{d+\ln\varepsilon^{-1}}=0. \end{equation*}

·指数收敛-(s,t)-弱易处理性 (EC-(s,t)-WT, s,t>0) 当且仅当

\begin{equation*} \lim\limits_{d+\varepsilon^{-1}\rightarrow \infty} \frac{\ln n^{L_\infty,\Lambda}_{X}(\varepsilon,d)}{d^{s}+(\ln\varepsilon^{-1})^{t}}=0. \end{equation*}

容易验证上述指数收敛易处理性概念具有如下的关系:

\begin{equation*} \text{EC-SPT}\Rightarrow \text{EC-PT} \Rightarrow \text{EC-QPT} \Rightarrow \text{EC-UWT}\Rightarrow \text{EC-WT}. \end{equation*}

如果把上述指数收敛易处理性定义中的 (1+\ln\varepsilon^{-1})\ln\varepsilon^{-1}换成 \varepsilon^{-1}, 相应地得到强多项式易处理性 (SPT)、多项式易处理性 (PT)、 拟多项式易处理性 (QPT)、一致弱易处理性 (UWT)、 弱易处理性 (WT) 和 (s,t)-弱易处理性 ((s,t)-WT).在本文中, 我们主要研究指数收敛易处理性.

3 多元 \mathbb{L}_{2}-逼近问题及关系

加权 Korobov 空间 H(K_{d,\textbf{a},\textbf{b}}) 中的多元 \mathbb{L}_{2}-逼近问题已经被广泛研究,参见文献 [2, 5, 19, 24, 26]. 该问题

\begin{equation*} \text{APP}_{d,2}:H(K_{d,\textbf{a},\textbf{b}})\rightarrow L_{2}([0,1]^{d}), \text{APP}_{d,2}(f)=f \end{equation*}

是由 H(K_{d,\textbf{a},\textbf{b}})L_{2}([0,1]^{d}) 的嵌入逼近.对于多元 \mathbb{L}_{2}-逼近问题, 我们仍然只需使用线性算法

\begin{equation*} A_{n,d}(f)=\sum_{k=1}^{n}\alpha_{k}L_{k}(f), \quad f\in H(K_{d,\textbf{a},\textbf{b}}), \end{equation*}

其中 L_{k}\in\Lambda^{\text{all}}=H^{*}(K_{d,\textbf{a},\textbf{b}})L_{k}\in\Lambda^{\text{std}}=\{L|L(g)=g(\textbf{x}), \forall \ g\in H(K_{d,\textbf{a},\textbf{b}}),\textbf{x}\in[0,1]^{d}\},1\leq k \leq n, \alpha_{k}L_{2}([0,1]^{d}) 中的元素.

\mathbb{L}_{\infty} 情形类似, 算法 A_{n,d} 在最坏框架下的误差定义为

\begin{equation*} e^{L_2}(A_{n,d},H(K_{d,\textbf{a},\textbf{b}})):=\mathop{ \sup\limits_{f\in H(K_{d,\textbf{a},\textbf{b}})}}_{\|f\|_{H(K_{d,\textbf{a},\textbf{b}})}\leq1} \| f-A_{n,d}(f) \|_{L_{2}([0,1]^{d})}, \end{equation*}

n 阶最小最坏框架误差定义为

e^{L_2,\Lambda}(n,d)=\inf_{A_{n,d}}e^{L_2}(A_{n,d},H(K_{d,\textbf{a},\textbf{b}})),

这里的下确界选自所有线性算法 A_{n,d}.n=0 时, 文献[2] 已得到初始误差

e^{L_2,\Lambda}(0,d)=1.

因此, 多元 \mathbb{L}_{2}-逼近问题的绝对误差和归一化误差标准没有区别.

对于 \varepsilon\in(0,1), d\in \mathbb{N} , \Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\},多元 \mathbb{L}_{2}-逼近的信息复杂度在绝对和归一化误差标准下的定义为

n^{L_2,\Lambda}(\varepsilon,d):= \min\{n\in \mathbb{N} :e^{L_2,\Lambda}(n,d)\leq \varepsilon \}.

由文献[7]知绝对误差标准下多元 \mathbb{L}_{2}-逼近的难度不比多元\mathbb{L}_{\infty}-逼近难, 即对于\Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\}

\begin{equation} n^{L_2,\Lambda}(\varepsilon,d)\leq n^{L_\infty,\Lambda}_{\text{abs}}(\varepsilon,d). \end{equation}
(3.1)

多元 \mathbb{L}_{2}-逼近问题的指数收敛易处理性概念 EC-SPT, EC-PT, EC-QPT, EC-UWT, EC-WT 和 EC-(s,t)-WT与多元 \mathbb{L}_{\infty}-逼近的类似, 只需要把定义中的 n^{L_\infty,\Lambda}_{X}(\varepsilon,d) 替换成n^{L_2,\Lambda}(\varepsilon,d) 即可. 关于 H(K_{d,\textbf{a},\textbf{b}}) 中的多元\mathbb{L}_{2}-逼近指数收敛易处理性结果已在文献 [2,5,19,24] 中得到研究.

注意到对于线性信息类 \Lambda^{\text{all}}, 其多元 \mathbb{L}_{2}-逼近和\mathbb{L}_{\infty}-逼近问题可以通过算子 W_{d} 的特征值得到很好地刻画, 其中W_{d}=\text{APP}_{d,2}^{*} \text{APP}_{d,2}:H(K_{d,\textbf{a},\textbf{b}})\longrightarrow H(K_{d,\textbf{a},\textbf{b}}),

W_{d}(f)=\int_{[0,1]^{d}}f(\textbf{t})K_{d,\textbf{a},\textbf{b}}(\cdot,\textbf{t})d\textbf{t}.

多元 \mathbb{L}_{2}-逼近的刻画结果参见文献[11,18],多元 \mathbb{L}_{\infty}-逼近的刻画结果已在文献[9] 中得到讨论.

事实上, 对于 \textbf{h}\in \mathbb{Z} ^{d}, 函数 e_{\textbf{h}} 定义为

e_{\textbf{h}}(\textbf{x})=(\omega_{\textbf{h}})^{1/2}\exp(2\pi {\rm i}\textbf{h}\cdot\textbf{x}).

此时, \{e_{\textbf{h}}\}_{\textbf{h}\in \mathbb{Z} ^{d}} 就是 Korobov 空间H(K_{d,\textbf{a},\textbf{b}}) 的完备正交基. 容易证明算子 W_{d} 的特征对为(\omega_{\textbf{h}},e_{\textbf{h}}), 即

W_{d}e_{\textbf{h}}=\omega_{\textbf{h}} e_{\textbf{h}}, \quad \forall \ \textbf{h}\in \mathbb{Z} ^{d}.

算子 W_d 非增重排后的特征值记为 \{\lambda_{d,k}\}_{k\in \mathbb{N} }, 即

\lambda_{d,1}\geq \lambda_{d,2}\geq \cdots \geq \lambda_{d,k}\geq\cdots\geq 0.

显然, \{\lambda_{d,k}\}_{k\in \mathbb{N} }=\{\omega_{\textbf{h}}\}_{\textbf{h}\in \mathbb{Z} ^{d}}\lambda_{d,1}=1.从而

e^{L_2,\Lambda^{\text{all}}}(n,d)=\lambda_{d,n+1}^{1/2}, \quad e^{L_\infty,\Lambda^{\text{all}}}(n,d)= \left(\sum\limits_{k=n+1}^{\infty} \lambda_{d,k} \right)^{1/2}.

\text{CRI}_{\text{abs}}=1, \quad \text{CRI}_{\text{nor}}=e^{L_\infty,\Lambda^{\text{all}}}(0,d)= \prod_{j=1}^{d} \left( 1+2\sum\limits_{h=1}^{\infty}\omega^{a_{j}h^{b_{j}}}\right)^{1/2},

X\in\{\text{abs},\text{nor} \} 时,

\begin{equation*} n^{L_2,\Lambda^{\text{all}}}(\varepsilon,d):= \min\{n\in \mathbb{N} : \lambda_{d,n+1} \leq \varepsilon^{2} \}, \end{equation*}
\begin{equation*} n^{L_\infty,\Lambda^{\text{all}}}_{X}(\varepsilon,d):= \min\left\{n\in \mathbb{N} : \sum\limits_{k=n+1}^{\infty} \lambda_{d,k} \leq \varepsilon^{2}\text{CRI}^{2}_{X} \right\}. \end{equation*}

另外, 对于多元 \mathbb{L}_{2}-逼近和 \mathbb{L}_{\infty}-逼近, 其 n 阶最小误差都由同一个算法

A_{n,d}^{*}(f)=\sum_{k=1}^{n} \langle f,\eta_{d,k}\rangle_{H(K_{d,\textbf{a},\textbf{b}})}\eta_{d,k}

获得, 这里 \eta_{d,k} 表示相应于特征值 \lambda_{d,k} 的特征函数[9].

4 主要结论

对于 \Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\} 及绝对和归一化误差标准,文献[7] 最早开始研究最坏框架下定义在 H(K_{d,\textbf{a},\textbf{b}})中的多元 \mathbb{L}_{\infty}-逼近问题的指数收敛易处理性. 该文献主要研究EC-SPT, EC-PT 和 EC-(1,t)-WT(t\geq 1), 最终通过权参数得到匹配的充分必要条件.最近, 文献[22] 进一步研究了多元 \mathbb{L}_{p}-逼近问题的指数收敛易处理性,其中 1\leq p\leq \infty. 对于线性信息类 \Lambda^{\text{all}} 及绝对和归一化误差标准,他们得到了 EC-SPT, EC-PT, EC-UWT 和 EC-WT 成立的条件. 下面给出多元\mathbb{L}_{\infty}-逼近问题的一些指数收敛易处理性结果.

对于 \Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\}及绝对和归一化误差标准[7],

·EC-SPT 成立当且仅当 EC-PT 成立当且仅当

\begin{equation*} \sum\limits_{j=1}^{\infty}\frac{1}{b_{j}}<\infty,\quad \liminf_{j\rightarrow\infty}\frac{\ln a_{j}}{j}>0. \end{equation*}

·EC-(1,t)-WT(t\geq 1) 成立当且仅当

\begin{equation*} \lim_{j\rightarrow\infty}a_{j}=\infty. \end{equation*}

对于 \Lambda^{\text{all}}及绝对和归一化误差标准[22]

·EC-UWT 成立当且仅当

\begin{equation*} \lim_{j\rightarrow\infty}\frac{\ln a_{j}}{\ln j}=\infty. \end{equation*}

对于 EC-SPT, EC-PT 和 EC-(1,t)-WT(t\geq 1),容易发现它们的结果在 \Lambda^{\text{all}}\Lambda^{\text{std}} 之间以及绝对误差和归一化误差标准之间并没有什么区别.但对于 EC-UWT, 并不知道其在 \Lambda^{\text{all}}\Lambda^{\text{std}}下的条件是否一致. 本文中, 我们对该公开的问题给出肯定的回答. 另外, 本文也研究EC-QPT 并得到 \Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\}和绝对误差标准下的充分必要条件, 该条件是否也对归一化误差标准成立仍然没有得到解决.我们采用类似文献[7] 中的构造性算法和方法得到 EC-QPT 和 EC-UWT的充分必要条件, 算法使用的采样点是具有不同网格尺寸的正则网格.这些正则网格之前在文献 [2,3] 中被用来解决多元逼近和积分问题.

定理4.1 考虑定义在 Korobov 空间 H(K_{d,\textbf{a},\textbf{b}}) 中的多元\mathbb{L}_{\infty}-逼近问题, 其中序列 \textbf{a}\textbf{b}满足 (2.1) 式.对于 \Lambda\in\{\Lambda^{\text{all}},\Lambda^{\text{std}}\} 和绝对误差标准, EC-QPT 成立当且仅当

\begin{equation} B^{*}:=\sup_{d\in\mathbb{N}}\frac{\sum_{j=1}^{d}b_{j}^{-1}}{1+\ln d}<\infty, \quad \quad \alpha_{*}:=\liminf_{j\rightarrow\infty}\frac{(1+\ln j)\ln a_{j}}{j}>0. \end{equation}
(4.1)

如果上式成立, 则 EC-QPT 的指数满足

t^{*}\in\left[\max\left(B^{*},\frac{\ln 2}{\alpha_{*}} \right), B^{*}+\frac{C}{\alpha_{*}}\right],

其中 C 为大于零的常数. 特别地, 若 \alpha_{*}=\inftyt^{*}=B^{*}.

定理4.2 考虑定义在 Korobov 空间 H(K_{d,\textbf{a},\textbf{b}}) 中的多元\mathbb{L}_{\infty}-逼近问题, 其中序列 \textbf{a}\textbf{b}满足(2.1) 式.对于 \Lambda^{\text{std}} 和绝对误差或归一化误差标准, EC-UWT 成立当且仅当

\begin{equation} A:=\lim_{j\rightarrow\infty}\frac{\ln a_{j}}{\ln j}=\infty. \end{equation}
(4.2)

5 标准信息类 \Lambda^{\text{std}} 下的预备知识

对于固定的 \omega\in(0,1), K_{d,\textbf{a},\textbf{b}} 为 (2.3)式定义的解析 Korobov 核, 且 \textbf{a},\textbf{b} 满足 (2.1) 式.由 (2.2) 式知再生核 Hilbert 空间 H(K_{d,\textbf{a},\textbf{b}})为单变量再生核 Hilbert 空间 H(K_{1,a_{j},b_{j}}) 的张量积, 其中 K_{1,a_{j},b_{j}}, j=1,\cdots,d为其再生核, 即

H(K_{d,\textbf{a},\textbf{b}})=H(K_{1,a_{1},b_{1}})\otimes H(K_{1,a_{2},b_{2}}) \otimes \cdots \otimes H(K_{1,a_{d},b_{d}}).

在证明我们的结果之前, 首先介绍关于标准信息类 \Lambda^{\text{std}} 的一些预备知识.文献 [23] 提供了一种基于函数值的特殊线性逼近算法, 利用该算法可得到误差界.对于给定的点集 P=\{\textbf{x}_{1}, \cdots,\textbf{x}_{n} \} 和函数值\{f(\textbf{x}_{1}), \cdots,f(\textbf{x}_{n}) \}, 样条 \sigma 定义为

\sigma(f;P):=\text{argmin}\left\{\|g\|_{H(K_{d,\textbf{a},\textbf{b}})} : g\in H(K_{d,\textbf{a},\textbf{b}}),g(\textbf{x}_{k})=f(\textbf{x}_{k}),k=1,2,\cdots,n \right\}.

于是可以利用样条 \sigma(f;P)\mathbb{L}_\infty 范数下逼近 f\in H(K_{d,\textbf{a},\textbf{b}}).而且, 该逼近算法使用的是 d 维网格 g_{n,d},此种正则网格在文献 [2,3,6,7] 中已被广泛讨论和研究.下面给出它们的定义.

给定 d\in \mathbb{N} , 具有网格尺寸 m_{1},\cdots,m_{d}\in \mathbb{N} 的正则网格 g_{n,d} 定义为如下的点集

g_{n,d}=\{(k_{1}/m_{1},\cdots,k_{d}/m_{d}) :k_{j}=0,1,\cdots,m_{j}-1, \ j=1,2,\cdots,d\},

其中 n=\prod_{j=1}^{d} m_{j}g_{n,d} 的基数. 用 g_{n,d}^{\bot} 表示 g_{n,d} 的对偶集合, 即

g_{n,d}^{\bot}=\{\textbf{l}=(l_1,\cdots,l_d)\in \mathbb{Z} ^{d}: l_{j}\equiv 0\ (\text{mod}\ m_{j}), \ j=1,2,\cdots,d\}.

对于正则网格 g_{n,d}, 集合 \mathbb{Z} ^{d} 可以表示成 g_{n,d}^{\bot} 和集合

V_{n}= \mathbb{Z} ^{d}\cap\prod_{j=1}^{d} \left(-\frac{m_j}{2}, \frac{m_j}{2}\right]

的直和, 即

\mathbb{Z} ^{d}=V_{n}\oplus g_{n,d}^{\bot} =\{\textbf{v}+\textbf{l}: \textbf{v} \in V_{n}, \textbf{l}\in g_{n,d}^{\bot}\}.

给定取自 g_{n,d} 的点 \textbf{x}_{1}, \cdots,\textbf{x}_{n},由文献 [23] 知样条 \sigma(f;P) 可以由基函数\phi_{k}, k=1,2,\cdots,n 表示, 这里 \phi_{k}K_{d,\textbf{a},\textbf{b}}(\cdot,\textbf{x}_{r}) 的线性组合. 更精确地说

\begin{aligned} \sigma(f; g_{n,d})(\textbf{x})=\sum\limits_{k=1}^{n}f(\textbf{x}_{k})\phi_{k}(\textbf{x}), \end{aligned}
(5.1)
\phi_{k}(\textbf{x})=\sum\limits_{r=1}^{n} K_{d,\textbf{a},\textbf{b}}(\textbf{x},\textbf{x}_{r})\xi_{r,k},

其中 \{\xi_{r,k}\} 是通过由 Kronecker delta 函数 \delta 表示的一个条件给出, 即

\delta_{j,k}=\sum\limits_{r=1}^{n} K_{d,\textbf{a},\textbf{b}}(\textbf{x}_{j},\textbf{x}_{r})\xi_{r,k}.

按照文献 [23,3.1 节] 类似的步骤, 文献 [7] 利用上面的正则网格最终得到如下估计

\begin{aligned} \left[ e^{L_\infty,\Lambda^{\text{std}}}(\sigma(\cdot; g_{n,d}),H(K_{d,\textbf{a},\textbf{b}})) \right]^{2} \leq 4\sum_{\textbf{h}\notin V_{n}}\omega_{\textbf{h}} = 4\sum_{\textbf{v}\in V_{n}}\sum_{\textbf{l}\in g_{n,d}^{\bot}\backslash \{\textbf{0}\}}\omega_{\textbf{v+l}}. \end{aligned}
(5.2)

6 定理 4.1 的证明

对于指数易处理性概念 EC-QPT, 容易验证当 d\in \mathbb{N} , \varepsilon\in(0,1) 时,

\begin{aligned} \exp\left(t(1+\ln d)(1+\ln(1+\ln\varepsilon^{-1}))\right)=(ed)^{t(1+\ln(1+\ln\varepsilon^{-1}))} =(e(1+\ln\varepsilon^{-1}))^{t(1+\ln d)}, \end{aligned}

其中 e=\exp{(1)}. 在下面的证明过程中将会用到上面的等价式子.

首先证明由 EC-QPT 成立得到条件 (4.1).由关系式 (2.5) 和 (3.1) 知只需考虑多元\mathbb{L}_{2}-逼近问题和线性信息类 \Lambda^{\text{all}}.因此由文献 [5] 得 (4.1) 式成立且 EC-QPT 的指数 t^{*}满足

t^{*}\geq\max\left(B^{*},\frac{\ln 2}{\alpha_{*}} \right).

下面证明条件 (4.1) 可以推出绝对误差标准下多元\mathbb{L}_{\infty}-逼近问题的 EC-QPT 成立.由关系式 (2.5) 知只需考虑标准信息类 \Lambda^{\text{std}}. 为了得到想要的结果,需要使用 (5.1) 式给出的算法 \sigma(\cdot; g_{n,d}), 其中采样点 \{\textbf{x}_{k}\}_{k=1}^{n}来自于正则网格 g_{n,d}. 而且, g_{n,d} 的网格尺寸 m_{1}, m_{2},\cdots,m_{d}

\begin{equation} m_{j}=2\left\lceil\left(\frac{x(M)}{a_{j}^{\beta}}\right) ^{1/b_{j}}\right\rceil-1, \quad j=1,2,\cdots,d, \end{equation}
(6.1)

这里 M>1, \beta\in(0,1),

x(M):=\frac{\ln M}{\ln\omega^{-1}}>0.

由 (6.1) 式知 m_{j}\geq1m_{j} 总为奇数. 另外, 容易证得如果 a_{j}\geq(x(M))^{1/\beta}, 则 m_{j}=1. 因为 \alpha_{*}>0, 这意味着对于每个 \tau\in (0,\alpha_{*}), 总存在正整数 N^{*}_{\tau} 使得 j\geq N^{*}_{\tau} 时,

\begin{equation} a_{j}\geq\exp\left(\frac{\tau j}{1+\ln j}\right). \end{equation}
(6.2)

对于 \ln x(M)\geq\beta\tau, 即 M\geq \omega^{-e^{\beta\tau}}, 由 (6.2) 式知, 若

j\geq N^{*}_{\beta,\tau, M}: =\max\left(N_{\tau}^{*},J(M)\right),

m_{j}=1, 其中 J(M) 为非线性方程

\begin{equation*} \frac{\ln x(M)}{\beta\tau}=\frac{J(M)}{1+\ln J(M)} \end{equation*}

的解. 因为 y\geq 1y/(1+\ln y) 是递增的, 故该方程的解是唯一的.而且由文献 [5]可知 M\rightarrow\infty 时,

\begin{equation*} J(M)=\frac{1+o(1)}{\beta\tau}[\ln\ln M]\ln\ln\ln M. \end{equation*}

由 (5.2) 式知

e_{n,d}^{2}:=\left[ e^{L_\infty,\Lambda^{\text{std}}}(\sigma(\cdot; g_{n,d}),H(K_{d,\textbf{a},\textbf{b}})) \right]^{2}\leq 4\sum_{\textbf{v}\in V_{n}}\sum_{\textbf{l}\in g_{n,d}^{\bot}\backslash \{\textbf{0}\}}\omega_{\textbf{v+l}}.

类似文献 [命题 1]的证明, 最终得到

\begin{equation} e_{n,d}^{2}\leq 4n \left[-1+\prod_{j=1}^{d}\left(1+2A_{1}\omega^{a_{j}[(m_{j}+1)/2]^{b_{j}}} \right)\right], \end{equation}
(6.3)

其中 A_{1}:=\sum_{l=1}^{\infty}\omega^{a_{1}(l^{b^{*}}-1)}<\infty.

首先估计不等式 (6.3) 中右边的因式. 利用 x\geq 0\ln (1+x)\leq x可得

\ln\left[\prod_{j=1}^{d}\left(1+2A_{1}\omega^{a_{j}[(m_{j}+1)/2]^{b_{j}}}\right) \right] \leq 2A_{1}\sum_{j=1}^{d}\omega^{a_{j}[(m_{j}+1)/2]^{b_{j}}}=:\gamma.

由 (6.1) 式 m_{j} 的定义得a_{j}[(m_{j}+1)/2]^{b_{j}}\geq a_{j}^{1-\beta}x(M)= a_{j}^{1-\beta}(\ln M)/\ln\omega^{-1}.因此

\omega^{a_{j}[(m_{j}+1)/2]^{b_{j}}}\leq \omega^{a_{j}^{1-\beta}(\ln M)/\ln\omega^{-1}}=\left(\frac{1}{M}\right)^{a_{j}^{1-\beta}}.

不失一般性, 假设 M\geq 2. 因为 j\leq N_{\tau}^{*}-1, a_{j}\geq a_{1},j\geq N_{\tau}^{*}a_{j}\geq\exp\left(\frac{\tau j}{1+\ln j}\right), 从而

\begin{equation} \gamma\leq 2A_{1}\left(\frac{N_{\tau}^{*}-1}{M^{a_{1}^{1-\beta}}} +\sum_{j=N_{\tau}^{*}}^{\infty}\left(\frac{1}{M}\right) ^{\exp\left(\frac{(1-\beta)\tau j}{1+\ln j}\right)}\right) \leq\frac{C_{\beta,\tau}}{M^{\min(a_{1}^{1-\beta},1)}}, \end{equation}
(6.4)

其中

C_{\beta,\tau}:=2A_{1}\left(N_{\tau}^{*}-1+ \sum_{j=N_{\tau}^{*}}^{\infty}\left(\frac{1}{2}\right)^{\exp\left(\frac{(1-\beta)\tau j}{1+\ln j}\right)-1} \right) <\infty.

由 (6.4) 式容易验证当 M\geq C_{\beta,\tau}^{1/\min(a_{1}^{1-\beta},1)}\gamma\leq 1.而且, 利用函数的凸性可以证得对所有 \gamma\in[0,1], -1+\exp(\gamma)\leq (e-1)\gamma.因此, 当 M\geq C_{\beta,\tau}^{1/\min(a_{1}^{1-\beta},1)} 时,

\begin{equation*} \begin{aligned} -1+\prod_{j=1}^{d}\left(1+2A_{1}\omega^{a_{j}[(m_{j}+1)/2]^{b_{j}}}\right) &\leq-1+\exp\left(2A_{1}\sum_{j=1}^{d}\omega^{a_{j}[(m_{j}+1)/2]^{b_{j}}} \right)\\ &=-1+\exp{(\gamma)}\leq (e-1)\gamma\\ &\leq\frac{(e-1)C_{\beta,\tau}}{M^{\min(a_{1}^{1-\beta},1)}}. \end{aligned} \end{equation*}

从而

\begin{equation} e_{n,d}^{2}\leq 4n \frac{(e-1)C_{\beta,\tau}}{M^{\min(a_{1}^{1-\beta},1)}}. \end{equation}
(6.5)

下面估计使用算法 \sigma(\cdot; g_{n,d}) 时所需要的函数值数量 n.x(M)\geq a_{1}^{\beta}, 即 M\geq \omega^{-a_{1}^{\beta}} 时,

\begin{matrix}\label{nn} n&=\prod_{j=1}^{d}m_{j}=\prod_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})}m_{j}\\ &\leq \prod_{j=1}^{N_{\beta,\tau,M}^{*}} \left(1+2\left(\frac{x(M)}{a_{j}^{\beta}}\right)^{1/b_{j}}\right) \leq \prod_{j=1}^{N_{\beta,\tau,M}^{*}}\left(3\left(\frac{x(M)}{a_{1}^{\beta}}\right)^{1/b_{j}}\right). \end{matrix}
(6.6)

因为 N^{*}_{\beta,\tau, M}\leq N_{\tau}^{*}+J(M), 从而由 (6.6) 式得

\begin{equation} \begin{aligned} n\leq 3^{N^{*}_{\beta,\tau, M}}\left(\frac{1}{a_{1}^{\beta}} \right)^{\frac{N^{*}_{\beta,\tau, M}}{b_{*}}}\left(\frac{\ln M}{\ln\omega^{-1}}\right)^{\frac{N^{*}_{\beta,\tau, M}}{b_{*}}}\leq \left(\frac{3^{b_{*}}}{a_{1}^{\beta}} \cdot \frac{\ln M}{\ln\omega^{-1}} \right)^{\frac{N_{\tau}^{*}+J(M)}{b_{*}}}. \end{aligned} \end{equation}
(6.7)

M\geq\max(2,C_{\beta,\tau}^{1/\min(a_{1}^{1-\beta},1)},\omega^{-e^{\tau}},\omega^{-a_{1}^{\beta}}) 时, 结合 (6.5) 和 (6.7) 式得

e_{n, d}^{2} \leq \frac{1}{M^{\min \left(a_{1}^{1-\beta}, 1\right)}}\left[4(\mathrm{e}-1) C_{\beta, \tau}\left(\frac{3^{b_{*}}}{a_{1}^{\beta}} \cdot \frac{\ln M}{\ln \omega^{-1}}\right)^{\frac{N_{\tau}^{*}+J(M)}{b_{*}}}\right] \leq \frac{D_{\beta \tau}}{M^{\min \left(a_{1}^{1-\beta}, 1\right) / 2}},

其中

D_{\beta\tau}:=\sup_{y\geq C_{\beta,\tau}^{1/\min(a_{1}^{1-\beta},1)}} \left[\frac{4(e-1)C_{\beta\tau}}{y^{\min(a_{1}^{1-\beta},1)/2}} \left( \frac{3^{b_{*}}}{a_{1}^{\beta}} \cdot \frac{\ln y}{\ln\omega^{-1}} \right) ^{\frac{N_{\tau}^{*}+J(y)}{b_{*}}} \right].

因为 y\rightarrow\inftyJ(y)=\frac{1+o(1)}{\beta\tau}[\ln\ln y]\ln\ln\ln y,由此可以得到 D_{\beta,\tau}<+\infty. 现在设

\begin{equation} M=\max\left(2,C_{\beta,\tau}^{1/\min(a_{1}^{1-\beta},1)},\omega^{-e^{\tau}},\omega^{-a_{1}^{\beta}}, D_{\beta,\tau}^{2/\min(a_{1}^{1-\beta},1)}\varepsilon^{-4/\min(a_{1}^{1-\beta},1)} \right), \end{equation}
(6.8)

从而

e_{n,d}=e^{L_\infty,\Lambda^{\text{std}}}(\sigma(\cdot; g_{n,d}),H(K_{d,\textbf{a},\textbf{b}})) \leq \varepsilon.

另外, 由 (6.8) 式知 \varepsilon\rightarrow 0\ln M=\mathcal{O}(\ln\varepsilon^{-1}).因此, 对于足够小的 \varepsilon 存在 C_{1}>1 使得

\begin{equation} \ln M\leq C_{1}\ln\varepsilon^{-1}. \end{equation}
(6.9)

由信息复杂度的定义和 (6.9) 式可知, 对于足够小的 \varepsilon

\begin{matrix}\label{nn33} n_{\text{abs}}^{L_\infty,\Lambda^{\text{std}}}(\varepsilon,d)&\leq n=\prod_{j=1}^{d}m_{j}=\prod_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})}m_{j} \leq \prod_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})} \left(1+2\left(\frac{\ln M}{a_{j}^{\beta}\ln\omega^{-1}}\right)^{1/b_{j}}\right)\\ &\leq 3^{\min(d,N_{\beta,\tau,M}^{*})}\prod_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})} \left(\frac{C_{1}\ln \varepsilon^{-1}}{a_{j}^{\beta}\ln\omega^{-1}} \right)^{1/b_{j}}\\ &\leq 3^{\min(d,N_{\beta,\tau,M}^{*})}C_{d,M}[e \ln \varepsilon^{-1}]^ {\sum_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})}b_{j}^{-1}}, \end{matrix}
(6.10)

其中

C_{d,M}:=\prod_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})} \left(\frac{C_{1}}{a_{j}^{\beta}e\ln\omega^{-1}}\right)^{1/b_{j}}.

由 (6.2) 式知 \lim a_{j}=\infty, 这意味着 C_{d,M} 中仅有有限个项大于 1.从而

\begin{equation} C_{d,M}\leq C_{2}:=\sup_{d\in \mathbb{N},M\in(1,\infty)}\prod_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})} \left(\frac{C_{1}}{a_{j}^{\beta}e\ln\omega^{-1}}\right)^{1/b_{j}}<\infty. \end{equation}
(6.11)

而且, 由假设条件 B^{*}<\infty

\begin{equation} \sum_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})}b_{j}^{-1}= \frac{\sum_{j=1}^{\min{(d,N_{\beta,\tau,M}^{*}})}b_{j}^{-1}}{1+\ln d} (1+\ln d)\leq B^{*}(1+\ln d). \end{equation}
(6.12)

现在结合 (6.11) 与 (6.12) 式, 可将 (6.10) 式变为

\begin{equation} n_{\text{abs}}^{L_\infty,\Lambda^{\text{std}}}(\varepsilon,d) \leq C_{2}\cdot 3^{\min(d,N_{\beta,\tau,M}^{*})} [e (1+\ln \varepsilon^{-1})]^{B^{*}(1+\ln d)}. \end{equation}
(6.13)

事实上, 我们可以使用类似文献 [5] 的方法去分析 (6.13) 式中的因式\rho:=3^{\min(d,N_{\beta,\tau,M}^{*})} 并最终得到

\rho\leq\max\left\{[\textrm{e} d]^{(\beta\tau)^{-1}(1+2\eta)(\ln 3) \left(1+\ln(1+\ln M)\right)}, [\textrm{e}(1+\ln M)]^{(\beta\tau)^{-1}(1+2\eta)(\ln 3)(1+\ln d) } \right\},

上式可等价地写成

\begin{equation} \rho\leq \exp\left(\left[\frac{1+2\eta}{\beta\tau}\ln 3\right](1+\ln d)(1+\ln(1+\ln M))\right). \end{equation}
(6.14)

因为对于足够小的 \varepsilon 存在 C_{1}>1 使得 \ln M\leq C_{1}\ln\varepsilon^{-1},故由 (6.14) 式得

\begin{matrix}\label{rhoineq2} \rho&\leq \exp\left(\left[\frac{1+2\eta}{\beta\tau}\ln 3\right](1+\ln d)(1+\ln(C_{1}+C_{1}\ln\varepsilon^{-1}))\right)\\ &\leq \exp\left(C_{3}\left[\frac{1+2\eta}{\beta\tau}\ln 3\right](1+\ln d)(1+\ln(1+\ln\varepsilon^{-1}))\right), \end{matrix}
(6.15)

其中 C_{3}:=1+\ln C_{1}.

由 (6.13) 和 (6.15) 式知 EC-QPT 成立, 且t\leq B^{*}+C_{3}(\beta\tau)^{-1}(1+2\eta)\ln 3.由于 \tau 可任意趋于 \alpha_{*}, \eta 可任意小,且 \beta 可任意趋于 1, 故得 EC-QPT 的指数满足t^{*}\leq B^{*}+\frac{C_{3}\ln 3}{\alpha_{*}}.

7 定理 4.2 的证明

首先证明由 EC-UWT 成立得到条件 (4.2).由 (2.5) 式可知只需考虑多元\mathbb{L}_{\infty}-逼近问题和线性信息类 \Lambda^{\text{all}}.因此, 由文献[22]可知条件 (4.2) 成立.

下面证明当使用标准信息类 \Lambda^{\text{std}} 时,条件 (4.2) 可以推出多元 \mathbb{L}_{\infty}-逼近问题的 EC-UWT 成立.由关系式 (2.6) 知, 只需证明绝对误差标准下的 EC-UWT 成立.需要强调的是我们使用的算法和正则网格与定理 4.1 的证明过程中一样.由于 A=\infty, 故对任意 \delta>0 存在正整数 N_{\delta}^{*}, 使得 j\geq N_{\delta}^{*} 时,

\begin{equation} a_{j}\geq j^{\delta}. \end{equation}
(7.1)

注意到 a_{j}\geq(x(M))^{1/\beta}m_{j}=1, 因此由 (7.1) 式可得当

j\geq N_{\beta,\delta,M}^{*}:=\max\left(N_{\delta}^{*}, (x(M))^{\frac{1}{\beta\delta}} \right)

时, m_{j}=1.

类似于定理 4.1 的证明方法, 得到

\begin{equation} e_{n,d}^{2}\leq 4n\left(-1+\exp(\gamma)\right), \end{equation}
(7.2)

其中

\begin{equation} \gamma=2A_{1}\sum_{j=1}^{d}\omega^{a_{j}[(m_{j}+1)/2]^{b_{j}}} \leq 2A_{1}\sum_{j=1}^{d}\left(\frac{1}{M}\right)^{a_{j}^{1-\beta}}. \end{equation}
(7.3)

不失一般性, 假设 M\geq 2. 因为 j\leq N_{\delta}^{*}-1a_{j}\geq a_{1},j\geq N_{\delta}^{*}a_{j}>j^{\delta}, 因此由 (7.3) 式得

\begin{equation} \gamma\leq 2A_{1}\left(\frac{N_{\delta}^{*}-1}{M^{a_{1}^{1-\beta}}} +\sum_{j=N_{\delta}^{*}}^{\infty}\left(\frac{1}{M}\right)^{j^{\delta(1-\beta)}}\right) \leq\frac{C_{\beta,\delta}}{M^{\min(a_{1}^{1-\beta},1)}}, \end{equation}
(7.4)

其中

C_{\beta,\delta}:=2A_{1}\left(N_{\delta}^{*}-1+ \sum_{j=N_{\delta}^{*}}^{\infty}\left(\frac{1}{2}\right)^{j^{\delta(1-\beta)}-1} \right) <\infty.

由于 M\geq C_{\beta,\delta}^{1/\min(a_{1}^{1-\beta},1)}\gamma\leq 1,故结合 (7.2) 和 (7.4) 式可得

\begin{equation} e_{n,d}^{2}\leq 4n(e-1)\gamma \leq 4n\frac{C_{\beta,\delta}(e-1)}{M^{\min(a_{1}^{1-\beta},1)}}, \end{equation}
(7.5)

这里第一个不等式利用了 \gamma\in[0,1]-1+\exp(\gamma)\leq (e-1)\gamma 的结论.

现在估计使用算法 \sigma(\cdot; g_{n,d}) 时所需要的函数值数量 n.因为 N_{\beta,\delta,M}^{*}\leq N_{\delta}^{*}+(x(M))^{\frac{1}{\beta\delta}},故由 (6.6) 式得

\begin{equation} \begin{aligned} n\leq 3^{N^{*}_{\beta,\tau, M}}\left(\frac{1}{a_{1}^{\beta}} \right)^{\frac{N^{*}_{\beta,\tau, M}}{b_{*}}}\left(\frac{\ln M}{\ln\omega^{-1}}\right)^{\frac{N^{*}_{\beta,\tau, M}}{b_{*}}}\leq \left(\frac{3^{b_{*}}}{a_{1}^{\beta}} \cdot \frac{\ln M}{\ln\omega^{-1}} \right)^{\frac{ N_{\delta}^{*}+\left(\frac{\ln M}{\ln\omega^{-1}}\right)^{\frac{1}{\beta\delta}}}{b_{*}}}. \end{aligned} \end{equation}
(7.6)

因此, 当 M\geq\max(2,C_{\beta,\delta}^{1/\min(a_{1}^{1-\beta},1)},\omega^{-a_{1}^{\beta}}) 时,结合 (7.5) 和 (7.6) 式得

\begin{aligned} e_{n, d}^{2} & \leq \frac{1}{M^{\min \left(a_{1}^{1-\beta}, 1\right)}}\left(4 n C_{\beta, \delta}(\mathrm{e}-1)\right) \\ & \leq \frac{1}{M^{\min \left(a_{1}^{1-\beta}, 1\right)}}\left[4 C_{\beta, \delta}(\mathrm{e}-1)\left(\frac{3^{b_{*}}}{a_{1}^{\beta}} \cdot \frac{\ln M}{\ln \omega^{-1}}\right)^{\frac{N_{\delta}^{*}+\left(\frac{\ln M}{\ln \omega^{-1}}\right)^{\frac{1}{\beta \delta}}}{b_{*}}}\right] \\ & \leq \frac{D_{\beta, \delta}}{M^{\min \left(a_{1}^{1-\beta}, 1\right) / 2}}, \end{aligned}

其中

\begin{equation*} D_{\beta,\delta}:=\sup_{y\geq C_{\beta,\tau}^{1/\min(a_{1}^{1-\beta},1)}} \left( \frac{4C_{\beta,\delta}(e-1) }{y^{\min(a_{1}^{1-\beta},1)/2}} \left(\frac{3^{b_{*}}}{a_{1}^{\beta}} \frac{\ln y}{\ln\omega^{-1}}\right)^{\frac{N_{\delta}^{*} +\left(\frac{\ln y}{\ln\omega^{-1}}\right)^{\frac{1}{\beta\delta}}}{b_{*}}} \right) <\infty, \end{equation*}

\delta>1/\beta.若取

\begin{equation} M=\max\left(2,C_{\beta,\delta}^{1/\min(a_{1}^{1-\beta},1)},\omega^{-a_{1}^{\beta}}, D_{\beta,\delta}^{2/\min(a_{1}^{1-\beta},1)}\varepsilon^{-4/\min(a_{1}^{1-\beta},1)}\right), \end{equation}
(7.7)

e_{n,d}=e^{L_\infty,\Lambda^{\text{std}}}(\sigma(\cdot; g_{n,d}),H(K_{d,\textbf{a},\textbf{b}})) \leq \varepsilon.

由 (7.7) 式知 \varepsilon\rightarrow 0\ln M=\mathcal{O}(\ln\varepsilon^{-1}).从而根据信息复杂度的定义和 (7.6) 式知

n_{\text{abs}}^{L_\infty,\Lambda^{\text{std}}}(\varepsilon,d)\leq n\leq \left(\frac{3^{b_{*}}}{a_{1}^{\beta}} \cdot \frac{\ln M}{\ln\omega^{-1}} \right)^{\frac{ N_{\delta}^{*}+\left(\frac{\ln M}{\ln\omega^{-1}}\right)^{\frac{1}{\beta\delta}}}{b_{*}}} =\mathcal{O} \left[\left(\ln \varepsilon^{-1}\right)^{(\ln \varepsilon^{-1})^{\frac{1}{\beta\delta}}}\right],

其中 \mathcal{O} 仅依赖于 \beta\delta.因此,

\ln n_{\text{abs}}^{L_\infty,\Lambda^{\text{std}}}(\varepsilon,d) =\mathcal{O}\left((\ln\varepsilon^{-1})^{\frac{1}{\beta\delta}}\ln\ln\varepsilon^{-1}\right).

对任意的 s,t\in(0,\infty), 如果选取 \delta>\max(1/\beta,1/(\beta t)), 则

\lim_{d+\varepsilon^{-1}\rightarrow\infty} \frac{\ln n_{\text{abs}}^{L_\infty,\Lambda^{\text{std}}}(\varepsilon,d)}{d^{s} +(\ln\varepsilon^{-1})^{t}}=0.

上面式子意味着 EC-UWT 成立.

参考文献

Chen J.

EC-(t_1,t_2)-tractability of approximation in weighted Korobov spaces in the worst case setting

Journal of Complexity, 2022, 73: 101680

[本文引用: 2]

Dick J, Kritzer P, Pillichshammer F, Woźniakowski H.

Approximation of analytic functions in Korobovspaces

Journal of Complexity, 2014, 30: 2-28

[本文引用: 9]

Dick J, Larcher G, Pillichshammer F, Woźniakowski H.

Exponential convergence and tractability of multivariateintegration for Korobov spaces

Math Comp, 2011, 80: 905-930

[本文引用: 5]

Gnewuch M, Woźniakowski H.

Quasi-polynomial tractability

Journal of Complexity, 2011, 27: 312-330

[本文引用: 1]

Irrgeher C, Kritzer P, Pillichshammer F, Woźniakowski H.

Tractability of multivariate approximation definedover Hilbert spaces with exponential weights

J Approx Theory, 2016, 207: 301-338

[本文引用: 9]

Kritzer P, Pillichshammer F, Woźniakowski H.

Multivariate integration of infinitely many times differentiablefunctions in weighted Korobov spaces

Math Comp, 2014, 83: 1189-1206

[本文引用: 3]

Kritzer P, Pillichshammer F, Woźniakowski H.

\mathbb{L}_\infty-approximation in Korobov spaces with exponentialweights

Journal of Complexity, 2017, 41: 102-125

[本文引用: 10]

Kritzer P, Pillichshammer F, Woźniakowski H.

Exponential tractability of linear weighted tensor productproblems in the worst-case setting for arbitrary linear functionals

Journal of Complexity, 2020, 61: 101501

[本文引用: 1]

Kuo F Y, Wasilkowski G W, Woźniakowski H.

Multivariate \mathbb{L}_\infty-approximation in the worst case settingover reproducing kernel Hilbert spaces

J Approx Theory, 2008, 152: 135-160

[本文引用: 2]

Liu Y P, Zhang J.

EC-tractability of approximation problems in function spaces defined over products ofsimplices

Journal of Complexity, 2019, 55: 101411

[本文引用: 2]

Novak E, Woźniakowski H. Tractability of Multivariate Problems: Volume I: Linear Information. Helsinki:European Mathematical Society, 2008

[本文引用: 3]

Novak E, Woźniakowski H. Tractability of Multivariate Problems: Volume II: Standard Information forFunctionals. Helsinki: European Mathematical Society, 2010

[本文引用: 1]

Novak E, Woźniakowski H. Tractability of Multivariate Problems: Volume III: Standard Information forOperators. Helsinki: European Mathematical Society, 2012

[本文引用: 1]

Siedlecki P.

Uniform weak tractability

Journal of Complexity, 2013, 29: 438-453

[本文引用: 1]

Siedlecki P, Weimar M.

Notes on (s,t)-weak tractability: A refined classification of problems with(sub)exponential information complexity

J Approx Theory, 2015, 200: 227-258

[本文引用: 1]

Sloan I H, Woźniakowski H.

When are Quasi-Monte Carlo algorithms efficient for high-dimensional integrals?

Journal of Complexity, 1998, 14: 1-33

[本文引用: 1]

Sloan I H, Woźniakowski H.

Multivariate approximation for analytic functions with Gaussian kernels

Journal of Complexity, 2018, 45: 1-21

[本文引用: 2]

Traub J, Wasilkowski G, Woźniakowski H. Information-Based Complexity. New York: Academic Press,1988

[本文引用: 3]

Wang H.

A note about EC-(s,t)-weak tractability of multivariate approximation with analytic Korobovkernels

Journal of Complexity, 2019, 55: 101412

[本文引用: 4]

Woźniakowski H.

Tractability and strong tractability of linear multivariate problems

Journal of Complexity, 1994, 10: 96-128

[本文引用: 1]

Woźniakowski H.

Tractability and strong tractability of multivariate tensor product problems

J ComputInform, 1994, 4: 1-19

[本文引用: 1]

Xu G.

EC-tractability of \mathbb{L}_p-approximation in Korobov spaces with exponential weights

J Approx Theory, 2020, 249: 105309

[本文引用: 6]

Zeng X, Kritzer P, Hickernell F J.

Spline methods using integration lattice and digital nets

Constr Approx, 2009, 30: 529-555

[本文引用: 3]

Zhang J.

A note on EC-tractability of multivariate approximation in weighted Korobov spaces for thestandard information class

Journal of Complexity, 2021, 67: 101573

[本文引用: 5]

Zhang J, Liu Y P.

Quasi-Monte Carlo tractability of integration problem in function spaces defined overproducts of balls

Int J Wavelets Multiresolut Inf Process, 2019, 17(6): 1950043

[本文引用: 1]

Zhang J, Liu Y P.

EC-tractability of multivariate approximation in Hermite spaces for the standard informationclass

Int J Wavelets Multiresolut Inf Process, 2022, 20(6): 2250029

[本文引用: 3]

/