DifferentialEvolution.py 文件源码

def NPV_prob(rho,g,eff,H,delta_t,p,Q,mean,std,a,b):
  t = np.arange(0,100000,0.5)
  PV = simps(np.exp(-t*r)*rho*g*eff*H*delta_t*p/1000*EV_flow(mean,std,Q), t)
  NPVal = PV - cost(a, b, Q)
  return -1*NPVal