补充资料：数论网格求积分法

    　　高维数值积分数论方法研究开始于20世纪50年代末，其理论基础是数论中的一致分布论。命U_s表示 s维单位立方体。假定是U_s上定义的函数,并假定存在且其绝对值以C为界。命 是U_s中具有偏差D(n)的点集。所谓数论方法就是用被积函数在p(k) (1≤k≤n)上值的算术平均作为U_s上定积分的近似值，而误差由下面的公式给出：
　　
　　J(??，p(k))就是由点集p(k)(1≤k≤n)定义的一个求积公式。因此寻求U_s上最佳求积公式的问题即等价于寻求U_s上最佳偏差的点集的问题。从计算方法的观点看，不仅要求点集p(k)(1≤k≤n)的偏差小，而且要求p(k)的形式简单，易于计算。
　　
　　① 科罗博夫-劳卡方法　命p表示素数,a=(α₁,α₂,...,α_s)表示整数向量，科罗博夫和E.劳卡证明了,对于任意p，皆存在a，使点集有偏差。也就是说用点集Q(k)(1≤k≤p)构造的求积公式有误差。对于p求出a的计算量为O(p²)次初等运算。因此当p较大时,算出a来很困难。
　　
　　② 分圆域方法　分圆域是一个次代数数域。利用 的独立单位组可得它的一个适合于
　　的单位列n_l(l＝1,2,...),其中表示n_l的共轭数。如果使则得点集
　　用这一点集构造的求积公式的误差为
　　
　　 式中ε为任意正数。算出n_l、h_jl(1≤j≤s-1)的计算量为O(logn_l)。因此算出n_l和没有困难，但缺点是误差略为偏大些。
　　
　　当2≤s≤18时，上述的p、a、n_l和h都已汇编成表，可供查阅。
　　
　　数论方法得到的求积公式的误差主阶均与维数无关，所以当s较大时，用数论方法近似计算U_s上的定积分比较合算。
　　
　　

参考书目
　　　华罗庚、王元著：《数论在近似分析中的应用》，科学出版社，北京，1978。
　　

说明：补充资料仅用于学习参考，请勿用于其它任何用途。