Ex1=60;

Ex2=1.7;

En1 =6;

En2=0.5;

H1 =0.1;

H2=0.01;

N=10000;

X1 (1:N)=zeros(1，N);

X2 (1∶N)=zeros(1，N);

En11 (1∶N)=zeros(1，N);

En22 (1∶N)=zeros(1，N);

Y（1∶N）=zeros （1，N）;

for i=1∶N

X1（1，i）=normrnd （Ex1，En1 ，1）;

X2（1，i）=normrnd （Ex2，En2 ，1）;

En11（1，i）=normrnd （En1，H1 ，1）;

En22（1，i）=normrnd （En2，H2 ，1）;

Y（1，i）=exp（-（X1（1，i）-Ex1）2/（2×En11（1，i）2）-（X21，i）-Ex2）2（2×En22（1，i））;

plot3（X1，X2，Y，'.'，'MarkerEdgeColor'，'k'，'MarkerSize'，4）;

xlabel（'体重 '，'fontsize'，25）;

ylabel（'身高 '，'fontsize'，25）;

zlabel（'隶属度 '，'fontsize'，25）;

title（'重约 60 kg，高约 1.7 m 的人 '，25）;

grid on;

end