2022复化积分法复化梯形求积复化Simpson公式变步长求积法MATLAB编程实验报告
《2022复化积分法复化梯形求积复化Simpson公式变步长求积法MATLAB编程实验报告》由会员分享,可在线阅读,更多相关《2022复化积分法复化梯形求积复化Simpson公式变步长求积法MATLAB编程实验报告(4页珍藏版)》请在装配图网上搜索。
1、复化积分法(复化梯形求积,复化Simpson公式,变步长求积法)MATLAB编程实验报告一、 问题描述:编写函数实现复化积分法。二、 实验环节(过程):(一) 复化求积法(1) 复化梯形求积:用复化梯形求积公式求解function f=Tn(a,b,n,y)syms t;h=(b-a)/n;f=0;for k=1:n+1 x(k)=a+(k-1)*h z(k)=subs(y,t,x(k);end for i=2:n f=f+z(i); end q=subs(y,t,a); if y=sin(t)/t&a=0 q=1; end p=subs(y,t,b); T=h/2*(q+p+2*f); T=
2、vpa(T,7)clc,clear;syms t;a=0;b=1;y=sin(t)/t;n=8;Tn(a,b,n,y);(2) 复化Simpson公式:用复化Simpson公式求解function f=simpson(a,b,n,y)syms t;h=(b-a)/n;f=0;l=0;for k=1:n+1 x(k)=a+(k-1)*h w(k)=0.5*h+x(k) z(k)=subs(y,t,x(k);end for i=2:n f=f+z(i); end for i=1:n l=l+w(i); end q=subs(y,t,a); if y=sin(t)/t&a=0 q=1; end p=
3、subs(y,t,b); T=h/2*(q+p+2*f);T=vpa(T,7)clc,clear;syms t;a=1;b=2;y=exp(1/t);n=5;simpson(a,b,n,y);(3) 变步长求积法:以课本例4.5为例function f=TN(a,b,y,R0)syms t;T=;f=0;q=subs(y,t,a);if y=sin(t)/t&a=0 q=1;endp=subs(y,t,b);T(1)=(b-a)/2*(q+p);i=2;n=i-1;h=(b-a)/n;z1=a+h/2;z2=subs(y,t,z1);T(2)=T(1)/2+h/2*z2;while (T(i)-T(i-1)/3)R0 i=i+1 n=i-1; n=2(n-1) h=(b-a)/n; f=0; for k=1:n x(k)=a+h*(k-1); w(k)=x(k)+h/2; z(k)=subs(y,t,w(k); f=f+z(k); end T(i)=T(i-1)/2+h/2*f if (T(i)-T(i-1)/3)=R0 break; endendtl=T(i)clc,clear;format long;syms t;y=sin(t)/t;a=0;b=1;R0=0.5*10(-3);TN(a,b,y,R0);结论如下:三、 结论:复化梯形求积:复化Simpson公式变步长求积法
- 温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 七年级数学:平面直角坐标系中不规则图形面积的计算课件
- 地理环境的差异性课件——曲靖市麒麟区第一中学-郭金涛1
- 高考数学文科一轮复习(北京卷B)ppt课件102双曲线及其性质
- 剖宫产瘢痕妊娠课件
- 七年级《走进初中-享受学习》主题班会ppt课件
- 2020高端楼盘营销渠道拓客课件
- 七年级道德与法治下册第四单元走进法治天地第九课法律在我们身边第1框生活需要法律ppt课件2新人教版
- 七年级生物血管和心脏课件
- 高考数学一轮复习讲义-第一章-1.2-命题及其关系、充分条件与必要条件ppt课件
- 2021公司元旦跨年预售活动方案课件
- 高考数学一轮复习-第十八章-第2讲-极坐标与参数方程ppt课件-文
- 高考数学一轮复习-第二章第四节-幂函数与二次函数ppt课件-理-(广东专用)
- 北欧家具设计课件
- 高考数学一轮复习5.35特殊数列求和ppt课件理
- 高考数学浙江专用二轮培优ppt课件专题六第1讲排列组合二项式定理