c语言数值分析实验报告,数值分析实验报告之数值积分.doc

数学与计算科学学院

实 验 报 告

实验项目名称 数值积分

所属课程名称 数值方法B

实 验 类 型 验证

实 验 日 期 2013.10.21

班 级

学 号

姓 名

成 绩

一、实验概述：【实验目的】

1．熟悉C语言与MATLAB的编程；

2．学会使用梯形公式、辛普森公式、复化梯形公式、复化辛普森公式求积分的方法；

3．比较各方法的精度；

4．用编程软件写出上述四个公式，并实例化。

5．此外，本实验还附加了cotes公式以及复化cotes公式的C语言源程序。

【实验原理】

1．梯形公式：

2．辛普森公式：

3．复化梯形公式：

4．复化辛普森公式：

其中。

【实验环境】

1.硬件环境：

HP

Microsoft

76481-640-8834005-23929

HP Corporation

Intel(R) Core(TM)

I5-2400 CPU @ 3.10GHz

3.09GHz,3.16GB的内存

2.软件环境：

Microsoft Windows XP

Professional

版本 2002

Service Pack 3

二、实验内容：【实验方案】

1．用复合求积公式计算、、，并比较各方法的精度；

2.分别讨论当区间n等分，当n=8，10, 100, 1000, 10000时比较n取值不同时对数值精度的影响的结果

【实验过程】(实验步骤、记录、数据、分析)

1.

用MATLAB中的int求积函数(源程序详见附录1)得到该积分的准确值为x*=0.111571776

梯形公式：0.1000000000

精度：一位有效数字

辛普森公式：0.1117647059

精度：三位有效数字

该问C语言编程详见附录1-4/1-5

(1)当n=8时，结果为：

Tn=0.1114023545

Sn=0.1115718133

(2)当n=10时，结果为：

Tn=0.1114633808

Sn=0.1115717910

(3)当n=100时，结果为：

Tn=0.1115706923

Sn=0.1115717757

表1 各积分方法的比较

等分段数n值

备注

复化梯形公式

复合辛普生公式

n=8

积分值

0.1114023545

0.1115718133

精度

三位有效数字

七位有效数字

n=10

积分值

0.1114633808

0.1115717910

精度

三位有效数字

七位有效数字

n=100

积分值

0.1115706923

0.1115717757

精度

五位有效数字

八位有效数字

n=1000

积分值

0.1115717648

0.1115717757

精度

七位有效数字

八位有效数字

n=10000

积分值

0.1115717755

0.1115717757

精度

八位有效数字

八位有效数字

2.

用MATLAB中的int求积函数得到该积分的准确值为x*=1.852112521

梯形公式：0.897530048810325

精度：一位有效数字

辛普森公式：1.135538476301203

精度：一位有效数字

该问MATLAB编程详见附录1-6/1-7

(1)当n=8时，结果为：

Tn=1.337355373803341

Sn=1.537749799432108

(2)当n=10时，结果为：

Tn=1.391477663937609

Sn=1.570881071946639

(3)当n=100时，结果为：

Tn=1.706124330945894

Sn=1.763114638290155

表2 各积分方法的比较

等分段数n值

备注

复化梯形公式

复合辛普生公式

n=8

积分值

1.337355373803341

1.537749799432108

精度

一位有效数字

二位有效数字

n=10

积分值

1.391477663937609

1.570881071946639

精度

一位有效数字

二位有效数字

n=100

积分值

1.706124330945894

1.763114638290155

精度

二位有效数字

三位有效数字

n=10