0

我正在尝试为我的日震学课程制定一个情节，并且该问题提供了一个分段函数来描述恒星中“流体”的动力学，就好像它是一回事，如果它是另一回事。我一遍又一遍地收到这个'Mul' object cannot be interpreted as an integer，但我正在处理实数中的数字，而不仅仅是整数集。我不知道如何解决这个问题，需要指导。代码如下。

import sympy as sy
from sympy import *
from sympy.physics.units import Unit
import numpy as np
import sys
import math
import scipy as sp
from scipy import special

phi = Symbol('phi', Variable = True)
x = Symbol('x', Variable = True, Real = True)
t = Symbol('t', Variable = True, Real = True)
xi = Symbol('xi', Function = True)
Solar_Radius = Symbol('R', Constant = True, unit = "meters")
Sound_Speed = Symbol('c', Constant = True, unit = "meters per second", Real = True)
gamma = Symbol('gamma', Constant = True)
gravity = Symbol('g', Constant = True, unit = "meters per second per second")

Solar_Radius = 6.963 * 10 ** 6
gamma = 5/3
g = 274.8265625336
gas_constant = 8201.25
c = 8.1 * 10 ** 3

for t in range(0,x/c):
    xi[x,t] = 0
for t in range(x/c,00):
    xi[x,t] = (1/2)*sy.exp(gamma*g*x/(2*c**2))*mpmath.besselj(0, (gamma*g/(2*c)*sy.sqrt(t**2 - ((x/c)**2))),derivative = 0)

完整追溯：

---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-50-3506376f1686> in <module>()
----> 1 for t in range(0,x/c):
      2     xi[x,t] = 0
      3 for t in range(x/c,00):
      4     xi[x,t] = (1/2)*sy.exp(gamma*g*x/(2*c**2))*mpmath.besselj(0, (gamma*g/(2*c)*sy.sqrt(t**2 - ((x/c)**2))),derivative = 0)

TypeError: 'Mul' object cannot be interpreted as an integer
4

1 回答 1

2

这里有很多问题：

  • 您传递给 Symbol 的关键字参数（常量、变量、单位、实数）都不是 SymPy 识别的东西。唯一接近的是real，它应该是小写的（如Symbol('x', real=True)）。其余的什么都不做。如果你想要单位，你应该在sympy.physics.units. 无需指定符号是常量还是变量。

  • 您已将Solar_Radius和重新定义gamma为数字。这意味着这些变量的符号定义毫无意义。

  • 如果您使用的是 Python 2，请确保包含from __future__ import division在文件的顶部，否则类似1/2and之类的5/3内容将被整数除法截断（这在 Python 3 中不是问题）。

  • range(0, x/c)没有意义。range创建一个数字列表，例如range(0, 3)-> [0, 1, 2]。但x/c不是一个数字，它是一个象征性的表达。

  • 另外，xi[x, t] = ...没有意义。xi是一个符号，它不允许索引，当然也不允许赋值。

  • 不要将数字（math、mpmath、numpy、scipy）函数与 SymPy 函数混合使用。它们不适用于符号表达式。您应该只使用 SymPy 函数。如果您创建一个符号表达式并希望将其转换为数字表达式（例如，用于绘图），请使用lambdify.

你在这里想要的是Piecewise. 语法是Piecewise((expr, cond), (expr, cond), ..., (expr, True)), where是当为真expr时使用的表达式（是“否则”条件）。cond(expr, True)

对于你的例子，我相信你想要

expr = Piecewise((0, t < x/c), (sy.exp(gamma*g*x/(2*c**2))*sy.besselj(0, (gamma*g/(2*c)*sy.sqrt(t**2 - (x/c)**2)))/2, t >= x/c))

如果要将其转换为 and 中的数值函数x，请t使用

xi = lambdify((x, t), expr)
于 2016-10-04T18:52:26.227 回答