计算方法c语言版实验1,计算方法上机实验1 西安石油大学.doc

计算方法上机实验1 西安石油大学

《计算方法》实验报告

实验一 非线性方程的迭代数值解法

二级学院： 计算机学院

专 业： 计算机科学与技术

指导教师： 爨莹

班级学号： 201007010226

姓 名： 张忠良

实验一 非线性方程的迭代数值解法

1、实验目的

1) 熟悉用牛顿法解非线性方程的过程；熟悉用弦截法求解非线性方程的过程

2) 编程实现牛顿法、弦截法求非线性方程的根。

2、实验要求

1)认真分析题目的条件和要求，复习相关的理论知识，选择适当的解决方案和算法；

2)编写上机实验程序，作好上机前的准备工作；

3、实验内容

1)用牛顿法求解的根，取初始值为10。

2) 用弦截法求解数学方程。

5、原理：

牛顿迭代法：

是用于求方程或方程组近似根的一种常用的算法设计方法。设方程为f(x)=0，用某种数学方法导出等价的形式 x(n+1) = x(n)–f(x(n))/f’(x(n)).然后按以下步骤执行：

(1) 选一个方程的近似根，赋给变量x1;

(2) 将x0的值保存于变量x1，然后计算g(x1)，并将结果存于变量x0;

(3) 当x0与x1的差的绝对值还小于指定的精度要求时，重复步骤(2)的计算。

若方程有根，并且用上述方法计算出来的近似根序列收敛，则按上述方法求得的x0就

认为是方程的根。

弦截法：

弦截法求方程的根是一种的基本方法，在计算机编程中常用。他的思路是这样的：任取两个数，判断这两个数的函数值，如果函数值是同号，换两个数再试，直到两个数x1，x2对应的函数值为异号时为止，这时方程的解肯定在这两个数x1，x2之间。连接这两点所对应的函数值，连线与x轴的交点为新的x，若f(x)与f(x1)同号，则把x当作新的x1，将新的x1与x2连接，如此循环……如果f(x)与f(x1)异号，则把把x当作新的x2，将x1与新的x2连接，循环……的根的基本思想是：将非线性方程f(x)=0逐步转化为线性方程来求解，即是依次用切线代替曲线，用线性函数的零点作为函数f(x)=0的近似值。

2)用弦截法求解根的基本思想是：

依次用弦线代替曲线，用线性函数的零点作为函数零点的近似值。

7、对应程序：

牛顿迭代法求解

#include

#include

double f1(double x){return pow(x,3)-155;}

double f2(double x){return 3*pow(x,2);}

void main()

{

int i=1;

double x0=10,x=0;

printf("X的初始值为： %.5lf\n\n",x0);

for(i=1;i<50;i++)

{

x=x0-f1(x0)/f2(x0);

if(fabs(x-x0)<1e-5) break;

x0=x;

printf("第 %d 次迭代得数为:X%d= %.5lf\n",i,i,x0);

}

printf("故 X=%.5lf\n",x0);

}

=============================================================

弦截法求解

求解数学方程的根0.002*9.81-1.4*e-5*X1.5-1.15*e-5*X2=0

/*求解数学方程的隔根区间

f(X)=0.002*9.81-1.4*e-5*X1.5-1.15*e-5*X2=0

f'(X)=-1.4*e-5*(1.5)*X0.5-1.15*e-5*(2)X ?? 0

f''(X)=-1.4*e-5*(1.5*0.5)X(-0.5)-1.15*e-5*(2)<0*/

#include

#include

#define N 0.5

double f(double x)

{ double f;

f=0.002*9.81-1.4*0.00001*pow(x,1.5)-1.15*0.00001*pow(x,2);

return f;

}

void main()

{

int i;

double Initial=36;

double a,b,fa,fb;

printf("求解根所在的隔根区间:\n");

printf("选取初始值为:%.2f,步长为:%.2f,列表可知:\n",Initial,N