数值分析+复合梯形公式实验报告.docVIP

数学与信息工程学院 实 验 报 告 课程名称： 数值分析 实 验 室： 实验台号： 班 级： 姓 名： 实验日期： 2012 年 3 月 14 日 实验名称 复合梯形公式 实验目的 和 要 求 1. 掌握复合梯形公式的基本思想。2. 编程实现用复合梯形公式求积分。3. 熟悉matlab软件的使用。 实验内容和步骤： 实验内容： 根据matlab编写复合梯形公式求积分代码 实验步骤： （1）开启软件平台——MATLAB，编程； （2）根据数值解法步骤编写M文件； ① function t=trapz(f,a,b,n) h=(b-a)/n; sum=0; for k=1:n-1 x=a+k*h; sum=sum+feval(f,x); end t=h/2*(feval(f,a)+feval(f,b)+2*sum); ② 执行M文件 trapz(f,a,b,n) （3）观察运行结果 实验数据记录： 数据1：trapz(inline(sqrt(x)),0.5,1,2) ans = 0.4299 数据2：trapz(inline(sqrt(x)),0.5,1,3) ans = 0.4305 数据3：trapz(inline(sqrt(x)),0.5,1,4) ans = 0.4307 数据4：trapz(inline(sqrt(x)),0.5,1,100) ans = 0.4310 数据5： trapz(inline(sqrt(x)),0.5,1,1000) ans = 0.4310 实验结果分析： 不同的数据代入运行结果不同，n越大数据越准确，n大到一定时，数据不再变化。 M文件保存的时候要与函数文件名一致； M文件中f(a)要写成feval(f,a)； 执行文件时函数前要用inline，否则运行不了； 成绩评定 签字： 年 月 日