我有两个点(x1和x2),希望在给定的步数中生成正态分布。x1和x2之间x值的y值之和为1。对于实际问题:
我对Python相当陌生,不知道为什么下面的代码会产生预期的结果,但是比PHP中的同一个程序慢了大约100倍。大约有2000个x1-x2对,每对大约有5个阶跃值。在
我试着用Cython编译,使用了多处理,但它只改进了2倍,仍然比PHP慢50倍。有什么建议至少能提高PHP的性能?在from scipy.stats import norm
import numpy as np
import time
# Calculates normal distribution
def calculate_dist(x1, x2, steps, slope):
points = []
range = np.linspace(x1, x2, steps+2)
for x in range:
y = norm.pdf(x, x1+((x2-x1)/2), slope)
points.append([x, y])
sum = np.array(points).sum(axis=0)[1]
norm_points = []
for point in points:
norm_points.append([point[0], point[1]/sum])
return norm_points
start = time.time()
for i in range(0, 2000):
for j in range(10, 15):
calculate_dist(0, 1, j, 0.15)
print(time.time() - start) # Around 15 seconds or so
编辑,PHP代码:
^{pr2}$
编辑2,分析每一行并发现规范.pdf()占用了98%的时间,因此找到了一个自定义的normpdf函数并对其进行了定义,现在的时间大约是0.67秒,这是相当快的,但仍然比PHP慢10倍左右。我还认为重新定义公共函数违背了Python的简单性思想?!在
自定义函数(源是其他Stackoverflow答案):from math import sqrt, pi, exp
def normpdf(x, mu, sigma):
u = (x-mu)/abs(sigma)
y = (1/(sqrt(2*pi)*abs(sigma)))*exp(-u*u/2)
return y