问题标签 [symbolic-computation]

我正在尝试解决许多变量的符号方程组。 方程的结构会根据 AM 的条目而改变，因此随便求解这些方程不会很好。这段代码符合我的预期，但是我需要解决大约 20,000 个案例，所以它太慢了。有什么方法可以加快速度（可能使用矢量化？）。

编辑：为澄清起见，将来我计划将 AM 的条目作为参数。这些条目会有所不同，但只有 1 或 0 个值。

我正在尝试使用 sympy 对以下形式的广义有理函数进行拉普拉斯逆变换：D/(A*s**2 + B*s + C)。

从数学分析中知道答案是：

但是 sympy 不会产生解决方案，并且代码将卡在无限 CPU 进程中，没有特定的增益。但是像这样放置eq：

sympy 将产生如下形式的方程：

这不是我需要的。关于如何使 sympy 屈服于这样一个人类可读的答案的任何建议？

我正在研究将 GF(2^N) 转换为 GF(2) 的多项式系统。在转换过程中，我需要进行大量模 2 计算，例如矩阵求和和乘法。

那么有没有办法在 Sympy 中使用模 2 算术对矩阵进行求和和乘法？

以下两段最大值代码应该是等价的：

和：

在第一种情况下，maxima 知道将表达式计算为 8。在第二种情况下，它不会简化并且给我留下 x[i] 的符号求和，因为 i = 1 到 2。

为什么 maxima 无法将第二个表达式识别为与第一个相同？

没有对答案进行硬编码（如第一段代码），在评估稍后指定上限的总和时，如何获得最大值以完全评估索引变量的总和？

我在适合 Ryacas 符号计算的类中mat1定义了两个矩阵：mat2Sym

和

现在的问题是如何通过 Ryacas 包对这两个矩阵进行符号乘法？

先感谢您。

我正在制作一个程序，该程序接受表达式的输入并将其集成并将结果打印给用户。我正在使用 SymbolicC++，但似乎无法弄清楚如何输入符号......我该怎么做 - 我正在为我的 IDE 使用 Visual Studio

我在 Sage 中使用以下地图：

对于x=0和f(x)=1对于x!=0 ，其计算结果为f(x) = 0 ；

在符号结果中，sgn(x)^2、sgn(x)^4和sgn(x)^8等被视为不相等，即使它们对于x的所有值都是相等的。有没有办法可以替代类似的东西：

对于这些关系的所有出现，以及x的所有符号值？

我可以为每个符号创建一个新的替换规则，例如

等等，但这似乎很难控制。

我的目标是使用 编写多维泰勒近似sympy，其中

• 使用尽可能多的内置代码，
• 计算两个变量的给定函数的截断泰勒近似
• 返回没有Big-O-remainder 项的结果，例如在sin(x)=x - x**3/6 + O(x**4).

这是我到目前为止所尝试的：

方法一

天真地，一个人可以series为每个变量组合两次命令，不幸的是，这不起作用，因为这个例子显示了函数sin(x*cos(y))：

方法二

基于这篇文章，我首先写了一个 1D 泰勒近似：

用一维示例检查它工作正常

但是，如果我知道在xand中进行两个扩展的链式调用y，则 sympy 挂断

有人知道如何解决这个问题或以另一种方式实现它吗？

http://www.doc.ic.ac.uk/~dsl11/klee-doxygen/overview.html

我正在尝试为二进制文件建模内存，所以我想看看 KLEE 如何处理它。但是，我无法理解概述中的解释。

我可以得到一个更简单的解释，并且可能是一个简单的例子来说明如何做到这一点？

我尝试使用这些命令

和

，但它们都产生

f(x) = 1

（是的..对于每个实数 x 包括 0）

问题是如何避免在命令“sym”中进行预评估，或者存在另一种处理此问题的方法。

非常感谢你！

2014 年 5 月 21 日更新：

让我稍微描述一下这个问题。

考虑

f(x) = x/x

和

g(x) = 1

很明显，域的f和分别g是R-{0}和R。

自动简化sym/syms可能会导致丢失一些信息。