《地球物理计算方法》第2章数值积分 2.pptx

《地球物理计算方法》;1、机械求积 2、Newton-cotes积分公式 3、复化求积方法 4、Romberg加速算法 5、Gauss积分公式 6、数值微分;问题的提出;A f (x );方法： 根据代数精度定义：求积具有2n-1 次代数精度，令 f (x) = 1, x, x2, …, x2n-1 代入机械求积公式等式准确成立，即建立 等式。如两点高斯公式;得到一个非线性方程组，;;;2. 高斯点的基本特性;设xk(k=1,2,..n)是求积公式的高斯点，作多项式w(x);定理：节点xk(k=1,2,..n)是高斯点的充分必要条件是，多项式w(x) 与一切=n-1的多项式pn-1(x)正交。 ?对任意的多项式pn-1(x)成立； ?最简单的多项式为1,x,x2…xn-1，所以根据下列积分等式，可 以求出各高斯点xk;3.勒让德多项式;13;高斯点 令勒让德多项式为零， 即为方程求根； 根据定义： 这些根为高斯点；;高斯积分公式的特点：;1、机械求积 2、Newton-cotes积分公式 3、Gauss积分公式 4、复化求积方法 5、Romberg加速算法 6、数值微分;由导数定义,差商近似导数,得到数值微分公式;差商比较示意图;f (a) ??.;2、中点方法的加速;G ?h / 2? ? f ??a? ? ?h /2?;G ?h? ? G ?h / 2? ? 3 ?h /2?