求助matlab solve 与fzero结果完全不一致。

解这个方程-87/2/x+1/2*exp(x)*(87+200*x)/x-782/5=0，用fzero命令：fzero(char(fangcheng),0.3)
可以算出正确值 0.1010
可是用solve命令：solve(fangcheng,x)，无论怎么样都得到0.的结果是怎么回事？

方法1:

>> syms x

>> solve(-87/2/x+1/2*exp(x)*(87+200*x)/x-782/5==0, x)
ans =
0.10099792968574978894955689083098
--------------------------

方法2:
fangcheng ='-87/2/x+1/2*exp(x)*(87+200*x)/x-782/5'
>> fzero(fangcheng, 0.3)
ans =
    0.1010
>> solve(fangcheng, 0.3)
ans =
0.10099792968574978894955689083098

这个跟你设置的初值有关系

多设置些初始值试一试。根本原因就是，这两个functions做optimization找local minimum 的方法是不同的，所以找点方式不同，会有这些问题，尤其是ill-conditioned problem。

mike68097 发表于 2014-10-21 10:29
方法1:>> syms x>> solve(-87/2/x+1/2*exp(x)*(87+200*x)/x-782/5==0, x)
ans =
0.100997929685749788949 ...

还是不行啊，用到solve时两个方法都不行，方法一：



方法二：

mike68097 发表于 2014-10-21 10:29
方法1:>> syms x>> solve(-87/2/x+1/2*exp(x)*(87+200*x)/x-782/5==0, x)
ans =
0.100997929685749788949 ...

我用的matlab7.0和win7的系统，是不是matlab坏了？

谢谢各位，问题已解决，这个方程没有解析解，而solve给出的是符号解析解，对于无解析解的方程效果很差。下面是一段百度百科的文字：“在MATLAB中，solve函数主要是用来求解代数方程（多项式方程）的符号解析解。也能解一些简单其他方程的数值解，不过对于解其他方程的能力很弱，此时求出的解往往是不精确或不完整的。注意可能得到的只是部分的结果，并不是全部解。”

不懂你的Matlab的设置是否有问题, 看起来小数点以下都被舍去或没显示, 也许你要试试以下:

format short
fangcheng ='-87/2/x+1/2*exp(x)*(87+200*x)/x-782/5'
solve(fangcheng, 0.3)

mike68097 发表于 2014-10-24 12:38
不懂你的Matlab的设置是否有问题, 看起来小数点以下都被舍去或没显示, 也许你要试试以下:

format short
...

还是一样的0.
你的matlab能得出正确的结果吗？

对了, matlab可以设定求解精度/次数等设定, 你可以试试以下, 看看设定的精度是否足够?  如果还不行, 那大哥你就重新安装吧, 免的忙来忙去一场空!

-------------------------------------
>> options=optimset('fsolve');可以试试去查看options的值, 并于必要时更改:

>> options.TolFun  
ans= 1.0000e-06

>> options.TolX
ans= 1.0000e-06

>>options.MaxIter
ans= 400

>> options.Algorithm
ans = trust-region-dogleg
------------------------------------------
如果还不行, 建议你换套matlab试验看看 https://bbs.pinggu.org/thread-3155044-1-1.html (matlab8.3).
这结果应该就会跟之前讲的一样了.

>> solve(fangcheng, 0.3)
ans =
0.10099792968574978894955689083098