问题标签 [symbolic-computation]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
2 回答
2462 浏览

r - R中矩阵的符号逆

如何在 R 中找到矩阵的符号逆;例如:

我知道有一个名为“Ryacas”的包,它是“yacas”的接口,但我无法应用它来进行此类计算。'yacas' 是一个用于对数学表达式进行符号操作的程序。请参阅链接了解更多详情。

谢谢你

0 投票
2 回答
494 浏览

f# - C# 库重载 ^ 运算符。如何改用**?

Symbolism库重载了算术运算符。虽然它是用 C# 编写的,但我可以从 F# 中使用它:

输出:

但是,它也使权力超载^。这当然不适用于 F#。

作为解决方法的一步,我为权力导出了一个方法组:

输出:

如何重载**以使用Aux.Pow方法组?

我可以做这样的事情:

这确实适用于MathObject价值观:

但是Aux.Pow也超载int

欢迎任何建议!

0 投票
1 回答
4699 浏览

python - 处理 sympy 积分返回的分段方程

在 sympy 中,我有一个返回 Piecewise 对象的积分,例如

我想只提取这个分段方程的第一个分支,换句话说,我希望能够做一些itest.parts(0)简单的事情sqrt(pi)*z/2。我似乎找不到任何方法来做到这一点,但也许我在文档中使用了错误的搜索词。有任何想法吗?

编辑

戳了一下,我设法发现如果我这样做,itest.args[0][0]我可以提取这个表达式。然而,这似乎有点像黑客。有更好的方法吗?

0 投票
1 回答
474 浏览

matlab - subs(sinc(K), K, 0) 其中 K 是符号将返回 NaN。不应该是1吗?

MATLABsinc(0)将按应有的方式返回 1。但是sinc(K)一些K用零值代替的符号将返回NaN.

以下代码说明了上述内容:

我可以sinc在符号情况下强制返回 1(事先不知道 的值K)吗?

0 投票
0 回答
809 浏览

python - 大三角函数的积分

在 SymPy 中,集成如下所示的大三角函数的推荐方法是什么?数值积分给出-27151.225.

此类函数在此处详述的问题中生成。

0 投票
1 回答
159 浏览

python - Sympy,是否可以强制一个类留在表达式的左侧或右侧

这个答案中,解释了如何避免args在 SymPy 类中进行排序,例如MulAdd等等。

对于像这里解释的这样一个新创建的类,它可以在乘以 a 时转到右侧sympy.core.numbers.Floatsympy.core.numbers.Pi或者sympy.core.numbers.Integer,例如,给出:

原始表达式0.在评估时给出,而新表达式给出D(x).

为了让这个微分算子正常工作,它必须留在左边:

是否有任何隐藏参数(例如_op_priority更改__mul__()优先级)告诉 SymPy 更靠左或靠右的类型?

0 投票
0 回答
1037 浏览

matlab - 符号计算在matlab中不起作用

我很惊讶为什么遵循符号表达式会出错

我检查了 syms 命令是否有问题,但它工作正常

工作没有任何问题,但我认为这条线有问题

请帮助 mt 修复它,也请查看此站点

0 投票
4 回答
686 浏览

python - 使用先前的术语 sympy 生成序列

我想生成一个序列,以便先前生成的元素包含在下一个元素中,我不确定如何执行此操作。

即生成列表,使其项目为:

其中 x 只是一个 Sympy 符号

[x,(x)*(x+1),(x)*(x+1)*(x+2)]

而不是[x,x+1,x+2]

我在想类似的东西

k.append(k*(K+o))

但我不断收到类型错误

非常感谢任何帮助!

0 投票
2 回答
601 浏览

python - 通过 SAGE 积分 cos(x)*cos(2x)*...*cos(mx)

我要找到 $I_m=\int_0^{2\pi} \prod_{k=1}^m cos(kx){}dx$,其中 $m=1,2,3\ldots$

简单的 SAGE 代码:

输出:

如您所见,数字答案是正确的,但是对于 $m=1,2,\ldots,7$ 而言,integrate(...) 的结果是正确的,然后存在一些错误。我们可以打印不定积分:

和输出:

所以对于 $m=7$ 的答案是正确的,通过 WolframAlpha 与 Indef_I_7比较 ,对于 $m=8$ 的答案是不正确的Indef_I_8 via WolframAlpha 应该有 Indef_I_8=$\frac{7x}{128}+\ldots$ 并且没有 $\ sin(x)$, $\sin(3x)$ 的总和,只有 $\sin(2k)$ for $k=1,2,3,\ldots 18$

对不起体积计算!问题是 - 这是符号集成中的错误,我说得对吗?

0 投票
2 回答
295 浏览

matlab - MATLAB符号系数问题

我正在使用 MATLAB 的符号工具箱,但在提取导数系数时遇到了一些问题。也许 MATLAB 不能做我正在寻找的东西。无论如何,重现我遇到的问题的代码如下所示:

上述代码部分中的注释突出显示了哪些有效,哪些不按预期工作。具体来说,输出是:

有谁知道为什么会发生这种情况或我是否做错了什么?