1 clear all;
 2 close all;
 3 clc;
 4 
 5 randn('seed',0);
 6 %%一维高斯函数
 7 mu=0;
 8 sigma=1;
 9 x=-6:0.1:6;
10 y=normpdf(x,mu,sigma);
11 plot(x,y);
12 figure;
13 
14 %%二维或多维高斯函数
15 mu=[0 0];
16 sigma=[0.3 0;0 0.35];
17 [x y]=meshgrid(linspace(-8,8,80)',linspace(-8,8,80)');
18 X=[x(:) y(:)];
19 z=mvnpdf(X,mu,sigma);
20 surf(x,y,reshape(z,80,80));
21 hold on;
22 %再生成一个
23 mu=[4 0];
24 sigma=[1.2 0;0 1.85];
25 [x y]=meshgrid(linspace(-8,8,80)',linspace(-8,8,80)');
26 X=[x(:) y(:)];
27 z=mvnpdf(X,mu,sigma);
28 surf(x,y,reshape(z,80,80));