R3=200;
F=zeros(2*R3+1,2*R3+1);
sigma=R3/3;
r=0;
for x=1:2*R3+1
for y=1:2*R3+1
r=(x-R3)*(x-R3)+(y-R3)*(y-R3);
F(x,y)=exp(-r/(2*sigma*sigma));
K=K+F(x,y);
end
end
F=zeros(2*R3+1,2*R3+1);
sigma=R3/3;
r=0;
for x=1:2*R3+1
for y=1:2*R3+1
r=(x-R3)*(x-R3)+(y-R3)*(y-R3);
F(x,y)=exp(-r/(2*sigma*sigma));
K=K+F(x,y);
end
end
浙公网安备 33010602011771号