如何正确计算以下等式？

使用python 2.7我试图计算下列相当困难的方程组。如何正确计算以下等式？

我已经成功地实现了前两个，但我与第三挣扎。这里是我的尝试，

def pi_tau(r): 
    return (1 - (1 - r)**(t + 1))/(2 - r) 

def mean_field(r): 
    return 1 - (1 - r)**t 

def pi_neighbour(r): 
    inside = [] 
    for kappa_dash in range(0, kappa - 1): 

     Binomial_Coefficient = (math.factorial(kappa - 1))/(math.factorial(kappa - 1 - kappa_dash)*math.factorial(kappa_dash))   

     top = ((mean_field(r)*pi_tau(r))**kappa_dash)*(1 - mean_field(r)*pi_tau(r))**(kappa - 1 - kappa_dash) 

     bottom = kappa_dash + 1 

    fraction = top/bottom 
    inside.append(kappa_dash) 
    inside[kappa_dash] = inside[kappa_dash] + fraction*Binomial_Coefficient 

return pi_tau*inside

我再尝试调用这个函数

# set parameters 
r = 0.15 
kappa = 2.1 
T = 10 

ppp_t = [] 
mmm_f = [] 
nnn_t = [] 

for t in range(0, T): 
    ppp_t.append(pi_tau(r)) 
    mmm_f.append(mean_field(r)) 
    nnn_t.append(pi_neighbour(r))

我收到以下错误消息

--------------------------------------------------------------------------- 

    TypeError         Traceback (most recent call last) 
    <ipython-input-41-9bcf462306f0> in <module>() 


    6  ppp_t.append(pi_tau(r)) 
     7  mmm_f.append(mean_field(r)) 
----> 8  nnn_t.append(pi_neighbour(r)) 

<ipython-input-39-d9acdd7490f9> in pi_neighbour(r) 
    19   inside[kappa_dash] = inside[kappa_dash] + fraction*Binomial_Coefficient 
    20 
---> 21  return pi_tau*inside 

TypeError: can't multiply sequence by non-int of type 'function'

我正在寻找关于如何实现的任何意见第三个功能和改进我的方法。

来源

2017-06-13 RedPen

你可以把标题和问题文本中的公式的名称，这样别人搜索可能会发现它们更容易。标题，例如，“在Python中计算名称1，名称，名称3”或其他内容。 –

'pi_tau'是一个函数，而不是一个数字。你打算叫'pi_tau'吗？ – Carcigenicate

这里有一些奇怪的事情：

你计算的顶部，底部和bionomial都在for循环，但不要总结在for环;和
你乘以pi_tau功能与列表？
```
return pi_tau*inside 
```
您使用range(0, kappa - 1)，但上界的range(..)是独家

然而，你做的事情挺难的。你可以用下面的办法：

from math import factorial 

def pi_neighbour(r): 
    sum = 0.0 
    for kappa_dash in range(0,kappa): 
     bin_coeff = factorial(kappa-1)/(factorial(kappa-1-kappa_dash)*factorial(kappa_dash))  
     top = ((mean_field(r)*pi_tau(r))**kappa_dash)*(1-mean_field(r)*pi_tau(r))**(kappa-1-kappa_dash) 
     sum +=bin_coeff*top/(kappa_dask+1) 
    return pi_tau(r)*sum

来源

2017-06-13 21:59:04

我相信你的意思做的是在回溯的线21以下：

return pi_tau(r)*inside

你忘了打电话的功能，因此它正试图通过一个功能乘以一个整数，而不是功能的返回值。

来源

2017-06-13 21:54:17

的错误是，你乘以一个列表，这是inside和一个整数，即return of pi_tau(r)

做到这一点，而不是（一list comprehension）

试这样的： return [i*pi_tau(r) for i in inside]

来源

2017-06-13 22:08:31 rafaeljneves

如何正确计算以下等式？

回答

相关问题