非线性方程求根实验报告.pdfVIP

计算机学院实验报告

课程：现代数值分析

实验名称：非线性方程求根

专业：网络工程

指导老师：

班级：1403班

学号：141110155

学生：

实验一非线性方程求根

一、实验目的

-x

1、用各种方法求方程f(x)=x-e=0在自变量区间[0,1]上的一个实

-5

根，精度为ε=10；

2、熟悉并理解二分法、简单迭代法和牛顿法的应用。

二、实验题目

-x

1、用二分法方程f(x)=x-e=0在自变量区间[0,1]上的一个实根，

-5

取定精度为ε=10。

-x

2、用简单迭代法方程f(x)=x-e=0在自变量区间[0,1]上的一个实

-5

根，取定精度为ε=10。

-x

3、用牛顿法方程f(x)=x-e=0在自变量区间[0,1]上的一个实根，取

-5

定精度为ε=10。

三、实验原理

1、二分法：对连续函数方程f(x)=0在[a，b]内求根。将所给区间二

分，在分点x=(b-a)/2判断是否f(x)=0；若是，则有根x=(b-a)/2。

否则，继续判断是否f（a）·f（x）0,若是，则令b=x,否则令a=x。

否则令a=x。重复此过程直至求出方程f(x)=0在[a,b]中的近似根为

止。

2、迭代法：将方程f(x)=0等价变换为x=ψ（x）形式，并建立相应

的迭代x=ψ（x）。

k+1

3、牛顿法：若已知方程的一个近似根x，则函数在点x附近可用一

00

阶泰勒多项式P(x)=f(x)+f’(x)(x-x)来近似，因此方程f(x)=0

1000

可近似表示为f(x)+f’(x)(x-x)=0设f’(x)≠0,则x=x-(f(x)/f’

000000

(x))。取x作为原方程新的近似根x，然后将x作为x代入上式。

0110

迭代为：xx-(f(x)/f’(x))。

k+1=0kk

四、实验内容

实验设备：MATLAB7.1软件、win7系统

实验步骤：

1、二分法

第一步：在MATLAB7.1软件，建立一个实现二分法的MATLAB函数文

件mbisec.m如下：

function[x,k]=mbisec(f,a,b,ep)

%用途：用二分法求非线性方程f(x)=0有根区间[a,b]中的一个根

%格式：x=mbisec(f,a,b,ep)f为函数表达式，

%a,b为区间左右端点，ep为精度，x,k分别返回近似根和二分次

数

x=(a+b)/2.0;k=0;

whileabs(feval(f,x))ep|(b-aep)//循环条件的限制

iffeval(f,x)