一个并不复杂的方程，solve为何死循环？

翻身仗 · 发表于 2019-9-3 10:38:43

您需要登录才可以下载或查看，没有帐号？立即注册

x

代码如下：
ta = 4.05756866;
syms x
fc = exp( -( (-log(x))^ta )^( 1/ta ) ) == 0.999;
sf = solve(fc,x)

点击运行之后，一直是运行状态，没有结果，死循环，什么原因呢？

言深深 · 发表于 2019-9-3 14:02:25

你这个函数是不是跟ta没关系啊？
对函数fc两边取自然对数，得到的就是log(x)=log(0.999)？check一下你的函数，手动推导一下看看

翻身仗 · 发表于 2019-9-3 22:07:03

言深深发表于 2019-9-3 14:02
你这个函数是不是跟ta没关系啊？
对函数fc两边取自然对数，得到的就是log(x)=log(0.999)？check一下你的函 ...

您说的对，推倒后就是您给出的那个方程。
这只是原方程参数取极值时的特殊情况，但就是这种最简单的情况，solve 也会陷入死循环

言深深 · 发表于 2019-9-4 07:15:36

翻身仗发表于 2019-9-3 22:07
您说的对，推倒后就是您给出的那个方程。
这只是原方程参数取极值时的特殊情况，但就是这种最简单的情况 ...

既然与x无关，那么求解的时候是不是自然有点问题？再者你是如何确定解方程陷入的是循环？逐步输出x看看，应该是迭代的。

翻身仗 · 发表于 2019-9-4 08:18:32

言深深发表于 2019-9-4 07:15
既然与x无关，那么求解的时候是不是自然有点问题？再者你是如何确定解方程陷入的是循环？逐步输出x看看， ...

解方程最长运行时间超过半小时，我手动停止了，
确实不能百分百确定是死循环，但是那么久的运行时间，现实中不方便操作

言深深 · 发表于 2019-9-4 12:00:10

翻身仗发表于 2019-9-4 08:18
解方程最长运行时间超过半小时，我手动停止了，
确实不能百分百确定是死循环，但是那么久的运行时间，现 ...

这种方程，可以手动编写程序去计算

翻身仗 · 发表于 2019-9-4 21:55:51

言深深发表于 2019-9-4 12:00
这种方程，可以手动编写程序去计算

感谢您的指教

		自动登录	找回密码
密码			立即注册

[源程序] 一个并不复杂的方程，solve为何死循环？