第3章矩阵、数组和符号运算 解非线性方程组的函数命令 fsolve 其调用格式为: ◆x= fsolve(fun,x0):解非线性方程组最简单的调用格式 该式中除两个输入参数外,其余输入输出参数都可以缺省; ◆x, fval, exitflag, output, jacob}= fsolve(funx0, options,P1,P2…):解非 线性方程组最完整的调用格式 x0是表示零点数是猜测值的向量; options是优化迭代所采用参数的结构数组。 P1和P2是向函数fun传递的参数; x和fval是输出参数,所求零点的自变量值和函数值; output是输出此命令所用的计算方法、迭代次数等信息。 jaob是函数在x处的 jacobian
解非线性方程组的函数命令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。 第3章 矩阵、数组和符号运算
第3章矩阵、数组和符号运算 求方程组x1-07mx1-02csx2=0的根。 x2-0.7cos xl+0.2sin x2=0 首先编制函数文件fcm function y=fc(x y(1)=x(1)-0.7*sin(x(1)-02*sin(x(2) y(2)=x(2)-0.7*c0s(X(1)+0.2*sin(x(2) output y=y1)y(2 firstorderopt: 5.7877e-009 然后用 fsolve求解 iterations: 5 >x, fval, exitflag, output, jacob]=fsolve(fc,[1, 1.,D) unc Count: 16 %,1为初值 iterations: 4 algorithm: [1x43 char 0.33670.5553 Jacob fval (1,1)0.393 1.0e-008* (2,1)0.2313 0.20290.5242 (1,2)-0.1700 exitflag (2,2)1.1700
第3章 矩阵、数组和符号运算 求方程组 的根。 首先编制函数文件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 − + = − − = 2 0.7cos 1 0.2sin 2 0 1 0.7sin 1 0.2cos 2 0 x x x x x x 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