注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

MATLAB 干涉条纹  

2015-04-06 09:56:44|  分类: 编程 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
%##########################%
%名称:等厚干涉模拟                       %
%作者:程序猿__int64Ago                %
%时间:2011-12-16                        %
%##########################%
clear                                                %清内存
lamda=632.8e-9;                             %设定入射光波长
n=1;                                                %缝隙折射率
R=5;                                                %柱透镜半径
ym=5e-3;                                         %设定屏幕的范围,根据实际情况调整
xs=ym;                                             %x方向无变化,单一值
N=1000;                                          %条纹柔和度
ys=linspace(-ym,ym,N);                      %y轴上坐标离散
for i=1:N                                          %依次求干涉光强
    hm=R-sqrt(R^2-ys(i)^2);                %柱透镜与底面间隙
    delta=2*n*hm+lamda/2;              %光程差
    phi=2*pi*delta/lamda;                  %相位差
    I(i,:)=4*cos(phi/2).^2;                  %干涉光强
end                                                  %求光强结束
M=255;                                           %设定灰度等级
II=(I/4.0)*M;                                    %设定黑白条纹对比度
subplot(1,2,1);                                  %分割figure1×2,并且在第一块上作图
image(xs,ys,II);                                  %画干涉条纹
colormap(gray(M));                             %二维光强分布灰度图
title('光强分布灰度图')                        %列出标题
subplot(1,2,2);                                  %分割figure1×2,并且在第二块上作图
plot(I,ys);                                          %一维光强分布曲线图
title('光强分布曲线图')                        %列出标题

--------------------------------
%##########################%
%名称:多缝夫朗禾费衍射模拟         %
%作者:程序猿__int64Ago                %
%时间:2011-12-16                        %
%##########################%
clear                                                %清内存
lamda=500e-9;                                %设定入射光波长
a=10e-6;                                         %缝宽
d=25e-6;                                         %光栏常数
N0=6;                                             %缝数
D=5;                                                %缝到像面的距离(这个随意,几乎不影响)
ym=1*lamda*D/a;                            %设定屏幕的范围,根据实际情况调整
xs=ym;                                             %x方向无变化,单一值
N=1000;                                          %条纹柔和度(不能太大)
ys=linspace(-ym,ym,N);                      %y轴上坐标离散
for i=1:N                                          %依次求衍涉光强
    sintheta=ys(i)/D;                           %以下由书上公式
    alpha=pi*a*sintheta/lamda;
    beta=pi*d*sintheta/lamda;       
    I(i,:)=(sin(alpha)./alpha).^2.*...     
        (sin(N0*beta)./sin(beta)).^2;    %上一行续行
end                                                 %求光强结束
M=255;                                           %设定灰度等级
II=(I/max(I))*M;                                 %设定黑白条纹对比度
subplot(1,2,1);                                  %分割figure1×2,并且在第一块上作图
image(xs,ys,II);                                  %画干涉条纹
colormap(gray(M));                             %二维光强分布灰度图
title('光强分布灰度图')                        %列出标题
subplot(1,2,2);                                  %分割figure1×2,并且在第二块上作图
plot(I,ys);                                          %一维光强分布曲线图
title('光强分布曲线图')                        %列出标题

http://blog.csdn.net/int64ago/article/details/7076015
  评论这张
 
阅读(1283)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016