复化积分法(复化梯形求积-复化Simpson公式-变步长求积法)MATLAB编程实验报告_第1页
复化积分法(复化梯形求积-复化Simpson公式-变步长求积法)MATLAB编程实验报告_第2页
复化积分法(复化梯形求积-复化Simpson公式-变步长求积法)MATLAB编程实验报告_第3页
复化积分法(复化梯形求积-复化Simpson公式-变步长求积法)MATLAB编程实验报告_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、精选优质文档-倾情为你奉上复化积分法（复化梯形求积，复化Simpson公式，变步长求积法）MATLAB编程实验报告一、 问题描述：编写函数实现复化积分法。二、 实验步骤（过程）：（1） 复化求积法（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=

2、subs(y,t,b); T=h/2*(q+p+2*f); T=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

3、); if y='sin(t)/t'&&a=0 q=1; end p=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

4、=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

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论