|
%矩孔衍射的解析计算2 K/ a$ C- _# Y1 ^8 \# o
clear;
. o9 ?% a$ k. h# |; SR=0.1;
2 ?3 Y- i# L1 e% `2 I! n$ X& Jlambda=1.064e-3;
9 h, H7 f7 a3 T R: Fk=2*pi/lambda;* J# d+ j5 I+ |" Q, i! l9 B) D
z=1.0e3;) o- c/ p2 u t1 V/ I
xmax=8*1.22*lambda/2/R*z;# h& |+ L8 J0 c3 v- u
x=linspace(-xmax,xmax,200);
. J( h9 u" I0 b' T* my=x;1 t$ J6 C4 `; ]
[x,y]=meshgrid(x,y);$ t0 \% Z7 o1 }" h. y5 `
IF=sin(k*x*R/z).ˆ2.*sin(k*y*R/z).ˆ2./(k*x/2/z).ˆ2./(k*y/2/z).ˆ2/lambdaˆ2/zˆ2;
% Q# z I' ?( m2 Y. Q6 W2 tsurf(x,y,IF.ˆ(1/2));
h& j6 ]0 u0 V- {/ o- B& \4 |' Ecolormap(’hot’)
, t9 v- h, u- R/ y# H) } Haxis equal
1 x5 V4 h9 H8 B, Sshading interp;
$ o" Y+ l5 _/ d" p2 f) E$ | |
|