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

数学与信息工程学院

实

课程名称：

实 验 室：

实验台号：

班 级：

姓 名：

实验日期： 验 报 告 数值分析 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

实验结果分析：

1. 不同的数据代入运行结果不同，n 越大数据越准确，n 大到一定时，数据不 再变化。 2. M 文件保存的时候要与函数文件名一致； 3. M 文件中 f(a)要写成 feval(f,a); 4. 执行文件时函数前要用 inline,否则运行不了；

成绩评定

签字：

年

月

日