数值分析方法.docVIP

 《数值分析方法》 论文 学生姓名： 专 业： 电子科学与技术 年 级： 2008 学 号： 080712110026 2011年 12月 12 日 一、内容 给定积分 和 ,分别用下列方法计算积分值要求准确到 ,并比较分析计算时间. 变步长梯形法; 变步长 Simpson 法; 3) Romberg 方法. 二、原理与方法 变步长梯形法： 输入 被积函数f(x),积分区间端点a,b和允许误差. 输出 复合梯形积分值 步1 步2 步3 反复执行步4→步10 步4 步5 反复执行步6→步7 步6 步7 若x≥b，则退出本层循环 步8 步9 步10 若e≤,则退出循环 步11 步12 输出 变步长Simpson法： 步1 定义被积函数f(x) 步2 输入积分区间端点a,b和允许误差 步3 开始计时：tic 步4 调用基于变步长Simpson积分函数[S,n]=quad(‘f’,a,b,tol) 步5 结束计时：t2=toc Romberg积分法： 输入 被积函数f(x),积分区间端点a,b和允许误差 输出 Romberg积分值 步骤 调用Romberg积分程序 三、调试过程及实验结 变步长梯形法： ① function y=f(x) y=exp(x); a=1;b=3;tol=0.000005; I=Vsm(f,a,b,tol) t1 = 0.0103 I = 17.3673 ② function y=f(x) y=1/x; a=1;b=3;tol=0.000005; I=Vsm(f,a,b,tol) t1 = 0.0049 I = 1.0986 变步长Simpson法： ① function y=f(x) y=exp(x); a=1;b=3;tol=0.00005; y=Simpson(f,a,b,tol) t2 = 4.2152e-004 y = 17.3673 ② function y=f(x) y=1./x; a=1;b=3;tol=0.00005; y=Simpson(f,a,b,tol) t2 = 3.9313e-004 y = 1.09861 Romberg积分法： ① function y=f(x) y=exp(x); a=1;b=3;tol=0.000005; I=Romberg(f,a,b,tol) t3 = 1.8141e-004 I = 17.3673 ② function y=f(x) y=1/x; a=1;b=3;tol=0.000005; I=Romberg(f,a,b,tol) t3 = 1.7178e-004 I = 1.0986 四、总结 1.由以上结果可以看出：Romberg积分法所用时间最短，效率最高。 2.变步长Simpson积分法，采用的调用已知函数quad，然而具体的程序还很欠缺。对于几种积分方法，可以得出它们都是自适应方法，本实验中采取的是等距节点，事实上，如果被积函数在某处变化很激烈，在该处就应当设置较多的节点；反之只要设置较少的节点，自适应求积公式的目标，就是设法根据被积函数的性质确定求积节点的合理分布以获得最优的计算效果。 1