cl;
m=31;
n=31;
img=zeros(m+1,n+1);
img=double(img);
pi=3.1415926;
sigma=10;
for i=-(m/2):m/2
for j=-(n/2):n/2
img(i+m/2+1,j+n/2+1)=(1/(2*pi*sigma*sigma))*exp(-(i*i+j*j)/(2*sigma*sigma));
end
end
img=mat2gray(img);
imshow(img);
imwrite(img,'pic.bmp');
%一个更好的实现,2012.12.18修改
h=5;
w=2;
[x y]=meshgrid(-w:w,-h:h);
sigma=5;
img = (1/(2*pi*sigma^2))*exp(-(x.^2+y.^2)/(2*sigma^2));
转载于:https://www.cnblogs.com/tiandsp/archive/2012/02/26/2368533.html
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)