0

我在 Mathematica 中为 f(x,t) 创建了一个原始方程:

uExact[x_, t_] := (1 + I*t^2)* Cos[2 * Pi *x];

f[x_, t_] := 
  4 * Pi^2 * (1 + I * t^2) * Cos[2* Pi * x] + 2*I*t*Cos[2 Pi*x] - 
   I * g * log[(t^4 + 1) * Cos^2[2*Pi*x]];

我已将其转换为以下代码:

def u_exact(x, t):

    return (1 + I*t**2)* cmath.cos(2 * cmath.pi *x)
    


#check!!!
def f(x,t):
    #1return 4 * cmath.pi**2 * (1 + I * t**2) * cmath.cos(2* cmath.pi * x) + 2*I*t*cmath.cos(2 * cmath.pi*x) - I * g * cmath.log((t**4 + 1) * cmath.cos**2(2*cmath.pi*x))
    

    #alternative

    p1 = 4 * pow(cmath.pi,2) * (1 + I * pow(t,2)) * pow(cmath.cos,2)* (2*cmath.pi * x)
    p2 = 2*I*t*cmath.cos(2 * cmath.pi*x)
    p3 = - I * g * cmath.log(pow(t,4) + 1) * pow(cmath.cos,2) * (2*cmath.pi*x)

    return (p1 + p2 + p3)

在 #1 和 #alternative 的情况下,我不断收到相同的错误:** 或 pow() 不支持的操作数类型:'builtin_function_or_method'

可能是什么问题呢?我错过了什么吗?

4

2 回答 2

1

cmath.cos是一个函数。您正在以函数为基础应用 pow 函数。

>>> pow(cmath.cos, 2)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
TypeError: unsupported operand type(s) for ** or pow(): 'builtin_function_or_method' and 'int'
>>> pow(cmath.cos(1), 2)
(0.2919265817264289+0j)
于 2021-01-05T12:00:51.727 回答
0

在您的两种选择中,您都试图将幂运算应用于函数(math.cos):

cmath.cos**2
pow(cmath.cos,2)

cmath.cos是函数,如果你想应用幂函数,你需要有一个浮点数,所以你需要调用函数:(cmath.cos(value)替换value为你需要的任何东西)。

于 2021-01-05T12:02:25.473 回答