问题标签 [lambdify]

我在这里用solve() 代码 成功求解了一个未知变量：

现在我试图从 Pandas 数据框中获得相同的结果，但没有成功。

我收到错误 - NameError: name 'y' is not defined
这很奇怪，因为 'y' 是我正在寻找的。

错误在于lambdify 函数。由于某种原因，它说语法错误。这就是确切地说：

文件“D:\Anaconda\lib\site-packages\sympy\utilities\lambdify.py”，第 434 行，在 lambdify func = eval(lstr, namespace) 中。

程序在用户输入容差后立即崩溃。这是一个牛顿法程序（或者至少是我的尝试）。任何有关如何改进此代码的建议也将不胜感激。

我输入f的dfdx内容如下：

和

分别。至于v0and eps，我分别输入1and 0.000001。无论我输入什么，程序都会以上述错误终止（假设它使用 x 作为变量，否则它会显示未定义，仅此而已）。

我正在用invert_realfrom反转一个函数sympy.solvers.solveset，因为由于某种原因既不solve也solveset不能这样做。结果是Intersection并且似乎是正确的。

我现在想用它来进行数值计算。当我尝试它.subs时效果很好，但是当我尝试时lambdify我得到

NameError：名称“交叉点”未定义

我已经发现我需要Intersection在modules=. lambdify但我找不到工作来源。

这是一个代码示例：

我正在对一个 sympy 分段函数进行lambdifying 尝试做这样的事情：

尽管

我明白了

看起来，（lambdified ？） - 在返回具有真实条件的表达式之前，分段计算所有表达式。有没有解决的办法？

这段代码可以得到 y 的导数。这很简单dx = 2 * x

现在我想获得 for 的dx值x = 2。

显然，dx = 2 * 2 = 4当x = 2

但是我怎样才能用 python 代码实现这一点呢？

谢谢你的帮助！

I am investigating the suitability of SymPy for some of my projects and have come across an issue with the interaction between lambdify and IndexedBase.

In short, my applications make heavy use of functions that employ doubled summed array structures. I need to be able to compute both the function and the first through 3rd derivative of the function with respect to an array element.

My questions are thus:

1. How do I make lambdify work with diff?
2. How do I make a lambdify function where I can specify the index with respect to which I want the derivative?
3. How do I extend the above to a second order derivative with a different index (i.e. second order with respect to index i and j)?

Simplified example:

Execution of this code yields the error:

I have a loss function L(u1,...,un) that takes the form

where

I want to numerically minimize L with scipy, hence I need to lambdify the expression.

However at this point in time lambdify does not natively support translating integrals. With some tricks one can get it to work with single parametric integrals, see Lambdify A Parametric Integral. However I don't see how the proposed solution could possibly be extended to this generalised case.

One idea that in principle should work is the following:

Take the computational graph defining L. Recursively, starting from the leaves, replace each symbolic operation with the corresponding numerical operation, expressed as a lambda function. However this would lead to an immense nesting of lambda function, which I suspect has a very bad influence on performance.

Ideally we would want to arrive at the same result as one would by hand crafting:

MWE: (using code from old thread)

我试图单独计算泰勒展开的元素，但没有得到我应该得到的结果。逼近的函数是 x**321，围绕 x=1 的泰勒展开式的前三个元素应该是： 1 + 321(x-1) + 51360(x-1)**2 出于某种原因，代码与第二个术语相关联的不起作用。请参阅下面的代码。

我获得运行此代码的打印是

我也使用了 .evalf 和 .lambdify 但结果是一样的。我不明白错误来自哪里。

我显然只是从语言开始，所以谢谢你的帮助。

我正在编写一个简短的程序，它应该找到函数的实部和虚部都为零的值。我不明白为什么我在运行程序后得到“无法将表达式转换为浮点数”。（请原谅我写代码时的乱七八糟！）我切断了符号a11-a88的定义，以节省您的阅读时间，但它们都是A类型cmath.exp(b x), A1*cmath.exp(1.0j *b1*x) 或 1.0j*A2*cmath.exp(1.0j*b2*x)。我一直使用 cmath 函数而不是数学（cmath.exp 不是 exp，cmath.sqrt 不是 sqrt）。

编辑：错误的完整追溯

我正在python中实现数据处理流程。我正在尝试尽可能多地使用符号计算 (sympy和numpy) 以获得与代码一致的清晰文档。因此，当我尝试获取点积并将其用于真实矩阵时（通过lambdify），我得到了其他东西：

给

链接到代码的实时版本

有没有人遇到过类似的问题？是否有已知的解决方法？

先感谢您。