计算方法 7.数值积分计算方法

这一节我们来学习如何估算积分：

这个Ak的计算方法其实就是把1,x,x^2都带进去…… 

求积公式是怎么来的？：

 这个就直接带公式。

既视感满满。

这波图穷匕见，先用插值法把f(x)换了，再求积分。换成多项式了我还不要用公式好好求，我还要用刚刚推出来的估算公式算。

Ak的值很好确认，Ln(x)=lk(x)*f(xk)，f(xk)只是常数，积分就和lk(x)有关。

 

理论看起来很复杂，算起来还是蛮简单的。

这下来等距了，真的既视感满满：

这个地方注意是对t积分。

例子：

研究一下这个公式的稳定性：

余项之前不是讲了吗？

既然项一多就会不稳定，那我就分开来球：

把那么多点分成2点一组，3点一组来求。（注意1/2的不算中间的点） 

隆贝哥积分法：

看起来贼复杂，至今没有搞懂原理。只能这么记：T0,i就是把原来的区间分成2^i份之后用复化梯形公式（以直代曲）求出来的公式求出来的结果，然后就用公式一步一步推到Tm,0。等到Ti-Ti-1小于要求的精度了就行。

gauss积分法：

之前讲过，两项的插值式阶数极限是2*1+1，不过要特定的点。

 例子：

Legendre是1,x,x^2……函数族的正交化。

这样可以强行转化为-1到1的积分，用同样的系数和Xk（记得这两个怎么算出来的吗？和F(x)没有关系，只是通过代数精度算出来的。）