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

复化积分法（复化梯形求积，复化Simpson公式，变步长求积法）MATLAB编程实验报告 问题描述： 编写函数实现复化积分法。 实验步骤（过程）： 复化求积法 复化梯形求积：用复化梯形求积公式求解 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 /ta 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 0;b 1; y sin t /t; n 8; Tn a,b,n,y ; 复化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 /ta 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 ; 变步长求积法：以书本例4.5为例 function [f] TN a,b,y,R0 syms t; T []; f 0; q subs y,t,a ; if y sin t /ta 0 q 1; end p 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; end end tl 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公式 变步长求积法