第八章代数方程和常微分方程求解.pptVIP

第8章 代数方程和常微分方程求解;8.1 代数方程求解;对于二维方程组图解，其解就是两条函数曲线的交点所对应的坐标数值。如果只有1个交点（或切点），则表示该方程组有1个解；如果有2个交点，则表示该方程组有2个解；如果没有交点，则表示该方程没有解。 例8-1 用图解法求解二维联立方程。 a=-2;b=2; % 定义横轴区间 ezplot(x^2+y^2-1.69,[a,b]); axis equal; % 控制坐标轴比例相等 hold on;grid on; ezplot(2.4*x^3-y+1.5,[a,b]); line([a,b],[0,0]);line([0,0],[b,a]); xlabel(\bf x);ylabel(\bf y); title(\bf 二维代数方程组的图解法) ;gtext(\bf f_1=x^2+y^2-1.3^2); gtext(\bf f_2=2.4x^3-y+1.5);;8.1.2 代数方程的解析解 求非线性方程或方程组解析解的函数调用格式： X=solve(fun,x) 其中，fun是符号方程的函数表达式，x是自变量，X是解析解。 应当指出，函数solve(fun,x)也可以用于求线性方组的解析解。 例8-2 求非线性解方程组解析解 ;% 二维非线性方程组的解析解 syms a b x y; f1=x^2-x*y-a; f2=y^2-x*y+b; disp( 二维非线性方程组的解析解：) [X,Y]=solve(f1,f2,x,y) M文件运行结果： 二维非线性方程组的解析解： x = a/(a-b)^(1/2) -a/(a-b)^(1/2) Y = 1/(a-b)^(1/2)*b -1/(a-b)^(1/2)*b;8.1.3 线性方程组的数值解 最简便方法是使用矩阵左除或是矩阵求逆的方法，求解线性方程组AX=b。 X= A\b X=inv(A)*b 其中，A是方程组的系数矩阵，b是常数向量，X是解析解。 例8-3 求线性方程组的数值解;% 线性方程组的数值解 AA=[1,1,1,1;1,2,-1,4;2,-3,-1,-5;3,1,2,11]; bb=[5;-2;-2;0]; % 线性方程组常数向量 disp( 采用矩阵左除求出线性方程组的解：) xx=AA\bb disp( 采用矩阵求逆求出线性方程组的解：) zx=inv(AA)*bb disp( 计算残量：) r=AA*zx-bb disp( 计算残量的模：) R=norm(r);M文件运行结果： 采用矩阵左除或矩阵求逆求出线性方程组的解： xx (zx)= 1.0000 2.0000 3.0000 -1.0000 计算残量： r = 1.0e-014 * 0.0888 0.2220 -0.4441 0.1776 计算残量的模： R = 5.3475e-015;8.1.4 非线性方程的数值解 1、一维非线性方程 对于一维非线性方程求解，可以看作是单变量的极小化问题，通过不断缩小搜索区间来逼近一维问题的真解。因此，可以使用一维非线性方程优化解函数来求解。其调用格式是： [x,fx,flag]=fzero(fun,x0) 其中，输入参数中：fun是非线性方程的函数表达式；x0是根的初值； 输出参数中：x是非线性方程的数值解；fx是数值解的函数值；返回参数flag0时，表示求解成功，否则求解出现问题。 函数fzero所使用的算法为二分法、secant法和逆二次插值法的组合。;例8-4 求解一维非线性方程 % 求解单变量x非线性方程 x0=0.1; % 解的初值 [xz,fz,flag]=fzero(atan(x)+exp(x),x0); disp( 求解成功性判断参数：), flag disp( 非线性方程的解：),xz disp( 非线性方程解的函数值：),fz M文件运行结果： 求解成功性判断参数： flag = 1 非线性方程的解： xz = -0.6066 非线性方程解的函数值： fz = -1.1102e-016;2、多维非线性方程组 对于多维非线性方程组 使用多维非线性方程组优化解函数求解的格式： [x,fval,flag]=fsolve(fun,x0) 其中，输入参数中fun是非线性方程组的向量函数表达式；x0是根的初值；输出参数中x是非线性方程（组）的数值解；fval是数值解的函数值；返回参数flag0时，表示求解成功. 函数f