数值分析上机(C++版)教程.docxVIP

数值分析论文 NUMERICAL ANALYSIS(THESIS) 题 目数值分析论文学生姓名指导教师张鸿雁学 院专业班级 2016年12月  目录 一 二分法与牛顿迭代法 1 二 拉格朗日插值法 6 三 最小二乘法 10 四 复合辛普森求积公式 15 五 改进欧拉公式 18 六 列主元消去法 21 七 不动点迭代法 26  PAGE \* MERGEFORMAT32 一 二分法与牛顿迭代法 一 问题简介 在有气体参加的恒容反应体系里，反应体系的总压会随着反应的进行而改变。通过反应的平衡常数可以求出反应后气体的总量，而气体分子的分压与其所占比率是成正比的。 这类问题中关键是要求出反应后各分子的分压，各分子分压之和即为总压，其实质即为求方程的根。例如： 在298K下，化学反应2OF2=O2+2F2的平衡常数为0.410atm，如在298K下将OF2通入容器，当t=0时为1atm，问最后总压是多少？计算精度为10-3。 二 数学模型 假设是理想气体，可知 2OF2=O2+2F2 设氧的分压为p，平衡时有 1-2p p p 平衡时有 4p3(1-2p)2=0.410 整理得 4p3-1.640p2+1.64p-0.410=0 函数关系式为 fp=4p3-1.640p2+1.64p-0.410=0 三 算法选择与算法过程 由计算得 f(0.2)=－0.1156，f(0.3)=0.0424 因此，有根区间为[0.2，0.3] 用求单根的二分法计算时： 可令a=0.2，b=0.3；最后计算得氧气分压：p=0.274609atm，总压为： P=3p+(1-2p)=1.274609atm 用牛顿迭代法计算时： 令 fx=4x3-1.640x2+1.64x-0.410 则 fx=12x2-3.280x+1.64 迭代公式为：xk+1=xk- fxfx 令x0=0.2进行迭代，最后计算得氧气分压为：p=0.274901，则总压为： P=3p+(1-2p)=1.274901 二分法的计算算法： #include iostream #includemath.h using namespace std; double a=0.2,b=0.3,e=0.001,n; double f(float c) { n=4*c*c*c-1.6400000*c*c+1.64*c-0.410; return n; } void main() { double c,m,p,a=0.2,b=0.3; for(p=a;b-a=0.001;) { c=(a+b)/2; m=f(c); if(m0) {a=c;} else if(m0) {b=c;} else{p=c;break;} p=(a+b)/2; } cout p=pendl; cout\n a=aendl; cout\n b=bendl; } 牛顿迭代法的算法： #includeiostream #includestring #includecmath using namespace std; double f(double x) { double m,n; m=4*x*x*x-1.640*x*x+1.64*x-0.410; n=12*x*x-3.280*x+1.64; return x-m/n; } void newton(double x,double d) { double a=x; double b=f(a); int k; //记录循环的次数 for(k=1;fabs(a-b)d;k++) { a=b; b=f(a); if(k100) { cout迭代失败，该函数可能不收敛!endl; } } cout a=aendl; cout\n b=bendl; cout\n k=kendl; return; } int main(){ cout请输入初始值x0和要求得结果的精度：; double x,d; cinxd; newton(x,d); return 0; } 四 数值实验过程 二分法通过VC6.0程序运行的结果如下图所示： 牛顿迭代法通过VC6.0的计算