文档介绍:一般的代数方程函数 solve 用于求解一般代数方程的根,假定 S为符号表达式,命令 solve (S) 求解表达式等于 0的根,也可以再输入一个参数指定未知数。例: syms a b c x S=a * x^2+b * x+c ; solve(S ) ans = [ 1/2/a * (-b+(b^2-4 *a* c)^(1/2))] [ 1/2/a * (-b-(b^2-4 *a* c)^(1/2))] b= solve(S,b ) b = -(a * x^2+c)/x 线性方程组线性方程组的求解问题可以表述为:给定两个矩阵 A和B,求解满足方程 AX=B 或 XA=B 的矩阵 X。方程 AX=B 的解用 X=A\B 或 X=inv (A) *B表示;方程 XA=B 的解用 X=B/A 或 X=B * inv (A) 表示。不过斜杠和反斜杠运算符计算更准确,占用内存更小,算得更快。线性微分方程函数 dsolve 用于线性常微分方程(组)的符号求解。在方程中用大写字母 D表示一次微分,D2 , D3 分别表示二阶、三阶微分,符号 D2y 相当于 y关于 t的二阶导数。函数 dsolve 的输出方式格式说明 y= dsolve (‘ Dyt =y0 *y’ ) 一个方程,一个输出参数[ u,v ]= dsolve (‘ Du =v’,’ Dv =u ’) 两个方程,两个输出参数 S= dsolve (‘ Df=g’,’ Dg =h’,’ Dh =-2 * f ‘)方程组的解以 结构数组的形式输出例 1 求21u dt du??的通解. 解输入命令: dsolve('Du=1+u^2','t')结果: u = tg(t-c ) 例 2 求微分方程的特解.?????????? 15 )0(',0)0( 0 29 4 2 2yy y dx dy dx yd解输入命令: y=dsolve('D2y+4 * Dy+29 * y=0','y(0)=0,Dy(0)=15','x') 结果为: y =3e -2x sin (5x) 例3 求微分方程组的通解.??????????????????zyxdt dz zyxdt dy zyxdt dx244 354 332解输入命令: [x,y,z]= dsolve('Dx =2 * x-3 * y+3 * z','Dy =4 * x-5 * y+3 * z','Dz =4 * x-4 * y+2 * z', 't') ; x=simple(x) % 将x化简 y=simple(y) z=simple(z) 结果为: x = (c 1 -c 2 +c 3 +c 2e -3t -c 3e -3t )e 2t y = -c 1e -4t +c 2e -4t +c 2e -3t -c 3e -3t +(c 1 -c 2 +c 3 )e 2t z = (-c 1e -4t +c 2e -4t +c 1 -c 2 +c 3 )e 2t 非线性微分方程[t, x]=so