牛顿-柯特斯数值积分公式及其MATLAB的实现

(2009-11-20 23:46:18)

至于Newton-Cotes数值积分的具体表达、推导、说明和附件,下面给只是给出Newton-Cotes数值积分的Matlab实现代码

复化Newton-Cotes数值积分公式

function y=mulNewtonCotes(fun,a,b,m,n)

% 复化Newton-Cotes数值积分公式,即在每个子区间上使用Newton-Cotes公式,然后求和 % 参数说明

% fun,积分函数的句柄,必须能够接受矢量输入

% a,积分下限

% b,积分上限

% m,将区间[a,b]等分的子区间数量

% n,采用的Newton-Cotes公式的阶数,必须满足n<8,否则积分没法保证稳定性 % (1)n=1,即复化梯形公式

% (2)n=2,即复化辛普森公式

% (3)n=4,即复化科特斯公式

%

% Example

% fun=@(x)sin(x).*cos(x)

% mulNewtonCotes(fun,0,2,10,4)

xk=linspace(a,b,m+1);

for i=1:m

s(i)=NewtonCotes(fun,xk(i),xk(i+1),n);

end

y=sum(s);

牛顿-科特斯数值积分公式

function [y,Ck,Ak]=NewtonCotes(fun,a,b,n)

% y=NewtonCotes(fun,a,b,n)

% 牛顿-科特斯数值积分公式

%

% 参数说明:

% fun,积分表达式,这里有两种选择

% (1)积分函数句柄,必须能够接受矢量输入,比如fun=@(x)sin(x).*cos(x)

% (2)x,y坐标的离散点,第一列为x,第二列为y,必须等距,且节点的个数小于9,比如:fun=[1:8;sin(1:8)]'

% 如果fun的表采用第二种方式,那么只需要输入第一个参数即可,否则还要输入a,b,n三

复化科特斯公式matlab_牛顿相关推荐

  1. 复化科特斯公式matlab_matlab实现复化NewtonCotes公式求积分的程序应用和代码

    命令窗口输入被积函数 2 1 2 0 t t e dt  . 输入应为:. Step2 :执行函数.输入形式为 mymulNewtonCotes(ft,a,b,m,n) : 其中 ft -被积函数, ...

  2. 复化辛普森公式n_复化辛普森公式及龙贝格方法求解积分

    一.实验目的及题目 1. 实验目的: (1) 学习用复化辛普森公式及龙贝格方法求解积分并掌握这种方法. (2)了解这些辛普森公式及龙贝格方法的概念,参考课本写出用复化辛普森算法以及龙贝格方法计算目标题 ...

  3. 【数值分析实验MATLAB】数值积分:梯形公式、辛普森公式、复化梯形公式、复化辛普森公式、龙贝格算法、自适应求积方法

    (1)Matlab代码: f1: if x==0y=1; else%y=exp(-x^2);y=(sin(x))/x; end 梯形公式: a=0; b=1; tx=(b-a)/2*(f1(a)+f1 ...

  4. 自适应复化辛普森公式求积算法(C语言实现)

    自适应复化辛普森公式求积算法(C语言实现) 利用复化辛普森公式求积分自适应步骤 基于C语言实现的代码 利用复化辛普森公式求积分自适应步骤 h为步长,a为积分下限,b为积分上限,f为积分函数,n为划分的 ...

  5. 【计算方法数值分析】复化梯形公式、复化辛普森公式和龙贝格数值积分

    [计算方法数值分析]复化梯形公式.复化辛普森公式和龙贝格数值积分 1. 复化梯形公式 %复化梯形公式 function t=agui_trapz(fname,d2fname,a,b,e) %fname ...

  6. MATLAB 复化梯形公式、复化Simpson公式

    f.m文件: function f = f(x) f = (x^2)*sin(x); end 复化梯形公式 Tn.m文件: function Tn=Tn(n) % n代表区间数a = -2; % 区间 ...

  7. C++实现复化辛普森公式求积分算法

    1. 算法原理简介 步1 将积分区间 [a,b] 分成 n 等分,分点为 xk=a+kh(k=0,1,⋯,n),其中 h=(b-a)/n. 步2 记区间 [xk,x(k+1)] 的中点为 x(k+1/ ...

  8. 数值分析复化梯形公式matlab,数值分析复化梯形公式,复化Simpson公式MATLAB程序

    <数值分析复化梯形公式,复化Simpson公式MATLAB程序>由会员分享,可在线阅读,更多相关<数值分析复化梯形公式,复化Simpson公式MATLAB程序(1页珍藏版)>请 ...

  9. Python实现梯形公式 、辛普森公式、复化梯形公式、复化辛普森公式

    数值分析:梯形公式 .辛普森公式.复化梯形公式.复化辛普森公式 Python实现梯形公式 .辛普森公式.复化梯形公式.复化辛普森公式 数值求积公式概念 梯形公式与辛普森公式 梯形公式与辛普森公式的余项 ...

  10. matlab中复化辛普森公式函数,复化梯形公式,辛普森公式的matlab程序

    复化梯形公式与辛普森公式的matlab程序 [程序代码] cclc; disp('1.复化梯形公式求解'); disp('2.simpson公式求解'); disp('请进行选择:'); c=inpu ...

最新文章

  1. Java项目:精品酒店管理系统(java+SSM+mysql+maven+tomcat)
  2. 国外发明的10大仿生机械,这才是真正的黑科技!
  3. CGContextRef学习笔记
  4. LINQ体验(6)——LINQ to SQL语句之Join和Order By
  5. 前端学习(2564):如何触发组件的更新
  6. CLR via C# 边读边想 03 - 本地程序集和强命名程序集
  7. 公钥 私钥_公钥 私钥 签名 验签 说的啥?
  8. 指派问题:匈牙利算法
  9. 三年内人人有FIL,FIL 世界零撸板块引发全球流量狂潮!
  10. Neo4j Desktop(Neo4j 桌面)安装插件
  11. 如果我们现在要做一个库存盘点的系统,你怎么设计表结构,设计哪些表
  12. Apache HBase JavaAPI 查询数据
  13. asp.net word操作
  14. Agile Office 365 免积分下载
  15. 公交充值之——北京一卡通过期延期操作
  16. 偏振器件传输矩阵matlab编程,偏振器件的琼斯矩阵
  17. 离散数学/组合数学:利用容斥原理,求满射函数的个数,注意它容斥原理求得的是”非满射函数的个数”,以在A到B形成的函数的总个数为“全集“求该全集中”非满射函数的个数”的补集才是目标个数。
  18. 山东青岛:青岛海科展让青岛这座“会展之都”更有“看头”
  19. 淮阴工学院计算机专业历年录取分数线,淮阴工学院录取分数线2021是多少分(附历年录取分数线)...
  20. python怎么爬取豆瓣首页_Python3 爬虫(二) -- 爬取豆瓣首页图片

热门文章

  1. 微信 h5(公众号)简单登录
  2. Golang 流媒体音视频网络传输开源项目-LAL
  3. 软考—程序设计语言—有限自动机、正规式、编译程序基础、语法树、解释程序基础
  4. 语音芯片c语言程序,语音芯片4004C语言.doc
  5. linux crontab每小时运行一次
  6. 我是一个人喜欢写博客的人
  7. 项目管理工具ShowDoc的部署
  8. excel转置怎么操作_EXCEL转置的方法介绍,这种函数80%的人没用过,教你如何转置996...
  9. CSS - 选择器优先级介绍
  10. 微信也能开门?微信开门解决方案来了