计算方法上机实验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