第二章 基于matlab的科学计算—非线性方程（组）.pdfVIP

科学计算—理论、方法 及其基于MATLAB 的实现与分析 解非线性方程（组） （一）直接法 二分法： 设方程   在区间 上有唯一解，并且     ，如方程 f x  0 a, b f a f b  0 f x  x3 2.3x2  x sin x 0.3  0   （1） 首先要确定适当的包含根的区间，这可以依据闭区间上连续函数的介值 f 1  1 sin1  0 f 2  2 sin 2 0.9  0 定理来确定，例如，   ，   ，所以 方程 （1）至少有一个实根属于区间[1，2] ，图 1 表明区间[1，2] 中只含 有一个根，显然方程 （1）的根不易直接求得。在区间[-1，0]、[0，1] 和[1，2]的情形，如下图１所示 例１ plotNL_fun01.m plotNL_fun01 3 2 The Image of f(x)=x -2.3*x +x*sin(x)+0.3 1 0.5 0 -0.5 -1 -1.5 -2 -2.5 -1 -0.5 0 0.5 1 1.5 2 　图１ clear x=-1:0.05:2; f=x.^3-2.3*x.^2+x.*sin(x)+0.3; plot(x,f,r,x,0*x,k) title(The Image of f(x)=x^3-2.3*x^2+x*sin(x)+0.3) xlabel(\fontsize {12} \fontname {宋体} 图１) axis square *   二分法的求根过程：用 表示方程 在区间 上的根，对于给 x f x  0 a, b a b 定的精度要求  0 ，取区间a, b的中点x1  ，并按下式进行判断： 2    * f x1 0 x1 x