1
私はいくつかの確率密度関数を持っている:分割統合確率密度
T = 10000
tmin = 0
tmax = 10**20
t = np.linspace(tmin, tmax, T)
time = np.asarray(t) #this line may be redundant
for j in range(T):
timedep_PD[j]= probdensity_func(x,time[j],initial_state)
私はxの2つの異なる領域の上にそれを統合したいです。私は2つの空間領域にtimedep_PD
配列を分割するために、次の試み、その後、統合を進め:
step = abs(xmin - xmax)/T
l1 = int(np.floor((abs(ab - xmin)* T)/abs(xmin - xmax)))
l2 = int(np.floor((abs(bd - ab)* T)/abs(xmin - xmax)))
#For spatial region 1
R1 = np.empty([l1])
R1 = x[:l1]
for i in range(T):
Pd1[i] = Pd[i][:l1]
#For spatial region 2
Pd2 = np.empty([T,l2])
R2 = np.empty([l2])
R2 = x[l1:l1+l2]
for i in range(T):
Pd2[i] = Pd[i][l1:l1+l2]
#Integrating over each spatial region
for i in range(T):
P[0][i] = np.trapz(Pd1[i],R1)
P[1][i] = np.trapz(Pd2[i],R2)
は、2つの空間領域に確率密度関数を分割して、統合について移動する簡単/より明確な方法はあります各時間領域の各空間領域内で?