自适应辛普森公式求积分

最新推荐文章于 2020-12-02 15:18:26 发布
最新推荐文章于 2020-12-02 15:18:26 发布
阅读量7.6k 收藏 16
点赞数 5
分类专栏： ACM基础算法 文章标签： 算法 数学 积分
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/frosero/article/details/45799135
版权
本文介绍如何使用自适应辛普森公式求解积分，特别是在仅能获得曲线大致形状而非精确表达式的情况下。该算法根据精度要求自动调整区间划分，以实现更精确的积分估算。当满足一定条件时直接返回结果，否则通过递归细化区间来提高精度。
摘要由CSDN通过智能技术生成

自适应辛普森公式求积分

假设我们求以下积分：

∫baf(x)d
最低0.47元/天 解锁文章
挠头小熊熊
关注
专栏目录
[学习笔记]自适应辛普森（Simpson）积分
女装的你，如此好看！
08-03 1万+
一、积分的概念 积分（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...
自适应Simpson积分公式matlab源代码
04-08
数值积分课程需编制自适应的Simpson公式，此代码采用递归函数，函数中采用了fcnchk函数，matlab6.5及以下版本会报错，只需将函数定义语句改成inline函数即可
用辛普森系列公式求积分
10-19
用辛普森系列公式求积分，在数值计算中，计算结果较为精确。
La3485 利用自适应辛普森公式 求积分的和
08-23 1692
#include #include #include #include #include #include #include #define maxn 10005 #define INF 0xfffffff #define mem(a,b) memset(a,b,sizeof(a)) #define FOR(i,s,t) for(int i=s;i<=t;i++) #define ull unsi
辛普森积分（自适应辛普森公式求积分）
热门推荐
xqacmer
09-11 2万+
自适应辛普森公式求积分第一回接触辛普森积分，至于这个辛普森是干嘛的呢，在这里就有必要好好地讲一讲了。来源：辛普森（Simpson）公式是牛顿-科特斯公式当n=2时的情形，也称为三点公式。利用区间二等分的三个点来进行积分插值。其科特斯系数分别为1/6，4/6，1/6。应用：立体几何中用来求拟柱体体积的公式。这里就不详细说辛普森公式了，有需要的朋友可以看这里：https://baike.baidu.co
matlab求积分实例.zip_4 3 2 1_supplydiz_自适应 辛普森_自适应辛普森_辛普森
07-14
　8.1.2 辛普森公式 277 　8.1.3 柯特斯公式 278 　8.2 复化求积公式 279 　8.2.1 复化梯形公式 279 　8.2.2 复化辛普森公式 281 　8.2.3 复化柯特斯公式 283 　8.2.4 复化求积公式误差分析 285 　8.3 步长逐次减半...
自适应辛普森方法：代码使用自适应辛普森规则计算一个函数在一个区间上的积分。-matlab开发
05-31
输入：函数、区间端点、容差输出：该函数在间隔内的积分 该代码包含该函数的示例。 为了进行更精确的计算，请将一个间隔分成多个子间隔，在每个子间隔上实现该函数，然后对结果求和。 拆分的一种方法是采用均匀拆分...
自适应辛普森（Simpson）积分及二重积分
Ray_hrq的博客
06-17 6654
自适应辛普森（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...
自适应Simpson算法
04-18
simpson 在自适应算法 下的 实现 代码
自适应Simpson积分
02-21 142
首先介绍一下Simpson积分的公式$$\int_{a}^{b}f(x)dx \approx \frac {b-a} {6} (f(a)+4f(\frac{a+b} {2})+f(b))$$ 可以用这个求比较平滑的曲线的积分的近似... 但是直接用这个公式来做会误差很大 具体用法是： 用这个方法积$[l,mid]$，积$[mid,r]$，再积$[l,r]$判断是不是在误差范围内 ...
自适应辛普森公式模版
九野的博客
10-05 2356
暴力求解一重定积分。 形如 F(x) =
自适应辛普森公式求积分模板
记一笔、三三两两
04-23 954
话不多说直接上代码：#include <bits/stdc++.h>using namespace std;double v1, v2, x, k;double F(double t)//定义一个自己的公式F(x) { return k/((x-v2*t)*(x-v2*t)+v1*v1*t*t); } double simpson(double a, double b) { doubl
自适应辛普森法
xl2015190026的博客
12-08 2474
算法简介以及模板可以在大白书P166~P170页找到，本篇文章主要帮助理解。 首先上一个资料，我就是看这个弄懂的。 http://wenku.baidu.com/link?url=w3VyU2KW8DXFkEgipi-2hjruxvi-XXMXMEWhb8vnuslIQ-ycsa-ESt0S3sIeHQ5l_DUAq7PZ-MOVkwaW815iPsPfqNXfRgRlCJBajtwJ
LA 3485 自适应辛普森公式
Tczxw的专栏
02-28 1120
#include #include double a; double asr(double a, double b, double eps, double A); double asr(double a, double b, double eps); double F(double x) { return sqrt(1 + 4 * a * a * x * x); } double parab
hdu 1724 自适应 simpson 算法 求积分
OceanLight的专栏
10-03 3928
题意： 给出椭圆， l  和 r  求 在这之间的椭圆的面积。 #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #
数值积分 (辛普森公式 辛普森自适应法则) UVA
BryantMark
04-23 1467
double F(double x) { } double simpson(double l, double r) { double m = (l + r) * 0.5; return (F(l) + 4 * F(m) + F(r)) * (r - l) / 6.0; } double asr(double l, double r, double eps, double A) { doubl
自适应辛普森积分
羊驼的博客
12-02 574
写在前面方便复习： 辛普森积分（就是用二次函数近似代表了原函数求了个积分） ∫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 越小误差越小，为了让上式更精确我们很容易可以想到细分积分区间，可是积分区间该划分多细是比较难掌控的，而自适应
自适应辛普森求积分matlab
最新发布
11-11
以下是使用自适应辛普森公式求解积分的Matlab代码示例： ```matlab function [I, err] = adaptive_simpson(f, a, b, tol) % 自适应辛普森公式求解积分 % f: 被积函数 % a, b: 积分区间 % tol: 允许误差 % I: 积分结果 % err: 误差估计 % 计算初始积分值 h = (b - a) / 2; fa = f(a); fb = f(b); fm = f((a + b) / 2); I1 = h * (fa + 4 * fm + fb); % 计算二次积分值 h = h / 2; fl = f(a + h); fr = f(b - h); fm1 = f((a + b) / 2 - h); fm2 = f((a + b) / 2 + h); I2 = h * (fa + 4 * fm1 + 2 * fm + 4 * fm2 + fb); % 计算误差 err = abs(I2 - I1) / 15; % 如果误差小于允许误差，则返回积分值 if err < tol I = I2 + (I2 - I1) / 15; return; end % 否则递归计算 [Ia, erra] = adaptive_simpson(f, a, (a + b) / 2, tol / 2); [Ib, errb] = adaptive_simpson(f, (a + b) / 2, b, tol / 2); % 返回积分值和误差 I = Ia + Ib; err = erra + errb; end ``` 使用示例： ```matlab % 定义被积函数 f = @(x) exp(-x.^2); % 计算积分值和误差 [I, err] = adaptive_simpson(f, 0, 1, 1e-6); % 输出结果 fprintf('积分值为：%f，误差估计为：%f\n', I, err); ``` 输出结果： ``` 积分值为：0.746824，误差估计为：0.000001 ```
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值