matlab中方程根的近似计算要点.docxVIP

实验一 方程根的近似计算 一、问题 求非线性方程的根 二、实验目的 1、学会使用 matlab 中内部函数 roots、solve、 fsolve、fzero 求解方程，并用之 解决实际问题。 4、熟悉 Matlab 的编程思路，尤其是函数式 M 文件的编写方法。 三、预备知识 方程求根是初等数学的重要内容之一，也是科学和工程中经常碰到的数值 计算问题。它的一般形式是求方程 f(x)=0 的根。如果有 x* 使得 f(x*)=0 ，则称 x* 为 f(x)=0 的根，或函数 f(x) 的零点。并非所有的方程都能求出精确解或解析 解。理论上已经证明，用代数方法可以求出不超过 3 次的代数方程的解析解， 但对于次数大于等于 5 的代数方程，没有代数求根方法，即它的根不能用方程 系数的解析式表示。至于超越方程，通常很难求出其解析解。不存在解析解的 方程就需要结合具体方程（函数）的性质，使用作图法或数值法求出近似解。 而计算机的发展和普及又为这些方法提供了广阔的发展前景，使之成为科学和 工程中最实用的方法之一。下面介绍几种常见的求近似根的方法。 求方程近似解的简单方法 1.1 图形方法 —放大法求根 图形的方法是分析方程根的性态最简洁的方法。不过，不要总是想得到根 的精确值。这些值虽然粗糙但直观，多少个根，在何范围，一目了然。并且还 可以借助图形局部放大功能，将根定位得更加准确一些。 例 1.1 求方程 x5+2x2+4=0 的所有根及其大致分布范围。 解 （ 1）画出函数 f(x)=x 5+2x2+4 的图形，确定方程的实数根的大致范围。为此， 在 matlab 命令窗中输入 clf ezplot x-x, grid on hold on ezplot(x^5+2*x^2+4,[-2*pi,2*pi]) x5 +2 x2+4 6000 4000 2000 0 -2000 -4000 -6000 -6 -4 -2 0 2 4 6 x 1-1 函数 f(x)=x5+2x2+4 的图形 clf x=-2*pi:0.1:2*pi; y1=zeros(size(x)); y2= x.^5+2*x.^2+4; plot(x,y1,x,y2) grid on axis tight title(x^5+2x^2+4) xlabel(x) x5 +2x2+4 8000 6000 4000 2000 0 -2000 -4000 -6000 -8000 -6 -4 -2 0 2 4 6 x 从图 1-1 可见，它有一个实数根，大致分布在 -2 与 2 之间。 （ 2）将作图范围不断缩小，用放大法可得到精度越来越高的根的近似值。在 matlab 命令窗中先后键入 subplot(2,2,1) ezplot x-x, grid on, hold on, ezplot(x^5+2*x^2+4,[-2,2]) subplot(2,2,2) ezplot x-x, grid on, hold on, ezplot(x^5+2*x^2+4,[-2,-1]) subplot(2,2,3) ezplot x-x, grid on, hold on, ezplot(x^5+2*x^2+4,[-1.6,-1.5]) subplot(2,2,4) ezplot x-x, grid on, hold on, ezplot(x^5+2*x^2+4,[-1.55,-1.54]) x5 +2 x2+4 x5+2 x2+4 30 20 0 10 0 -10 -10 -20 -2 -1 0 1 2 -2 -1.5 -1 x x x5 +2 x2+4 x5+2 x2+4 1 0.1 0 0 -1 -0.1 -1.6 -1.55 -1.5 -1.55 -1.545 -1.54 x x 图 1-2 放大法求函数 f(x)=x 5+2x2+4 的根 由图 1-2 可知，方程的根在 -1.545 与 -1.54 之间。 1.2 数值方法 非线性方程 f(x)=0 求根的方法有区间法和迭代法两大类，二分法、弦位法 是区间法，简单迭代法和牛顿迭代法及其变形是迭代法，这里只给出二分法、 简单迭代法和牛顿迭代法的构造过程。 （1）根的隔离与二分法 根的隔离思想来源于连续函数的零点定理： 若函数 f(x) 在闭区间 [a,b] 上连续，且 f(a)f(b)0 ，则方程 f(x)=0 在(a,b)内至 少有一根 x* 。 二分法是最简单的求根方法，它是利用连续函数的零点定理，将含根区间 逐次减半缩小，取区间的中点构造收敛点列 {x n} 来逼近根 x *。 用该方法求 f(x)=0 的近似解可分两步做： 第一步，确定根的近似位置或大致范围，即确定一个区间 [a,b],使所求根是 位于这个区间内的唯一实