登录论坛

查看完整版本 : [求助]计算问题


sohu820528
2008-05-28, 23:46
syms x y z
eq1='1/x+log(10000-z)-log(y)-((10000-z)/y)^x*log((10000-z)/y-((10000-z)/y)^x*log((10000-z)/y))=0';
eq2='-x/y+(x/y)*((10000-z)/y)^x+(x-y)*((10000-z)/y)^x=0';
eq3='-(x-1)/(10000-z)+(x/y)*((10000-z)/y)^(x-1)+(x/y)*((10000-z)/y)^(x-1)=0';
fsolve(eq1,eq2,eq3)
x=ans.x
y=ans.y
z=ans.z
以上运算显示如下:
Warning: Warning, solutions may have been lost
Warning: Explicit solution could not be found.
请高手指教

fanxing39
2008-05-29, 01:13
你用 help fsolve 看下 这个函数的用法,你好像用的不对哦

sohu820528
2008-05-29, 09:46
你用 help fsolve 看下 这个函数的用法,你好像用的不对哦

syms x y z
eq1='1/x+log(10000-z)-log(y)-((10000-z)/y)^x*log((10000-z)/y-((10000-z)/y)^x*log((10000-z)/y))=0';
eq2='-x/y+(x/y)*((10000-z)/y)^x+(x-y)*((10000-z)/y)^x=0';
eq3='-(x-1)/(10000-z)+(x/y)*((10000-z)/y)^(x-1)+(x/y)*((10000-z)/y)^(x-1)=0';
solve(eq1,eq2,eq3)
x=ans.x
y=ans.y
z=ans.z
以上运算显示如下:
Warning: Warning, solutions may have been lost
Warning: Explicit solution could not be found.
请高手指教
是solve不是fsolve 写错了

fanxing39
2008-05-29, 13:16
报错不是程序错的意思,是用solve 无法求解出这道题目!
solve 解一些简单的方程还行,难点的就没办法了
你的这个问题是 多元非线性方程组 ,你找下其他的求解命令 再试下!