你好,代码供参考,参数都可以调:
clc;clear
r = linspace(0,10,100);
theta = linspace(0,2*pi,100);
[R,T] = meshgrid(r,theta);
rho = 0.5;
sigma1 = 3;
sigma2 = 3;
mu1 = 10;
mu2 = 10;
x = mu1 + R.*cos(T);
y = mu2 + R.*sin(T);
f = 1/(2*pi*sigma1*sigma2*sqrt(1-rho^2))*...
exp(-1/(2*(1-rho^2))*((x-mu1).^2/sigma1^2-...
2*rho*(x-mu1).*(y-mu2)/(sigma1*sigma2)+(y-mu2).^2/sigma2^2));
mesh(x,y,f)
效果:
有帮助还望给个采纳哟