-
0 votesanswersviews
使用Sympy在IS-LM模型中推导出均衡
我不是编码员,对于我的经济课,我们必须使用python . 其中一个问题要求我们使用sympy得出IS-LM模型的均衡Y和r . 任何帮助将非常感激 . 对于IS-LM模型,IS方面看起来像:Y = CIG,其中Y是实际GDP,C是消费,I是计划投资支出,G是政府采购C = C0 cY,其中C0是一次性消费率收入(Yt)为0,cY是边际消费倾向I = I0 - br LM侧看起来像:(M / ... -
8 votesanswersviews
在SymPy中反转矩阵?
我想知道如何在Python中使用SymPy创建矩阵并计算其逆? 例如,对于此符号矩阵: -
4 votesanswersviews
在同情中得到矩阵乘法的元素方程
我有2个矩阵,其中第一个是稀疏的整数系数 . import sympy A = sympy.eye(2) A.row_op(1, lambda v, j: v + 2*A[0, j]) 第二个是象征性的,我在它们之间执行操作: M = MatrixSymbol('M', 2, 1) X = A * M + A.col(1) 现在,我想要的是获得元素方程式: X_{0,0} = A_{0,0} ... -
2 votesanswersviews
将Sympy矩阵替换为多项式
我有一个Sympy矩阵A和一个多项式表达式P,我想计算P(A) . 这是一个例子: x = Symbol('x') P = x**2 - 3*x + 5 A = Matrix([ [1,3], [-1,2] ]) P.subs(x, A) 我希望Sympy计算 A**2 - 3*A + 5*eye(2) (在这个例子中,结果是零矩阵) . 但是这失败并显示错误消息: AttributeError... -
3 votesanswersviews
如何用SymPy中的向量替换矩阵的对角元素?
我有一个矢量 X ,我这样创建: from sympy import * x1 = Symbol('x1') x2 = Symbol('x2') x3 = Symbol('x3') X = Matrix([x1, x2, x3]) 然后我还有一个矩阵 myMat ,它只包含一个: myMat = ones(3, 3) Matrix([ [1, 1, 1], [1, 1, 1], [1, 1... -
1 votesanswersviews
Matrix在SymPy中定义了未定义的函数
我正在寻找指定依赖于变量的矩阵量的可能性 . 对于如下工作的标量,使用未定义的函数: from sympy import * x = Function('f')(t) diff(x,t) 对于像Matrix这样的符号 x = MatrixSymbol('x',3,3) 我找不到相应的东西 . 有 i,j = Symbols('i j') x = FunctionMatrix(6,1,Lambd... -
7 votesanswersviews
如何在SymPy中扩展矩阵表达式?
在SymPy中,我试图执行矩阵乘法并在之后展开它 . 但是,SymPy似乎不支持矩阵表达式的扩展 . 例如,这是矩阵的4阶Runge-Kutta(RK4): from sympy import init_session init_session() from sympy import * A = MatrixSymbol('A', 3, 3) x = MatrixSymbol('x', 3, 1... -
3 votesanswersviews
如何使用Sympy导出Matrix元素
给定矩阵和向量的乘积 A.v 对于形状(m,n)和v的dim n,其中m和n是符号,我需要相对于矩阵元素计算导数 . 我还没有找到使用正确向量的方法,所以我从2 MatrixSymbol 开始: n, m = symbols('n m') j = tensor.Idx('j') i = tensor.Idx('i') l = tensor.Idx('l') h = tensor.Idx('h')... -
17 votesanswersviews
SymPy - 任意数量的符号
我正在编写一个函数来解决任意数量的联立方程 . 方程的数量由函数的一个参数设置,每个方程由许多符号构成 - 与方程一样多的符号 . 这意味着我不能简单地对方程式进行硬编码,甚至不需要将方程组合在一起所需的符号;该函数需要能够处理任意数量的方程 . 所以,我的问题是,如何生成符号列表? 我有一个可能的解决方案,但我的直觉告诉我,它不会非常有效 . 如果有更好的方法,请告诉我 . 我是SymPy的新手... -
3 votesanswersviews
使用Sympy,python2.7求解方程组
我想解决一个方程组 . 但我希望能够精确地将值变为“获得”,并将其作为“什么”的函数 . 为了更好地理解,我从here开始,我修改了一个例子: import sympy as sp x, y, z = sp.symbols('x, y, z') rho, sigma, beta = sp.symbols('rho, sigma, beta') f1 = sigma * (y - x) f2 = x... -
1 votesanswersviews
在Python中解决符号线性方程给出了意想不到的答案
首先,请原谅即将到来的长壁测试 . 我需要解决由坐标变换引起的以下方程组: 给定v_n,v_m,x_p,z_p,y_p,s,d,解决下面的系统v_yp,v_zp: I: v_n = v_yp * (1/(s*d)) + v_xp * y_p * (1/(s*d^2)) II: v_m = v_zp * (1/(s*d)) + v_xp * z_p * (1/(s*d^2)) III: v_... -
2 votesanswersviews
如何用sympy分析和数值地求解一阶线性微分方程?
如何在_609045中解决像这样的简单线性微分方程? y' + p(t)y = q(t) 我希望以两种方式解决它:如果可能的话,象征性地(分析地),如果sympy可以得出积分因子等,还可以用数字方式进行,以便可以比较两者 . 怎么能在 sympy 完成?是 sympy.mpmath.odefun 正确的地方吗? -
0 votesanswersviews
非线性方程组可以在Maple中求解,但不能在同一时间解决
我很新见,并享受pythonic语法 . 但是,我遇到了一个问题,无法通过sympy解决,但可以在Maple中轻松解决 . 所以,我有以下系统: 0.0165 * exp( -2.0405-0.33*b-0.5*n)+0.031 * exp(-4.164-0.62*b-0.5*n)=k*p 0.025 * exp( -2.0405-0.33*b-0.5*n) +0.025 * exp(-4.164... -
3 votesanswersviews
使用scipy.integrate.odeint时出现TypeError
我正试图用 scipy.integrate.odeint 求解一组耦合微分方程 . 但是,当我尝试运行该程序时,我收到以下错误: TypeError: Cannot cast array data from dtype('O') to dtype('float64') according to the rule 'safe' odepack.error: Result from function ... -
7 votesanswersviews
在SymPy中自动填充矩阵元素
有没有办法在SymPy中隐式定义符号矩阵的元素遵循规则,例如: symbol 后跟矩阵中的 subindices (或数字对) 例如,我想定义一个名为 M 的3 x 2矩阵,我希望SymPy自动创建它并将其填充为: M = [ M_11 M_12] [ M_21 M_22] [ M_31 M_32] 如果无法隐式执行此操作,那么显式执行此操作的最简单方法是什么(例如循环)? -
1 votesanswersviews
是否有相当于python / scipy中的Mathematicas RootApproximant函数,numpy,sympy?
我在python或python lib中寻找相当于Mathematicas RootApproximant函数:https://reference.wolfram.com/language/ref/RootApproximant.html 基本上,该函数找到数字的代数根,例如, RootApproximant[1.414213] -> sqrt(2) 谢谢! -
3 votesanswersviews
使用sympy和Mathematica对Modified Bessel函数进行符号积分
我想用 sympy 解决一个我知道Mathematica可以解决的定积分 . 在Mathematica中,以下行 Integrate[z^2 (BesselI[0, z^2] - BesselI[1, z^2]) Exp[-z^2], {z, 0, x}] 产量 1/3 x^3 HypergeometricPFQ[{1/2,3/2},{1,5/2},-2 x^2]-1/10 x^5 Hyperg... -
1 votesanswersviews
Python(SymPy):如何获得与Mathematica相同的结果?
我使用 Mathematica 和 Python 计算了以下内容 . Mathematica 使用以下代码 f[x_] = a*b/(a - b)^2*Exp[-r*x] (Exp[-b*x] - Exp[-a*x]) (a*Exp[-b*x] - b*Exp[-a*x]) Assuming[{a > 0, b > 0, r > 0}, Integrate[f[x], {x, ... -
1 votesanswersviews
使用sympy / pydy限制/控制输出scipy的odeint
我正在使用PyDy为三体系统创建ODE系统: right_hand_side = pydy.codegen.generate_ode_function( forcing_vector, model.coordinate_symbols, model.speed_symbols, model.constant_symbols, m... -
2 votesanswersviews
在python中使用sympy包将值替换为10x10矩阵
我正在使用符号包sympy以10个变量存储10X10反对称矩阵 . 然后在每个迭代步骤中,我将数值替换为矩阵的条目 . 但是,替换值需要很长时间,因为符号表达式非常复杂 . 每次替换大约需要100秒,我需要至少10,000次迭代才能获得有意义的结果 . 所以,我想知道是否有任何快速的方法来取代 Value 观 . 现在,我只是将值替换为矩阵的顶部对角线一半,这使得替换的数量减半,但它仍然需要太长时... -
3 votesanswersviews
使用sympy计算矩阵的伪逆
我应该如何使用sympy计算矩阵的伪逆(不使用numpy,因为矩阵具有符号常数,我也想要符号中的逆) . 正常 inv() 对于sympy中的非方矩阵不起作用 . 例如,如果 M = Matrix(2,3, [1,2,3,4,5,6]), pinv(M) 应该给 -0.9444 0.4444 -0.1111 0.1111 0.7222 -0.2222 -
5 votesanswersviews
Sympy“全球”替代
我有很多符号表达式,我可能会发现其中一个系数为零 . 我想,也许是因为我习惯了mathematica,以下是有道理的: from sympy import Symbol x = Symbol('x') y = Symbol('y') f = x + y x = 0 f 令人惊讶的是,返回的是 x + y . 除了在每个等式上明确调用"subs"之外,有没有办法让 f 只返回... -
0 votesanswersviews
使用字典的Sympy矩阵替换
我知道this post但我没有看到答案 . 我想将数值代入一个sympy矩阵 . 更详细地说,我确实 Build 了一个系统,在某些情况下我想插入基于numpy数组的数值 . 我认为可以使用词典 . 但要么我误解了这个概念,要么它不是该语言的一个特征 . 我有一个矩阵: import sympy as sp cp = sp.MatrixSymbol('cp',3,3) C = sp.Matrix... -
1 votesanswersviews
写一个sympy矩阵到python文件
我有一个自动生成方程式的脚本,并将它们添加到列表中 . 例如: eq = [a*b, a + b*c -d, c**2 + a, d*a - 2] 用a,b,c,d表示所有符号 . 然后它将列表转换为sympy矩阵并计算雅可比矩阵 eq = sympy.Matrix(eq) jacobi = eq.jacobian([a, b, c, d]) 我想将这个jacobian保存在python文件中... -
2 votesanswersviews
如何在sympy中创建符号变量的向量
如何在sympy中创建符号变量的向量?我想做点什么 x, x1, x2, x3 = symbols ('x x1 x2 x3') A = [x+x1,x+x2,x+x3] B = A * Transpose(A) print (B) A 是符号变量数组 . 我查看了同情documentation,但无法弄清楚 . (Python 2.7.6,sympy 0.7.4.1) 更新: 我想做点什么 x... -
2 votesanswersviews
在同义词中对矩阵进行逐列(或逐行)操作
是否有办法在同义词中对矩阵进行逐列(或逐行)操作?例如,将矩阵的每一列除以其范数,或将矩阵的每一行乘以其范数? -
0 votesanswersviews
如果它存在于SymPy中,则找到奇异矩阵方程的答案
假设我们有矩阵方程 A*x = b ,其中 A 是正方形但是奇异的 . 然后一般来说,方程式要么没有解决方案要么没有多少解决方案 . 如果确实有解决方案,我想使用SymPy以自由变量的形式象征性地找到它们的形式 . 我知道如果我们能找到一个特定的解决方案 x ,那么通过将 A 的零空间中的向量添加到 x 上,我们可以找到其余的,所以问题就减少到找到一个解决方案 . 在Mathematica中,函数... -
-1 votesanswersviews
如何将Numpy矩阵转换为Sympy
我构建了一个随机生成的指定大小的矩阵,该部分工作得很好 . 给它一个行和列大小,并从0到100繁殖一个整数的矩阵 . 最近我尝试对一个numpy矩阵执行一个sympy操作,python一直在崩溃我 . 我很快就知道,来自sympy的操作无法在numpy矩阵上运行 . 因此,我研究了如何将numpy转换为sympy,但更多时候我只是使用lambdify找到了numpy的同情 . 我想知道我是否可以使... -
3 votesanswersviews
如何将一个sympy多项式转换为一个列表?
我正在尝试使用sympy来自动化多项式拟合问题的解决方案 . 例如,假设我定义了一个多项式函数: #!/usr/bin/python from sympy import * x, xa, xb, a0, a1, a2, a3 , a4 = symbols(('x', 'xa', 'xb', 'a0', 'a1', 'a2', 'a3', 'a4')) def p4(x): re... -
5 votesanswersviews
有限域上的SymPy多项式
import sympy as S F = S.FiniteField(101) 当我调用 f = S.poly(y ** 2 - x ** 3 - x - 1,F) 时,我收到以下错误: 'FiniteField'对象没有属性'is_commutative' 但有限域是定义的可交换的!所以我不确定这个错误应该是什么意思! 有没有人遇到过这个?如何在有限域上声明多项式?