SymPy，使用已知模式或子表达式的简化/替换
我有以下表达式: from sympy import pi, sin, cos, var, simplify var('j,u,v,w,vt,wt,a2,t,phi') u0 = v*a2*sin(pi*j/2 + pi*j*t*phi**(-1)/2) + pi*vt*a2*cos(pi*j/2 + pi*j*t*phi**(-1)/2)*j*phi**(-1)/2 + pi*w*a2*
..
symbolic-math相关内容
如何使用sympy.physics.quantum换向器?
我想研究一些换向器操作，并在sympy中找到此工具.它似乎按预期工作(但实际上该文档不存在，或者至少我发现的很少，但是请参阅下面的Dalton Bentley的评论)，但是我遇到了以下问题. from sympy.physics.quantum import Commutator as Comm from sympy.physics.quantum import Operator A = O
..
用于数学字母数字符号的RegEx
我一直在尝试删除这样的粗体字符: 𝐒𝐎𝐍𝐘𝐀 - New song 您看到的是没有标签的粗体文本. 我曾经尝试过，假设角色可能是一种表情符号，但是我没有任何运气. $regex_emoticons = '/[\x{1F600}-\x{1F64F}]/u'; $clear_string = preg_replace($regex_emoticons, '', $
..
SymPy和复数的平方根
当使用solve计算二次方程式的根时，SymPy返回可以简化的表达式，但我无法简化它们.一个最小的示例如下所示: from sympy import * sqrt(-24-70*I) 在这里，SymPy仅返回sqrt(-24-70*I)，而Mathematica或Maple会以5-7*I等效. 我知道有两个平方根，但是这种行为需要例如SymPy从中返回相当复杂的解决方案 z
..
Sympy:删除多项式中的高阶项
使用Sympy，假设我们有一个表达式f，它是符号"x"(以及其他可能的符号)的多项式. 我想知道是否存在一种有效的方法来删除所有大于f的整数n的f中的所有项. 作为一个特例，我有一个非常复杂的函数，但我只想将x的术语保持在2阶以内.有效的方法是什么? 一个明显的，不是很有效的方法是对于每个小于n的m，取m个导数并将x设置为0，以获得x ^ m的系数.我们以这种方式获得每个系数，然
..
千里马:将矩阵转换为列表
我通过以下方式将列表转换为Maxima中的矩阵: DataL : [ [1,2], [2,4], [3,6], [4,8] ]; DataM: apply('matrix,DataL); 怎么做呢?如何将给定的矩阵 DataM 转换为列表 DataL ? 解决方案 我知道比赛已经很晚了，但就其价值而言，有一种更简单的方法. my_matrix : matrix ([a, b
..
有没有一种矢量化的方法来计算sympy中的梯度?
如何计算sympy中多元函数的(符号)梯度? 很明显，我可以分别计算每个变量的导数，但是有矢量化运算可以做到吗? 例如 m=sympy.Matrix(sympy.symbols('a b c d')) 现在，对于i = 0..3，我可以这样做: sympy.diff(np.sum(m*m.T),m[i]) 这将起作用，但是我宁愿做类似的事情: sympy.di
..
处理Sympy积分返回的分段方程
在sympy中，我有一个积分，该积分返回一个逐段对象，例如 In [2]: from sympy.abc import x,y,z In [3]: test = exp(-x**2/z**2) In [4]: itest = integrate(test,(x,0,oo)) In [5]: itest Out[5]: ⎧ ___
..
如何用sympy求解矩阵方程?
在sympy中，给出了矩阵方程 M * x + N * y = 0 (或更复杂..) 如何解决x的问题? (M，N =矩阵，x，y =向量) 我尝试使用普通符号进行此操作，但显然失败了.使用MatrixSymbol效果不佳.有什么方法可以做到吗，或者sympy无法做到? 解决方案 正如MRocklin指出的那样，MatrixExpressions还不支持此功能，但是非可交
..
Clojure的计算机代数
短版: 我对某些Clojure代码感兴趣，这将使我能够指定x的变换(例如，排列，旋转)，函数f(x)的值不变，这样我就可以高效地生成满足r的x序列= f(x). Clojure的计算机代数是否有发展? 对于(一个简单的)示例 (defn #^{:domain #{3 4 7} :range #{0,1,2} :invariance-group :full
..
R中方程组的求解
可以使用Ryacas库在R中以符号方式求解方程式. 例如 library(Ryacas) yacas("Solve(x/(1+x) == a, x)") 给予 expression(list(x == a/(1 - a))) 有人知道如何(象征性地)求解方程组吗? 谢谢. 解决方案 好吧，我使用出色的python库 sympy 进行符号计算. 使用符号，直接
..
C#库，用于代数简化和求解
网络上有很多代数求解器和简化器(例如，algebra.com上的代名词).但是，我正在寻找可以作为更大项目的一部分插入C#的东西(我正在制作自己的计算器，但显然我会征求许可等). 理想情况下，我会使用以下代码: String s = MathLib.Simplify("5x*(500/x^2*(sqrt(3)/4)+1)+2x^2+(sqrt(3)/2)*x^2"); "s"将简化
..
当Mathematica中的r为符号时，如何找到函数的rth导数?
我有一个功能f(t)=2/(2-t).在不使用Mathematica的情况下，很难在t = 0处获得rth导数(即2^(-r)*r!).在Mathematica计算的情况下，当r = 4时，我可以得到第r个导数，如下所示:D[2/(2-t), {t, 4}].但是，当r为ANY整数时，如何在Mathematica中的t = 0处获得rth导数?我尝试使用此表达式，但未按预期工作: Simpl
..
从表达式中提取完整的参数和函数名称(仅限标准数学函数)
对于类似 的表达式 a1 := sin(1+a/b); a2 := log(1+ a*b); a3 := abs(a^2+b); 如何获取各个函数的表达式.例如 someCommand(a1) # should get an output 1+a/b someCommand(a2) # should get an output 1+a*b someCommand(a3) # s
..
加快SymPy中符号行列式的计算
我有一个4x4矩阵A，在每个条目中都有相当长但很简单的符号表达式.涉及大约30个不同的符号. “简单"是指仅使用加/减，乘/除和整数幂来组合这些符号. “长"是指如果我打印出矩阵，它将覆盖三到四个屏幕. 我需要这个矩阵的行列式.或者，更具体地说，我知道行列式是一个特定符号中的四阶多项式，并且我需要该多项式的系数. A.det()在运行数小时后不会终止，因此我需要一种不同的方法.有任何想法吗?
..
函数，SymPy和SciPy矩阵上的数值积分
从我的SymPy输出中，我得到如下所示的矩阵，必须将其集成为2D格式.目前，我正在按元素进行操作，如下所示.此方法有效，但对于我的实际情况(其中A及其功能要大得多(请参见下面的编辑))，它变得太慢(对于sympy.mpmath.quad和scipy.integrate.dblquad而言): from sympy import Matrix, sin, cos import sympy im
..
评估符号功能
我想找到cos(5).为什么此表达式无效: syms x f=sin(x) disp(diff(f)(5)) 错误是 Line: 3 Column: 12 Indexing with parentheses '()' must appear as the last operation of a valid indexing expression. 解决方案 您的错误与符号变量
..
如何从包含特定变量MAPLE的符号表达式中获得部分表达式?
我有一个如下的符号表达式 y1 = (1/a)-(b/a^2)+x*a*b-x/b y2 = a*b+a*x+b*sqrt(x) 现在我需要获取具有特定术语的部分表达式.像 xFunction(y1, x) # should return x*a*b-x/b xFunction(y2,x) # should return a*x+b*sqrt(x) 任何建议或想法都非常有用
..
在sympy中分解polys
我正在做一个非常简单的概率计算，即从A-Z集中获取X，Y，Z的子集(具有相应的概率x，y，z). 由于公式很繁琐，为了处理它们，我正在尝试简化 (或 collect 或 factor -我不知道确切的定义)使用符号表示这些多项式表达式. 因此..(具有从具有相应概率x，y，z的A-Z集合中获取X，Y，Z的子集的非常简单的概率计算表达式) import sympy as sp x
..
使用SymPy将符号表达式转换为Python函数
我有一个相当大的符号函数，该函数针对循环中参数的不同值进行求值.在每次迭代中，在找到函数的表达式之后，将导出偏导数.像这样: from sympy import diff, symbols,exp def lagrange_eqs(a): x,y,z= symbols('x y z') FUNC=x**2-2*x*y**2+z+a*exp(z) d_lgrng_1=
..