MATLAB光学仿真——楔板

5 用平行光 点光源 准直镜 照射如图 7 所示不同形状楔板 观察干涉图样 1 程序代码如下 Clear 清除原有变量 Lambda 500 1e 9 设定波长 theta 1 1e 6 设定倾角 ni 500 微元个数 ds linspace 0 0 0005 ni 竖直方向分割 for k 1 ni y k ds k sin theta 水平方向的对应坐标 Delta 2 ds k Lambda 2 对应点的光程差 Phi 2 pi Delta Lambda 对应点的相位差 B k 4 cos Phi 2 2 光强与相位差相联系 end figure gcf 以下为形成图像过程 将光强与图像辉度 相 对应 NCLevels 250 Br B 4 0 NCLevels image 0 y Br colormap gray NCLevels title 二维强度分布 2 程序代码如下 clear 清除原有变量 Lambda 500 1e 9 设置波长 r 5 设置介质的尺寸 t 3 ni 1000 分割份数 y linspace 0 005 0 005 ni y方向上进行分割 for k 1 ni h k t r sqrt r 2 y k 2 对应点的介质厚度 Delta 2 h k Lambda 2 对应点的光程差 Phi 2 pi Delta Lambda 对应点的相位差 B k 4 cos Phi 2 2 对应点的光强 end figure gcf 以下为形成图像过程 将光强与图像辉度相 对应 NCLevels 250 Br B 4 0 NCLevels image 0 y Br colormap gray NCLevels title 二维强度分布 3 程序代码如下 clear 清除原有变量 lambda 500 1e 9 R 5 H 2 设置波长和介质尺寸 ni 2000 设置分割份数 x linspace 0 005 0 005 ni 在x y方向上进行分割 y linspace 0 005 0 005 ni for i 1 ni for j 1 ni r x i 2 y j 2 点距介质中心的直线距离 delta 2 H 2 R sqrt R 2 r lambda 2 相应点的光程差 phi 2 pi delta lambda 相应点的相位差 B i j 4 cos phi 2 2 相应点的光强 end end figure gcf 以下为形成图像过程 将光强与 图像辉度相对应 NCLevels 250 Br B 4 0 NCLevels image x y Br colormap gray NCLevels title 二维强度分布 4 程序代码如下 clear lambda 500 1e 9 R 5 H 2 ni 2000 T 0 005 omega 2 pi T a2 axes Position 0 1 0 1 图片在figure中的位置 前两个代表x方向的 坐标 后两个是y方向的坐标 x linspace T T ni y linspace T T ni for i 1 ni for j 1 ni r sqrt x i 2 y j 2 delta 2 H sin omega r 1e 6 phi 2 pi delta lambda B i j 4 cos phi 2 2 end end figure gcf NCL