我想计算三维中两个正态分布(高斯)之间的体积。
我得到了两种不同的正态分布:1)一种高斯分布,适用于x,y,z数据(红色),2)一种模拟高斯分布(蓝色)。我可以绘制两个函数,但我不知道如何计算它们之间的差异。
所以我想要的是红色和蓝色的音量差。这有一个物理背景:卷对应于电源,我想知道电源的区别。
我考虑过整合,但我对此并不熟悉。
谢谢您。
from matplotlib import pyplot;
from pylab import genfromtxt;
import matplotlib.pyplot as plt
import numpy as np
from mpl_toolkits.mplot3d import Axes3D
########### FITTED GAUSSIAN (red) ##############
# Load file into mat0
mat0 = genfromtxt("0005.map");
#PLOT FIGURE
fig = plt.figure(figsize=(20,10))
ax = plt.axes(projection='3d')
#Define Gaussian function
def twoD_Gauss((x,y),amplitude,x0,y0,sigma_x,sigma_y,offset):
x0=float(x0)
y0=float(y0)
return offset + amplitude*np.exp(-(((x-x0)**(2)/(2*sigma_x**(2))) + ((