-
2019-12-05 08:48:27
为什么要引入数值积分
由于Gauss型求积公式属于数值积分的内容,学东西总要知道它的来龙去脉,下面我简单介绍一下为什么要引入数值积分
给定函数 f ( x ) ∈ C [ a , b ] f(x)\in C[a,b] f(x)∈C[a,b],考虑积分
I ( f ) = ∫ a b f ( x ) d x I(f)=\int_{a}^{b} f(x) dx I(f)=∫abf(x)dx
的计算问题,从数学分析中知道,当已知 f ( x ) f(x) f(x)的原函数为 F ( x ) F(x) F(x)时,由牛顿-莱布尼兹公式,有
I ( f ) = ∫ a b f ( x ) d x = F ( b ) − F ( a ) I(f)=\int_{a}^{b}f(x)dx=F(b)-F(a) I(f)=∫abf(x)dx=F(b)−F(a)
然而,在实际计算中,被积函数的 f ( x ) f(x) f(x)的原函数经常无法用初等函数表示,但过于复杂。还有时, f ( x ) f(x) f(x)只在一些离散点上给出。在这样的情况下,就有必要借助数值方法来求 I ( f ) I(f) I(f)的近似值。Gauss型求积公式的定义
从正交多项式理论可知,在区间[a,b]上,对给定的权函数 ρ ( x ) \rho(x) ρ(x),存在正交多项式系 ω k ( x ) k = 0 ∞ {\omega_{k}(x)}_{k=0}^{\infty} ωk(x)k=0∞并且可以将其构造出来。进一步,已经证明了 ω k ( x ) \omega_{k}(x) ωk(x)在[a,b]上恰有k个相异的根,取 ω n + 1 ( x ) \omega_{n+1}(x) ωn+1(x)的n+1个根为求积结点,构造插值型求积公式,将得到具有2n+1次代数精度的求积公式,称如此构造出来的求积公式为Gauss型求积公式即:
如果求积结点 x 0 , x 1 , . . . . , x n x_{0},x_{1},....,x_{n} x0,x1,....,xn,使插值型求积公式 ∫ − 1 1 f ( x ) d x ≈ ∑ k = 0 n A k f ( x k ) \int_{-1}^{1}f(x)dx\approx\sum_{k=0}^{n}A_kf(x_k) ∫−11f(x)dx≈∑k=0nAkf(xk)的代数精度为2n+1,则称该求积公式为Gauss型求积公式,称这些求积结点为Gauss点。
特别地,Gauss求积公式当[a,b]=[-1,1]时,取 ρ ( x ) \rho(x) ρ(x)=1,Gauss型求积公式称为Gauss求积公式,此时的求积结点
x k ( k = 1 , 2 , 3.... ) x_k(k=1,2,3....) xk(k=1,2,3....)为n次Legendre多项式的根。
且 A k = ∫ − 1 1 ω ( x ) ( x − x k ) ω ′ ( x k ) d x A_k=\int_{-1}^{1} \frac{\omega(x)}{(x-x_k)\omega'(x_k)}dx Ak=∫−11(x−xk)ω′(xk)ω(x)dx
而Legendre多项式为:
打代码太南了,我手写偷下懒。Gauss公式的应用
Gauss公式主要应用为二点Gauss公式和三点Gauss公式
二点Gauss公式:当n=2时,有 x 1 = − 1 3 , x 2 = 1 3 , A 1 = A 2 = 1 x_1=-\sqrt\frac{1}{3},x_2=\sqrt\frac{1}{3},A_1=A_2=1 x1=−31,x2=31,A1=A2=1,求积公式为 G 2 ( f ) = f ( − 1 3 ) + f ( 1 3 ) G_2(f)=f(-\sqrt\frac{1}{3})+f(\sqrt\frac{1}{3}) G2(f)=f(−31)+f(31)
三点Gauss求积公式:当n=3时,有 x 1 = − 3 5 , x 2 = 0 , x 3 = 3 5 , A 1 = A 3 = 5 9 , A 2 = 8 9 x_1=-\sqrt\frac{3}{5},x_2=0,x_3=\sqrt\frac{3}{5},A_1=A_3=\frac{5}{9},A_2=\frac{8}{9} x1=−53,x2=0,x3=53,A1=A3=95,A2=98,求积公式为 G 3 ( f ) = 5 9 f ( − 3 3 ) + 8 9 f ( 0 ) + 5 9 f ( 3 5 ) G_3(f)=\frac{5}{9}f(-\sqrt\frac{3}{3})+\frac{8}{9}f(0)+\frac{5}{9}f(\sqrt\frac{3}{5}) G3(f)=95f(−33)+98f(0)+95f(53)
为了计算[a,b]上的定积分 ∫ a b f ( x ) d x \int_{a}^{b} f(x) dx ∫abf(x)dx可以通过自变量的变换 x = 1 2 ( a + b + t ∗ ( b − a ) ) x=\frac{1}{2}(a+b+t*(b-a)) x=21(a+b+t∗(b−a))将其转换为[-1,1]上的积分,然后使用Gauss积分公式计算其近似值。
例:
用两点Gauss公式计算 ∫ 0 1 s i n x x d x \int_{0}^{1}\frac{sinx}{x}dx ∫01xsinxdx
解:变换令x=0.5(t+1)
∫ 0 1 s i n x x d x = ∫ − 1 1 s i n 0.5 ( t + 1 ) t + 1 d x \int_{0}^{1}\frac{sinx}{x}dx=\int_{-1}^{1}\frac{sin0.5(t+1)}{t+1}dx ∫01xsinxdx=∫−11t+1sin0.5(t+1)dx
取 t 0 = − 1 3 , t 1 = 1 3 t_0=\frac{-1}{\sqrt 3},t_1=\frac{1}{\sqrt 3} t0=3−1,t1=31
∫ 0 1 s i n x x d x = 0.5 [ s i n 0.5 ( t 0 + 1 ) 0.5 ( t 0 + 1 + s i n 0.5 ( t 1 + 1 ) 0.5 ( t 1 + 1 ] \int_{0}^{1}\frac{sinx}{x}dx=0.5[\frac{sin0.5(t_0+1)}{0.5(t_0+1}+\frac{sin0.5(t_1+1)}{0.5(t_1+1}] ∫01xsinxdx=0.5[0.5(t0+1sin0.5(t0+1)+0.5(t1+1sin0.5(t1+1)]Matlab程序
function [w,p] = Gauss_point_1D(n,a,b) % Gauss quarature point on [-1,1] if n == 2 w = [1,1]; p = [-1/sqrt(3),1/sqrt(3)]; elseif n == 4 w = [0.3478548451,0.3478548451,0.6521451549,0.6521451549]; p = [0.8611363116,-0.8611363116,0.3399810436,-0.3399810436]; elseif n == 8 w = [0.1012285363,0.1012285363,0.2223810345,0.2223810345,0.3137066459,0.3137066459,0.3626837834,0.3626837834]; p = [0.9602898565,-0.9602898565,0.7966664774,-0.7966664774,0.5255324099,-0.5255324099,0.1834346425,-0.1834346425]; end % Gauss quarature point on [a,b] w = 0.5*(b-a)*w; p = 0.5*(b-a)*p+0.5*(b+a); function q = Gauss_quadrature(fun,a,b,n) % Gauss quadrature [w,p] = Gauss_point_1D(n,a,b); q = sum(w.*fun(p));
更多相关内容 -
高斯积分matlab程序,高斯求积公式matlab程序,matlab
2021-09-10 17:54:43计算在同一区域上的二重高斯-勒让德积分。 -
两个复化的Gauss-Legendre型求积公式及其误差分析 (2011年)
2021-05-10 16:15:04对两点和三点Gauss-Legendre公式进行复化,建立了两个新的数值积分公式,并分析了它们的积分误差和收敛阶。数值例子表明,我们的方法是高效的。 -
两点的高斯求积分公式自编写亲测可用
2019-12-29 14:27:38两点的高斯求积分公式自编写亲测可用-高斯积分公式-自编写-matlab书写-用于两点高斯求积直接运行.m文件即可 -
《数值分析》-- 高斯求积公式
2021-12-08 21:36:31适当选择节点,可使公式的精度最高达到2n+1,这就是本节所要介绍的高斯求积公式。...包括高斯型求积公式的一般理论,其中重点是高斯型求积公式和高斯点。主要讲解了高斯-勒让德求积公式,高斯-切比雪夫求积公式。文章目录
概述
- 问题
那么,在节点个数一定的情况下,是否可以在[a,b]上自由选择节点的位置,使求积公式的精度提得更高?
- 例题
一、高斯型求积公式的一般理论
- 一般理论
1.1 高斯型求积公式和高斯点
1.2 高斯点的特征
- 利用正交多项式构造高斯求积公式
习题
二、常用的高斯求积公式⭐
2.1 高斯-勒让德求积公式( Gauss-Legendre )
- 概念
勒让德多项式
- Legendre 多项式族:
低阶Legendre多项式
- 高斯-勒让德求积公式(G-L求积公式)⭐
- 总结
习题
- 一般有界区间[a, b]上的高斯-勒让德求积公式
(G-L求积公式)
目的转换区间到 [-1,1] :
习题
2.2 高斯-切比雪夫求积公式( Gauss-Chebyshe)
- 截断误差
总结
-
例题
用待定系数法构造高斯求积公式:
- 问题
为什么2点Gauss公式有应该有三次代数精度?
一般n+1个节点的求积公式的代数精度最高为2n+1次
- 详细过程
- 详细过程
-
例题
利用正交多项式构造高斯求积公式:
-
例题
- 公式回忆
- 复合梯形公式
三个点: T 2 T2 T2
五个点: T 4 T4 T4 - Simpson公式
- 三个点的高斯-勒让德求积公式
-
例题
- 公式
- 两点高斯-勒让德求积公式
-
例题
- 问题
-
数值计算方法(二)——复化求积公式
2019-10-20 15:16:53将区间n等分,每个小区间分别用梯形公式(两个节点)求积分,化简得到如上公式。 /** *@name Cotes:复化梯形公式 *@param1 below:区间下限 *@param2 upper:区间上限 *@param3 n:划分子区间的个数 **/ double...数值积分
直接利用某些节点上的函数值计算积分值,将积分求值问题归结为函数值的计算,可概述为如下公式
复化梯形公式
将区间n等分,每个小区间分别用梯形公式(两个节点)求积分,化简得到如上公式。/** *@name Cotes:复化梯形公式 *@param1 below:区间下限 *@param2 upper:区间上限 *@param3 n:划分子区间的个数 **/ double trapezoid(double below,double upper,int n) { double h=(upper-below)/n; double x=below; double S=0; for (int i=1;i<n;i++) { //计算内部各节点的函数值 1~(n-1) x+=h; S+=function(x); } return (h/2)*(2*S+function(below)+function(upper)); }
复化辛普生公式
同样将区间n等分,每个小区间分别用辛普生公式(两个节点加上中点)求积分,化简得到如上公式。/** *@name Cotes:复化辛普生法 *@param1 below:区间下限 *@param2 upper:区间上限 *@param3 n:划分子区间的个数 **/ double Simpos(double below,double upper,int n) { double h=(upper-below)/n; double S1=0; double S2=function(below+h/2); double x=below; for (int i=1;i<n;i++) { //计算各节点处的函数值 1~(n-1) x+=h; S1+=function(x); } x=below+h/2; for (int j=1;j<n;j++) { //计算各区间中点处的函数值 0.5~(n-0.5) 共n项 x+=h; S2+=function(x); } return (h/6)*(function(below)+function(upper)+2*S1+4*S2); }
复化柯特斯公式
同样将区间n等分,每个小区间分别用柯特斯公式(两个节点加上三个四等分点)求积分,化简得到如上公式。/** *@name Cotes:复化柯特斯公式 *@param1 below:区间下限 *@param2 upper:区间上限 *@param3 n:划分子区间的个数 **/ double Cotes(double below,double upper,int n) { double h=(upper-below)/n; double S1=0; double S2=0; double S3=0; double S4=0; double x=below; x=(below+h/4); for (int i=0;i<n;i++) { //计算各区间1/4处的值 1/4~(n-3/4)共n项 S1+=function(x); x+=h; } x=(below+h/2); for (int i=0;i<n;i++) { //计算各区间1/2处的值 1/2~(n-1/2)共n项 S2+=function(x); x+=h; } x=below+3*h/4; for (int i=0;i<n;i++) { //计算各区间3/4处的值 3/4~(n-1/4)共n项 S3+=function(x); x+=h; } x=below+h; for (int i=1;i<n;i++) { //计算各区间3/4处的值 1~(n-1)共n-1项 S4+=function(x); x+=h; } return (h/90)*(7*function(below)+32*S1+12*S2+32*S3+14*S4+7*function(upper)); }
-
26.高斯求积公式
2020-06-13 23:04:31 -
机械求积公式
2021-05-01 18:22:31数值积分与微分 代数精度 基本定义,只要满足 ∫abXm+1dx≠∑k=0nAkXkm+1 \int_a^bX^{m+1}dx\neq\sum_{k=0}^nA_kX_k^{m+1} ∫abXm+1dx=k=0∑nAkXkm+1 在m次以内都是相等的那么就称...插值型求积公式就是利 -
常见高斯型求积公式简介
2021-10-11 16:36:11一维区间上的求积公式 区间 [a,b][a,b][a,b] 上带权函数 ρ(x)\rho(x)ρ(x) 的插值型求积公式的一般形式为 I(f)=∫abρ(x)f(x)dx≈∑k=0nAkf(xk)(1) I(f) = \int^b_a \rho(x) f(x) dx \approx \sum_{k=0}^nA_k f(x_k... -
科学网—数值分析--低阶求积公式matlab程序 - 殷春武的博文
2021-04-21 04:39:22西北工业大学自动化学院 Email: yincwxa2013@mail.nwpu.edu.cn%% All rights reservedclearclcformat longsyms xf=input('输入被积函数f=')a=input('输入积分区间[a b]=')%%%%%%%%%%%%%%%%%%%%%梯形积分公式('梯形... -
数值计算实验报告---复合求积公式(梯形、抛物线、龙贝格)、导数求值(三点、四点、五点公式)
2021-04-23 19:35:15----------------------个人作业,如果有后辈的作业习题一致,可以参考学习,一起交流,请勿直接copy··复合抛物线公式:··龙贝格公式:四、实验内容------1实验题目1中所用到的三种算法的matlab实现代码具体如下... -
数值积分 (二)| 求积公式的代数精度 + 构造方法
2020-06-19 09:19:073. 求积公式的代数精度 求积公式的精度可以用余项R来表示,即: ∫abf(x)dx=∑k=0nAif(xi)+R \int_a^bf(x)dx=\sum_{k=0}^nA_if(x_i)+R ∫abf(x)dx=k=0∑nAif(xi)+R 一般地,R越小,表明求积公式的精度越高。... -
高斯型求积公式
2018-07-06 16:12:29下面讲高斯型求积公式的由来以及一般情况,比较枯燥,但若有志于弄懂高斯型求积公式的还是可以参考一下的。但对于使用者,做题型的要求来看,可快进到第二篇,或者直接看题。知道了高斯点,再求积求积系数,岂不太... -
Gauss-Legendre:使用 Gauss-Legendre 求积法对积分进行数值计算。-matlab开发
2021-05-31 04:55:22这个 zip 有 3 个功能: -polegende 生成legendre polinomyal。 -flege 执行积分变量更改以将任何区间 [a,b] 更改为 [-1,1],以便可以正确评估积分。 -gausslegende numericy 计算函数 [f] 的积分 -
Guass-Legendre(高斯-勒让德)求积方法 | Guass型求积公式 + Legendre多项式
2020-06-29 11:28:061. Guass型求积公式 定义1:在区间[a,b][a,b][a,b]内,如果由节点x0,x1,⋯ ,xnx_0,x_1,\cdots,x_nx0,x1,⋯,xn构造的插值型求积公式 ∫abf(x)dx≅∑k=0nAkf(xk) \int_a^b f(x)dx \cong \sum_{k=0}^nA_k f(x_k)... -
数值分析(10):数值积分之Gauss型求积公式
2021-05-23 15:37:38Gauss型求积公式2.1 Gauss型求积公式的定义2.2 Gauss点的性质2.3 构造Gauss型求积公式2.4 Gauss型求积公式的余项2.5 Gauss型求积公式的稳定性与收敛性2.6 Gauss-Legendre求积公式2.7 Gauss-Chebyshev求积公式 ... -
Gauss-Legendre求积公式-高斯型积分公式
2020-07-26 22:42:52网上找到的高斯勒让德求积公式一般都是最高4个插值点,正好这次做项目需要高一点的精度就写了这个程序。最高可以到8个插值点,并且插值点个数可以选择。 刚开始学习写程序,感觉这个" if "写的有点啰嗦但是自己也不... -
牛顿科特斯求积公式.ppt
2021-04-26 16:13:37n=4时的Cotes求积公式 按Newton-Cotes系数公式可以计算出 由此可得Cotes求积公式: a b 计算方法 余项公式为: 计算方法 §7.1 牛顿-科特斯求积公式 我们知道,若函数f(x)在区间[a,b]上连续且其原函数为F(x),则可用... -
高斯-勒让德求积公式及Matlab实现
2018-03-08 18:03:26微分实际上是求一函数的导数,而积分是已知一函数的导数,求这一函数。所以,微分与积分互为逆运算。实际上,积分还可以分为两部分。第一种,是单纯的积分,也就是已知导数求原函数,称为不定积分。相对而言,另一种... -
【Python学习笔记】6:用Gauss-Legendre求积公式近似求积分值
2017-05-04 21:50:52用高斯公式求积分的近似值,精确度是非常高的,一般用几个点就可以得到很不错的近似值。这里用了三点高斯积分和五点高斯积分。 至于这些点怎么找,实际上它们是勒让德多项式的零点,因为这个我没学,老师直接给出了... -
梯形求积公式 和 复合梯形求积公式 Matlab 实现
2017-11-04 18:22:52梯形求积公式 和 复合梯形求积公式 Matlab 实现梯形求积公式 仅使用区间两点x1,f(x1),x2,f(x2)x_1,f(x_1),x_2,f(x_2) 组成的梯形面积S代替∫x2x1f(x)dx \int_{x_1}^{x_2} f(x)dx 的近似方法 ∫x2x1f(x)dx≈S=... -
数值积分 (一)| 基本思想 + 一般求积公式
2020-06-19 09:16:49在Newton-Leibnitz公式应用受到限制时,可以用数值方法构造计算公式,即所谓的数值积分方法来解决定积分的计算问题。 表格函数,如某个微分方程的解不能用解析方法解出,而用数值方法给出的解是表格函数 被积函数f... -
Gauss型(Gaussian quadrature)求积公式和方法
2021-01-08 10:52:54The Gauss-Legendre quadrature rule is not typically used for integrable functions with endpoint singularities.(端点奇点) (1)基本概念 注:P0没有根(与x轴无交点),...,P2有2个根(与x轴有两个交点)... -
牛顿—柯特斯(Newton-Cotes)求积公式
2018-07-05 17:36:40看此博文需要首先知道插值型求积公式,见博文:插值型求积公式看插值型求积公式,又需要看拉格朗日插值法,见博文:拉格朗日插值法......————————————————————————————————————... -
matlab 高斯求积 简单
2013-12-31 16:13:02简单的高斯求积,用matlab编写,易懂易用 -
复合辛普森求积公式原理_辛普森法则:公式及其原理
2020-08-16 05:55:11复合辛普森求积公式原理Simpson's rule is a method for numerical integration. In other words, it's the numerical approximation of definite integrals. 辛普森法则是一种数值积分方法。 换句话说,它是定积分... -
Gauss型求积公式
2017-01-14 11:35:39高斯求积公式是变步长数值积分的一种,基本形式是计算[-1,1]上的定积分。理论证明对于 n个节点的上述求积公式,最高有 2n - 1 次的代数精度,高斯公式就是使得上述公式具有 2n - 1次代数精度的积分公式。 -
复化Cotes求积公式 | matlab
2020-07-10 15:26:59function c=cotes(a, b, TOL) % a为积分下限 % b为积分上限 % TOL 为误差容限 % cotes_f()为积分函数 n = 1; h = b - a; delt = 1;... % 如果两次计算的误差大于给定误差则进入循环 k = k + 1; h -
python实现数值分析之龙贝格求积公式
2018-10-24 16:34:23复合梯形公式的提出: 1.首先,什么是梯形公式: 梯形公式表明:f(x)在[a,b]两点之间的积分(面积...但高次的缺陷是当次数大于8次,求积公式就会不稳定。因此,我们用于数值积分的牛顿-科特斯公式通常是一次的梯形... -
数学----向量点积公式推导
2019-03-10 09:52:01向量的点积有两种形式的定义,代数定义和几何定义。 一 几何定义: 向量点积:a·b=|a||b|cosα 注意:该定义只对2维3维空间有效。 二 代数定义: 设二维空间内有两个向量和,定义它们的数量积(又叫内积、点...