基于Matlab的迈克尔逊干涉仪仿真模拟

clear;close all;clc; m=1; % 定义单位米 cm=1e-2*m; % 定义单位厘米 mm=1e-3*m; % 定义单位毫米 nm=1e-9*m; % 定义单位纳米 rad=1; % 定义弧度单位 mrad=1e-3*rad; lambda=632.8*nm; % 定义波长 size=10*mm; N=500; R=3*mm; % 定义孔径 z1=3*cm; % 定义距离Z1 z4=5*cm; % 定义距离Z4 RBS=0.5; % 定义分光比 ty=0.0*mrad; % 定义反射镜2的y向倾斜角度 tx=5*mrad; % 定义反射镜2的x向倾斜角度 f=40*cm; % 定义透镜焦距 z3=7*cm+lambda/4; % 定义距离Z3 z2=13*cm; % 定义距离Z2 %透镜光束会聚模拟 F=LPBegin(size,lambda,N); F=LPGaussAperture(R,0,0,1,F); % 高斯光束模拟 F=LPLens(f,0,0,F); % 透镜 %传播至分束镜 F=LPForvard(z1,F); %经过分束镜后传播至反射镜2 F2=LPIntAttenuator(1-RBS,F); F2=LPForvard(z3,F2); %反射镜2引入倾斜，并反射回至分束镜 F2=LPTilt(tx,ty,F2); F2=LPForvard(z3,F2); F2=LPIntAttenuator(RBS,F2); %第二束光(参考光模拟)传播至反射镜1 F10=LPIntAttenuator(RBS,F); F1=LPForvard(z2*2,F10); F1=LPIntAttenuator(1-RBS,F1); %将两个光束重新组合并传播到屏幕： F=LPBeamMix(F1,F2); F=LPInterpol(size/3,N,0,0,0,1,F); F=LPForvard(z4,F); I=LPIntensity(1,F); %在屏幕上绘制强度（位图bitmap） figure(1);image(I*80);colormap(gray);axis off;axis square;