数值分析上机c 版.pdfVIP

..

数值分析论文

NUMERICALANALYSIS(THESIS)

题目数值分析论文

学生

指导教师鸿雁

学院

专业班级

2021年12月

目录

一二分法与牛顿迭代法1

二拉格朗日插值法6

三最小二乘法10

四复合辛普森求积公式15

五改进欧拉公式18

六列主元消去法21

七不动点迭代法26

.v.

..

一二分法与牛顿迭代法

一问题简介

在有气体参加的恒容反响体系里，反响体系的总压会随着反响的进展而改变。通过

反响的平衡常数可以求出反响后气体的总量，而气体分子的分压与其所占比率是成正比

的。

这类问题中关键是要求出反响后各分子的分压，各分子分压之和即为总压，其实质

即为求方程的根。例如：

在298K下，化学反响2OF=O+2F的平衡常数为0.410atm，如在298K下将OF

2222

-3

通入容器，当t=0时为1atm，问最后总压是多少.计算精度为10。

二数学模型

假设是理想气体，可知

2OF=O+2F

222

设氧的分压为p，平衡时有1-2ppp

平衡时有

32

整理得4p-1640p+164p-0410=0

函数关系式为

三算法选择与算法过程

由计算得f(0.2)=－0.1156，f(0.3)=0.0424

因此，有根区间为[0.2，0.3]

用求单根的二分法计算时：

可令a=0.2，b=0.3；最后计算得氧气分压：p=0.274609atm，总压为：

.v.

..

P=3p+(1-2p)=1.274609atm

用牛顿迭代法计算时：

令

那么

迭代公式为：

令进展迭代，最后计算得氧气分压为：p=0.274901，那么总压为：

P=3p+(1-2p)=1.274901

二分法的计算算法：

*includeiostream

*includemath.h

usingnamespacestd;

doublea=0.2,b=0.3,e=0.001,n;

doublef(floatc)