(第次符号计算内容多)矩阵数组和符号运算.ppt

第2章 矩阵、数组和符号运算 用 fzero 命令求解函数 的零点 （ 1）建立函数 f(x)的 M 文件。 function y=fun1(x) y=x.^4-4*x-5; （ 2）建立水平横轴的 M 文件。 function y=fun2(x) y=0; （ 3）用作图法估计函数零点位置。 fplot(fun1,[-5,5],r) hold on fplot(fun2,[-5,5],r) （ 4）用 zoom 和 ginput 命令获得零点的初始近似值 在程序中输入下列命令，可得到函数的局部放大图及鼠标操作线。 zoom on %局部放大命令 [tt]=ginput(1) [yy]=ginput(1) %用鼠标获取 2 个零点猜测值 zoom off %恢复原来图形大小 显示所得零点初始猜测值，结果为： tt = -0.9838 0.0001 yy = 1.8762 0.0091 （5）用函数 fzero 命令求函数的精确零点 [x,fval,exitflag]=fzero(’fun1’,tt(1),[]) %靠近 tt(1)点处的精确零点 [x,fval,exitflag]=fzero(’fun1’,yy(1),[]) %靠近 yy(1)点处的精确零点 结果为： Zero found near tt. x = -1 fval = 0 exitflag = 1 第2章 矩阵、数组和符号运算 Zero found near yy. x = 1.8812 fval = -6.2172e-015 exitflag = 1 解非线性方程组的函数命令fsolve 其调用格式为： ◆ x=fsolve(fun,x0) ：解非线性方程组最简单的调用格式。该式中除两个输入参数外，其余输入输出参数都可以缺省； ◆[x,fval,exitflag,output,jacob]=fsolve(fun,x0,options,P1,P2...)：解非线性方程组最完整的调用格式。 x0 是表示零点数是猜测值的向量； options 是 优 化 迭 代 所 采 用 参 数 的 结 构 数 组 。 P1 和 P2是向函数 fun 传递的参数； x 和 fval是输出参数，所求零点的自变量值和函数值； output是输出此命令所用的计算方法、迭代次数等信息。 jacob是函数在 x 处的 jacobian。 第2章 矩阵、数组和符号运算 第2章 矩阵、数组和符号运算 求方程组 的根。 首先编制函数文件fc.m function y=fc(x) y(1)=x(1)-0.7*sin(x(1))-0.2*sin(x(2)); y(2)=x(2)-0.7*cos(x(1))+0.2*sin(x(2)); y=[y(1),y(2)]; 然后用fsolve求解 [x,fval,exitflag,output,jacob]=fsolve(‘fc’, [1.,1.],[]) %[1.,1.]为初值 x = 0.3367 0.5553 fval = 1.0e-008 * 0.2029 0.5242 exitflag = 1 output = firstorderopt: 5.7877e-009 iterations: 5 funcCount: 16 cgiterations: 4 algorithm: [1x43 char] jacob = (1,1) 0.3393 (2,1) 0.2313 (1,2) -0.1700 (2,2) 1.1700 * * 第2章 矩阵、数组和符号运算 掌握内容： （1）了解 MATLAB 6.0 的符号变量，掌握 MATLAB 符号表达式、符号矩阵的两种创建方法。 （2）掌握 M