问题标签 [symbolic-computation]

Please consider this example. I would like to solve x^3 - 2x > 0. I try the following commands:

and Matlab returns this

ans = solve([0.0 < x^3 - 2.0*x], [x])

which is not what I expect. Therefore I use

which returns

ans = Dom::Interval(2^(1/2), Inf) Dom::Interval(-2^(1/2), 0)

Can someone explain that why solve works successfully only in the second case?

以下命令

和

返回结果

求解([0.0 < (8.0*exp(-1.0/x^2))/x^4 - (2.0*exp(-1.0/x^2))/x^2 - (6.0*log(x^2) *exp(-1.0/x^2))/x^4 + (4.0*log(x^2)*exp(-1.0/x^2))/x^6], [x == RD_NINF..RD_INF ])

这不是我所期望的。

第一个问题：是否可以强制 Matlabsolve返回所有解决方案的集合？

（这与这个问题有关。）此外，当我尝试求解方程时

返回

答案=
-1.5056100417680902125994180096313

我不满意，因为没有返回所有解决方案（它们大约是从 WolframAlpha 获得的 -1.5056、1.5056、-0.5663 和 0.5663）。

我知道vpasolve通过一些初步的猜测可以处理这个问题。但是，我不知道我通常如何找到初始猜测值来获得所有解决方案，这是我的第二个问题。

欢迎其他解决这些问题的解决方案或建议。

考虑

命令

返回

我尝试

哪个返回

和

分别。

这真的是一个错误numeric::solve吗？

或者我只是错过了命令中的一些要点。

欢迎所有意见和建议。谢谢你。

我有以下无限系列需要评估。目前，根据输入参数（kappa、sigmav、thetav），这个函数需要很长时间来计算。有没有可能加快速度？耗时的表达式是 symsum(..., 0, Inf)

PS。这是非中心卡方分布的 pdf ，有一些小的修改，因此我不能使用ncx2pdf. 我已经尝试过https://stackoverflow.com/a/15966126/321749这只会导致无穷大。

代码是

这产生了一个情节

对于 [-2,2] 中的每个 x，v(x) 不应该为零吗？

考虑使用 定义的矩阵Symbolic Toolbox：

我们获得：

我们怎样才能只提取这个矩阵的第二行？在这种情况下，索引方案是否不同（f(X,Y)(2,:)）？

考虑以下类定义：

类上定义的操作行为如下：

根据标题，我的问题是：为什么 Python 运行时将运算符的输出显式强制为布尔类型，而不是其他预期返回布尔值（ 、、、等）in的比较运算符的输出？==<=>

编辑：文档分别对in运算符和丰富的比较进行了以下说明：

• 对于定义__contains__()方法的用户定义类，x in y当且仅当y.__contains__(x)为真时为真。

• 按照惯例，成功比较返回 False 和 True。但是，这些方法可以返回任何值，因此如果在布尔上下文中使用比较运算符（例如，在 if 语句的条件中），Python 将对值调用 bool() 来确定结果是真还是假.

问题仍然存在：为什么需要区分？

我有这两个方程，其中x, y和z是变量，p1, p2和p3是参数。MATLAB 可以找到x y z基于的参数/符号解决方案p1, p2, p3吗？

• 2(x−p1)+2(xy−p3)y = 0
• 2(y−p2)+2(xy−p3)x = 0

我将以下非线性系统传递给 Matlab：

并使用 syms 象征性地找到 x 和 y 的解决方案，但我得到了：

和

我不明白 z 来自哪里？什么是z？

必须存在一种良好的惯用方式来表达 Haskell 中类型级别的一般计算。我能想到的就是这个（非法的）OO 模仿。

理想情况下，我想保持开放，所以这种方法对我没有吸引力。我要求太多了吗？