数值实验2.docVIP

数值分析实验报告 实验序号：02　　 实验项目名称：数值积分与数值微分 学　　号 201141210230 姓　　名 谢冬冬 专业、班 信计（2）班 实验地点 K7-405 指导教师 汪玉霞 时 间 2013.10.23 一．实验目的 （1）由函数的个节点处函数值得出复化梯形积分公式； （2）由函数的个节点处函数值得出复化Simpson积分公式； （3）函数的变步长积分公式； （4）熟悉Romberg求积方法的原理与计算流程 二、实验原理 （1） （2） （3） （4） 三、实验设备要求与实验原理 1、支持Intel Pentium Ⅲ及其以上CPU，内存256MB以上、硬盘1GB以上容量的微机；软件配有Windows98/2000/XP操作系统及MATLAB软件； 2、实验过程中，务必分析实验结果，按要求写出实验报告。（建议同时网上提交电子版实验报告：wangty6895@126.com） 四、实验内容与步骤 第一个程序:（复化梯形公式） 用复化梯形公式计算y=1/(1+x^2),在x∈[0,1]上的积分的近似值。 function I=xiedongdong1(x,y) n=length(x);h=(x(n)-x(1))/(n-1);a=[1/2 ones(1,n-2) 1/2]; I=h*sum(a.*y); end 运行结果 将区间分成10个小区间 x=0:0.1:1;y=1./(1+x.^2);I=fuhetixing(x,y) I = 0.784981497226790 pi/4 ans = 0.785398163397448 WUCHA=pi/4-I WUCHA = 4.166661706586083e-004 第二个程序：（复化Simpson公式） (1)用复化Simpson公式计算y=1/(1+x^2),在x∈[0,1]上的积分的近似值。 将区间二等分。 预先定义被积函数 function y=bjhs(x) y=1/(1+x^2); end function s=xiedongdong2simpson(a,b,n) if nargin3 n=100; end banh=(b-a)/(2*n);s1=0;s2=0; for k=1:n %计算分数节点函数值的和 x=a+banh*(2*k-1);s1=s1+bjhs(x); end for k=1:n-1 %计算区间内部整数节点处的函数值的和 x=a+banh*2*k;s2=s2+bjhs(x); end s=banh*(bjhs(a)+bjhs(b)+4*s1+2*s2)/3; 运行结果： s=xiedongdong2simpson(0,1,2) s = 0.785392156862745 将区间分成10个小区间并计算误差 s=xiedongdong2simpson(0,1,10) s = 0.785398163242446 wucha=pi/4-s wucha = 1.550019002038994e-010 （2）用复化Simpson公式计算y=sin(x),在x∈[0, π]上的积分的近似值。 将区间10等分。 function y=bjhs(x) y=sin(x); end s=xiedongdong2simpson(0,pi,10) s = 2.000006784441801 第三个程序：（变步长公式） function I=xiedongdong3bianbuchang(fun,a,b,ep) %fun为被积函数的解析式 if nargin4;ep=1e-5;end N=1;h=b-a;T=h*(feval(fun,a)+feval(fun,b))/2%即为T1 while 1 %一直执行 h=h/2;I=T/2; for k=1:N I=I+h*feval(fun,a+(2*k-1)*h);%公式4.17，K=1时即T2 end if abs(I-T)ep%判断计算停止跳出while的条件 break; end N=2*N;%步长减半 T=I; end 运行结果 用变步长积分公式计算fun=1/(1+x^2),在x∈[0,1]上的积分的近似值。 fun=inline(1/(1+x^2));I=xiedongdong3bianbuchang(fun,0,1,1e-5) I = 0.785395620265938 用变步长积分公式计算fun=sin（x）,在x∈[0, π]上的积分的近似值 fun=inline(sin(x));I