在数学分析、优化理论以及计算几何等多个领域,凸函数扮演着核心角色。如何有效地用简单的函数(例如分段线性函数)来逼近复杂的凸函数,是一个基础且重要的问题。这不仅关系到函数的数值表示和计算效率,也深刻联系着凸几何中的一些基本概念。一个经典的范例是凸体的多面体逼近问题:一个凸体可以被一个顶点数给定的外接多面体以何种精度逼近?Gruber等人的开创性工作表明,这种逼近误差的渐近行为与凸体的仿射表面积紧密相关。仿射表面积是凸几何中的一个基本不变量,它度量了凸体边界的“曲率”分布,在Brunn-Minkowski理论中具有核心地位。
自然的问题是,对于定义在Rn上的凸函数,是否存在类似的逼近理论?具体来说,如果我们用一个由有限个(例如m个)仿射函数取最大值构成的、且从下方逼近原凸函数u的凸函数lm来逼近u,那么其加权Lp误差△p(u, P(m)c, ω) = min {∫dom u(u(x) - lm(x))pω(x, u(x)) dx : lm∈ P(m)c} 随着m的增加会如何衰减?其衰减速率和极限常数是否同样揭示了函数u本身的某种内在几何特征?这正是发表在《Advances in Applied Mathematics》上的这项研究旨在回答的核心问题。
为了回答这些问题,研究人员开展了一项理论分析研究。他们考虑定义在Rn的紧子集X上的、足够光滑(例如C+2类)的凸函数u,以及一个连续正权函数ω: Rn+1→ R。研究的目标是刻画当逼近函数lm属于由至多m个仿射函数最大值构成的、且满足lm(x) ≤ u(x)的函数类P(m)c时,加权Lp误差△p(u, P(m)c, ω)在m → +∞时的精确渐近行为。
本研究主要运用了凸分析、几何测度论以及渐近分析的方法。关键技术包括利用Zador定理关于点集在给定测度下最佳逼近的渐近结果,对凸函数的Monge-Ampère测度MA(u; ·)及其与函数Hessian行列式det(D2u)的关系进行分析,以及通过Legendre变换将函数空间Convld(Rn)与其对偶空间ConvMA,h(Rn; R)联系起来。研究还涉及对函数epi-收敛(epi-convergence)和τ**-收敛性的讨论,以确保逼近过程的良好定义。
主要结果
定理1:加权Lp逼近误差的渐近公式
研究表明,对于满足一定正则性条件(如属于C+2类)的凸函数u,其加权Lp逼近误差满足以下精确渐近公式:
limm→+∞m2p/n△p(u, P(m)c, ω) = (δp,n/ 2p) [ ∫dom u(det(D2u(x)))p/(n+2p)ω(x, u(x))n/(n+2p)dx ](n+2p)/n。
其中,常数δp,n是Zador定理中出现的常数,仅依赖于维度n和指数p。这一结果的核心在于,逼近误差的渐近行为完全由函数u的Hessian行列式(即其Monge-Ampère测度的密度)与权函数ω的一个特定加权积分所决定。
推论2:非加权情形的特例
当权函数恒为1(即ω ≡ 1)时,定理1退化为一个更简洁的形式:
limm→+∞m2p/n△p(u, P(m)c, 1) = (δp,n/ 2p) [ ∫dom u(det(D2u(x)))p/(n+2p)dx ](n+2p)/n。
这直接揭示了逼近误差与函数u的“仿射几何”本质之间的联系。
对偶框架与仿射表面积的联系
研究进一步在一个对偶框架下阐述了上述结果。通过Legendre变换,凸函数u与其共轭函数u**之间建立了对应关系。在此对偶框架下,定义了一个泛函Zζ**(v) = ∫Rnζ̃(det(D2v(x))) dx,其中v属于一个特定的凸函数类ConvMA,h(Rn; R)。论文指出,当p=1且权函数取为ω(x,t) = e-t时,定理1中的渐近常数项恰好与加权泛函仿射表面积(weighted functional affine surface area)Z(v) = ∫Rn(det(D2v(x)))1/(n+2)e-n v(x)/(n+2)dx 相关。这表明,本文的逼近误差渐近结果在本质上推广和统一了凸体仿射表面积的经典概念到凸函数的情形。
定理7与定理8:局部化与一致性估计
为了证明主要定理,论文进行了细致的局部化分析。定理7和定理8表明,对于定义在紧集X上的凸函数v,可以在其Monge-Ampère测度的支集supp(MA(v; ·))上考虑逼近问题,并且得到的渐近结果与全局情况一致。这通过构造函数v的局部二次逼近qx(y) = y · D2v(x) y,并利用函数的光滑性保证在足够小的邻域内,qx与在某个固定点ai的二次型qai是“可比”的(即存在常数λ>1使得(1/λ)qai(y) ≤ qx(y) ≤ λ qai(y))。最终,通过将定义域X分割为有限个这样的“一致性”区域Ji,并将Zador定理应用于每个区域上的逼近问题,再综合各区域的结果,从而完成了主要定理的证明。
研究结论与意义
本研究成功地将凸体多面体逼近的经典理论推广到了凸函数的范畴。主要结论是,用从下方逼近的、最多由m个仿射函数最大值构成的凸函数来逼近一个光滑凸函数u时,其加权Lp误差以m-2p/n的速率衰减至零,并且该衰减速率前的常数因子由函数u的Monge-Ampère测度(通过其Hessian行列式表示)的一个特定加权积分完全确定。
这项研究的重要意义体现在多个层面。首先,它在理论上建立了凸函数最佳逼近与凸几何中核心概念——仿射表面积——之间的深刻联系,为理解凸函数的几何结构提供了一个新的视角。当p=1时,该常数泛函恰好对应于凸函数图像的(加权)仿射表面积,这为仿射表面积概念提供了一个自然的函数论解释。其次,该结果为凸函数逼近的算法设计提供了理论极限,指出了任何基于最大仿射函数逼近的方案,其误差衰减速率不可能优于m-2p/n,并且最优逼近器的构造应与函数局部曲率的分布密切相关。此外,文中引入的对偶框架(τ**-收敛)和分析方法(如局部二次逼近、epi-收敛)为处理更广泛的凸优化和变分问题提供了有力的工具。
总之,这项工作不仅深化了我们对凸函数逼近本质的理解,而且架起了函数逼近论、凸几何和Monge-Ampère方程理论之间的桥梁,为后续在这些交叉领域的深入研究奠定了坚实的基础。