自适应辛普森

很多时候，我们会面临一些求积分的问题，无论是直接给你函数，让你想办法求解积分，还是对于一些计算几何问题，无法直接推导积分，我们都可以用这种方法来求一段区域的积分，积分的相关基础概念这里就不再赘述，今天主要就是说明他的大致原理，和他的用法。一个普通的思路就是，把要求的积分区间分解成一个个小的区间，然后累加在一起，但是这样做的缺点就是，如果为了保证精度，每个区间都会很小，所以这样就很容易会TLE。至于为啥要弄15， 这里有个论文给了解释，感兴趣可以了解，不感兴趣，存个板子直接用就好了。这里在介绍一个优化技巧。

简述： 对于普通的二次函数有f(x)=ax2+bx+cf(x) = ax ^ 2 + bx + cf(x)=ax2+bx+c，原函数F(x)=ax33+bx22+cxF(x) = \frac{a x^3}{3} + \frac{bx ^2}{2} + cxF(x)=3ax3​+2bx2​+cx。 有 ⎰lrf(x)=F(r)−F(l)=a3(r−l)3+b2(r−l)2+c(r−l)=(r−l)6(2a(l2+r2+lr)+3b(l+r)+c)=(r−l)6((al2+bl+c)+(ar2+br+c)+4(

一、积分的概念 积分（integral）的几何意义是函数的曲线上 xxx 的一段区间与 xxx 轴围成的曲边梯形的面积： xxx 的区间为 [a,b][a,b][a,b] ，那么上图阴影面积为： ∫baf(x)dx∫abf(x)dx\int_a^bf(x)\text{d}x 计算方法一：分割成无穷多个小区间。 ∫baf(x)dx=limn→∞∑i=1nb−anf(a+b−ani)∫a...

数值分析——自适应辛普森积分 自适应辛普森积分法是利用二次函数对被积函数进行插值的一种模拟积分法，常用于估计积分的数值，利用递归自动适应所需精度。 公式推导 基本思想公式： ∫abf(x)dx≈∫ab(Ax2+Bx+C)dx \int_{a}^{b} f(x) dx \approx \int_{a}^{b} (Ax^2+Bx+C)dx ∫ab​f(x)dx≈∫ab​(Ax2+Bx+C)dx 然后我们对这个二次函数积分进行求解： ∫ab(Ax2+Bx+C)dx=(b−a)6(Aa2+Ba+C+Ab2+Bb+C

何为自适应Simpson法 自适应辛普森算法(Adaptive Simpson′s rule)是一类近似算法(Approximation algorithm)，主要用于在信息计算时求较难反导的函数的定积分。其思想是利用二次函数曲线来不断**拟合(Overfitting)所求曲线（显然这传统的直接用矩形或直边梯形作为微元更精准），而所谓的Adapative（自适应）**则是用于优化时间复杂度的方法。 ----wikipedia 在计算几何里，显然可在特殊情

写在前面方便复习： 辛普森积分（就是用二次函数近似代表了原函数求了个积分） ∫abf(x)dx≈b−a6(f(a)+f(b)+4f(a+b2))\int_a^bf(x)dx \approx \frac {b-a} 6(f(a)+f(b)+4f(\frac{a+b} 2))∫ab​f(x)dx≈6b−a​(f(a)+f(b)+4f(2a+b​)) 自适应辛普森积分 很显然上式中的 b−ab-ab−a 越小误差越小，为了让上式更精确我们很容易可以想到细分积分区间，可是积分区间该划分多细是比较难掌控的，而自适应

辛普森积分是数值积分的一种，是中点公式和梯形公式的改进。 ——这是我在APIO2017时听某位大神讲课所了解的。 下面就讲讲辛普森积分到底是啥玩意儿。假定我们要求如下定积分 ∫baf(x)dx\int_a^bf(x)dx 略微懂一点微积分知识的都知道，对于一个黎曼可积的函数，我们要求其在某个闭区间上的定积分，要先求该函数的不定积分，即先求原函数。就是找到一个函数F(x)F(x)，使得F

一.啥子是自适应辛普森法呢？ 简而言之，就是一种用二次函数来逼近被积函数。 把求原来函数的积分换成求二次函数的积分的一种近似求积分的方法 二.自适应辛普森法的推导 什么？你想看推倒过程？ 小小年纪怎么净想着推倒？ 莫得推倒，这就是simpson公式（嘿嘿） 那这玩应和自适应辛普森法又有啥子关系呢？ 所谓自适应 就是我们可以以二分的方式，把函数划分成一个个小的区间进行求和，使得函数划...

自适应辛普森公式求积分第一回接触辛普森积分，至于这个辛普森是干嘛的呢，在这里就有必要好好地讲一讲了。来源：辛普森（Simpson）公式是牛顿-科特斯公式当n=2时的情形，也称为三点公式。利用区间二等分的三个点来进行积分插值。其科特斯系数分别为1/6，4/6，1/6。应用：立体几何中用来求拟柱体体积的公式。这里就不详细说辛普森公式了，有需要的朋友可以看这里：https://baike.baidu.co

先来看一道题： 给你三条直线(x=a,x=b,y=0)和一个函数 求该函数图像与这三条直线围成的面积。 【输入】 第一行n,a,b。 第二行n+1个整数A0~An。 【输出】 一行一个整数，即这四条先围成的面积，保留两位小数。 【样例输入1】 2 0 1 1 0 0 【样例输出1】 0.33 【样例输入2】 3 -1 0 1 0 0 0 【样例输出2】 -

一、概念： 辛普森积分是数值积分的一种，是中点公式和梯形公式的改进。 二、解析： 假定我们要求如下定积分： 略微懂一点微积分知识的都知道，对于一个黎曼可积的函数，我们要求其在某个闭区间上的定积分，要先求该函数的不定积分，即先求原函数。就是找到一个函数 F(x)，使得 F′(x)=f(x)，然后根据牛顿—莱布尼茨公式去搞。即： 而有时候原函数并不好求，比如要求原的函数长得......

各种数值积分方法总结（龙贝格积分、高斯积分等）常用的数值积分方法牛顿-科茨（Newton-Cotes）积分公式梯形法则 ，辛普森法则，布尔法则，牛顿-科茨积分公式，复合梯形积分复合辛普森积分等

有时候我们需要计算一个函数的定积分，粗略上可以使用估算的方法。如图所示，将原本的曲线粗略地看成一个梯形。这个方法叫。也叫做一阶其中所谓一阶，指的就是n=1的情况。最理想的情况就是把这个图像分割成 无数 个梯形，便可求出对应的定积分。但是在实际操作的情况下，梯形法则为了保证速度无法取极多的点，这样照成梯形法则误差较大。分割成无限个梯形其实就等效于因此我们将考虑更高阶的公式，本文将要介绍的便是二阶即将函数近似看成一条抛物线。显然一阶牛顿-柯特斯闭型积分公式需要在首尾取两个点方可得到f(x)的解析式。

试计算积分结果保留至小数点后 6 位。数据保证计算过程中分母不为 0 且积分能够收敛。一行，包含 6 个实数 a,b,c,d,L,R。一行，积分值，保留至小数点后 6 位。2.732937a,b,c,d∈[−10,10]，−100≤L<R≤100 且 R−L≥1。

数值微积分

120701101数值积分 复合梯形法辛普森1/3法CONTENTSCONTENTS1. 复合梯形法1.1 代码1.2 计算结果1.3 图像结果2. 辛普森1/3法2.1 代码2.2 计算结果2.3 图像结果3. 源码编辑器1. 复合梯形法1.1 代码# -*- coding: utf-8 -*- 1.2 计算结果运行代码即可得到结果，该区间内的数值积分为1.3 图像结果复合梯形2. 辛普森1/3...

自适应辛普森（Simpson）积分及二重积分 概述 辛普森积分法利用拉格朗日插值法来近似拟合原函数，然后通过对拟合函数的积分来简化近似。 近似公式： ∫abf(x)d⁡x≈b−a6(f(a)+4f(a+b2)+f(b)\int_{a}^{b}{f(x)}\operatorname d{x} \approx \frac{b-a}6(f(a)+4f(\frac{a+b}2)+f(b)∫ab​f(x...

一种计算定积分的数值算法。所谓数值算法就是给不出解析式时，可以通过数值模拟的方式得到近似值，在一定精度内是对的已知被积函数，但是原函数求不出来，又想求定积分，就可以上自适应辛普森法。原理大概是有个数学家发明的辛普森公式，可以用二次函数来逼近其他函数，那么我们在一个区间定积分时，可以拆到更小区间，如果在一个小区间内，辛普森公式的拟合误差小于一个阈值，我们就用拟合出来的二次函数的积分，作为被积函数在这个小区间内的积分。如果拟合精度不够，继续拆分小区间，向下递归。