稀疏网格

来自testwiki
跳转到导航 跳转到搜索

稀疏网格是表示、积分或插值高函数的数值计算技术。最初是由俄罗斯数学家Sergey A. Smolyak (Lazar Lyusternik的学生)基于稀疏张量积构造发展。高效实现此类网格的计算机算法后来由Michael GriebelChristoph Zenger 开发。

维度诅咒

表示多维函数的标准方式是采用张量或完全网格。故用于存储、运算的基函数或节点的数量与维数指数增加。即使以今天的计算能力,也不可能处理超过 4 或 5 维的函数。Template:来源请求

维度诅咒可以表示为使用Nl个格点进行l阶积分积分误差。若函数的正则性为r,即r次可微,维数为d,则

|El|=O(Nlrd)

Smolyak求积法则

Smolyak 发现了基于单变量求积规则Q(1)的计算上更为高效的多维函数积分方法。对d维函数f,Smolyak积分Q(d)一个函数的可以写成具有张量积的递归公式:

Ql(d)f=(i=1l(Qi(1)Qi1(1))Qli+1(d1))f

Q的下标是离散化的水平,我们不妨令一维i阶的积分要对O(2i)个点求值。[1]正则性为r的函数的误差估计是:

|El|=O(Nlr(logNl)(d1)(r+1))

延伸阅读

外部链接