下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上执行函数为mymulNewtonCotes.m1、使用方法:Step1:在MATLAB命令窗口输入被积函数。输入应为:。Step2:执行函数。输入形式为mymulNewtonCotes(ft,a,b,m,n);其中ft被积函数,此体重,已经在第一步赋值;a积分下限,本题中为0;b积分上限,本题中为1;m将区间a,b等分的子区间数量,本题可选为10;n采用的Newton-Cotes公式的阶数,必须满足n8,否则积分没法保证稳定性。当n=1时,即为复化梯形公式;n=2时,即为复化复化辛普森公式。 所以,分别输入mymulNewtonCotes(ft,0,1,10,1)和m
2、ymulNewtonCotes(ft,0,1,10,2)就可以得到两种方法的积分计算结果。2、计算结果而根据积分运算,可得:说明复化梯形和复化辛普森公式计算出的结果基本一致,与实际结果相符。3、程序代码function yy = mymulNewtonCotes(ft,a,b,m,n)% 复化Newton-Cotes数值积分公式,即在每个子区间上使用Newton-Cotes公式,然后求和,% 参考的输入形式为mymulNewtonCotes(ft,0,1,10,2)% 参数说明:% ft被积函数,此题中ft=(t)t.*exp(t2/2)% a积分下限% b积分上限% m将区间a,b等分的子区
3、间数量% n采用的Newton-Cotes公式的阶数,必须满足n8,否则积分没法保证稳定性% (1)n=1时为复化梯形公式% (2)n=2时为复化辛普森公式xx = linspace(a,b,m+1);for l = 1:m s(l) = myNewtonCotes(ft,xx(l),xx(l+1),n);endyy = sum(s);function y,Ck,Ak = myNewtonCotes(ft,a,b,n)% 牛顿-科特斯数值积分公式% Ck科特斯系数% Ak求积系数% y牛顿-科特斯数值积分结果xk = linspace(a,b,n+1);for j = 1:n+1 ff(j) = ft(xk(j);end% 计算科特斯系数for i=1:n+1 k=i-1; Ck(i)=(-1)(n-k)/factorial(k)/factorial(n-k)/n*quadl(t)intfun(t,n,k),0,n);end% 计算求积系数Ak=(b-a)*Ck;% 求和算积分y=Ak*ff;function f=intfun(t,n,k)%
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 美容、护肤、护发产品相关行业投资规划报告
- 加气柱相关项目投资计划书范本
- 空调过滤网相关项目投资计划书
- 3-5-Dimethylphenol-d6-生命科学试剂-MCE
- 2-3-Dipalmitoyl-sn-glycero-1-phosphocholine-99-D-DPPC-99-生命科学试剂-MCE
- 人音版 音乐五年级上册第1课 朝夕 渔舟唱晚 教案
- 专项素养巩固训练卷(六) 电学综合计算教学设计 2024-2025学年教科版物理九年级上册
- 统编海上日出课件
- 准备单元 倒立的小丑 教学设计 -2024-2025学年科学三年级上册大象版
- 灾情统计工作汇报
- 光电子学-光电探测器
- 中学日语社团方案
- 胃肠穿孔护理常规课件
- 中学生普法知识小案例
- 工程竣工资料目录
- 五年级举一反三-第1周-平均数
- 人教A版数学高中选择性必修一《3.1.1椭圆及其标准方程》教学设计
- 肺部感染并呼吸衰竭病人疑难病例讨论
- 四年级心理健康教育教案-第三课 胸可以比海宽 通用版
- 手动啤机安全操作规程
- 高中语文-司马穰苴列传教学课件设计
评论
0/150
提交评论