菲涅尔公式工程光学.docx
《菲涅尔公式工程光学.docx》由会员分享,可在线阅读,更多相关《菲涅尔公式工程光学.docx(12页珍藏版)》请在冰豆网上搜索。
菲涅尔公式工程光学
实验名称菲涅尔公式的认识
一、实验目的:
加深理解菲涅尔公式,对给出的反射波或折射波与入射波振幅的相对变化进行分析,以及对相位变化进行分析。
二、实验原理:
任一方位振动的光矢量E都可以分解成互相垂直的两个分量
称平行于入射面振动的分量为光矢量的p分量,记为EP。
称垂直于入射面振动的分量为光矢量的s分量,记为ES。
1.菲涅耳公式:
表示反射波、折射波与入射波的振幅和位相关系。
(1)S波(垂直于入射面分量)的菲涅耳公式
——S波的振幅反射系数
——S波的振幅透射系数
(2)P波(平行于入射面分量)的菲涅耳公式
——P波的振幅反射系数
——P波的振幅透射系数
2.光从光疏介质入射到光密介质(如空气射向玻璃)
当时,即垂直入射时,都不为零,表示存在反射波和折射波。
当时,即掠入射时,即没有折射光波。
、
随
的增大而减小;
随
的增大而增大,直到等于1;
值在
时,有
,即反射光波中没有p波,只有s波,产生全偏振现象。
3.光从光密介质入射到光疏介质(n2/n1<1时)
当时,即垂直入射时,都不为零,表示存在反射波和折射波。
当(θc为θ2=90度时对应的θ1)时,
表示发生全反射现象。
、
都大于1,且随
的增大而增大。
4.当平面波在接近正入射或掠入射下从光疏介质与光密介质的分界面反射时,反射光的电矢量相对于入射光的电矢量产生了的相位突变(半波损失:
反射时损失了半个波长)。
实验内容
根据菲涅尔公式,分别画出从光疏介质入射到光密介质和从光密介质入射到光疏介质时,反射波和折射波与入射波振幅的相对变化关系以及相位变化关系,并进行讨论。
所用的程序:
n1=input('n1=');%输入折射率n1
n2=input('n2=');%输入折射率n2
theta_i=0:
0.001:
pi/2;%定义入射的变换范围
z=(n1/n2)*sin(theta_i);%计算折射角的正弦值和范围
E=z<=1;%取x中小于1的元素,E与x同尺寸,x中小于1的元素对应下标为1,否则为0
Y=z(E);%Y只能取E为1的下标对应的x的元素,Y的尺寸只是x中小于的元素的个数
F=z>1;%取x中大于1的元素,为了得到z大于1的元素的尺寸
x=[Y,ones(size(z(F)))];
theta_t=-j*log(j*z+(ones(size(z))-z.^2).^0.5);%计算折射角,并用log的方式表示
r_pa=(n2*cos(theta_i)-n1*cos(theta_t))./(n2*cos(theta_i)+n1*cos(theta_t));%p波的振幅反射系数
t_pa=(2*n1*cos(theta_i))./(n2*cos(theta_i)+n1*cos(theta_t));%p波振幅透射系数
r_pe=(n1*cos(theta_i)-n2*cos(theta_t))./(n1*cos(theta_i)+n2*cos(theta_t));%s波的振幅反射系数
t_pe=(2*n1*cos(theta_i))./(n1*cos(theta_i)+n2*cos(theta_t));%s波的振幅透射系数
[NM]=min(abs(r_pa));%取n1的p波振幅反射系数的模的最小值
theta_i=theta_i/pi*180;%入射角弧度值
theta_cri=asin(n2/n1)*180/pi;%入射角反正弦值
figure
(1)
plot(theta_i,real(t_pa),'-',theta_i,real(t_pe),':
',theta_i,real(r_pa),'-',theta_i,real(r_pe),'--',theta_i(M),0,'o');%绘制入射角和反射角的p波、s波的振幅反射系数和振幅反射系数的模随入射角变化而变化的图
m1=min([real(t_pa)real(t_pe)real(r_pa)real(r_pe)]);%取入射角和反射角的p波、s波的振幅反射系数和振幅反射系数模的最小值
M1=max([real(t_pa)real(t_pe)real(r_pa)real(r_pe)]);%取入射角和反射角的p波、s波的振幅反射系数和振幅反射系数模的最大值
legend('t_p_a','t_p_e','r_p_a','r_p_e')
text(theta_i(M),0.01*(M1-m1),['\leftarrow']);
text(theta_i(M),0.018*(M1-m1),['Brewsterangle=',num2str(theta_i(M))]);%加坐标轴的名称
ifn1>=n2
text(theta_cri,-0.1*(M1-m1),['Critialangle=',num2str(theta_cri)]);
end
xlabel('incidentangle');axis([0120m1*1.1M1*1.1]);title('Realpartofcoefficents')
gridon
figure
(2)
plot(theta_i,imag(t_pa),'-',theta_i,imag(t_pe),':
',theta_i,imag(r_pa),'-',theta_i,imag(r_pe),'--',theta_i(M),0,'o');%绘制入射角和反射角的p波、s波的振幅反射系数和振幅反射系数的虚部随入射角变化而变化的图
m1=min([imag(t_pa)imag(t_pe)imag(r_pa)imag(r_pe)]);%取入射角和反射角的p波、s波的振幅反射系数和振幅反射系数虚部的最小值
M1=max([imag(t_pa)imag(t_pe)imag(r_pa)imag(r_pe)]);%取入射角和反射角的p波、s波的振幅反射系数和振幅反射系数虚部的最大值
m11=min([abs(imag(t_pa))abs(imag(t_pe))abs(imag(r_pa))abs(imag(r_pe))]);%取入射角和反射角的p波、s波的振幅反射系数和振幅反射系数的绝对值的最小值
M11=max([abs(imag(t_pa))abs(imag(t_pe))abs(imag(r_pa))abs(imag(r_pe))]);%取入射角和反射角的p波、s波的振幅反射系数和振幅反射系数的绝对值的最大值
legend('t_p_a','t_p_e','r_p_a','r_p_e')
text(theta_i(M),0.01*(M11-m11+0.2),['\leftarrow']);
text(theta_i(M),0.03*(M11-m11+0.2),['Brewsterangle=',num2str(theta_i(M))]);
ifn1>=n2
text(theta_cri,-0.1*(M1-m1),['Critialangle=',num2str(theta_cri)]);
end
xlabel('incidentangle');axis([0120m1*1.1-0.1M1*1.1+0.1]);title('Realpartofcoefficents')
gridon
figure(3)
plot(theta_i,abs(t_pa),'-',theta_i,abs(t_pe),':
',theta_i,abs(r_pa),'-',theta_i,abs(r_pe),'--',theta_i(M),0,'o');%绘制入射角和反射角的p波、s波的振幅反射系数和振幅反射系数的绝对值随入射角变化而变化的图
m1=min([abs(t_pa)abs(t_pe)abs(r_pa)abs(r_pe)]);%取入射角和反射角的p波、s波振幅反射系数和振幅发射系数的绝对值的最小值
M1=max([abs(t_pa)abs(t_pe)abs(r_pa)abs(r_pe)]);%取入射角和反射角的p波、s波振幅反射系数和振幅发射系数的绝对值的最大值
legend('t_p_a','t_p_e','r_p_a','r_p_e')
text(theta_i(M),0.01*(M1-m1),['\leftarrow']);
text(theta_i(M),0.03*(M1-m1),['Brewsterangle=',num2str(theta_i(M))]);
ifn1>=n2
text(theta_cri,0.1*(M1-m1),['Critialangle=',num2str(theta_cri)]);
end
xlabel('incidentangle');axis([0120m1*1.1M1*1.1]);title('Magnitudeofcoefficents')
gridon
figure(4)
plot(theta_i,real(j*180/pi*log(t_pa)),'-',theta_i,real(j*180/pi*log(t_pe)),':
',theta_i,real(j*180/pi*log(r_pa)),'-',theta_i,real(j*180/pi*log(r_pe)),'--',theta_i(M),0,'o');%绘制入射角和反射角的p波、s波的振幅反射系数和振幅反射系数的模随入射角变化而变化的图,用log表示
m1=real(min([-j*180/pi*log(r_pa)-j*180/pi*log(r_pe)]));%取入射角和反射角的p波、s波的振幅反射系数和振幅反射系数模的最小值,并用log方式表示
M1=max([-j*180/pi*log(r_pa)-j*180/pi*log(r_pe)]);%取入射角和反射角的p波、s波的振幅反射系数和振幅反射系数模的最大值,并用log表示
legend('t_p_a','t_p_e','r_p_a','r_p_e')
text(theta_i(M),0.015*(M1-m1),['\leftarrow']);
text(theta_i(M),0.03*(M1-m1),['Brewsterangle=',num2str(theta_i(M))]);
ifn1>=n2
text(theta_cri,-0.1*(M1-m1),['Critialangle=',num2str(theta_cri)]);
end
xlabel('incidentangle');axis([0120-185185]);title('Phaseangleofcoefficents')
gridon
当:
n1=1.5
n2=1
表示光密介质入射到光疏介质的情况,当
=0时,
、
都不为0,表示存在反射波和折射波,,当
时,
=
=1,表示发生全反射现象。
绘制入射角和反射角的p波、s波的振幅反射系数虚部随入射角变化而变化的图,随角度的变化,反射系数的虚部由0变为非0再变为0,即是有一开始没有发生反射,后来发生反射和折射,最后没有反射。
绘制入射角和反射角的p波、s波的振幅反射系数和振幅反射系数的绝对值随入射角变化而变化的图,一开始p光折射大于反射,达到全反射后折射系数变为1
绘制入射角和反射角的p波、s波的振幅反射系数和振幅反射系数的模随入射角变化而变化的图,用log表示
当