第六章非线性方程求根-席答案.ppt

* * * * * * * * * * * * * * * * * * * * 五、割线法 在牛顿法中用差分近似来计算 ，也可以用最 新计算的两个 f 值来近似。 设定 开始迭代， 可以取得任意 小，迭代到满足容许误差为止。 六、非线性方程组的求解方法 1、一般非线性方程组及其向量表示法 2、不动点迭代法 与单个非线性方程迭代法类似,先化为等价的方程组 称为迭代向量函数 由此就可以建立一个迭代格式: 一般迭代法的收敛条件与单个方程迭代法的收敛条 件很类似。 “不动点迭代法” 3、Newton迭代法 设非线性方程组 若知道解的初值，解可以写为： 其中 是初值， 是未知的校正值。 关于 的一阶Taylor展开为： 矩阵形式为： 设令其等于0 其中 是Jacobian矩阵。 Jacobian矩阵的程序： 例：syms x1 x2; f1=x1*x2+x2^2; f2=x1^2+x2; Df=jacobian([f1,f2],[x1,x2]) 偏导数也可以用差分近似来计算： 是可选的任意小值。 Df = [ x2, x1+2*x2] [ 2*x1, 1] 已知非线性方程组如下： 3x12-cos(x2x3)-1/2=0 x1-81(x2+0.1)2+sin(x3)+1.06=0 exp(-x1x2)+20*x3+(10π-3)/3=0 求解，要求精度达到0.00001。 首先建立向量值函数 function f=fun5(x); %定义非线性方程组如下 %变量x1 x2 x3 %函数f1 f2 f3 syms x1 x2 x3 f1=3*x1-cos(x2*x3)-1/2; f2=x1^2-81*(x2+0.1)^2+sin(x3)+1.06; f3=exp(-x1*x2)+20*x3+(10*pi-3)/3; f=[f1 f2 f3]; 建立函数dfun来求方程组的Jacobian矩阵 function df=dfun(x); f=fun5(x); df=[diff(f,x1);diff(f,x2);diff(f,x3)]; %diff求函数的导数，后面会介绍 df=conj(df); 求解非线性方程组的newton_non.m function x=newton_non(x0,eps,n) con=0;%其中x0为迭代初值,eps为精度要求,n为最大%迭代步数,con用来记录结果是否收敛 for ii=1:n ff=subs(fun5(x0),{x1,x2,x3},{x0(1), x0(2), x0(3)}); df=subs(dfun(x0),{x1,x2,x3},{x0(1), x0(2), x0(3)}); %subs—元素替换运算 运行程序 在matlab中输入以下内容 x=newton_non([0.1 0.1 -0.1],0.00001,20) 结果： x= 0.5000 0.0000 -0.5236 x=x0-ff/df; if norm(x-x0)eps con=1;break;end x0=x; end * * * * * * * * * * * * * * * * * * * * * 第六章 非线性方程求根 非线性方程 的解称为 的根或零点。 本章主要研究单个非线性方程求根的一些方法，如图解法、二分法、牛顿迭代法、割线法等，另外简单介绍一下求解非线性方程组的逐次代换法和牛顿迭代法。 MATLAB命令 例1：解方程 解：在MATLAB里输入命令： X=solve(8*x^9+17*x^3-3*x=-1) 可解出9个解。 例2：解方程 解： X=solve(sin(x)=0) 可惜的是运行后只输出一个根x=0。 缺点：用solve命令不能求出周期函数所对应的全 部根。 例3：解非线性方程组 解：在MATLAB命令窗口输入命令： E1=sym(x^x-4=0);E2=sym(2*x*y+x=1); [x,y]=solve( E1,E2) x1=double(x),y1=double(y) 或 syms x y [x,y]=s