复化积分法(复化梯形求积-复化Simpson公式-变步长求积法)MATLAB编程实验报告.doc
复化积分法(复化梯形求积,复化Simpson公式,变步长求积法)MATLAB编程实验报告一、 问题描述编写函数实现复化积分法。二、 实验步骤(过程)(1) 复化求积法(1) 复化梯形求积用复化梯形求积公式求解function fTna,b,n,ysyms t;hb-a/n;f0;for k1n1 xkak-1*h zksubsy,t,xk;end for i2n ffzi; end qsubsy,t,a; if ysint/ta0 q1; end psubsy,t,b; Th/2*qp2*f; TvpaT,7clc,clear;syms t;a0;b1;ysint/t;n8;Tna,b,n,y;(2) 复化Simpson公式用复化Simpson公式求解function fsimpsona,b,n,ysyms t;hb-a/n;f0;l0;for k1n1 xkak-1*h wk0.5*hxk zksubsy,t,xk;end for i2n ffzi; end for i1n llwi; end qsubsy,t,a; if ysint/ta0 q1; end psubsy,t,b; Th/2*qp2*f;TvpaT,7clc,clear;syms t;a1;b2;yexp1/t;n5;simpsona,b,n,y;(3) 变步长求积法以书本例4.5为例function fTNa,b,y,R0syms t;T;f0;qsubsy,t,a;if ysint/ta0 q1;endpsubsy,t,b;T1b-a/2*qp;i2;ni-1;hb-a/n;z1ah/2;z2subsy,t,z1;T2T1/2h/2*z2;while Ti-Ti-1/3R0 ii1 ni-1; n2n-1 hb-a/n; f0; for k1n xkah*k-1; wkxkh/2; zksubsy,t,wk; ffzk; end TiTi-1/2h/2*f if Ti-Ti-1/3R0 break; endendtlTiclc,clear;at long;syms t;ysint/t;a0;b1;R00.5*10-3;TNa,b,y,R0;结论如下三、 结论复化梯形求积复化Simpson公式变步长求积法