matlab对符号方程组的求解（线性方程组以及非线性方程组）

下面举一个栗子吧，下面是三个变量的符号非线性方程组，对于它的求解，我们可以用入戏的

syms ro1 m1 s1 ro0 m0 s0 p1 p0 r
y1=sym('ro1*m1*s1=ro0*m0*s0');
y2=sym('ro1*m1^2*s1+p1*s1=ro0*m0^2*s0+p0*s0');
y3=sym('2*r*p1/((r-1)*ro1)+m1^2=2*r*p0/((r-1)*ro0)+m0^2');
[m0,m1,ro1]=solve(y1,y2,y3,m0,m1,ro1);
pretty(m0(1))%%m0可能为两个结果，所以我们要分两次去把m0的结果规范化，使人比较容易懂的那种。
pretty(m0(2))
pretty(m1(1))
pretty(m1(2))
pretty(ro1)%%经过solve函数得出的ro1值，可以知道，ro1的结果只有一个。

下面再给出一个简单的栗子吧。

syms a b c
y1=sym('a*b+2=c');
y2=sym('a+b=c');
[a,b]=solve(y1,y2,a,b)%输出符号解
c=2;
subs(a)%输出c为2时候的数值解
subs(b)%输出c为2时候的数值解