好文档就是一把金锄头!
欢迎来到金锄头文库![会员中心]
电子文档交易市场
安卓APP | ios版本
电子文档交易市场
安卓APP | ios版本

复化积分法(复化梯形求积,复化Simpson公式,变步长求积法)MATLAB编程实验报告.doc

5页
  • 卖家[上传人]:m****
  • 文档编号:410828195
  • 上传时间:2024-01-07
  • 文档格式:DOC
  • 文档大小:29KB
  • / 5 举报 版权申诉 马上下载
  • 文本预览
  • 下载提示
  • 常见问题
    • 复化积分法(复化梯形求积,复化Simpson公式,变步长求积法)MATLAB编程实验报告一、问题描述:编写函数实现复化积分法二实验步骤(过程):(一)复化求积法(1)复化梯形求积:用复化梯形求积公式求解J1聾dx0xfunction[f]二Tn(a,b,n,y)symst;h=(b-a)/n;f=0;fork=1:n+1x(k)=a+(k-1)*hz(k)二subs(y,t,x(k));endfori=2:nf=f+z(i);endq=subs(y,t,a);ify=='sin(t)/t'&&a==0q=1;endp二subs(y,t,b);T二h/2*(q+p+2*f);T=vpa(T,7)clc,clear;symst;a=0;b=1;y=sin(t)/t;n=8;Tn(a,b,n,y);(2)复化Simpson公式:用复化Simpson公式求解€2exdxifunction[f]二simpson(a,b,n,y)symst;h=(b-a)/n;f=0;l=0;fork=1:n+1x(k)=a+(k-1)*hw(k)=0.5*h+x(k)z(k)=subs(y,t,x(k));endfori=2:nf二f+z(i);endfori二上nl=l+w(i);endq二subs(y,t,a);ify=='sin(t)/t'&&a==0q=1;endp=subs(y,t,b);T=h/2*(q+p+2*f);T=vpa(T,7)clc,clear;symst;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)。

      点击阅读更多内容
      关于金锄头网 - 版权申诉 - 免责声明 - 诚邀英才 - 联系我们
      手机版 | 川公网安备 51140202000112号 | 经营许可证(蜀ICP备13022795号)
      ©2008-2016 by Sichuan Goldhoe Inc. All Rights Reserved.