sympy相关内容

如何设置 sympy 以执行标准微分几何任务?

我是一名工程专业的学生.我需要做的几乎所有数学运算都是在 R2 或 R3 中进行的，并且涉及微分几何.当然，我真的很喜欢 sympy，因为它使我的计算可重用且易于展示. 我发现了什么: sympy 中的东西与我所知道的功能非常接近，即标量或向量值到标量或向量值的映射，具有名称并连接到表达式似乎是某种形式 functionname=sympy.Lambda(元组中的变量，表达式) 或 ..
发布时间：2021-09-03 20:15:14 其他开发

估计 SymPy 函数的进度

我发出了一个 sympy 命令来求解某个方程或其他方程.现在已经过去几天了，我不知道什么时候才能完成. 我可以让 sympy 记录对 .solvers.solve 的调用进度吗?如果不是，我如何估计 sympy 求解一组方程所需的最坏情况时间? 例如 import sympy, sympy.solvers从 sympy 导入 sqrta,c,d,e,f,x = tuple(map(s ..
发布时间：2021-09-03 20:15:11 Python

刻意简化分数指数

我有一个涉及分数指数的表达式，我想将其转化为可识别的多项式，以便 sympy 求解.如有必要，我可以使用 Rational 编写指数，但无法使其工作. 我能做什么? >>>从 sympy 导入 *>>>var('d x')(d, x)>>>(0.125567*(d + 0.04) - d**2.25*(2.51327*d + 6.72929)).subs(d,x**4)0.125567*x ..
发布时间：2021-09-03 20:15:05 其他开发

使所有符号在 sympy 表达式中可交换

假设你在一个 sympy 表达式中有许多非交换符号，比如 a, c = sympy.symbols('a c', 可交换=假)b = sympy.Symbol('b')expr = a * c + b * c 使表达式中的所有符号都可交换的首选方法是什么，例如，sympy.simplify(allcommutative(expr)) = c * (a + b)?> 在 this answer ..
发布时间：2021-09-03 20:15:02 Python

Sympy autowrap (cython):定义“助手"对于 sympy.Max, sympy.Heaviside

我有一个 sympy.Matrix，称为 J_sym，我想自动换行(最好使用 cython 后端)；相应的符号存储在列表 list_args 中. 然而，我遇到的问题是显然不支持某些 sympy 函数，在我的例子中特别是 sympy.Max 和 sympy.Heaviside. 具体来说，如果我尝试 J_num_autowrap = autowrap(J_sym, backend=" ..
发布时间：2021-09-03 20:14:58 其他开发

sympy - 有没有办法区分抽象变量?

我现在正在学习 SymPy，我想知道是否有一种方法可以通过一般形式的一个变量来区分函数. 考虑这个例子: 有一个向量，让我们把它的分量写成 x_1、x_2、x_3.这种向量的长度 r 将是 r = sqrt(x_1^2 + x_2^2 + x_3^2). 我想在 x_i 上区分这个向量，我没有指定 i，它应该给出类似 x_i/r^3 的东西. 可以在 SymPy 中做到这一 ..
发布时间：2021-09-03 20:14:55 Python

Sympy:使用符号表达式作为数值被积函数

我需要以符号方式操作一个函数，然后对该函数进行数值积分.如何在被积函数中正确使用我的表达式 f.如果这甚至是明智的做法，我该如何正确使用lambdify?非常感谢. from sympy import *导入 scipy.integrate 作为集成r = symbols('r') #定义符号f = diff(r*r) #进行符号操作def integrand(x): #定义要积分的函数retu ..
发布时间：2021-09-03 20:14:52 Python

如何使用 sympy.lambdify 和 Max 函数来替换 numpy.maximum 而不是 numpy.amax?

我正在尝试使用 sp.Max(x, 0) 对大型解析表达式进行 lambdify.我想使用 numpy 来向量化我的计算，所以 x 将是一个数组.我需要 x 和 0 的元素最大值.不过，默认情况下，sympy 将 sp.Max 更改为 np.amax.它沿轴找到最大值，这不是我需要的.lambdify 中的“modules"关键字不像我预期的那样工作.我试过了: 将 numpy 导入为 np将 ..
发布时间：2021-09-03 20:14:49 Python

sympy 将分数与变量分开

使用 sympy 如何将分数与变量分开 Mul(Fraction(3,5), Pow(K, Integer(2)))23·K────5到3 2─ K5 我知道这个简化版还不错，但是当我有很大的方程时，它会变得很乱 解决方案 我对漂亮打印或 LaTeX 打印不是很熟悉，但我设法想出了一些东西.将 UnevaluatedExpr 放在 Mul 的每个参数中: from sympy imp ..
发布时间：2021-09-03 20:14:47 其他开发

打印没有数值计算的替换表达式

是否可以在不计算的情况下打印替代的 SymPy 表达式?我想打印替换的快递和结果. 例如 x = Symbol('x')expr = x**2pprint(expr) # 打印表达式结果 = expr.subs({x:2})print(result) # 这个打印结果 4 如何打印“中间结果"，表达式2**2? 解决方案 你可以为此传入UnevaluatedExpr，如下图: ..
发布时间：2021-09-03 20:14:44 Python

使用 Sympy 集成到 Python 中

我目前正在使用 Sympy 来帮助我进行数学计算.现在，我正在尝试执行数值积分，但每次运行脚本时都会出现错误.这是脚本: from sympy import *cst = {'qe':1.60217646*10**-19,'m0':N(1.25663706*10**-6)}d = 3.6*10**-2l = 20.3*10**-2n = 217.0我 = 10.2# 循环圈数圈 = l/n;# ..
发布时间：2021-09-03 20:14:41 Python

sympy中的链式法则

我正在尝试使用 SymPy 进行一些基本练习.我想对直角坐标中的函数相对于极坐标中的半径参数进行象征性的二阶导数. 我想要一个很好的链式规则符号表达式，它可以计算可以计算的内容，而未评估不能进一步简化的内容. from sympy import *init_session()x, y, r, t = symbols('x y r t') # r (半径), t (角度θ)f, g = 符号 ..
发布时间：2021-09-03 20:14:38 Python

在 sympy.stats 中绘制混合分布

(此问题的要点 此处) 我想创建两个 Gamma 分布的混合，并绘制在给定范围内评估的结果. sympy.stats 似乎能够做到这一点，因为它能够计算混合物的期望值并从中提取样本.我对 sympy 很陌生，所以不确定在这种情况下是否有比我一直在使用的方法更喜欢的评估和绘图方法. %matplotlib 内联从 matplotlib 导入 pyplot 作为 plt从 sympy.s ..
发布时间：2021-09-03 20:14:35 Python

将每个元素映射到表达式

我有这个代码: import sympy将 numpy 导入为 np从 sympy.utilities.lambdify 导入 Lambdify从集合导入 OrderedDictarr = [np.array([ 1, 2]), np.array([ 5, 6])]a,b = sympy.symbols('a b')var = [a,b]expr = ['a+cos(b)', 'a+cos(b) ..
发布时间：2021-09-03 20:14:33 Python

如何求解形式幂级数中的微分方程?

我想先得到几个由微分方程隐式定义的形式幂级数的系数. 示例. import sympy as spsp.init_printing() # 数学如乳胶from IPython.display 导入显示z = sp.Symbol('z')F = sp.Function('F')(z)F_ = sp.Derivative(F, z)方程 = sp.Eq(F**2 + 1 - F_, 0)显示( ..
发布时间：2021-09-03 20:14:30 其他开发

负实数的立方根

我正在尝试绘制一个非常复杂的函数，即 log(x/(x-2))**Rational(1,3).我只使用实数.如果我尝试绘制它，sympy 只会绘制它的 x>2 部分. 我发现实际上复数开始起作用，例如，root(-8,3).n() 给出: 1.0+1.73205080756888i 这是合理的，即使这不是我想要的(因为我只对真实结果感兴趣). 阅读sympy › 原理根我发 ..
发布时间：2021-09-03 20:14:27 Python

Sympy:如何将乘积的对数简化为对数之和?

为了激发这个问题，sympy.concrete 有一些有效的工具来处理符号和.为了将这些工具应用于符号产品，必须取对数.但是，直接取对数并不会自动给出转换: import sympy as spsp.init_printing() # 将数学显示为乳胶z = sp.Symbol('z')j,k = sp.symbols('j,k')Prod = sp.Product( (z + sp.sqrt( ..
发布时间：2021-09-03 20:14:24 Python

Sympy 自定义函数

我注意到，没有像 sech(x) 和 csch(x) 这样的函数. 有什么办法可以快速的分别定义为 1/cosh(x) 和 1/sinh(x) 吗? 另外，我怎样才能使 sympy 将 arccos 视为 acos? 我正在使用解析器，因此正确解析了“cos(pi/2)".但我希望以类似的方式解析“arccos(pi/2)". 解决方案 您可以定义实现所需行为的函数. ..
发布时间：2021-09-03 20:14:21 Python

SymPy 相当于 mathematica 中的 HoldForm

在 Mathematica 中，可以防止系统在输入表达式时对其进行简化.语法如下: HoldForm[x/x] 是否可以用 SymPy 做类似的事情? 解决方案 以下方法可以达到类似的效果.可能还有其他可用的我不知道. import sympy as spx = sp.symbols('x')expr1 = x/xexpr2 = sp.S('x/x',evaluate=False)e ..
发布时间：2021-09-03 20:14:18 其他开发