形式陈述
“函数很光滑,所以谱方法很快”还缺少什么定量信息?对解析函数,决定几何衰减速度的是它能向复平面延伸多远。
设 , 为 周期函数,并在闭条带 的某个邻域有全纯延拓公理库全纯函数Holomorphic function · Complex analytic function在复平面开集的每一点都复可微的函数。,在该闭条带上满足 。其Fourier 系数公理库Fourier 级数Fourier series · 傅里叶级数把周期函数投影到整数频率的正交指数基上所得的离散频谱展开。满足
因此截取 的尾和满足 。如果只知道开条带内解析,则对任何更窄闭条带分别取有界常数。
区间上的类似几何对象是 Bernstein 椭圆:,。若 在该椭圆及内部的邻域解析,模长不超过 ,则Chebyshev 展开公理库Chebyshev 多项式与节点Chebyshev polynomials and nodes · Chebyshev nodes以余弦表示和端点聚集节点控制区间上一致逼近、插值条件性与高次振荡。 的 系数满足 ,次数 截断的误差不超过 。
这里陈述的是级数截断。节点插值还包含高模态折叠,但在同样解析条件下也有相应几何阶误差界,常数和具体表达需另行核对。
直觉
对正频率 ,将 Fourier 系数积分路径从实轴下移到 。由Cauchy 积分定理公理库Cauchy 积分定理与公式Cauchy integral theorem · Cauchy integral formula全纯函数在可缩闭路上的积分为零,并可由边界积分重建内部值与所有导数。,矩形边界上的积分为零;周期性使两侧竖边抵消,而 多出因子 ,积分剩余部分由 控制。负频率向上移。路径能移多远,取决于是否先遇到复奇点。
Chebyshev 情形令 。椭圆内外邻域的解析性,使 在包含闭环域 的开环域上全纯。Laurent 定理公理库Laurent 级数Laurent series · 洛朗级数全纯函数在指定环域上有唯一的双边幂级数展开,其系数由环域内圆周积分确定。给出唯一展开 ;由 及系数唯一性,有 。又由 上的系数积分公式, 对 成立。利用 ,得到 、()。几何尾和再给出上述截断误差界。
例子与边界
实轴之外的极点怎样决定速度
取 ,实区间 上没有奇点。最近实外极点在 ,对应椭圆参数 。令 ,可验证展开
故相邻非零系数之比恒为 ,每增加一个次数,系数约缩到原来的四分之一。 时所有 ,截断尾部恰为 ,给出可直接核对的误差。
一般有界解析定理只能对 直接应用,因为边界碰到极点时 无界。这个例子的精确系数仍显示极点位置控制的极限衰减率。
光滑与解析的边界
仅有有限阶导数时,常只能得到代数衰减;即使实轴上无限次可微,也未必在复邻域解析。间断函数的系数衰减更慢,并在跳跃附近出现 Gibbs 现象。应按实际正则性选择误差预期。
微分会给系数乘上 或更高次幂。解析函数的几何优势仍在,但导数误差常带次数的代数因子。到舍入误差平台后,继续增加次数可能因微分条件性而使误差反弹。
推论与应用
实际实验应同时看系数尾部、独立检验点误差和离散求解稳定性。系数几何下降支持解析分辨率判断;平台提示有限精度或噪声;代数尾部提示有限光滑度或边界失配。它们是诊断信号,应与方程及边界条件一同解释。
参考资料