首页/文章/ 详情

数值积分|高斯积分

7月前浏览4900

在区间[a,b]上,函数的积分可以近似用梯形公式表示:

 

如图a所示。这样当然会造成很大的误差。如果在区间内部找两个点,且通过这两个点的直线与区间端点构成的梯形面积最大限度地接近精确值,即图b中A1+A2=A3,这就是高斯积分的思路。

两点高斯积分公式

 

其中C0,C1为权系数。x0,x1为区间[-1,1]内的未知点(x0不等于x1),常识告诉我们,要想知道这四个未知量,就要解一个四元方程组。假设有四个函数 ,根据上面的公式可以得到联立解得

 
 
 

于是,两点高斯积分公式为

 

式(1)很好理解,就是一个矩形面积嘛。对于式(2),取一般的一次函数    ,如图c式(2)右边的积分就是梯形的面积,就是所熟知的上底加下底乘高除以二。还可以用梯形中位线表示

 

上式的意义是:一次函数的高斯积分需要一个高斯积分点即x=0的位置,确定的权重是2,积分点的函数值是f(0)。

对于式(3),取一般的二次函数    ,可以验证:

 
 
 

上式的意义是:二次函数的高斯积分需要两个高斯积分点    和    ,权重各为1,就可以计算积分了。

换句话说,对于线性函数,只要你告诉我f(0),我无需知道函数表达式,就可以算出你的积分了!同样,对于二次函数,只要你告诉我这俩函数值,我不需要知道函数的表达式,只要把这俩函数值    和    各乘以权重(都为1)相加即可算出积分值了.再来看三次函数    ,可以验证:

 
 

由此得到的规律是:四次,五次曲线有三个高斯积分点,六次曲线和七次曲线则需要四个高斯积分点,规律也是一样的。也就是说,n个高斯积分点可以计算2n-1次及以下的函数积分。

高斯积分点是强制使这种数值积分结果与前2n-1阶多项式的积分相等解出来的。比如你打算使用n个点,你还有n个未知权重。你就要使这种数值积分的结果等于对应的从0到2n-1的所有多项式项在区间内的积分结果。这样你就有一个2n阶的非线性方程组,解了它,就能获得积分点和权重值。


来源:数值分析与有限元编程
非线性
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-04-02
最近编辑:7月前
太白金星
本科 慢慢来
获赞 5粉丝 10文章 321课程 0
点赞
收藏
作者推荐

数值积分|泰勒(Taylor)公式求积分

泰勒(Taylor)公式大致可以叙述为:函数在一个点的邻域内的值可以用函数在该点的值及各阶导数值组成的无穷级数表示出来。ƒ(x)在x=a处的泰勒展开式为:注意,等号右边是无穷多项。特别地,当a=0时,又叫麦克劳林(Maclaurin)展开式ƒ(x)=e^x在x=0处分别展开几项展开多项式的函数图像与ƒ(x)=e^x对比ƒ(x)=cosx在x=0处展开多项式的函数图像与ƒ(x)=cosx对比可以看到,展开多项式项数越多,得到的图像和原函数越接近。实际应用中不可能无穷多项展开,故将其写成有限项Pn(x)和余项Rn(x)之和即ƒ(x)=Pn(x)+Rn(x)。例如c介于0和x之间。对于一个正常数M,有例如M=1时[算例]1.求积分要求误差小于0.001展开得x=1代入如果要求误差小于10^-6, 则保留前五项泰勒公式,应用于数学、物理领域,是一个用函数在某点的信息描述其附近取值的公式。如果函数足够平滑的话,在已知函数在某一点的各阶导数值的情况之下,泰勒公式可以用这些导数值做系数构建一个多项式来近似函数在这一点的邻域中的值。泰勒公式还给出了这个多项式和实际的函数值之间的偏差。来源:数值分析与有限元编程

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈