MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   MATLAB论坛 (https://www.labfans.com/bbs/forumdisplay.php?f=6)
-   -   [求助]计算问题 (https://www.labfans.com/bbs/showthread.php?t=3323)

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

[QUOTE=fanxing39;11522]你用 help fsolve 看下 这个函数的用法,你好像用的不对哦[/QUOTE]

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 解一些简单的方程还行,难点的就没办法了
你的这个问题是 多元非线性方程组 ,你找下其他的求解命令 再试下!


所有时间均为北京时间。现在的时间是 11:05

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.