Question

怎么用matlab算含有未知系数的方程?(即方程中未知量多于方程个数)高手进...帮忙啊

例如:有未知量x,y,a (a可以看为已知的系数)

有方程
a*x+y=10
8x-a*y^2=100

则x,y一定可以表达为a的多项式..

一下为我的m文件:
eq1=sym('a*x+y=10')
eq2=sym('8*x-a*y^2=100')
[x,y]=solve(eq1,eq2)

电脑得出适合的解:
x =

[ 1/2/a^3*(20*a^2+8+4*(20*a^2+4-25*a^3)^(1/2))]
[ 1/2/a^3*(20*a^2+8-4*(20*a^2+4-25*a^3)^(1/2))]

y =

[ -1/2/a^2*(20*a^2+8+4*(20*a^2+4-25*a^3)^(1/2))+10]
[ -1/2/a^2*(20*a^2+8-4*(20*a^2+4-25*a^3)^(1/2))+10]

上面仅仅是一个例子... 在含有3个未知数,两条方程的时候.. 把a看为未知数

然后我实际的问题是这个m文件:
eeq1=sym('m/(-a)=200/b')
eq2=sym('(c+2.11)^2+d^2=2.25')
eq3=sym('760/d=(1000-2*m)/(-c)')
eq4=sym('(a+1.47)^2+b^2=2.25')
[a,b,c,d]=solve(eq1,eq2,eq3,eq4)

基本的类型还是一样的a,b,c,d,m 5个未知量.. 4个方程... 把m看为已知系数
但是好像解不出..
有兴趣的,有matlab的可以crrl+c+v进去试试看啊....

不知道有没有办法解这个方程??

Answer 1

用[a,b,c,d]=solve(eeq1,eq2,eq3,eq4,'a,b,c,d')
结果：
a =

[ -.15000000000000000000000000000000e-1*m/(m^2+40000)*(98*m+100*(m^2+1584)^(1/2))]
[ -.15000000000000000000000000000000e-1*m/(m^2+40000)*(98*m+100*(m^2+1584)^(1/2))]
[ -.15000000000000000000000000000000e-1*m/(m^2+40000)*(98*m-100*(m^2+1584)^(1/2))]
[ -.15000000000000000000000000000000e-1*m/(m^2+40000)*(98*m-100*(m^2+1584)^(1/2))]

b =

[ 3.0000000000000000000000000000000/(m^2+40000)*(98*m+100*(m^2+1584)^(1/2))]
[ 3.0000000000000000000000000000000/(m^2+40000)*(98*m+100*(m^2+1584)^(1/2))]
[ 3.0000000000000000000000000000000/(m^2+40000)*(98*m-100*(m^2+1584)^(1/2))]
[ 3.0000000000000000000000000000000/(m^2+40000)*(98*m-100*(m^2+1584)^(1/2))]

c =

[ .13157894736842105263157894736842e-2/(25*m^2+9860000-25000*m)*(-40090*m+20045000+1900*(225*m^2-225000*m+24451676)^(1/2))*(-500.+m)]
[ .13157894736842105263157894736842e-2/(25*m^2+9860000-25000*m)*(-40090*m+20045000-1900*(225*m^2-225000*m+24451676)^(1/2))*(-500.+m)]
[ .13157894736842105263157894736842e-2/(25*m^2+9860000-25000*m)*(-40090*m+20045000+1900*(225*m^2-225000*m+24451676)^(1/2))*(-500.+m)]
[ .13157894736842105263157894736842e-2/(25*m^2+9860000-25000*m)*(-40090*m+20045000-1900*(225*m^2-225000*m+24451676)^(1/2))*(-500.+m)]

d =

[ 1/2/(25*m^2+9860000-25000*m)*(-40090*m+20045000+1900*(225*m^2-225000*m+24451676)^(1/2))]
[ 1/2/(25*m^2+9860000-25000*m)*(-40090*m+20045000-1900*(225*m^2-225000*m+24451676)^(1/2))]
[ 1/2/(25*m^2+9860000-25000*m)*(-40090*m+20045000+1900*(225*m^2-225000*m+24451676)^(1/2))]
[ 1/2/(25*m^2+9860000-25000*m)*(-40090*m+20045000-1900*(225*m^2-225000*m+24451676)^(1/2))]

Answer 2

楼上的可以将结构 simple 一下，这样看起来太难看了

一般的数学规模的方程都是可以解的，只是你的机器的配置的问题吧~