2020年复化梯形公式和复化Simpson公式.pptxVIP

学 海 无 涯 数值计算方法上机题目 3一、计算定积分的近似值：2?xe2 ? xe dx 1要求：（1）若用复化梯形公式和复化 Simpson 公式计算，要求误差限? ? 1 ?10?7 ，分别利用他们的余项估计对每种算法做出步长的事前估2计；分别利用复化梯形公式和复化 Simpson 公式计算定积分；将计算结果与精确解比较，并比较两种算法的计算量。1.复化梯形公式程序：程序 1（求f（x）的 n 阶导数：syms x f=x*exp(x) % 定 义 函 数 f（x） n=input( 输 入 所 求 导 数 阶 数 :) f2=diff(f,x,n) %求f(x)的n阶导数 结果1 输入 n=2 f2 =2*exp(x) + x*exp(x)程序 2：clc 学 海 无 涯 clear syms x % 定 义 自 变 量 x f=inline(x*exp(x),x) %定义函数f(x)=x*exp(x)，换函数时只需换该函数表达式 即可 f2=inline((2*exp(x) + x*exp(x)),x) %定义f(x)的二阶导数,输入程序1里求出的 f2即可。 f3=-(2*exp(x) + x*exp(x)) %因fminbnd（）函数求的是表达式的最小值，且要求表 达式带引号，故取负号，以便求最大值 e=5*10^(-8) %精度要求值 a=1 % 积 分 下 限 b=2 % 积 分 上 限 x1=fminbnd(f3,1,2) %求负的二阶导数的最小值点，也就是求二阶导数的最大值 点对应的x值 for n=2:1000000 % 求 等 分 数 n Rn=-(b-a)/12*((b-a)/n)^2*f2(x1) %计算余项 if abs(Rn)e % 用 余 项 进 行 判 断 break % 符 合 要 求 时 结 束 end end h=(b-a)/n % 求 h Tn1=0 for k=1:n-1 % 求 连 加 和 xk=a+k*h Tn1=Tn1+f(xk) end Tn=h/2*((f(a)+2*Tn1+f(b))) z=exp(2) R=Tn-z % 求 已 知 值 与 计 算 值 的 差 fprintf(用复化梯形算法计算的结果 Tn=) disp(Tn) fprintf(等分数 n=) disp(n) % 输 出 等 分 数 fprintf(已知值与计算值的误差 R=) disp(R) 输出结果显示：用复化梯形算法计算的结果 Tn= 7.3891等分数 n=7019已知值与计算值的误差 R= 2.8300e-008学 海 无 涯 2. Simpson 公式程序：程序 1：（求 f（x）的 n 阶导数）：syms x f=x*exp(x) % 定 义 函 数 f（x） n=input( 输 入 所 求 导 数 阶 数 :) f2=diff(f,x,n) %求f(x)的n阶导数 结果1 输入 n=4 f2 =4*exp(x) + x*exp(x)程序 2：clc clear syms x % 定 义 自 变 量 x f=inline(x*exp(x),x) %定义函数f(x)=x*exp(x)，换函数时只需换该函数表达式 即可 f2=inline((4*exp(x) + x*exp(x)),x) %定义f(x)的四阶导数，输入程序1里求 出的f2即可 f3=-(4*exp(x) + x*exp(x)) %因fminbnd（）函数求的是表达式的最小值，且要 求表达式带引号，故取负号，一边求最大值 e=5*10^(-8) % 精 度 要 求 值 a=1 % 积 分 下 限 b=2 % 积 分 上 限 x1=fminbnd(f3,1,2) %求负的四阶导数的最小值点，也就是求四阶导数的 最大值点对应的x值 for n=2:1000000 % 求 等 分 数 n Rn=-(b-a)/180*((b-a)/(2*n))^4*f2(x1) %计算余项 if abs(Rn)e % 用 余 项 进 行 判 断 break % 符 合 要 求 时 结 束 end end h=(b-a)/n % 求 h 学 海 无 涯 Sn1=0 Sn2=0 for k=0:n-1 % 求 两 组 连 加 和 xk=a+k*h xk1=xk+h/2 Sn1=Sn1+f(xk1) Sn2=Sn2+f(xk) end Sn=h/6*(f(a)+4*Sn1+2*(Sn2-f(a))+f(b)) %因Sn2多加了k=0时的值，故减去f(a) z=exp(2) R=Sn-z % 求 已 知 值 与 计 算 值 的 差 fprintf