实验6-干涉的Matlab模拟.docVIP

《MATLAB在光信息处理中的应用》课程上机实验（三峡大学2013年） 实验6 干涉的Matlab模拟 一、实验目的：掌握双缝干涉、牛顿环的matlab 模拟。 二、实验内容： 折射率 折射率n=1.4, 厚度为5 1、 D=1m;d=2mm ?当在一个缝的位置放置一个折射率n=1.4, 厚度为5mm的物质，不考虑损耗，请画出此时的双缝干涉的图样（与课件相似的图） clear all clc D=1; %设置双缝到光屏的距离（1000mm） Lambda = 550e-009; %设置光线波长(550nm) d=2e-3; N = 100; %定义变量 x = linspace(-5,5,N); %设置图像精度 y = linspace(-5,5,N); %设置图像精度 for i = 1:N %定义变量确定横向点的坐标 I(i) = 4*cos(pi/Lambda*(d*x(i)/D-0.002)).^2; %光强分布公式 end subplot(2,1,1); %给窗口分栏,定义光强分布曲线的位置 plot(x,I); %画出光强分布曲线 Xlabel(x/mm); %标出横坐标（单位mm） Ylabel(I(x)); %标出纵坐标单位 title(光强分布曲线); %标出标题 axis([-5 5 0 4]); %标出坐标分度值 A=255; %定义干涉图像灰度分度值 B=(I/4*A); %换算 subplot(2,1,2); %定义干涉图样位置 image(x,y,B); %画出干涉图像原形 colormap(gray(A)); %将图样转化为灰度图 Xlabel(x/mm); %标出横坐标（单位mm） Ylabel(y/mm); %标出纵坐标（单位mm） title(干涉图样); %标出标题 axis([-5 5 0 5]); %标出坐标分度值 2、请画出d, D, 波长分别变化（假如对不同波长折射率一样）时候的双缝干涉的强度分布的2维图形。 clear all clc D=1; %设置双缝到光屏的距离（1000mm） Lambda = 500e-009; %设置光线波长(500nm) N = 800; %定义变量 d=linspace(5e-3,6e-3,N); x =2 ; for i = 1:N %定义变量确定横向点的坐标 I(i) = 4*cos(pi*d(i)*x/(Lambda*D)).^2; %光强分布公式 end subplot(3,1,1); plot(d,I); %画出光强分布曲线 Xlabel(d/m); %标出横坐标（单位m） Ylabel(I(x)); %标出纵坐标单位 title(随d变化的光强分布曲线); %标出标题 axis([0.005 0.006 0 4]); %标出坐标分度值 D=linspace(1,2,N); d=2e-3; x =2 ; for i = 1:N %定义变量确定横向点的坐标 I(i) = 4*cos(pi*d*x/(Lambda*D(i))).^2; %光强分布公式 end subplot(3,1,2); plot(D,I); %画出光强分布曲线 Xlabel(D/m); %标出横坐标（单位m） Ylabel(I(x)); %标出纵坐标单位 title(随D变化的光强分布曲线); %标出标题 axis([1 1.5 0 4]); %标出坐标分度值 Lambda=linspace(400e-009,700e-009,N); D=1; x=2; d=2e-3; for i =