数学实验第三次讲稿.ppt

数学实验第三讲; 1) 掌握方程求解的三种解法：解析 法、数值解法以及图形表示解的方法; 2) 学会使用MATLAB软件求解析解、数值解 和图形解； 3） 通过范例学习怎样建立方程模型和 分析问题的思想。;一，建立方程;数学建模的一般步骤;二，方程求解;1，方程求解之解析方法;2，方程求解之图形放大法; 方程 f(x)=0 1）建立坐标系，画曲线f(x)； 2）观察曲线f(x)与x轴相交的交点； 3）将其中一个交点进行局部放大； 4）该交点的横坐标值就是方程的根。 ;例: 求??程 x5 +2x2 + 4 = 0 的一个根.; 逐次缩小区间，观察一个根在-1.55~-1.5之间。;3,方程求解之迭代法;解： 第一步 构造迭代函数： x=j (x) ;第二步 迭代 设定初值 x0=1， xn+1 =j (xn)，n =0，1，… 用 MATLAB 编程（died2.m文件);序号; 迭代函数j2(x)和j3(x)的选取是成功的。精确解为 x=1.8393。 并且选取函数j2(x)、j3(x)其收敛速度不一致，前者的速度快些！;当今最流行的迭代法是牛顿法以及由此改进的一些方法，比如拟牛顿法等。其基本的思想就是构造迭代格式是利用函数的导数，这类方法有收敛速度快，稳定性好等特点。对低维和高维情况都适合，也是当今一些软件均采用的方法。当然，因为需要函数的导数信息，所以自然对不可微的问题受到制约。该方法的迭代格式为： xk+1=xk-f(xk)/f’(xk);例:利用牛顿法求方程x3 -x2 -x-1 = 0的根. 方法:第一步,给出函数的导函数3x2-2x-1; 第二步,给出函数的迭代格式: xk+1=xk-f(xk)/f’(xk); 设置一定的精度要求,达到即终止. 定义函数m文件: function ff=mynewton(x) ff=(x^3 –x^2 -x-1 )/(3*x^2-2*x-1); 定义命令m文件: x0=1; x1=x0-mynewton (x0); while abs(x1-x0)>0.0001 x0=x1; x1=x0-mynewton (x0); end;三，解方程函数格式及例子;输出: -1.2131 -0.9017 + 0.5753i -0.9017 - 0.5753i -0.2694 + 0.9406i -0.2694 - 0.9406i 0.4168 + 0.8419i 0.4168 - 0.8419i 0.8608 + 0.3344i 0.8608 - 0.3344i;4、线性方程组：AX = b 其中A是m×n阶矩阵，b是m维向量。; solution to the following linear system of equations: You can formulate and solve the problem as A = [ 3 11 -2; 1 1 -2; 1 -1 1]; b = [ 7; 4; 19]; x = A\b x = 13.2188 -2.3438 3.4375;例: AX = b,;函数fzero;1、方程(组)， f1(x) = 0,…,fn(x) = 0, x = (x1,…,xn);输出:[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))];例2： 解方程： x3-2x2=x-1 解： s=solve('x^3-2*x^2=x-1') double(s); 例3 求解方程： tan(x)-sin(x)=0 ;输入： [x,y]=solve('x^2*y^2,x-(y/2)-b') 输出：x = [ 0] ， y =[ -2*b] [ 0]， [ -2*b] （符号解