问题标签 [sympy]

我遇到了一些问题，nsolve很难找到一些给出一些初步猜测的函数的解决方案。然后我想尝试 numpy/scipy 求解器。

这是一个使用 sympy 的程序，并且在提供此解决方案时效果很好：[0.0, -9.05567e-72, 9.42477, 3.14159]

这是同一个程序，但使用 numpy 兼容函数：

但是当我运行程序时：

我尝试给 x0 赋值 0，第二个程序（使用 numpy）工作，给出一个接近 0 的数值，但是从 pi/4 开始，它给出了错误消息。我在 numpy 中错过了什么吗？

I'm wondering if there's a way to treat array variables symbolically (something like sympy but for array instead of numerical variables). So that I can have array expressions such as

or something even more advanced such as nested array relations

By treating A, B symbolically, I can change the contents of A, B by using substitution , e.g.:

I can then add f1 and f2, etc.

如何使用类似上面的结构编写一个函数来返回它的参数的 y 截距？

编辑：我得到我的方程式的参考包含几个错误。我已经在这里修好了。解决方案现在可能真的有意义了！

当两层流体流过地形时，根据流体中流速和波速的相对大小，存在许多不同的解决方案。

这些被称为“超临界”、“亚临界”和“临界”（我在这里将前两个称为“超临界”）。

以下方程定义了 (h, U0) 参数空间中临界和超临界行为之间的边界线：

我想消除d_1c（即我不在乎它是什么）并在(h, U_0).

简化因素：

• 我只需要给定的答案 d_0
• 我不需要精确的解决方案，只需要解决方案曲线的轮廓，因此可以通过解析或数值方式解决。
• 我只想绘制区域 (h, U0) = (0,0) 到 (0.5, 1)。

我想使用 Enthought 发行版中提供的模块（numpy、scipy、sympy）来解决这个问题，但真的不知道从哪里开始。真正让我困惑的是变量 d1c 的消除。

以下是python中的方程式：

我期待一个具有许多解决方案分支的解决方案（并不总是物理的，但不要担心）并且看起来大致如下：

我该如何实施呢？

sympy 中的求解功能需要一些时间才能提出解决方案。是否可以通过降低所需的精度来加快速度（我真的不需要小数点后的 15 位数字！）？

我想编写一个 Boost-Python 程序来从用户那里获取一个符号 python 函数并在我的程序中评估它的导数。

例如，用户提供了一个 python 文件 (Function.py)，它定义了一个函数，如 F = sin(x)*cos(x)。

然后我想使用 Sympy 的符号微分能力访问 F'(x)（F(x) 的导数）。我不想使用数值微分。

有没有办法使用 Boost-Python 在 C++ 中使这样的函数 F'(x) 可访问。

可能重复：
SymPy - 任意数量的符号

我正在尝试创建一个通过迭代形成的方程组。每个方程生成：

summation(binomial(n,i) * (1-t) ** (n-i) * t ** i * P[i]) - V[j]

所有V[j]都作为输入给出。这个想法是使用周围创建的方程组P[i]来解决每个问题，P[i]但要这样做，它们必须是 sympy 符号对象，我不知道如何生成这样的列表。

是否可以有一个无名列表或生成一个符号列表？

如果问题需要澄清，请告诉我，我很难用语言表达。

提前致谢 ：）

现在我收到错误：应用程序中的服务器错误。无法导入已键入的名称 说明：执行当前 Web 请求期间发生未处理的异常。请查看堆栈跟踪以获取有关错误及其源自代码的位置的更多信息。异常详细信息：IronPython.Runtime.Exceptions.ImportException：无法导入名称类型的源错误：第 44 行：expr))；第 45 行：第 46 行：script.Execute(scope); 第 48 行：返回 scope.GetVariable("result");

我有一个带有可变数量参数的python函数：

我想针对每个参数自动生成表示 F 的导数的 N 个函数。

例如，我可以同时给出 F(x1) = sin(x1) 和 F(x1, x2) = sin(x1) * cos(x2) 并自动获得所有导数。

Edit2：如果函数 F 是 2 个变量（固定数量的参数），我可以使用

我对为什么以下代码的答案是 10 而不是 1 感到困惑。有人可以帮助我了解发生了什么lambdify或导致错误答案的原因吗？

为什么不是答案k1或 1，因为x = 0小于loc = 0.5？

但是，如果我这样做，它会返回正确的答案

我需要将其kfun作为 x 和 y 的函数，因为稍后，我将其用作积分参数的一部分。它最终也将取决于 y。

我在 Mac 10.6.x 上使用 python 2.6.8。