matlab编程实例2.doc

数值计算方法上机题目2 计算定积分的近似值： 要求： 若用复化梯形公式和复化Simpson公式计算，要求误差限，分别利用他们的余项估计对每种算法做出步长的事前估计； 分别利用复化梯形公式和复化Simpson公式计算定积分； 将计算结果与精确解比较，并比较两种算法的计算量。 （1）采用复化梯形公式对上述函数用余项估计区间划分个数n，步长h，近似解Tn，以及与真实值误差w的源程序如下 function Fuhuatixing clc f=input (请输入积分函数f:); a=input (请输入积分下限值a:); b=input (请输入积分上限值b:); eps=input (请输入计算误差eps:); syms x f1 f2; f1=diff(f); %求函数的一阶导数 f2=diff(f1); %求函数的二阶导数 f2=-f2; %以下为求二阶导数的最大值fmax ff=inline(sym(f2),x); [xmax,fmax,exitflag,output]=fminbnd(ff,a,b); fmax; n=1; t=1; while abs(t)eps %求满足精度的等分区间个数n n=n+1; t=(-(b-a)/12)*((b-a)/n)^2*fmax; end disp(区间等分数n) n disp(复化梯形公式确定步长h) h=(b-a)/n %求复化梯形公式的近似值Tn Tn=subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b); for k=1:n-1 x1=a+h*k; Tn=Tn+2*subs(sym(f),findsym(sym(f)),x1); end Tn=Tn*h/2; disp(复化梯形公式求得近似解Tn) Tn=vpa(Tn,10) disp(复化梯形公式求积分的误差) z=int(sym(f),x,a,b); %求函数积分的真实值z w=Tn-z; %求复化梯形公式求值的误差w w=vpa(w,8) 运行结果： 请输入积分函数f:x*exp(x) 请输入积分下限值a:1 请输入积分上限值b:2 请输入计算误差eps:0.5*10^-7 区间等分数n n = 7019 复化梯形公式确定步长h h = 1.4247e-004 复化梯形公式求得近似解Tn Tn = 7.389056127 复化梯形公式求积分的误差 w = 0.000000028299571 （2）采用复化Simpson公式对上述函数用余项估计区间划分个数n，步长h，近似解Sn，以及与真实值误差w的源程序如下 function Fuhuasimpson clc f=input (请输入积分函数f:); a=input (请输入积分下限值a:); b=input (请输入积分上限值b:); eps=input (请输入计算误差eps:); syms x f1 f2; f1=diff(f); %求函数的一阶导数 f1=diff(f1); %求函数的二阶导数 f1=diff(f1); %求函数的三阶导数 f2=diff(f1); %求函数的四阶导数 f2=-f2; %以下为求二阶导数的最大值fmax ff=inline(sym(f2),x); [xmax,fmax,exitflag,output]=fminbnd(ff,a,b); fmax; n=1; t=1; while abs(t)eps %求满足精度的等分区间个数n n=n+1; t=(-(b-a)/2880)*((b-a)/n)^4*fmax; end disp(区间等分数n) n disp(复化simpsion公式确定步长h) h=(b-a)/n %求复化simpsion公式的近似值Tn Sn=subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b); for k=0:n-1 x1=a+h*(k+0.5); Sn=Sn+4*subs(sym(f),findsym(sym(f)),x1); end for k=1:n-1 x1=a+h*k; Sn=Sn+2*subs(sym(f),findsym(sym(f)),x1); end Sn=Sn*h/6; disp(复化simpsion公式求得近似解Sn) Sn=vpa(