计算机在材料科学中的应用上机,计算机在材料科学中的应用上机指导书-.doc

计算机在材料科学中的应用上机指导书-

《计算机在材料科学中的应用》上机指导书

山东科技大学材料学院

目录

实验一 数值计算程序设计 4学时 3

实验二 最小二乘法曲线拟合及回归分析 4学时 13

实验三 相图热力学计算 4学时 16

实验四 互联网在材料科学中的应用 2学时 18

实验一 数值计算程序设计 4学时

一、实验目的及任务要求

1、实验目的：

1)熟悉程序设计语言，了解简单的程序设计方法；

2)熟悉迭代法求方程根及牛顿法求方程的算法设计；

4)熟悉雅可比迭代法及其算法设计；

5)熟悉追赶法解线性方程组及其算法设计。

2、任务要求：

1)熟悉机房环境；

2)理解算法原理及程序编译过程。

二、实验操作过程

(一)方程求根 1学时

1、熟悉C语言编译软件界面

回顾C语言程序的编写方法，熟悉C编译软件的特点。

2、熟悉算法

1)二分法求方程根基本思路：

假设方程

f(x)=0

在区间[a，b]上有且仅有一根x*。

(1)首先取区间的中点x0，则

x0=(a+b) [a,b]对分为两个区间[a,x0]与[x0,b]；

(2)判断：f(a)*f(x0)>0：若不等式成立，则方程的根在[x0,b]，若不等式不成立，则方程的根在[a,x0]，根所在的新区间的端点仍然命名为［a,b］；

(3)判定根是否达到预定的精度要求

|xk+1-xk|

2)牛顿法求根

牛顿迭代法求根的方法是，将已知方程式f(x)=0看作函数式y＝f(x)。对函数f(x)在xk处进行一阶泰勒展开,即可得:

f(x)=f(xk)+f'(xk)(x-xk)+…

函数f(x)=0时x值即为方程的根，即函数曲线与x轴的交点值。

f(x)=f(xk)+f'(xk)(x-xk)＝0

x=xk-f(xk)/f'(xk)

设x=xk+1,则可得牛顿迭代求根的计算公式为：xk+1=xk-xk-f(xk)/f'(xk)

例如：已知方程式 f(x)=x3+4x2-10=0

求得导函数 f'(x)=3x2+8x

因此可得牛顿迭代式为：xk+1=xk-

几何解释：

x1=x0-

牛顿迭代法求根计算步骤：

1)输入已知条件x0、ε、N：

x0—初始迭代点，可判断其迭代收敛性；

ε――计算精度，一般取ε=0.00001

N――迭代预定最大次数，即当迭代次数超过N次之后，则主观认为迭代发散或失败，此时停止计算。一般取N=1000

2)计算函数值f(x0)与导数值f’(x0)

3)迭代求解：x1=x0-

4)精度判断：｜x1-x0|<=ε，条件成立时，则输出计算结果；条件不成立，则x0=x1,k=k+1，返回第2步重新计算。

3、程序设计

1)读下面的程序，并改写采用二分法求解方程x3-x-1=0在区间［0，2］中的根(精度为e=10-6)

二分法求一元n次方程的近似根:?

?#include

#include

main()

{ float x0,x1,x2,fx0,fx1,fx2;

do

{ printf("Enter x1,x2:");

scanf("%f,%f",&x1,&x2);

fx1=2*x1*x1*x1-4*x1*x1+3*x1-6;

fx2=2*x2*x2*x2-4*x2*x2+3*x2-6;

} while(fx1*fx2>0); /* 判断x1,x2 之间是否有根*/

do

{ x0=(x1+x2)/2; /*取中点*/

fx0=2*x0*x0*x0-4*x0*x0+3*x0-6;

if (fx0*fx1<0) /*判断x1,x0之间是否有根*/

{ x2=x0; /*是的话用x0来替换x2*/

fx2=fx0;

}

else { x1=x0;