3.2 非线性方程方程组求解.pptxVIP

MATLAB 程序设计MATLAB 程序设计主讲人：王佐才 “黄山青年学者教授”道路与桥梁工程系 副主任Department of Civil EngineeringHefei University of Technology mailto:zw3z8@zw3z8@; wangzuocai1127@163.com第三章 非线性方程组本章内容：非线性代数方程的求解常微分方程求解 (不讲）第三章 非线性方程组1. 单变量非线性方程求解用到的Matlab命令fzeroFZERO Single-variable nonlinear zero finding. X = FZERO(FUN,X0) tries to find a zero of the function FUN near X0, if X0 is a scalar. It first finds an interval containing X0 where the function values of the interval endpoints differ in sign, then searches that interval for a zero. FUN is a function handle. FUN accepts real scalar input X and returns a real scalar function value F, evaluated at X. The value X returned by FZERO is near a point where FUN changes sign (if FUN is continuous), or NaN if the search fails.第三章 非线性方程组例1： （1）To find the zero of cosine between 1 and 2 （求出cos(x)在区间[1 2]之间的零解）x = fzero(@cos,[1 2]) （2）To find a zero of the function f(x)?=?x3?–?2x?–?5（求出方程f(x)?=?x3?–?2x?–?5的零解）Solution：write an anonymous function f:f = @(x)x.^3-2*x-5;Then find the zero near 2:z = fzero(f,2) 第三章 非线性方程组2. 多变量非线性方程组求解用到的Matlab命令fsolveSolve system of nonlinear equationsEquationSolves a problem specified byF(x) = 0for x, where x is a vector and F(x) is a function that returns a vector value.Syntaxx = fsolve(fun,x0)x = fsolve(fun,x0,options)x = fsolve(problem)[x,fval] = fsolve(fun,x0)[x,fval,exitflag] = fsolve(...)[x,fval,exitflag,output] = fsolve(...)[x,fval,exitflag,output,jacobian] = fsolve(...)第三章 非线性方程组例2： 先创建一个函数文件function y=fun(x)y=[x(1)-0.5*sin(x(1))-0.3*cos(x(2)),x(2)-0.5*cos(x(1))+0.3*sin(x(2))];然后求解x0=[0.1,0.1];fsolve(@fun,x0,optimset(fsolve))第三章 非线性方程组例3： This example solves the system of two equations and two unknowns: Rewrite the equations in the form F(x)?=?0:??Start your search for a solution at x0 = [-5 -5].First, write a file that computes F, the values of the equations at x.function F = myfun(x) F = [2*x(1) - x(2) - exp(-x(1)); -x(1) + 2*x(2) - exp(-x(2))];Save this function file as myfun.m somewhere on your MATLAB p