这是我的代码和出现的问题:

R=0.0005;
lmda=500e-9;
f=1;
m=150;
xm=6000*lmda*f;
def=1e-5;
I0=1;
[x,y]=meshgrid(-xm:def:xm);
t=2*pi*R*sqrt(x^2+y^2)/(lmda*sqrt(x^2+y^2+f^2));
I=4*I0*besseli(1,t).^2/(t.^2);
imshow(x,y,I)
Warning: Matrix is singular to working precision.


出现的图片:



请大侠帮忙修改一下,我知道可能与自变量的范围有关,但是不知道怎么改进!摆脱了。。。