Question

我的数学问题在这里 - https://math.stackexchange.com/questions/2063507/solving-this-integral-involving-ei-function

这与种群动态有关，t =时间，N_t =时间t的种群，r =生长速度，K = cqrrying能力。

我的代码附在下面。

Python无法计算N_next的值，因为它在指数积分函数scipy.special.expi()内。我怎么能绕过这个呢？它现在说＆＃34; Scipy没有特殊属性＆＃34;但根据这个 - https://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.special.expi.html - 它应该是。

import math
import scipy
import numpy as np
import matplotlib
import matplotlib.pyplot as plt

t_f =100
N_0 = 10
t = []
N_t = [N_0,]
r = 2.5
K = 1000

for i in range(0,100):
    scipy.special.expi(r*N_next/K) = i*math.exp(r) + scipy.special.expi(r/K * N_t[i])
    N_t.append(N_next)
    t.append(i)

plt.plot(t,N_t)
plt.show

Answer 1

该错误是Python解释的奇特方式：您无法将表达式结果分配给函数调用

也许你在你的例子中混合了这一行？在这种情况下，请编辑您的帖子，以便我们为您提供帮助。

Answer 2

https://math.stackexchange.com/questions/2063507/solving-this-integral-involving-ei-function中的答案提供了一个公式，您希望为 t 的值序列进行评估。在构建这样的脚本时，我建议你采用自下而上的方法。这是我到目前为止所写的内容。

from numpy import exp
import matplotlib.pyplot as plt

K = 1000
r = 2.5

N_0 = 10
N_t = 50

def Ei(x):
    return x

def Ei_inv(x):
    return 1/x

def N(t):
    return K*Ei_inv(exp(r)*t+Ei(r*N_0/K))/r

t = [_ for _ in range(40)]
N = [N(_) for _ in t]

plt.plot(t,N)
plt.show()

我知道一些常数是什么;因此，我把它们放在代码中。
我不知道Ei或它的逆是什么，或者我可能在哪里找到这些功能;作为一项临时措施，我写了“假”功能来取代他们的位置。
现在我可以根据虚函数编写我实际想要评估的函数作为具有单个参数的另一个Python函数。那时我使用print (N(50))对我所拥有的内容进行了试用，以验证代码的各个部分是否在一起表现得很礼貌。
我添加了两行，一行用于定义 x ，或者是时间，用于绘图的值，第二行用于 y 或 N 价值观。

轮到你，用计算Ei及其倒数所需的任何东西替换虚函数。

Answer 3

所以，我明白了。嗯，有点。代码现在非常耗费计算量。但它至少可以运作。它在下面。

import math
from scipy import *
from scipy.special import expi
import numpy as np

t_f =100
N_0 = 10
t = [0,]
n = [N_0,]
r = 2.5
K = 100
N_t = [N_0,]

for i in range(0,100):
    a= i*math.exp(r) + expi(r/K * N_0)
    n.append(a)
    t.append(i+1)

g = np.arange(0,1000,0.0001)
d = []
e = []

for i in g:
    x = expi(i)
    d.append(x)
    e.append(i)    

for j in range(1,99):
    b = n[j]
    c = []

    for i in g:
        c.append(b)

    t = np.subtract(d,c)
    for i in range(len(t) - 1):
        if t[i] == 0. or t[i] * t[i + 1] < 0.:
            y = e[i]
            print(y)

    N_next = math.floor(y*K/r)
    N_t.append(N_next)

print(N_t)
plot(t,N_t)

变量内指数积分函数。 Python无法计算

3 个答案: