dsolve解微分方程

Matlab中solve函数主要是用来求解线性方程组的解析解或者精确解。对于得出的结果是符号变量，可以通过vpa()得出任意位数的数值解！

solve函数的语法定义主要有以下四种：

solve(eq)
solve(eq, var)
solve(eq1, eq2, …, eqn)
g = solve(eq1, eq2, …, eqn, var1, var2, …, varn)

eq代表方程，var代表的是变量。

例1：

syms a b c x; solve(‘a*x^2 +b*x + c’)

当没有指定变量的时候matlab默认求解的是关于x的一元二次方程的解，求解的结果为：

ans = -(b + (b^2 – 4*a*c)^(1/2))/(2*a) -(b – (b^2 – 4*a*c)^(1/2))/(2*a)

当指定变量为b的时候：

syms a b c x; solve(‘a*x^2 +b*x + c’,’b’)

求解的结果为：

ans = -(a*x^2 + c)/x

从上面的例子很容易理解语法1,2。

例2：对于方程组的情况

syms x; S = solve(‘x + y = 1′,’x – 11*y = 5′); S = [S.x S.y]

求解的结果为：

S = [ 4/3, -1/3]

例3：

syms a u v; A = solve(‘a*u^2 +v^2′, ‘u – v = 1′, ‘a^2 – 5*a + 6′)

的求解结果为

A = a: [4×1 sym] u: [4×1 sym] v: [4×1 sym] 对于查看具体的数值可以通过

Aa = A.a Au = A.u Av = A.v

命令来查看。

PS：对于solve求解的方程，默认的为eq=0，eq1=0，eq2=0….eqn=0;