码迷,mamicode.com
首页 > 编程语言 > 详细

【python-plt】二元正态密度函数图像

时间:2020-03-06 21:47:58      阅读:123      评论:0      收藏:0      [点我收藏+]

标签:sqrt   false   pre   def   lse   width   src   numpy   att   

import numpy as np
import matplotlib.pyplot as plt
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter

x=np.arange(-10,10,0.5)
y=x
X,Y=np.meshgrid(x,y)

def norm_d(x1,x2,mu1,mu2,s1,s2,rho):
    num=((1)/(2*np.pi*s1*s2*np.sqrt(1-rho**2)))
    A=((x1-mu1)**2)/(s1**2)
    B=2*rho*(((x1-mu1)*(x2-mu2))/(s1*s2))
    C=((x2-mu2)**2)/(s2**2)
    D=-1/(2*(1-rho**2))*(A-B+C)
    pdf=num*np.exp(D)
    return 2dnm

R=norm_d(X,Y,0,0,1,3,0.75)
fig = plt.figure()
ax = fig.gca(projection='3d')
surf = ax.plot_surface(X, Y, R, cmap=cm.coolwarm,linewidth=0, antialiased=False)
plt.show()

\(\sigma_1^2=1,\sigma_2^2=3,\rho=0.75\)

技术图片

\(\sigma_1^2=1,\sigma_2^2=1,\rho=0\)

技术图片

【python-plt】二元正态密度函数图像

标签:sqrt   false   pre   def   lse   width   src   numpy   att   

原文地址:https://www.cnblogs.com/rrrrraulista/p/12431004.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!