MATLAB辛普森法则

梯形规则和辛普森规则是牛顿-科特规则的特例,它们使用更高阶的函数进行数值积分。

让抛物线代表图形的曲线。

Y =αX2 +βX+γ………… ..equation 1

在此方法下,间隔-h≤x≤h的面积为

曲线穿过三个点(-h,y 0 ),(0,y 1 )和(h,y 2 )。然后,通过等式,我们得到:

现在,我们可以评估系数α,β,γ并以h,y 0 ,y 1和y 2表示等式2。

通过将等式3的(b)代入(a)和(c)并重新排列,我们得到

αX2-βh= Y 0 -Y 1 … ..equation4

αX2 +βH= Y 2 -Y 1 … ..equation5

将方程式4与方程式5相加得出

2αh2 = Y 0 -2y 1 + Y 2 …… equation6

通过代入等式2,我们得到

现在,我们可以将方程式8应用于任意曲线y = f(x)的连续段,其间隔为a≤x≤b,如图所示。

我们观察到抛物线可以通过曲线的两端和中点近似于宽度为2h的每个线段。因此,线段AB下的面积为

同样,BC段下的面积为

等等。将每个细分下的区域相加后,我们得到

由于每个线段的宽度为2h,因此要应用辛普森数值积分法则,细分的数量n必须为偶数。此限制不适用于数值积分的梯形规则。

方程11的值可从

0

相关文章:MATLAB-循环类型 MATLAB-循环类型 在某 […]...

讨论MATLAB 讨论MATLAB MATLAB是MathWorks开& […]...

MATLAB简介 什么是MATLAB? MATLAB是用于高性能Math 计算,可视化和编程环境的软件包。它提供了具有数百种内置 […]...

MATLAB中的工作区,变量和函数 MATLAB中的工作区,变量和函数 工作空间 工作区包含我们在MATLAB中工作时创建的所有变量。 每当我们为 […]...

MATLAB while循环 MATLAB while循环 while循环重复执行语句,而指定的语句为true。 句法:...

MATLAB 2D Bar() MATLAB Bar() 条形图是其中每个点都由垂直条或水平条表示的图。 句法...

MATLAB 2D Silhouette() MATLAB轮廓线() 等高线图显示了矩阵Z的等值线。 句法...

MATLAB特征值和特征向量 特征值和特征向量 方阵A的特征值和特征向量为标量λ和满足条件的非零向量v Av =λv 在这个等式中, A是一 […]...

matlab中复华simpson,MATLAB Simpson的规则相关推荐

  1. 数值分析龙贝格matlab,龙贝格matlab程序

    k>=15 [龙贝格求积算法 Matlab 主程序] function[t]=rbg(f,a,b,c) t=zeros(15,4); %定义龙贝格积分函数,f 为待积函数,a 与 b 为积 分上 ...

  2. 数值积分与数值微分MATLAB,MATLAB程序设计教程(8)——MATLAB数值积分与微分

    MATLAB程序设计教程(8)--MATLAB数值积分与微分 第8章MATLAB数值积分与微分 8.1  数值积分 8.2  数值微分 8.1数值积分 8.1.1  数值积分基本原理 求解定积分的数值 ...

  3. matlab不少项一阶差分,MATLAB在有限差分法中的应用

    第!"卷第!期 !(("年)月 桂林工学院学报 *+,-'./+01,2/2'2'3424,45+04567'+/+18#$%&!"'$&!.9:& ...

  4. matlab lscov c语言,Matlab 函数大全 - 文章教程

    一.一般函数命令(General) 1.常用信息(General information) help 联机帮助命令,在MATLAB命令窗口显示帮助主题(在help后加函数或命令的特殊字符,可得到具体命 ...

  5. matlab 符号运算 数值,MATLAB中的微积分运算(数值符号)

    显然这个函数是单词differential(微分)的简写,用于计算微分.实际上准确来说计算的是差商. 如果输入一个长度为n的一维向量,则该函数将会返回长度为n-1的向量,向量的值是原向量相邻元素的差, ...

  6. 用MATLAB编程正弦稳态相量图,matlab课程设计--利用MATLAB对线性电路正弦稳态特性分析...

    matlab课程设计--利用MATLAB对线性电路正弦稳态特性分析 课程设计任务书 学生姓名: 专业班级: 指导教师: 刘 新 华 工作单位:信息工程学院 题 目: 利用MATLAB对线性电路正弦稳态 ...

  7. matlab设计译码器,基于MATLAB的循环码编译码器设计与仿真.doc

    扳昂旨螺冈唉陨裤外狸尿恨铸伸隧刽搅必勒诚天腑皖漂豌鲁靳碑缆键兽峙棘陶宽槐撒层僧袁廖颤渐魄货鼎躬薛扬衍逮西兰迫依煤鲁虐渠惫平合啥昭并屿己笆坍痞庐披吏去凄嘛兄察突徊溅今箩直藩潦咙锨谓崇若制匹扮复淌颐糖嗅你 ...

  8. matlab引入外部数据,matlab引入外部数据

    (1)适用与 matlab 等非编译语言. (2)适用与 C 和 Fortran 等编译语言. B. 如何在外部程序中告诉 ANSYS 改变单元属性和载荷变化: 在从 ANSYS 中导入数据后.... ...

  9. matlab com组件 发布,matlab生成com组件

    (中国科学院声学研究所北京100084) 摘要:介绍了采用MaflabCOMBuilder制作COM组件的方法,通过支持COM方式的编译工具调用生成的COM组件,能简 单方便地实现Matlab利其他. ...

  10. 加权voronoi图 matlab,加权voronoi图matlab

    广义Voronoi图的快速生成算法_电力/水利_工程科技_专业资料.27 卷第 ... Voronoi图理论与应用新成... 3页 免费 Voronoi图理论与应用新成... 3页 免费 Vorono ...

最新文章

  1. delphi获取闲置时间
  2. Android Service介绍
  3. 使用Spring JUnit规则进行参数化集成测试
  4. 【原始图像】对数及对比度扩展变换
  5. android开发 解析 b5,张绍文android开发高手课读书笔记4-启动优化篇
  6. Octotree插件
  7. amaze ui教程
  8. 带通滤波器电路图设计(转)
  9. C3之text属性的补充
  10. 劳埃镜 matlab,光学习题及答案
  11. MTK FAQ最新资料,MTK_on_line_FAQ_SW_ALPS_Framework+-+Battery
  12. 优秀的流程图是如何制作的?简单的教程讲解
  13. Safe Browsing API
  14. 王垠:编程宗派之我见
  15. Python爬虫必备知识点--整理核心基础
  16. 并行传输和串行传输的区别
  17. 有中国电信手机一定要看。CTWAP和CTNET是什么意思?有什么区别?
  18. 什么是代码签名?代码签名的好处
  19. TCP套接口的FIN_WAIT_2状态超时
  20. SQL中的DML语句

热门文章

  1. 彻底理解样本方差为何除以n-1
  2. 相控阵天线有源驻波测试
  3. [JZOJ4940]前鬼后鬼的守护/[JZOJ4623]搬运干草捆
  4. linux系统怎么安装office软件,如何在Linux上安装Microsoft Office | MOS86
  5. realme 制作一键root
  6. 新高考(3+1+2)模式下高中考试考场自动安排(Excel+VBA)
  7. npm ERR! Command failed: git clone
  8. SPIR-V 研究:编译器基本原理(三) - Chomsky文法分类
  9. java 判断日期格式_java判断字符是否是日期格式
  10. Three things can't discuss with with people: Religion, Politics and The Great Pumpkin. Why?为什么不能谈论?