计算方法与(5-2) .ppt

* 西南交通大学电磁所 5.6 高斯求积公式 在Newton—Cotes公式中，取xi为等距节点，使得C i(n)的计算容易，也使整个N─C公式简便，我们由此得到的逐次等分区间的变步长方法、 Romberg方法均建立在等距节点基础上，是非常方便的。 但 xi 限定为等距节点，却可能使其代数精度受到限制. 对 ，其几何意义是曲边梯形aABb的面积，若以梯形公式计算I ( f )≈(b-a)[f(a)+f(b)]/2 . 即以a, b为节点(n=1)，用直线AB近似曲线。也即以梯形aABb的面积近似替代曲边梯形面积。 A B B ? A ? a b x1 x2 由于假定是等距节点，因此n=1只能以a,b为插值节点，若不限定为等距节点，即可选择两个节点，如以x1,x2为节点，做直线A?B?近似曲线,也即以梯形面积aA?B?b近似曲边梯形面积，更准确一些，这表明：我们可选择合适的节点（不限定为等距节点）有可能进一步提高求积公式的代数精度。 我们已知N—C公式是由n+1个点构造一个小于等于n次的插值多项式Ln(x)逼近f (x)，所以插值型求积公式 ： 的代数精度不会低于n次，即对f (x)=1,x,x2,…,xn 准确成立。 （这一点在误差估计式中f（n+1）( ? ) = 0，因此误差为0非常清楚）亦即插值型求积公式的代数精度至少为n次，即代数精度最低为n次，那么最高能达到多少次代数精度呢？ 若只取n个点，N-C 公式为： 则最低次代数精度为n－1次，最高为多少次？若在上述公式中，限定xi为等距节点，则只有Ai可选择，为确定Ai(i=1,2,…,n)，按待定系数法，可设上述公式对f(x)=1,x,x2,…,xn-1成立?构成n个方程的线性方程式求解出A1, A2,…,An 。 假如xi也可选择，则在上述公式中，有xi , Ai (i=1,2,…,n)共2n个参数待定，按待定系数法可设上述公式对f (x)=1,x,x2, …,xn-1, xn, xn+1, …,x2n-1成立?构成2n个方程的非线性方程组可确定x1, x2,…, xn及A1, A2,…, An。 的代数精度不会超过2n－1次。 这表明插值型求积公式： 证（反证法）：假定最高代数精度m＞2n-1。则对f ( x )为2n次多项式时，上述求积公式准确成立： 下面证明插值型求积公式的代数精度最高不超过 2n－1 次。 与上式结果相矛盾，所以最高次代数精度m≤2n－1。 上述推导和证明表明：可适当选择节点和系数，使求积公式代数精度达到最高次。 例 10: 试确定x1, x2, A1, A2使其具有3次代数精度。 解：设对f ( x )=1, x, x2, x3准确成立，可得： 解方程组，可得： 这里n=2,两个点代数精度能够达到3次,即2n－1=3,为Gauss型求积公式，其几何意义，如右图：以梯形OA?B?1，近似曲边梯形OAB1，称x1,x2为Gauss点。 A B B ? A ? o x1 x2 1 5.7 数值微分 对于给定的(xi,yi) (i=0,1,2,…,n)，要求f (x)在节点xi处的导数值f (k)(xi)，通常的做法是建立插值多项式Ln(x)，以Ln(x)近似f (x)，而以Ln(x)的的导数近似 f (x)的导数： 称为插值型的微分公式。 上述误差估计式中，尽管有了一个估计式，但要确定也比较困难，注意到在结点xi处，ω(x)为零，因此，对等距节点xi = x0+ih, (i =0,1,2,…,n)，我们可求在节点xi处的导数值。 1. 两点公式： * 西南交通大学电磁所