|
%矩孔衍射的解析计算
8 U4 J( r" F' ?: ^1 n* Hclear;
) @3 ?3 N8 {) XR=0.1;, ~: t* c" Q% O& l4 Y
lambda=1.064e-3;
/ \$ y/ H# ?1 F4 G* E- sk=2*pi/lambda;
2 q% v' F/ f" z' w$ x7 oz=1.0e3;5 l1 w: G. C$ m$ [
xmax=8*1.22*lambda/2/R*z;) O9 b% d4 x0 V+ X: i6 w! @
x=linspace(-xmax,xmax,200); |; F2 k6 D( L: e
y=x;% W) ~- e. [9 x3 R& H- K
[x,y]=meshgrid(x,y);
" Q5 Q& X* i& P3 ]; B3 k' H9 Y RIF=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;% j- C8 g; C; T! ?. ?+ {; i6 I( o
surf(x,y,IF.ˆ(1/2));+ e1 U/ L9 a5 ?6 |- k8 o
colormap(’hot’) P5 l8 J6 M; X( l) w( |3 F1 D# [2 V
axis equal
# j: I( A& c) ?6 ^4 Cshading interp;
\& O1 F, `2 z; E% o" U |
|