基于matlab干涉系统仿真 - 图文(5)

[xs, ys] = meshgrid(xs, ys); %生成网格采样点 n*n矩阵 r1=sqrt(xs.^2+ys.^2+D^2); %光屏上点（xs,ys）距光源1距离r1 E1=A1./r1.*exp(1i*r1*2*pi/lam);%点光源1在接受屏上复振幅E1 E2=A2; %平面波2在接受屏上复振幅E2 E=E1+E2; %复振幅叠加 I=abs(E).^2; %光强

pcolor(xs, ys, I); %生成干涉图样 shading flat; colormap gray

ys：-0.002—0.002

光源振幅：A1=A2=1；（可调）

z

② 干涉仿真图样：

2）斜入射（与接收屏法线方向60°）

①模型参数：

光线波长：lam=500nm； 点光源到接收屏距离：D=1mm； 接收屏范围：xs：- 0.00051—0.00051

②matlab代码： clear;

lam=500e-9; %设定波长lam（500纳米） D=10e-3; %点光源孔到接收屏距离 A1=1; %两光源均为单位幅度 A2=A1; xm=0.00051;

ym=xm; %接受屏的范围ym，xm n=1001;

ys：- 0.00051—0.00051

光源振幅：A1=A2=1；（可调）

xs=linspace(-xm,xm,n); %用线性采样法生成两个一位数组

%xs，ys（n为总点数）

ys=linspace(-ym,ym,n);

[xs, ys] = meshgrid(xs, ys); %生成网格采样点 n*n矩阵 r1=sqrt(xs.^2+ys.^2+D^2); %光屏上点（xs,ys）距光源1距离r1 E1=A1./r1.*exp(1i*r1*2*pi/lam); %点光源1在接受屏上复振幅E1 E2=A2*exp(-1i*xs*cos(pi/3)*2*pi/lam); %平面波2在接受屏上复振幅E2 E=E1+E2; %复振幅叠加 I=abs(E).^2; %光强

pcolor(xs, ys, I); %生成干涉图样 shading flat; colormap gray

干涉仿真图样：

改变平面光波入射方向：

zz

5、用平行光（点光源+准直镜）照射如图7所示不同形状楔板，观察干涉图样

不同形状楔板

1） 斜面楔板

当薄膜为空气时，其折射率n=1，而且上下表面BB与AA之间的夹角又很小，使光线几乎垂直入射，则C在BB表面上，光线11'和22'的光程差为：

（1）

式中，是因为光线由光疏媒质入射到光密媒质在AA界面上反射时，发生的半波损失引起的光程差。

式（1）只与薄膜的厚度d和光波波长有关。 当光程差：

（k=1,2,3……）即：

当光程差：

（k=1,2,3……）即：

时产生亮条纹。 时产生暗条纹；

从以上亮纹或暗纹公式易导出，从一个条纹过渡到另一个条纹，平板的厚度变化为间距，且

。对应光程差变化为，楔板的楔角 又

所以

。

。e为条纹

①matlab代码： clear;

lam=500e-9; %设定波长lam（500纳米） A1=1; %两光源均为单位幅度 A2=A1; %默认为相等 theta=pi/90; %楔板倾角theta n=1; %楔板折射率n xm=10e-5;

ym=xm; %接受屏的范围xm，ym n=1001;

xs=linspace(0,xm,n); %用线性采样法生成两个一维数组

%xs，ys（n为总点数）

ys=linspace(0,ym,n); [xs, ys] = meshgrid(xs, ys);

deta=2*n*xs*tan(theta); %光程差 phi=2*pi*deta/lam; %相位差

E1=A1; %平面波1在接受屏上复振幅E1 E2=A2.*exp(-1i*phi); %平面波2在接受屏上复振幅E2 E=E1+E2; %复振幅叠加 I=abs(E).^2; %光强

pcolor(xs, ys, I); %生成干涉图样 shading flat; colormap gray ③ 干涉图样