Matlab在基础物理学中的应用举例4.docx

Matlab在基础物理学中的应用举例(四)（光学：牛顿环的模拟）2011-02-13 13:30:18| 分类： Matlab Sky | 标签：matlab 光学 |字号大中小 订阅 下面是模拟牛顿环的图案，给出了一个用户交换操作界面的实现程序：close all;clear;figure(Position,90 164 873 483); % 设置图形窗口位置L=632.8;R=5;H=0.005; % 设置默认参数值：波长L、半径R和厚度Ha1=axes(Position,0.83,0.3,0.15,0.4); % 设置坐标轴位置，在其上绘制简易的原理图hold on;axis(0,1,0,1); % 设置坐标轴范围plot(0.25,0.25,0.5,0.8,r); % 画竖直线，入射光线plot(0.5,0.5,0.5,0.8,r); % 画竖直线，入射光线plot(0.75,0.75,0.5,0.8,r); % 画竖直线，入射光线fill(0.22,0.25,0.27,0.5,0.44,0.5,r); % 画出箭头fill(0.47,0.5,0.52,0.5,0.44,0.5,r); % 画出箭头fill(0.72,0.75,0.77,0.5,0.44,0.5,r); % 画出箭头z=1.8*exp(i*(linspace(-0.21,0.21,30)-pi/2)+2.1*i+0.5; % 生成平凸透镜的轮廓fill(real(z),imag(z),w); % 利用fill函数画平凸透镜的轮廓rectangle(Position,0.1,0.18,0.8,0.12); % 画玻璃板底面set(gca,xtick,ytick,box,on); % 删去坐标轴刻度a2=axes(Position,0.4,0.16,0.4,0.7); % 生成第二个坐标轴，在其上画牛顿环x,y=meshgrid(linspace(-0.005,0.005,200); % 生成坐标矩阵r2=(x.2+y.2); % 坐标矩阵各点到中心的距离Di=2*H+2*(R-sqrt(R2-r2)*1e9/L; % 计算光程差In=abs(cos(Di*pi*2); % 得到牛顿环上的光强数值cr=abs(L-560)/200; % 红色的分量值cg=1-cr; % 绿色的分量值cb=abs(L-600)/240; % 蓝色的分量值Ik(:,:,1)=In*cr; % 转化为RGB三基色表示颜色矩阵数值（红色）Ik(:,:,2)=In*cg; % 转化为RGB三基色表示颜色矩阵数值（绿色）Ik(:,:,3)=In*cb; % 转化为RGB三基色表示颜色矩阵数值（蓝色）Pc=imshow(Ik,); % 绘制缺省参数下的数值title(the pattern of Newtons rings,fontsize,12); % 添加图题% 实时显示波长的文本框Lt=uicontrol(gcf,style,text,. unit,normalized,position,0.06,0.86,0.21,0.06,. BackgroundColor,0.7*1,1,1,ForegroundColor,0.8,0.1,0.9,. string,wavelength: 632.8nm,fontsize,12,fontname,times new roman); s1=uicontrol(gcf,style,slider,. unit,normalized,position,0.06,0.76,0.21,0.04,. BackgroundColor,0.7*1,1,1,ForegroundColor,0.1,0.1,0.9,. SliderStep,0.01,0.01,value,(632.8-360)/400,.callback,L=get(s1,value)*400+360;,.set(Lt,string,wavelength: ,num2str(L/10),nm);,.Di=2*H+2*(R-sqrt(R2-r2)*1e9/L;,.In=abs(cos(Di*pi*2);cr=abs(L-560)/200;cg=1-cr;,.cb=abs(L-600)/240;Ik(:,:,1)=In*cr;Ik(:,:,2)=In*cg;,.Ik(:,:,3)=In*cb;set(Pc,CData,Ik);); % 利用滑动条改变波长数值uicontrol(gcf,style,text,. unit,normalized,position,0.04,0.81,0.08,0.04,. BackgroundColor,0.8*1,1,1,ForegroundColor,0.1,0.1,0.9,. string,360,fontsize,12,fontname,times new roman); % 显示最小波长数值uicontrol(gcf,style,text,. unit,normalized,position,0.22,0.81,0.08,0.04,. BackgroundColor,0.8*1,1,1,ForegroundColor,0.1,0.1,0.9,. string,760,fontsize,12,fontname,times new roman); % 显示最大波长数值Rt=uicontrol(gcf,style,text,. unit,normalized,position,0.06,0.66,0.23,0.06,. BackgroundColor,0.7*1,1,1,ForegroundColor,0.8,0.1,0.9,. string,radii:,fontsize,12,fontname,times new roman); %实时显示半径数值的文本框s2=uicontrol(gcf,style,slider,. unit,normalized,position,0.06,0.56,0.21,0.04,. BackgroundColor,0.7*1,1,1,ForegroundColor,0.1,0.1,0.9,. SliderStep,0.01,0.01,.callback,R=get(s2,value)*7+5;,.set(Rt,string,radii: 5m,num2str(R),m);,.Di=2*H+2*(R-sqrt(R2-r2)*1e9/L;,.In=abs(cos(Di*pi*2);cr=abs(L-560)/200;cg=1-cr;,.cb=abs(L-600)/240;Ik(:,:,1)=In*cr;Ik(:,:,2)=In*cg;,.Ik(:,:,3)=In*cb;set(Pc,CData,Ik);); % 通过滑动条改变半径数值uicontrol(gcf,style,text,. unit,normalized,position,0.04,0.61,0.08,0.04,. BackgroundColor,0.8*1,1,1,ForegroundColor,0.1,0.1,0.9,. string,5,fontsize,12,fontname,times new roman); % 显示最小半径uicontrol(gcf,style,text,. unit,normalized,position,0.22,0.61,0.08,0.04,. BackgroundColor,0.8*1,1,1,ForegroundColor,0.1,0.1,0.9,. string,12,fontsize,12,fontname,times new roman); % 显示最大半径% 实时显示厚度的文本框Ht=uicontrol(gcf,style,text,. unit,normalized,position,0.06,0.46,0.23,0.06,. BackgroundColor,0.7*1,1,1,ForegroundColor,0.8,0.1,0.9,. string,thickness: 5nm,fontsize,12,fontname,times new roman); s3=uicontrol(gcf,style,slider,. unit,normalized,position,0.06,0.36,0.21,0.04,. BackgroundColor,0.7*1,1,1,ForegroundColor,0.1,0.1,0.9,. SliderStep,0.01,0.01,value,0.05,.callback,H=get(s3,value)*0.01;,.set(Ht,string,thickness: ,num2str(H),nm);,.Di=2*H+2*(R-sqrt(R2-r2)*1e9/L;,.In=abs(cos(Di*pi*2);cr=abs(L-560)/200;cg=1-cr;,.cb=abs(L-600)/240;Ik(:,:,1)=In*cr;Ik(:,:,2)=In*cg;,.Ik(:,:,3)=In*cb;set(Pc,CData,Ik);); % 通过滑动条改变厚度数值uicontrol(gcf,style,text,. unit,normalized,position,0.04,0.41,0.08,0.04,. BackgroundColor,0.8*1,1,1,ForegroundColor,0.1,0.1,0.9,. string,0,fontsize,12,fontname,times new