matlab常用解方程及方程组函数资料.doc

  r=roots(c) 注意:c为一维向量,者返回指定多项式的所有根(包括复根),poly和roots是互为反运算,还有就是roots只能求解多项式的解  还有下面几个函数poly2sym、sym2poly、eig  >>syms x >>y=x^5+3*x^3+3; >>c=sym2poly(y);%求解多项式系数 >>r=roots(c); >>poly(r)    [r, p, k] = residue(b,a) >>b = [ 5 3 -2 7] >>a = [-4 0 8 3] >>[r, p, k] = residue(b,a)  (组)——使用最多的  g = solve(eq1,eq2,...,eqn,var1,var2,...,varn)  注意:eqn和varn可以是符号表达式,也可以是字符串表达式,但是使用符号表达式时不能有“=”号,假如说varn没有给出,使用findsym函数找出默认的求解变量。返回的g是一个结构体,以varn为字段。由于符号求解的局限性,好多情况下可能得到空矩阵,此时只能用数值解法  解方程A=solve('a*x^2 + b*x + c')  解方程组B=solve('a*u^2 + v^2', 'u - v = 1', 'a^2 - 5*a + 6')    [x,fval,exitflag,output]=fzero(fun,x0,options,p1,p2...)  fun是目标函数,可以是句柄(@)、inline函数或M文件名   x0是初值,可以是标量也可以是长度为2的向量,前者给定一个位置,后者是给定一个范围 options是优化参数,通过optimset设置,optimget获取,一般使用默认的就可以了,具体参照帮助 p1,p2...为需要传递的其它参数假如说(x/1446)^2+p/+(t/)*(log(1-x/1446)+(1-1/)*x/1446)=0的根,其中p,t是已知参数,但是每次都改变  那么目标函数如下三种书写格式,效果完全等效。注意参数列表中,未知数一定放第一位,其他参数放后面 (1)objfun=@(x,p,t)(x/1446).^2+p/+(t/).*(log(1-x/1446)+(1-1/).*x/1446); (2)objfun=inline('(x/1446).^2+p/+(t/).*(log(1-x/1446)+(1-1/).*x/1446)','x','p','t')  此时的调用格式如下 fzero(objfun,x0,options,p,t)%如果options使用的默认的话,那直接使用[],p和t就是我们需要传递的参数 fzero(@(x)objfun(x,p,t),x0,options)%这种格式与上面的等效 区别就是前者,将参数p和t作为fzero的参数进行传递,而后者是将p和t作为objfun的参数进行传递,没有本质区别  (3)function f=objfun(x,p,t)%以M文件格式书写目标函数    f=(x/1446).^2+