数值分析

数值分析上机报告

前言

随着计算机技术的高速发展，越来越多的科技工作者使用计算机进行科学研究和解决工程技术问题。数值分析（或计算方法）课程的内容是科学工程计算的必备知识，已经成为众多理工科大学生、研究生的必修课程，越来越受到重视。

由于工程实际中所遇到的数学模型求解过程迭代次数很多，计算量很大，所以需要借助很多编程软件来解决，得到一个满足误差限的解。本文所计算题目，均采用C++编程。在本文中使用C++编写了牛顿法、牛顿-Steffensen法方程求解的程序和雅格比法、高斯－赛德尔迭代法求解方程组的程序及Ru n ge-Kutt a4阶算法，并通过实例求解验证了其可行性，比较了求解同一种问题时不同方法之间的优缺性，其中包含解的精确度和解的收敛速度两个重要指标。

一 牛顿法和牛顿-Steffensen 法迭代求解的比较

1. 计算题目

分别用牛顿法，及基于牛顿算法下的Steffensen 加速法

(1) 求ln(x +sin x )=0的根。初值x0分别取0.1, 1,1.5, 2, 4进行计算。

(2) 求sin x =0的根。初值x0分别取1，1.4，1.6, 1.8，3进行计算。

分析其中遇到的现象与问题。

2. 计算过程和结果

1.对方程ln(x +sin x )=0，其导数有些复杂，我们可以对其进行变形，即求解x+sinx=1的解。使用牛顿法，令1sin )(-+=x x x f ，则x x f cos 1)(+='，直至51101||-+⨯<-k k x x 时，结束迭代；然后再使用基于牛顿法的Steffensen 加速法进

行计算，直至51101||-+⨯<-k k x x 时，结束迭代。其迭代结果与迭代次数如下表所

示（注N1为牛顿法迭代次数，N2为基于牛顿法Steffensen 加速法迭代次数）：

2.对方程sin x =0，使用牛顿法时，令x x f sin )(=，使用牛顿法计算，直至51101||-+⨯<-k k x x 时，结束迭代；然后依据Steffensen 加速法进行编程计算，直

至51101||-+⨯<-k k x x 时，结束迭代。其迭代结果与迭代次数如下表所示：

附：

第一问牛顿法

#include

#include

#define EPS 1.0e-5

int main()

{

double x=0.0;

double y=0.0;

int n=0;

cout<<"请输入初值X0"<

cin>>x;

do

{

y=x;

x=x-(x+sin(x)-1)/(1+cos(x));

n++;

}

while(fabs(x-y)>=EPS);

cout<<"X"<

cout<<"迭代次数为N="<

}

第一问牛顿-Steffensen法

#include

#include

#define EPS 1.0e-5

int main()

{

double x=0.0;

double y=0.0;

double z=0.0;

double a=0.0;

int n=0;

cout<<"请输入初值x0"<

cin>>x;

do

{

a=x;

y=x-(x+sin(x)-1)/(1+cos(x));

z=y-(y+sin(y)-1)/(1+cos(y));

x=x-(y-x)*(y-x)/(z-2*y+x);

n++;

}

while(fabs(x-a)>=EPS);

cout<<"x"<

cout<<"迭代次数N="<

}

第二问牛顿法

#include

#include

#define EPS 1.0e-5

int main()

{

double x=0.0;

double y=0.0;

int n=0;

cout<<"请输入初值X0"<

cin>>x;

do

{

y=x;

x=x-tan(x);

n++;

}

while(fabs(x-y)>=EPS);

cout<<"X"<

cout<<"迭代次数为N="<

}

第二问牛顿-Steffensen法

#include

#include

int main()

{

double x=0.0;

double y=0.0;

double z=0.0;

double a=0.0;

int n=0;

cout<<"请输入初值x0"<

cin>>x;

do

{

a=x;

y=x-tan(x);

z=y-tan(y);

x=x-(y-x)*(y-x)/(z-2*y+x);

n++;

}

while(fabs(x-a)>=1.0E-05);