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;
//解一元二次方程
>> syms a b c x
>> solve('a*x^2+b*x+c=0',x)
ans =
-1/2*(b-(b^2-4*a*c)^(1/2))/a
-1/2*(b+(b^2-4*a*c)^(1/2))/a
求表达式的值(符号表达式转换为数值表达式)
>> subs(s,'x1,x2',{x1,x2}) %是是符号表达式,里面有X1,X2两个变量