画二维正态分布函数图像matlab代码

function ewnormpdf(u1,u2,sigma1,sigma2,p) [x,y]=meshgrid(-5:0.1:5,-5:0.1:5); %X=[x,y];u=[0,0]; %B=[1,0;0,1]; %f=1/(2*pi*sqrt(det(B)))*exp(-1/2*(X-u)*inv(B)*(X-u)'); %f=(1/(2*pi))*exp((-1/2)*(x.^2+y.^2)); f=1/(2*pi*sigma1*sigma2*sqrt(1-p*p))*exp(-1/(2*(1-p*p))*(((x-u1).^2)./(sigma1*sigma1)-2*p*((x-u1)*(y-u2))./(sigma1*sigma2)+((y-u2).^2)./(sigma2*sigma2))); mesh(x,y,f)