高等数学上教学课件作者张圣勤黄勇林姜玉娟MATLAB实验ex55课件.pptVIP

机械工业出版社 Jixie gongye chubanshe 高等数学实验 MATLAB7.0实验五 代数方程的符号解 符号线性方程（组）的求解 非线性符号方程的求解 常微分方程的符号解 符号方程的求解 上海电机学院 张圣勤 符号方程求解的solve函数调用格式 solve(eq) 功能：对系统默认的符号变量求方程eq=0的根。 solve(eq,v) 功能：对指定变量v求解方程eq(v)=0的根。 [x1,x2,…xn]=solve(eq1,eq2,…eqn) 功能：对系统默认的一组符号变量求方程组eqi=0(i=1,2,…n)的根。 [v1,v2,…vn]=solve(eq1,eq2,…eqn,v1,v2,…vn) 功能：对指定的一组符号变量v1,v2,…vn求方程组eqi=0(i=1,2,…n)的根。 linsolve(A,B) 功能：求符号线性方程（组）AX=B的解。相当于X=sym(A)\sym(B) fsolve(f,x0) 功能：从x0开始搜索f=0的解。 fsolve(f,x0,options) 功能：根据指定的优化参数options从x0开始搜索f=0的解。 fsolve(f,x0,options,p1,p2…) 功能：优化参数option不是默认时，在p1,p2…条件下求f=0解。优化 参数option可取的值有0（默认）和1 [x,fv]=fsolve(f,x0,options,p1,p2…) 功能：优化参数option为默认时，在p1,p2…条件下求f=0解，并输出 根和目标函数值。 [x,fv,ex]=fsolve(f,x0,options,p1,p2…) 功能：优化参数option为默认时，在p1,p2…条件下求f=0解，并输出根 和目标函数值，并通过exitflag返回函数的退出状态。 [x,fv,ex,out]=fsolve(f,x0,options,p1,p2…) 功能：优化参数option为默认时，在p1,p2…条件下求f=0解，并给出 优化信息。 [x,fv,ex,out,jac]=fsolve(f,x0,options,p1,p2…) 功能：优化参数option为默认时，在p1,p2…条件下求f=0解，输出 值为x处的jacobian函数。 [x1,x2,…]=dsolve(eq1,eq2,…,cond1,cond2…,v) 功能：在初始条件为cond1,cond2…时求微分方程组eq1,eq2,…对指定 变量v的特解。 [x1,x2,…]=dsolve(eq1,eq2,…,cond1,cond2…,v) 功能：同[x1,x2,…]=dsolve(eq1,eq2,…,cond1,cond2…,v) 一、代数方程的符号解 求代数方程的符号解是通过调用solve函数实现的。用solve 函数求解一个代数方程时的调用格式一般是： solve（代数方程,未知变量）或x=solve（代数方程,未知变量） 当未知变量为系统默认变量时，未知变量的输入可以省略。 当求解由n个代数方程组成的方程组时调用的格式是： [未知变量组]=solve(代数方程组，未知变量组) 未知变量组中的各变量之间，代数方程组的各方程之间用逗号分隔，如果各未知变量是由系统默认的，则未知变量组的输入可以省略。 例1、求解高次符号方程 和方程 对y的解。 syms x y z a b %定义符号变量 solve(x^4-3*a*x^2+4*b) %求解高次方程 ans = 1/2*(6*a+2*(9*a^2-16*b)^(1/2))^(1/2) -1/2*(6*a+2*(9*a^2-16*b)^(1/2))^(1/2) 1/2*(6*a-2*(9*a^2-16*b)^(1/2))^(1/2) -1/2*(6*a-2*(9*a^2-16*b)^(1/2))^(1/2) solve(x^3+2*a*x*y-3*b*y^2,y) %对指定变量求解方程 ans = 1/6/b*(2*a+2*(a^2+3*b*x)^(1/2))*x 1/6/b*(2*a-2*(a^2+3*b*x)^(1/2))*x 例2、求解多元高次方程组 [x,y]=solve(x^3+2*x*y-3*y^2-2,x^3-3*x*y+y^2+5) %求解多元高次方程组 x = 1.80618931290919