0

我试图找到一个表达式的拉普拉斯逆，除了一个变量之外的所有变量在声明时都已经定义：

from numpy import *
import mpmath as mp
p0 = 1
E = 2
c= 3
L = 4
x = 2.5
t = linspace(1,5,10)
ulaplace = []

def U(s):
    return(c*p0*(-exp(L*s/c) + exp(s*(L + 2*x)/c))*exp(-s*x/c)/(E*s**2*(exp(2*L*s/c) + 1)))

for ti in t:
    ulaplace.append(mp.invertlaplace(U, ti, method='talbot'))

但我收到此错误：

Traceback (most recent call last):
  File "D:\TEMP\IDLEscripts\CompareAnalyticalSolutions2.py", line 46, in <module>
    ulaplace.append(mp.invertlaplace(U, ti, method='talbot'))
  File "C:\Python35\lib\site-packages\mpmath\calculus\inverselaplace.py", line 805, in invertlaplace
    fp = [f(p) for p in rule.p]
  File "C:\Python35\lib\site-packages\mpmath\calculus\inverselaplace.py", line 805, in <listcomp>
    fp = [f(p) for p in rule.p]
  File "D:\TEMP\IDLEscripts\CompareAnalyticalSolutions2.py", line 43, in U
    return(c*p0*(-exp(L*s/c) + exp(s*(L + 2*x)/c))*exp(-s*x/c)/(E*s**2*(exp(2*L*s/c) + 1)))
TypeError: attribute of type 'int' is not callable

我也尝试了文档网站lambda function建议的格式，但仍然出现同样的错误。

该mpmath.invertlaplace函数是否要求在定义时所有内容都以数字形式表示？我问是因为这有效：

>>> import mpmath as mp
>>> def F(s):
    return 1/s

>>> mp.invertlaplace(F,5, method = 'talbot')
mpf('1.0')

如果是这样，我需要能够规避这一点。对我来说，重点是玩弄其他变量，看看它们如何影响逆拉普拉斯算子。此外，有人会认为函数在传递给mpmath.

如果不是，那么这里到底发生了什么？

4

1 回答 1

1

好吧，我明白了。基本上mp.invertlaplace需要自己的函数只使用mpmath定义的函数。在我exp从numpy库中使用的原始问题中提供的代码中。exp(x)确实如此numpy.exp(x)。为了使代码正常工作，它需要mpmath.exp按如下方式调用函数：

def U(s):
    return -p0*mp.exp(s*x/c)/(E*s*(-s*mp.exp(L*s/c)/c - s*mp.exp(-L*s/c)/c)) + p0*mp.exp(-s*x/c)/(E*s*(-s*mp.exp(L*s/c)/c - s*mp.exp(-L*s/c)/c))

我没有在原始问题中提供的简化示例上测试上述内容，因为它是更通用脚本的子集。但是它应该可以工作，这似乎是问题的根源。

于 2017-08-04T08:02:23.593 回答