辛普森求积公式 和 复合辛普森求积公式 Matlab 实现
辛普森求积公式 和 复合辛普森求积公式 Matlab 实现
辛普森求积公式
利用区间[a,b]的端点及中点计算积分
∫x2x1f(x)dx≈b−a6×(f(a)+4×f(b−a2)−+f(b))\int_{x_1}^{x_2} f(x)dx \approx \frac{b-a}{6} \times (f(a) + 4 \times f(\frac{b-a}{2}) -+f(b))
复合辛普森求积公式
将求积区间[a,b]分为n个区间,每个区间步长为h( h=b−an h = \frac{b-a}{n})在每个区间求辛普森积分
∫x2x1f(x)dx≈h6∫(f(a)+4×∑k=0n−1f(xk+12)+2×∑k=1n−1f(xk)+f(b))\int_{x_1}^{x_2} f(x)dx \approx \frac{h}{6} \int (f(a) + 4 \times\sum_{k = 0}^{n-1}f(x_\frac{k+1}{2}) + 2 \times\sum_{k = 1}^{n-1}f(x_k) + f(b))
其中 xk+12=xk+h2 x_\frac{k+1}{2} = x_k + \frac{h}{2}
体会
∑k=0n−1f(xk+12)\sum_{k = 0}^{n-1}f(x_\frac{k+1}{2}) 这项是复合辛普森求积公式比复合梯形求积公式多出的式子
这项是辛普森求积公式多余出计算量也正是 复合更精确的原因
说明
上述公式是我用mathjax写的,如有错误请联系我修正
敬请指正
概述省略了部分推导过程,请查阅详细推导资料
Matlab 实现代码
辛普森求积公式
将该函数存为m文件
function res = Simpson(f,a,b)format long;if b < ac = b;b = a;a = c;endres = (b-a)/6 * (f(a) + 4*f((a+b)/2) + f(b));
end
调用下面语句测试函数
f = inline('sin(x)','x')
Trapezium(f,0,pi/2)
复合辛普森求积公式
将该函数存为m文件
function res = ComSimpson(f,n,a,b)format long;if b < ac = b;b = a;a = c;endh = (b-a)/n;d = f(a);for i = a+h:h:b-h d = d + (2 * f(i));endfor i = a+h/2:h:b-h/2 d = d + (4 * f(i));endd = d + f(b);res = (d * h / 6);
调用下面语句测试函数
f = inline('sin(x)','x')
ComSimpson(f,4,0,pi/2)
辛普森求积公式 和 复合辛普森求积公式 Matlab 实现相关推荐
- 复合辛普森求积公式原理_辛普森法则:公式及其原理
复合辛普森求积公式原理 Simpson's rule is a method for numerical integration. In other words, it's the numerical ...
- 抛物线求积公式求积分算法matlab,基于Matlab的数值积分公式问题.doc
基于Matlab的数值积分公式问题 数值分析 学 号: 学 生 姓 名 :教 师 : 教师 2数值积分算法介绍............................................. ...
- 抛物线求积公式求积分算法matlab,数值计算实验报告---复合求积公式(梯形、抛物线、龙贝格)、导数求值(三点、四点、五点公式)...
----------------------个人作业,如果有后辈的作业习题一致,可以参考学习,一起交流,请勿直接copy ··复合抛物线公式: ··龙贝格公式: 四.实验内容 ------1 实验题目 ...
- 复合梯形公式C语言程序,复合梯形公式、复合辛普森公式 matlab(示例代码)
1. 用1阶至4阶Newton-Cotes公式计算积分 程序: function I = NewtonCotes(f,a,b,type) % syms t; t=findsym(sym(f)); I= ...
- matlab 复合辛普森公式,matlab中如何用复合辛普森公式求二重积分?
function q=DblSimpson(f,a,A,b,B,m,n) if(m==1 && n==1) %辛普森公式 q=((B-b)*(A-a)/9)*(su ...
- matlab数值与符号运算
matla数值与符号运算 1. 多项式计算 主函数 clc clear close %多项式 x^4-12*x^3+25*x+116 p=[1 -12 0 25 116]%多项式求值函数 polyva ...
- 数值积分-求积公式余项,牛顿-柯特斯公式,辛普森公式,复合梯形公式,复合辛普森公式
文章目录 1.求积公式余项 1.1 定义 1.2 Python实现求积公式余项 2.牛顿-柯特斯公式 2.1 定义 2.2 Python实现牛顿-柯特斯公式 3.复合梯形公式 3.1 定义 3.2 P ...
- 数值计算大作业:数值积分(梯形、辛普森与龙贝格方法在Matlab实现)
作为研究生的入门课,数值计算的大作业算是所有研究生开学的重要编程作业. 针对数值积分的编程,我把梯形.辛普森与龙贝格方法在MATLAB中编写的程序放在文章最后了,需要的同学自取. PS:附录中的程序并 ...
- matlab用辛普森公式求积分_数值积分常用方法
数值积分的基本思想 由积分中值定理可知,在积分区间 内存在一点 ,成立 式的几何意义即为:底为 而高为 的矩形的面积恰等于所求曲边梯形的面积 .因此,要想求出 式左端积分,我们只需要知道三个值: 即可 ...
最新文章
- pandas 排序 给excel_懂Excel轻松入门Python数据分析包pandas(二十五):循环序列分组...
- 一步一步做一个linux文件传输软件(一)
- VS 调用外部DLL时提示“无法嵌入互操作类型”
- php代码expl,php – 参数号无效:参数未定义Explination
- ZT:Java代码编写的30条建议
- fiddler启动文件_如何成为Fiddler大神,请看这里!
- 30个免费网页设计模板
- paip.提升性能--多核编程中的java .net php c++最佳实践 v2.0 cah
- 数据分析思维框架与分析方法总结
- LIF-Seg 激光雷达与相机图像融合的三维激光雷达语义分割
- springboot的web进阶知识(2)
- 前后端开发的心得体会_web前端学习心得体会
- js实现简易打点计时器
- 文件上传控件-如何上传文件-文件夹上传
- XlsReadWriteII 版本声明操作要点
- C++课程设计快递业务管理系统实验报告
- php卡片式,卡片式设计的优点和不适用性
- 洛谷P1796 汤姆斯的天堂梦【DP-数字三角形模型】【橙】
- 人工智能在各个领域里的应用场景
- 小勇机器人如何绑定_小勇机器人app