sympy相关内容

如何使用 Matplotlib 处理渐近线/不连续性

在绘制具有不连续性/渐近线/奇异性/任何内容的图形时,是否有任何自动方法可以防止 Matplotlib 在“中断"处“加入点"?(请参阅下面的代码/图片). 我读到 Sage 有一个看起来不错的 [detect_poles] 工具,但我真的希望它与 Matplotlib 一起使用. 将 matplotlib.pyplot 导入为 plt将 numpy 导入为 np从 sympy 导入 symp ..
发布时间:2021-12-08 10:00:09 Python

如何使用 Python 求解一对非线性方程?

使用 Python 求解一对非线性方程的(最佳)方法是什么.(Numpy、Scipy 或 Sympy) 例如: x+y^2 = 4 e^x+ xy = 3 解决上述问题的代码片段会很棒 解决方案 数值求解,可以使用fsolve: http://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.f ..
发布时间:2021-12-08 09:09:52 Python

快速评估大量输入值的数学表达式(函数)

以下问题 计算字符串中的数学表达式 Python 中的方程解析 解析用户提供的数学公式的安全方法在 Python 中 根据 Python 中不安全的用户输入评估数学方程 和他们各自的答案让我思考如何解析单个数学表达式(一般来说,按照这个答案https://stackoverflow.com/a/594294/1672565) 由(或多或少受信任的)用户有效提供,用于来 ..
发布时间:2021-11-26 17:40:13 Python

错误:function() 至少需要 n 个参数(n 给定)

我正在尝试使用 SymPy 取余数,在这种情况下是余切函数.我有一个集成()函数: import sympy as sy将 numpy 导入为 npdef 积分(f,z,伽玛,t,下,上,精确=真):'''沿等高线 gamma(t) 对 f(z) 进行积分:[lower, upper] -->C输入:f - SymPy 表达式.应该代表一个从 C 到 C 的函数.z - SymPy 符号.应该是 ..
发布时间:2021-11-17 05:16:15 Python

求解两个方程时出现 KeyError

1st 我需要得到两个最长线长度的两个方程 我将长度与 eq 放在列表中,如下所示 [( length 1 , eq 1 ) ,.....] 反向排序列表 求两条最长直线的两个方程 运行以下代码时 from sympy import *导入数学类形状():def __init__(self,T):自我.T=Tself.Ax = self.T[0][0]self.Ay = self ..
发布时间:2021-09-03 20:17:48 Python

为 `x` 赋值不计算总和

我有以下问题: x=Symbol('x')书={1:2*x,2:3*x}x=2打印(书)>>>{1:2*x,2:3*x} 我曾希望它会打印 {1:4,2:6} 但是如果我设置 book={1:2*x,2:3*x}就在打印语句之前,我得到了想要的结果. 令人沮丧的是,如果我改为编写 book=book,这应该是相同的(对吗?),就在打印语句之前,我得到 {1:2*x,2:3*x} - ..
发布时间:2021-09-03 20:17:45 Python

使用 SymPy 包打印产品符号 (∏)

我想知道如何使用 Python 中的 SymPy 包打印这张图片中的符号. 解决方案 符号积是用 Product(expr, (index, low, high)) 创建的,其中 high与大多数 Python 代码不同,它包容.例如: P00 = IndexedBase("P^{00}")P02 = IndexedBase("P^{02}")i = Idx("i")t = 符号(“t" ..
发布时间:2021-09-03 20:17:39 Python

sympy 元素明智的向量划分

我有一个价格向量和一个数量向量.我想在 sympy 中象征性地划分这两个向量.我什至找不到矢量设施.所以我完全失去了对如何创建矢量符号以及如何操作它们的同情. 最终我想使用拉格朗日方法最大化的是 prod((x/p) **b) 受约束的 sum(x) = 1.我可以用标量做到这一点,但不能用矢量: import sympy as sp进口同情从 sympy.abc 导入,p1,p2,l,x ..
发布时间:2021-09-03 20:17:37 Python

Sympy 张量包中的 LeviCivita 张量

我想使用 LeviCivita 张量定义一个带有 sympy 张量模块的矢量积.我没有找到任何方法来做到这一点.有什么线索吗? 解决方案 参见 SymPy 文档.使用示例显示在该页面上. ..
发布时间:2021-09-03 20:17:34 Python

将方程改写为多项式

from sympy import *K, T, s = 符号('K T s')G = K/(1+s*T)方程式1 =方程式(G+1,0) 我想用sympy作为多项式重写方程Eq1:1+K+T*s==0 我该怎么做? 我花了几个小时搜索和尝试简化方法,但找不到优雅、简短的解决方案. SymPy 中的实际问题: from IPython.display 导入显示将 sympy ..
发布时间:2021-09-03 20:17:31 Python

为什么 sympy 不接受我的符号?

我的符号定义如下: import sympy类构造(对象):def __init__(self, value):self.value = 价值def __add__(self, 符号):返回构造(self.value+symbol.value)def __mul__(self,symbol):返回构造(self.value*symbol.value)a = 构造(2)b = 构造(10) (上 ..
发布时间:2021-09-03 20:17:28 Python

将 xy 转换为 x*y?

我正在开发一个项目,其中用户以 xy 或 2x 或 2ab 的形式输入输入,但 sympy 无法明白了,我想把它们分别转换成 x * y, 2 * x, 2 * a * b 解决方案 sympy 实际上能够理解这一点.您需要使用 parsing 来帮助它: from sympy.parsing.sympy_parser import parse_expr从 sympy.parsing.s ..
发布时间:2021-09-03 20:17:25 Python

如何使用 sympy 在 Python 中求解联立方程

我有以下两个等式: 40 = Vmax*5.041667 + (Vmax**2/Amax)*sympy.exp((-Amax/Vmax)*5.041667) - (Vmax**2/Amax)20 = Vmax*2.897628 + (Vmax**2/Amax)*sympy.exp((-Amax/Vmax)*2.897628) - (Vmax**2/Amax) 我在 sympy 中编写了以下代码 ..
发布时间:2021-09-03 20:17:22 Python

用 Sympy 求解线性方程组

为了评估某种弹性属性,我想使用 sympy 来可视化方程组.我使用以下代码: 将 numpy 导入为 np将 sympy 导入为 symb1=sym.Array([[Rational(-1/2),sqrt(3)/2,0],[-sqrt(3)/2,Rational(-1/2),0],[0,0],1]])西格玛=[]对于范围内的 i (0,3):对于范围内的 j (0,3):对于范围内的 k (0, ..
发布时间:2021-09-03 20:17:20 Python

Einstein-Convention 喜欢隐式求和与 sympy

给定一些简单的张量 A = sympy.tensor.array.Array() 和表示基变换的矩阵 T 和 T1=T.inv(),是否可以使用像 这样的符号 B[i,j,k] = T[i,a] * A[a,b,c] * T1[b,j] * T1[c,k] 计算转换后的张量? 原则上,似乎可以使用爱因斯坦求和约定,但我遇到了多个问题: 代码片段 from sympy 导入符 ..
发布时间:2021-09-03 20:17:17 其他开发

是否有等价于 Mathematica 的“用法"?功能?

在数学中,可以定义一个函数的标签,例如 f::usage = "f[x] 给出 (x - 1)(x + 1)"; 然后可以这样调用: ?f SymPy 中是否有等价物? 解决方案 您可以“猴子补丁"文档字符串如下: import sympyf = sympy.Function("f")f.__doc__ = "一个没有任何特殊属性的裸函数" 您可以以相同的方式检索文档字符串,例 ..
发布时间:2021-09-03 20:17:14 其他开发