专题实验四 方程近似解的求法.pptVIP

专题实验四：方程近似解的求法 一、代数方程近似求解 相关概念 （1）二分法 （2）迭代法 牛顿法迭代公式 二、方程求解 (1)、roots(p) %求多项式的根，期中p是多项式向量 例：求x3-x2+x-1=0的根 解：p=[1 -1,1,-1]; roots(p) (2)、solve(fun) %求方程fun=0的符号解，如果不能求得精 确的符号解，可以计算可变精度的数值解。 例：用solve命令求方程x9+x8+1=0的根 解：solve(x^9+x^8+1) 给出了方程的数值解（32位有效数字的符号量） (3)、solve(fun,var) %对指定变量var求方程fun=0的符号解 例：解方程ax2+bx+c=0 解：syms a b c x; f=a*x^2+b*x+c; solve(f) 结果：1/2/a*(-b+(b^2-4*a*c)^(1/2)) 1/2/a*(-b-(b^2-4*a*c)^(1/2)) 如果不指明变量，系统默认为x，也可指定自变量，比如指 定b为自变量。 solve(f,b) 结果：-(a*x^2+c)/x (4)、fsolve(fun,x0) %求非线性方程fun=0在 估计值x0附近的近似解 解：fsolve(x-exp(-x),0) 结果： 0.56714316503697 (5)、fzero(fun,x0) %求函数fun在x0附近的根 例：求方程x-10x+2=0在x0=0.5附近的根 解：y=@(x)x-10^x+2; fplot(y,[-1,1]) %注意不能用plot fzero(y,0.5) 结果：0.3758 三、 函数极值的求法 多元函数的最小值问题，在Matlab中有2个经常使用的函数： 1、fminsearch 2、fminunc 1、 fminsearch 例：求函数f(x,y)= -(x+y)+(x^2+y^2+1)在x=1,y=2 附近的最小值点。 解决步骤： 1、建立M文件，保存函数f； M文件内容为： function f1=f1(x) a=-(x(1)+x(2)); b=(x(1)^2+x(2)^2+1); f1=a+b; 2、调用fminsearch函数求最值. 在命令窗口中，输入： x0=[1,2]; [x,fval]=fminsearch(@f1,x0) 3、输出结果为： X = 0.5000 0.5000 fval =0.5000 2、 fminunc * * 实验目的： （1）掌握求方程近似解的二分法、迭代法和牛顿迭代法的算法原理，会用MATLAB语言编程实现以上方法。 （2）学会使用MATLAB中内部函数fzero、fsolve、roots求解方程。 解方程（代数方程）是最常见的数学问题之一，也是众多应用领域中不可避免的问题之一。 目前还没有一般的解析方法来求解非线性方程，但如果在任意给定的精度下，能够解出方程的近似解，则可以认为求解问题已基本解决，至少可以满足实际需要。 这里主要介绍一些有效的求解方程的数值方法：二分法，迭代法 和 切线法。同时要求大家学会如何利用MATLAB 来求方程的近似解. 如果 f(x) 是一次多项式，称上面的方程为线性方程；否则称之为非线性方程 线性方程 与 非线性方程 这里主要讨论非线性方程的数值求解。 求方程的近似解，可分两步来做： 1、确定根的大致范围，就是确定一个区间[a,b]，使所求的根是位于这个区间内的唯一实根，这一步工作称为根的隔离，区间[a,b]称为所求实根的隔离区间。为了确定根的隔离区间，可以先画出y=f(x)的图形，然后从图上定出它与x轴交点的大概位置。 2、以根的隔离区间的端点作为根的初始近似值，逐步改善根的近似值的精确度，直到求得满足精确度要求的根的近似值。完成这一步工作有多种方法，这里介绍二分法、迭代法和切线法。 将有根区间进行对分，判断出解在某个分段内，然后再对该段对分，依次类推，直到满足给定的精度为止 数学原理：介值定理 设 f(x) 在 [a, b] 上连续，且 f(a) f(b)0，则由介值定理可得，在 (