下面举一个栗子吧,下面是三个变量的符号非线性方程组,对于它的求解,我们可以用入戏的
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时候的数值解
版权声明:本文为qq_41823255原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。