MATLAB中的solve()函数 求代数方程的符号解析解
(2013-08-05 11:32:43)
标签:
matlab函数教育字符串 |
分类: matlab基础 |
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: [4x1 sym]
u: [4x1 sym]
v: [4x1
sym]
对于查看具体的数值可以通过
Aa = A.a Au = A.u Av = A.v 命令来查看。
PS:对于solve求解的代数方程,默认的为eq=0,eq1=0,eq2=0….eqn=0;