本文基于光学原理,用Matlab实现菲涅尔直边衍射现象仿真。
Project Code
% 作者:ZQJ% 日期:.5.18 星期二%***********************菲尼尔直边衍射仿真**************************clear,clc,close all;% 本例输入:Lambda = 632.8; L = 0.2;Lambda = input('输入光的波长(单位为nm):');Lambda = Lambda * 1e-9; %波长单位转为:mL = input('直边挡板到接收屏的距离(单位为m):'); x = -0.001:0.000003:0.001;w = sqrt(2/Lambda/L).*x;I0 = 1; %定义标准光强C = linspace(0,0,length(x)); %预定义数组S = linspace(0,0,length(x)); % 进行菲涅尔积分运算******************************for n=1:length(x)C(n) = integral(@cos_fun,0,w(n));%菲涅尔积分实部S(n) = integral(@sin_fun,0,w(n));%菲涅尔积分虚部I = I0/2.*((C+0.5).^2+(S+0.5).^2);%光强数值F = C + 1i.*S;end% 光强分布曲线************************************plot(x,I);xlabel('x'),ylabel('I(x)');title('Intensity');% 考纽卷线************************************figure, plot(F);title('Cornu Spiral');% 光强分布************************************figure;imagesc(x,max(x),I);colormap hot;title('Intensity Distribution');% 函数调用**********************************************function Z_Fresnel_integral = cos_fun(t)% 函数功能:菲涅尔积分余弦函数Z_Fresnel_integral = cos(0.5*pi*t.^2);endfunction Z_Fresnel_integral = sin_fun(t)% 函数功能:菲涅尔积分余弦函数Z_Fresnel_integral = sin(0.5*pi*t.^2);end
仿真结果图:
光强函数曲线图:考纽卷线图:
光强分布图:
专栏内容供作者本人或大家学习使用,多多指教~