工程光学等厚干涉Matlab仿真

1、第五题四个图均为等厚干涉。因此相位差，光强的式子均相同，只要更改光程差，并设置相关的参数即可。.二,2,、I =4cos 侦）1. 平面楔板:=2nh 2(2) 参数设置波长 =500 nm 楔角 =0.1rad 折射率n = 1(3) Matlab仿真程序2 -clc:clear all: close all:3 -Lambda=500?'9:4 -that a=1e-1;5 -n=l:6 -vect orl=linspace (-0. 0L 0. 0L 400);7 -vect or2=linspace Ol?Olj 400):8 -x,y=meshEridfvectorlj ve

2、ctor?):9 -Az+t an (theta):10 -Dell a=2*n*A+Lambda/2 :11 -Phi=2*pi*D?lta- /Lambda:12 -I=4*cos<Phi/2). n213 -mesh(x,yj I);14 - pcolor (Xj yj I) ;shading interp%设置波长*设置楔角藕设置折射率鬼设置*方向向璧空间鬼设置y方向向量空间%构成%y二维向崖空间%计算光程差瓶十箕相位差普计算光强%根据光强绘制干于条纹%在网格片内采用颜色插值处理 坂设置颜色为灰度15 -colormap gray;(4)仿真结果2. 柱面楔板(1) 光程差分析2

3、22722y x = R = y -、R - xA = y -h = . R2 -x2 -h:=2nhA -2(2) 参数设置波长 =500 nm平行平板距柱面圆心透镜距离h=0.5mm柱面透镜半径R=4.5mm折射率n = 1(3) Matlab仿真程序2 - clc;clear all; closs all;3 -Lanbda=500e-9;4 -R=4. 5;5 -h=0, 5;6 n=l;7 vectorl=linspace r-0. Olj 0. 01, 400);8 -vector2=linEpace(-0, 0L 0, 01, 400);9 一xjyleshgridtvector

4、lj vector2):10 -A=0T2r. '2). "0.5-h:11 DeltaF2*A*n+Lanbda/2 :12 Phi=2 *pa. /Lambda:13 I=4#cos(Phi/2), 214 -y, I);15 一 pcolcr(x, y, I):shading int erp%设置波长靳殳置楔角畏设置平行平板到柱面圄心的距离强设置折射率强设置x方向向重空间强设置y方向向重空间畏构成二维向量空间%计尊光程差强计算相位差贤计箕光强%根据光强绘制干涉条绞%在网格片内采用颜色插值处理%设置颜色为切度16 一 colormap gray;(4)仿真结果(2)参数设

5、置ho -2 r2R3. 球面楔板(1) 光程差分析222R = (R -h0) r =A = R - h0 - h:=2nhA -2(3)Matlab波长 =500 nm平行平板距球面球心透镜距离h=0.5mm球面透镜半径R=4.5 mm折射率n =1仿真程序-210 -11 -12 -13 -14 15 -16 17 -clc;clear all; close all:Lajnbda=500e9 :R=4.5;h=Q.5;n=l. 5vectorl=linspace (-0, 01 j 0, Olj 400); vector2=linspace (-0. 01 a Q, Olj 400)；

6、 be, yleshEridtvctcrlj vect or2); r=abs'x+i*y);A=R-r. *2/R/2-h'D e 11 2 *A*n+L ajub da/ 2 :Phi=2*pi#D?lta, /Lambda;I=4*cos(Phi/2)-"2neshI);pcolor(x, y, I):shading interp colormap gray:$设置波长%设置楔角%设置平行平板到球面球心的距离%设置折射翠%设置x方向向里空间%设置y方向向重空间%构成xj y_维向量空间吊将平面宣角坐标转换为极坐标第计算光程差%计算相位差计算光强第根据光强绘制干涉

7、条纹那在网格片内采用颜色插值处理贤将颜色设置为独度(4)仿真结果4.不规则柱面楔板(1) 光程差分析以柱面楔板干涉为基础，将以常数R为半径的柱面更改为 R = (100x)2+0.1x的二(2)次曲线，即可满足题目要求。参数设置(3 )波长 =500nm平行平板距柱面圆心透镜距离h=0.5mm折射率n =1Matlab仿真程序2 -tic:clear all: close all:3 -Lambda二5。0巳-9：4 -h=0- 5:5 -n=l;6 -vectorllinspace (-0» 0L 0* 01? 400)7 - vector2=luispace (-0» OL Q. Olj 400):S -x,y =nieshgridfvectorvectar2);9 -R=(IOO*k), "2+0.l*x;10 -A=(R'2-x,11 -Delta-2+A4n+Lajhbda/2 :12 -Phi=2ipi*Delta. /Lambda:13 -I=4*cos(Phi/2). *214 -nesh(K, yt I);15 - p co Lor (x, Yj I) ； shading interp16 -color