杨氏干涉的模拟matlab源码
% 杨氏干涉的模拟
close all
figure('position',[217 266 694 244]);
set(gcf,'DoubleBuffer','on');
axes('position',[0.12,0.08,0.6,0.8]);hold on;
rectangle('position',[-1,1.1,0.1,3],'FaceColor',[0.1,0.3,0.4]);
rectangle('position',[-1,-4,0.1,3],'FaceColor',[0.1,0.3,0.4]);
rectangle('position',[-1,-0.8,0.1,1.6],'FaceColor',[0.1,0.3,0.4]);
axis([-4,12,-5,5]);
ha=plot([-3.8,-3.8],[0.9,0.9],'r','linewidth',3);
hb=plot([-3.8,-3.8],[-0.9,-0.9],'r','linewidth',3);
title('Double slit interference \copyright zjliu')
L=-3.8;
rectangle('position',[11,-4,0.1,8],'FaceColor',[0,0,0]);
plot([-1,4.1],[-3.5,-3.5],'k');
plot([5.9,11],[-3.5,-3.5],'k');
plot([-0.7,-1,-0.7],[-3.3,-3.5,-3.7],'k');
plot([10.7,11,10.7],[-3.3,-3.5,-3.7],'k');
text(4.8,-3.5,'3m');
plot([-1.5,-1.5],[0.9,0.35],'k');
plot([-1.5,-1.5],[-0.9,-0.35],'k');
plot([-1.65,-1.5,-1.3],[0.6,0.9,0.6],'k');
plot([-1.6,-1.5,-1.3],-[0.6,0.9,0.6],'k');
text(-2.1,0,'2mm');
while L<=-0.95;
L=L+0.05;
set(ha,'xdata',[-3.8,L]);
set(hb,'xdata',[-3.8,L]);
pause(0.05);
end
xas=-0.95;yas=0.9;
xbs=-0.95;ybs=-0.9;
po=[-3.5:3.5];
ka=(po-yas)*i+(11+0.95);ka=ka./abs(ka);
kb=(po-ybs)*i+(11+0.95);kb=kb./abs(kb);
for n=1:8;
ah(n)=plot(xas,yas,'r');
bh(n)=plot(xbs,ybs,'r');
end
r=0;
za=xas+yas*i;zb=xbs+ybs*i;
Le=sqrt(11.95^2+3.5^2);
while r<=Le;
for n=1:8;
Ta=za+r*ka(n);
Ya=imag(Ta);
if abs(imag(Ta))>=abs(po(n));
Ya=po(n);
end
set(ah(n),'xdata',[xas,min(11,real(Ta))],'ydata',[yas,Ya]);
Tb=zb+r*kb(n);
Yb=imag(Tb);
if abs(imag(Tb))>=abs(po(n));
Yb=po(n);
end
set(bh(n),'xdata',[xbs,min(11,real(Tb))],'ydata',[ybs,Yb]);
end
r=r+0.05;
pause(0.05);
end
s=meshgrid(linspace(4,-4,300))';
De=abs((s-1)*i+3000)-abs((s+1)*i+3000);
lambda=0.6328e-3;
de_A=De/lambda*pi*2;
It=1-cos(de_A);
axes('position',[0.72,0.16,0.26,0.64]);
cc=cat(3,It/2,zeros(size(It)),zeros(size(It)));
imshow(cc,[]);
xlabel('\copyright zjliu')
相关文章
- MATLAB中readtimetable函数用法
- Matlab中functions函数的使用
- Matlab中句柄类析构函数
- [matlab] 8.蚁群算法解决TSP问题
- [matlab] 4.M函数
- 基于matlab的车牌识别
- 【MATLAB教程案例88】通过matlab模拟光的干涉现象
- 【MATLAB教程案例86】通过matlab实现lorenz混沌系统
- 【MATLAB教程案例43】初识matlab中三维数据处理相关函数——plot3,surf,mesh,isosurface,slice等
- 【MATLAB教程案例29】基于Baker映射和Logistic混沌序列的图像加解密matlab实现
- 【MATLAB教程案例17】基于NSGAII多目标优化算法的matlab仿真及应用
- 【MATLAB教程案例14】基于ACO蚁群优化算法的函数极值计算matlab仿真及其他应用
- 基于matlab的精馏塔作业模拟仿真
- 基于MATLAB的指纹识别算法仿真实现
- 飞机器减震器的matlab建模和仿真模拟
- 【RFID】基于MATLAB的RFID 系统的空中接口过程以及防碰撞算法仿真
- 【人脸识别】基于matlab从图像中检测人脸部分源码程序
- 使用matlab模拟、检验和估计泊松过程
- ML算法的matlab源码
- 图像压缩编码码matlab实现——JSjpeg
- 【Matlab小问题】matlab启动时出现Warning: Name is nonexistent or not a directory
- 基于Matlab的Robotics Toolbox工具箱的机器人仿真函数介绍(运动学)
- MATLAB安装机器人工具箱 Matlab_Robotic_Toolbox_v10.2
- 【Matlab 六自由度机器人】运动学逆解(附MATLAB机器人逆解代码)
- 【Matlab算法】粒子群算法求解一维非线性函数问题(附MATLAB代码)
- 【Matlab算法】MATLAB求解背包问题(附MATLAB代码)
- 【Matlab】系统预定义变量——元胞数组与结构体
- 用MATLAB一直画直线竟然得到了简单的禅绕画,论一直画直线的人有多无聊
- matlab常用函数——数学函数
- matlab常用函数——数据类型函数