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;

 

//解一元二次方程

>> 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两个变量