首页 文章
  • 0 votes
     answers
     views

    使用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 votes
     answers
     views

    在SymPy中反转矩阵?

    我想知道如何在Python中使用SymPy创建矩阵并计算其逆? 例如,对于此符号矩阵:
  • 4 votes
     answers
     views

    在同情中得到矩阵乘法的元素方程

    我有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 votes
     answers
     views

    将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 votes
     answers
     views

    如何用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 votes
     answers
     views

    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 votes
     answers
     views

    如何在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 votes
     answers
     views

    如何使用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 votes
     answers
     views

    SymPy - 任意数量的符号

    我正在编写一个函数来解决任意数量的联立方程 . 方程的数量由函数的一个参数设置,每个方程由许多符号构成 - 与方程一样多的符号 . 这意味着我不能简单地对方程式进行硬编码,甚至不需要将方程组合在一起所需的符号;该函数需要能够处理任意数量的方程 . 所以,我的问题是,如何生成符号列表? 我有一个可能的解决方案,但我的直觉告诉我,它不会非常有效 . 如果有更好的方法,请告诉我 . 我是SymPy的新手...
  • 3 votes
     answers
     views

    使用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 votes
     answers
     views

    在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 votes
     answers
     views

    如何用sympy分析和数值地求解一阶线性微分方程?

    如何在_609045中解决像这样的简单线性微分方程? y' + p(t)y = q(t) 我希望以两种方式解决它:如果可能的话,象征性地(分析地),如果sympy可以得出积分因子等,还可以用数字方式进行,以便可以比较两者 . 怎么能在 sympy 完成?是 sympy.mpmath.odefun 正确的地方吗?
  • 0 votes
     answers
     views

    非线性方程组可以在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 votes
     answers
     views

    使用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 votes
     answers
     views

    在SymPy中自动填充矩阵元素

    有没有办法在SymPy中隐式定义符号矩阵的元素遵循规则,例如: symbol 后跟矩阵中的 subindices (或数字对) 例如,我想定义一个名为 M 的3 x 2矩阵,我希望SymPy自动创建它并将其填充为: M = [ M_11 M_12] [ M_21 M_22] [ M_31 M_32] 如果无法隐式执行此操作,那么显式执行此操作的最简单方法是什么(例如循环)?
  • 1 votes
     answers
     views

    是否有相当于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 votes
     answers
     views

    使用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 votes
     answers
     views

    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 votes
     answers
     views

    使用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 votes
     answers
     views

    在python中使用sympy包将值替换为10x10矩阵

    我正在使用符号包sympy以10个变量存储10X10反对称矩阵 . 然后在每个迭代步骤中,我将数值替换为矩阵的条目 . 但是,替换值需要很长时间,因为符号表达式非常复杂 . 每次替换大约需要100秒,我需要至少10,000次迭代才能获得有意义的结果 . 所以,我想知道是否有任何快速的方法来取代 Value 观 . 现在,我只是将值替换为矩阵的顶部对角线一半,这使得替换的数量减半,但它仍然需要太长时...
  • 3 votes
     answers
     views

    使用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 votes
     answers
     views

    Sympy“全球”替代

    我有很多符号表达式,我可能会发现其中一个系数为零 . 我想,也许是因为我习惯了mathematica,以下是有道理的: from sympy import Symbol x = Symbol('x') y = Symbol('y') f = x + y x = 0 f 令人惊讶的是,返回的是 x + y . 除了在每个等式上明确调用"subs"之外,有没有办法让 f 只返回...
  • 0 votes
     answers
     views

    使用字典的Sympy矩阵替换

    我知道this post但我没有看到答案 . 我想将数值代入一个sympy矩阵 . 更详细地说,我确实 Build 了一个系统,在某些情况下我想插入基于numpy数组的数值 . 我认为可以使用词典 . 但要么我误解了这个概念,要么它不是该语言的一个特征 . 我有一个矩阵: import sympy as sp cp = sp.MatrixSymbol('cp',3,3) C = sp.Matrix...
  • 1 votes
     answers
     views

    写一个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 votes
     answers
     views

    如何在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 votes
     answers
     views

    在同义词中对矩阵进行逐列(或逐行)操作

    是否有办法在同义词中对矩阵进行逐列(或逐行)操作?例如,将矩阵的每一列除以其范数,或将矩阵的每一行乘以其范数?
  • 0 votes
     answers
     views

    如果它存在于SymPy中,则找到奇异矩阵方程的答案

    假设我们有矩阵方程 A*x = b ,其中 A 是正方形但是奇异的 . 然后一般来说,方程式要么没有解决方案要么没有多少解决方案 . 如果确实有解决方案,我想使用SymPy以自由变量的形式象征性地找到它们的形式 . 我知道如果我们能找到一个特定的解决方案 x ,那么通过将 A 的零空间中的向量添加到 x 上,我们可以找到其余的,所以问题就减少到找到一个解决方案 . 在Mathematica中,函数...
  • -1 votes
     answers
     views

    如何将Numpy矩阵转换为Sympy

    我构建了一个随机生成的指定大小的矩阵,该部分工作得很好 . 给它一个行和列大小,并从0到100繁殖一个整数的矩阵 . 最近我尝试对一个numpy矩阵执行一个sympy操作,python一直在崩溃我 . 我很快就知道,来自sympy的操作无法在numpy矩阵上运行 . 因此,我研究了如何将numpy转换为sympy,但更多时候我只是使用lambdify找到了numpy的同情 . 我想知道我是否可以使...
  • 3 votes
     answers
     views

    如何将一个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 votes
     answers
     views

    有限域上的SymPy多项式

    import sympy as S F = S.FiniteField(101) 当我调用 f = S.poly(y ** 2 - x ** 3 - x - 1,F) 时,我收到以下错误: 'FiniteField'对象没有属性'is_commutative' 但有限域是定义的可交换的!所以我不确定这个错误应该是什么意思! 有没有人遇到过这个?如何在有限域上声明多项式?

热门问题