登录论坛

查看完整版本 : [求助]请大家看看这个方程组用solve为什么解不出来啊


chinapenguin
2009-02-27, 00:05
我用solve函数,6个方程,6个未知数,为什么解不出来呢。请大家指点一下吧。我还试过fsolve,也不行。是不是需要数值方法呢?菜鸟一只,请多包涵。谢谢!

syms A b1 b2 g1 g2 R0

S1=g1*(log((R0+A)*exp(b1*16000)-A)-log(R0))/(b1*(R0+A))-1000; S2=g2*(log((R0+A)*exp(b2*44)-A)-log(R0))/(b2*(R0+A))-1000;

S3=g2*(log((R0+A)*exp(b2*22)-A)-log(R0))/(b2*(R0+A))+g1*(log((R0+A)*exp(b2*22+b1*3781)-A)-log((R0+A)*exp(b2*22)-A))/(b1*(R0+A))-1000;

S4=g1*(log((R0+A)*exp(b1*50)-A)-log(R0))/(b1*(R0+A))+g2*(log((R0+A)*exp(b1*50+b2*86)-A)-log((R0+A)*exp(b1*50)-A))/(b2*(R0+A))-1000;

S5=g1*(log((R0+A)*exp(b1*150)-A)-log(R0))/(b1*(R0+A))+g2*(log((R0+A)*exp(b1*150+b2*84)-A)-log((R0+A)*exp(b1*84)-A))/(b2*(R0+A))-1000;

S6=g1*(log((R0+A)*exp(b1*200)-A)-log(R0))/(b1*(R0+A))+g2*(log((R0+A)*exp(b1*200+b2*116)-A)-log((R0+A)*exp(b1*116)-A))/(b2*(R0+A))-1000;

result=solve(S1,S2, S3, S4, S5, S6,'A,b1,b2,g1,g2,R0');

mikewg
2009-03-05, 10:09
书写的格式不对,你需要加''的,看下help文件,你面有一些例子

chinapenguin
2009-03-09, 08:49
谢谢啊!我再去试试看!