gaussplot.m
上传用户:mozhenmi
上传日期:2008-02-18
资源大小:13k
文件大小:1k
源码类别:

其他小程序

开发平台:

Matlab

  1. function h=gaussplot(mu, Sigma, ncontours)
  2. % GAUSSDRAW Plot some contours of constant probability density of a 2D Gaussian centered on the mean
  3. % h=gaussdraw(mu, Sigma, ncontours)
  4. if ~exist('ncontours'), ncontours = -1; end
  5. % Modelled afer demgauss from netlab.
  6. ngrid = 50;
  7. sx = 1*Sigma(1,1);
  8. sy = 1*Sigma(2,2);
  9. Xmin = mu(1)-sx; Xmax = mu(1)+sx; 
  10. Ymin = mu(2)-sy; Ymax = mu(2)+sy; 
  11. Xvals = linspace(Xmin, Xmax, ngrid);
  12. Yvals = linspace(Ymin, Ymax, ngrid);
  13. [X1, X2] = meshgrid(Xvals, Yvals);
  14. XX = [X1(:), X2(:)];
  15. % the i'th row of XX is the (x,y) coord of the i'th point in the raster scan of the grid
  16. probs = gauss(mu, Sigma, XX);
  17. probs = reshape(probs, ngrid, ngrid);
  18. if ncontours==-1
  19.   [C,h]=contour(X1, X2, probs);
  20. else
  21.   [C,h]=contour(X1, X2, probs, ncontours, 'k');
  22. end