CN105929531B - 一种激光角反射器远场衍射强度的数值仿真方法 - Google Patents
一种激光角反射器远场衍射强度的数值仿真方法 Download PDFInfo
- Publication number
- CN105929531B CN105929531B CN201610545255.5A CN201610545255A CN105929531B CN 105929531 B CN105929531 B CN 105929531B CN 201610545255 A CN201610545255 A CN 201610545255A CN 105929531 B CN105929531 B CN 105929531B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- laser
- mtr
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Optical Elements Other Than Lenses (AREA)
Abstract
本发明涉及一种激光角反射器远场衍射强度数值仿真方法,属于光学设计与仿真领域。本发明以激光角反射器实际加工过程中的二面直角误差、面形误差和倒棱等误差参数、实际使用过程中的光束倾斜效应和激光偏振态效应等因素为输入条件,以激光角反射器的衍射区域、出射激光束的复振幅分布和简化的远场衍射强度模型为基础,通过入射激光束和接收屏的离散划分以及光线的追迹方法,仿真得到综合影响因素条件下激光角反射器远场衍射强度的数值分布。所采用的仿真方法充分考虑了影响激光角反射器远场强度分布的所有因素,扩充和完善了激光角反射器远场衍射强度的仿真范畴,特别适用于远距离条件下激光角反射器的使用性能指标的综合评估。
Description
技术领域
本发明属于光学设计与仿真领域,具体涉及一种激光角反射器远场衍射强度数值仿真方法,适用于激光测距合作目标的性能指标综合评价。
背景技术
激光角反射器是由三个相互垂直的直角反射面和一个底面组成的实心四面体,其几何结构见图1(a)所示。平面OBC、OAC和OAB为激光角反射器的直角反射面,分别以数字1、2和3来表示,平面ABC为激光角反射器的底面。通常情况下,激光角反射器底面采用圆形切割以便于实际装配,见图1(b)所示。当入射激光束从激光角反射器底面入射时,将先后经过底面的第一次折射、反射面的三次反射和底面的再次折射,以六种不同的反射次序:1→2→3、1→3→2、2→1→3、2→3→1、3→1→2、3→2→1,从底面出射。对于理想的激光角反射器而言,考虑到其特殊的几何结构,出射激光束将与入射激光束的方向平行,即理想激光角反射器具有定向反射特性,因此,通常其作为激光测距中的合作目标,用于增加测距仪接收机位置处的激光强度。
激光角反射器在使用过程中,其与激光测距机之间距离远大于激光角反射器的有效口径,因此,出射激光束必然受到衍射效应的制约,即在激光测距仪接收机位置处接收到的激光强度是激光角反射器远场衍射强度。由于激光角反射器在加工制造过程会出现二面直角误差和面形误差,若偏振入射激光束倾斜入射至激光角反射器底面时,则激光角反射器远场衍射强度的分布将受到激光角反射器几何参数、光束入射条件和接收区域有关。
建立如图2所示的坐标框架,包括激光角反射器棱坐标系和底面坐标系、接收屏坐标系,其中,激光角反射器棱坐标系的原点位于顶点处,其xl轴、yl轴和zl轴分别与激光角反射器直角边重合;底面坐标系的原点位于底面中心,zf轴与底面的法线方向平行,yf轴是yl轴在底面的投影,xf轴与yf轴和zf轴构成右手螺旋法则;接收屏坐标系的原点是入射激光束中心与接收屏的交点,zr轴与入射激光束方向平行,xr轴和yr轴位于接收屏范围内。因此,在距离激光角反射器底面d0位置处所接收到的远场衍射强度可以表示为:
其中,为入射激光束在激光角反射器底面上的入射角,λ为入射激光束的波长,k=2π/λ,Σ为激光角反射器的有效衍射区域,Ep(xf,yf)和Es(xf,yf)分别为出射激光束在激光角反射器底面上的p波和s波光场分量,r是激光角反射器底面上的点(xf,yf)与接收屏上的点(xr,yr)之间的距离。考虑到激光角反射器远场衍射强度的计算模型不存在解析解,并且其受到的影响因素较多,因此,构建完善的激光角反射器远场衍射强度的综合仿真方法是十分复杂和困难的。
近些年来,已经有一些关于激光角反射器远场衍射强度仿真的报道,如在垂直入射并忽略入射激光束偏振态条件下的远场衍射强度分析(聂辉,等.角锥棱镜的远场衍射特性.光学学报,23(12):1470-1474,2003)、在忽略面形误差和入射激光束偏振效应的远场衍射强度数值模拟(叶一东,等.角锥后向反射器的数值模拟研究.光学学报,23(4):485-490,2003)、在垂直入射并忽略激光角反射器加工误差的远场衍射强度仿真(Murphy T W,etal.Polarization and far-field diffraction patterns of total internalreflection corner cubes.Applied optics,52(2):117-126,2013.),这些仿真方法是在仅考虑影响激光角反射器远场衍射强度分布的部分因素前提下,利用相应的假设近似来获得激光角反射器远场衍射强度的分布规律。采用这些方法仅能仿真得到满足某些特定激光角反射器的几何参数或测量条件的远场衍射强度的仿真,无法实现激光角反射器实际使用过程中的远场衍射强度的真实仿真。
发明内容
本发明主要是解决现有仿真方法所存在的问题:提出一种综合考虑激光角反射器实际加工过程中的二面直角误差、面形误差和倒棱等误差参数,实际使用过程中的光束倾斜效应和激光偏振态效应等因素,通过入射激光束和接收屏的离散划分以及光线的追迹方法,仿真得到非理想条件下激光角反射器远场衍射强度的数值,以实现对激光角反射器实际使用性能的精确评估。
本发明的上述技术问题主要是通过下述技术方案得以解决的:
一种激光角反射器远场衍射强度的仿真方法,包括以下过程:首先根据激光角反射器远场衍射强度数值仿真的初始参数计算激光束在激光角反射器内部的光线方向;然后利用光线方向计算激光束与激光角反射器表面的交点坐标;其次根据各交点坐标的限制条件计算激光角反射器的衍射区域范围;再次利用激光角反射器加工误差和入射激光偏振态计算出射激光束的附加复振幅分布;最后根据衍射区域范围和出射激光束的附加复振幅分布的结果,计算接收屏位置处的激光角反射器远场衍射强度分布。
在上述的一种激光角反射器远场衍射强度的仿真方法,所述的激光角反射器远场衍射强度数值仿真的初始参数包括:
激光角反射器的几何参数,包括口径和倒棱宽度;
激光角反射器的加工误差,包括二面直角误差、反射面的面形误差和底面的面形误差;
激光角反射器材料的折射率和反射面膜层的折射率,若未反射面镀膜,则膜层折射率置为1;
入射激光束的入射条件,包括入射激光束的入射角和方位角,入射光场p波与s波分量的复振幅分布;
入射激光束的波长和用于激光角反射器面形检测的激光波长;
接收屏的大小及其与激光角反射器之间的距离;
激光角反射器底面和接收屏细分的点数。
在上述的一种激光角反射器远场衍射强度的仿真方法,计算激光束在激光角反射器内部的光线方向,包括以下子步骤:
步骤3.1,计算入射激光束的入射矢量R和折射光矢量A:
A=[A1,A2,A3]T=refr(-MR,N0);
式中,和θ分别为入射激光束在激光角反射器底面坐标系的入射角和方位角,N0为激光角反射器底面的法线方向,M=[M1;M2;M3], refr函数为光线折射定律的矢量形式,R1、R2和R3分别为入射光矢量在激光角反射器底面坐标系下的分量,A1、A2和A3分别为折射光矢量在激光角反射器反射面坐标系下的分量。
步骤3.2,计算激光束在激光角反射器内部经过6种不同反射次序后所对应的光线方向矢量Bm、Cm、Dm和Em,m=1、2、3…、6:
Bm=γ1A,m=1、2;Bm=γ2A,m=3、4;Bm=γ3A,m=5、6;
Cm=γ2Bm,m=1、6;Cm=γ3Bm,m=2、4;Cm=γ1Bm,m=3、5;
Dm=-A;Em=R;
其中,
在上述的一种激光角反射器远场衍射强度的仿真方法,计算激光束与激光角反射器表面的交点坐标,包括以下子步骤:
步骤4.1,将入射激光束细分为nx×ny根激光线,每根激光线在激光角反射器底面上的横坐标xfab和纵坐标yfab分别满足:
xfab=-r+2(a-1)r/nx,a=1、2、3…、nx;
yfab=-r+2(b-1)r/ny,b=1、2、3…、ny;
其中,r为激光角反射器底面的半径。
步骤4.2,利用激光角反射器底面坐标系和反射面坐标系之间的转换关系,获得每根激光线在反射面坐标系下的坐标(xlab,ylab,zlab):
[xlab,ylab,zlab]T=M[xfab,yfab,0]T+Mr,
式中,L为激光角反射器的棱长,通常情况下,
步骤4.3,根据入射激光束细化光线的起始坐标以及激光束在激光角反射器内部的矢量方向,计算激光束与激光角反射器底面的交点坐标:
[xomab,yomab,zomab]T=g(βk,Cm)·[xlab,ylab,zlab]T;
[xtmab,ytmab,ztmab]T=g(βj,Bm)·[xomab,yomab,zomab]T;
[xhmab,yhmab,zhmab]T=g(βi,A)·[xtmab,ytmab,ztmab]T;
式中,(xomab,yomab,zomab)、(xtmab,ytmab,ztmab)和(xhmab,yhmab,zhmab)分别表示激光束与激光角反射器反射面的交点坐标,(xcmab,ycmab,zcmab)表示激光束与激光角反射器底面的交点坐标。βi、βj、βk和β4表示系数矩阵,下标i、j和k的数值与光线反射次序有关,它们取值为:
当m=1、2时,i=1;当m=3、4时,i=2;当m=5、6时,i=3;
当m=3、5时,j=1;当m=1、6时,j=2;当m=2、4时,j=3;
当m=4、6时,k=1;当m=2、5时,k=2;当m=1、3时,k=3;
系数矩阵β1=[1,0,0],β2=[0,1,0],β3=[0,0,1],β4=[1,1,1],函数g(β,J)满足关系式:
其中,β和J表示函数的输入量,β=[βx,βy,βz],βx、βy和βz分别为β的分量;J=[Jx,Jy,Jz],Jx、Jy和Jz分别为J的分量。
在上述的一种激光角反射器远场衍射强度的仿真方法,计算激光角反射器的衍射区域范围,包括以下子步骤:
步骤5.1,激光束与反射面的交点必须限制在激光角反射器反射面范围内,则符合该条件的激光线序列号集合(Ra,Rb)必须满足以下形式:
(Ra,Rb)={(a,b)|T1∩T2∩T3}
其中,符号“∩”表示交集运算,集合T1、T2和T3表达式为:
式中,l表示激光角反射器的倒棱宽度。
步骤5.2,入射激光线和激光束与底面交点必须限制在激光角反射器底面范围内,则符合该条件的激光线序列号集合(Fa,Fb)必须满足:
步骤5.3,符合反射区域和底面区域限制要求的激光线所占的底面范围,即为激光角反射器的衍射区域(Ha,Hb),其表达形式为:
(Ha,Hb)=(Ra,Rb)∩(Fa,Fb)。
在上述的一种激光角反射器远场衍射强度的仿真方法,计算出射激光束的附加复振幅分布,包括以下子步骤:
步骤6.1,计算由激光角反射器反射面面形误差所导致的出射激光束附加光程ωrm(m=1、2、3...、6),它满足以下关系:
式中,ξi(i=1、2、3)表示第i个反射面的面形光圈值,nr为表示激光角反射器材料的折射率,λc为面形检测的激光波长,Pi(i=1、2、3)可以表示为:
其中,
步骤6.2,计算由激光角反射器底面面形误差所导致的出射激光束附加光程ωfm(m=1、2、3...、6),满足以下关系:
式中,ξf为激光角反射器底面的面形光圈值,α为入射激光束在底面的折射角,
步骤6.3,计算由二面直角误差所导致的出射激光束附加光程,包括以下子过程:
(6.3.1)计算存在二面直角误差时,经过激光角反射器反射面三次反射后的反射光束的附加矢量dFm,m=1、2、3…、6:
dF1=[2(A2δ1+A3δ2),2(A3δ3-A1δ1),-2(A2δ3+A1δ2)]T;
dF2=[2(A3δ2+A2δ1),-2(A3δ3+A1δ1),2(A2δ3-A1δ2)]T;
dF3=[2(A3δ2-A2δ1),2(A1δ1+A3δ3),-2(A1δ2+A2δ3)]T;
dF4=-dF2;dF5=-dF3;dF6=-dF1,
式中,δ1、δ2和δ3为激光角反射器三个反射面两两之间的二面直角误差。
(6.3.2)计算反射光束的附加矢量dFm经过激光角反射器底面折射后,出射激光束的附加矢量dRm:
dRm=[dRmx,dRmy,dRmz]T=refr(dFm,-N0),m=1、2、3…、6。
(6.3.3)计算由二面直角误差导致的出射激光束附加光程ωam(m=1、2、3...、6),满足以下关系:
步骤6.4,计算由入射激光偏振态所导致的出射激光束的附加复振幅分布,包括以下子过程:
(6.4.1)计算入射激光束经过一次折射和三次反射后,在激光角反射器底面和反射面上的p波和s波的矢量方向,满足以下关系式:
式中,(pr,sr)和(pfm,sfm)分别为折射光矢量A在底面和第一次反射面上的p波和s波的矢量方向,(pem,sem)和(ptm,stm)分别为反射光矢量Bm在第一次反射面和第二次反射面上的p波和s波的矢量方向,(pum,sum)和(pvm,svm)分别为反射光矢量Cm在第二次反射面和第三次反射面上的p波和s波的矢量方向,(pgm,sgm)和(phm,shm)分别为反射光矢量Dm在第三次反射面和底面上的p波和s波的矢量方向,下标m为正整数,m=1、2、3…、6,符号表示向量的叉乘运算,N4为激光出射时底面的法线方向,N4=-N0,N1、N2和N3为反射面的法线方向,它们满足以下形式:
N1=[-1,0,0]T,m=1、2;N1=[0,-1,0]T,m=3、4;N1=[0,0,-1]T,m=5、6;
N2=[-1,0,0]T,m=3、5;N2=[0,-1,0]T,m=1、6;N2=[0,0,-1]T,m=2、4;
N3=[-1,0,0]T,m=4、6;N3=[0,-1,0]T,m=2、5;N3=[0,0,-1]T,m=1、3。
(6.4.2)计算折射光和反射光在激光角反射器表面的转换矩阵Γt(t=1、2、3、4),满足以下关系式:
式中,符号“⊙”表示矢量的点乘运算。
(6.4.3)计算折射光和反射光的p波和s波分量的反射系数矩阵和折射系数矩阵,满足以下关系式:
式中,Tf和Tl分别为第一次折射和最后一次折射的折射系数矩阵,Rf、Re和Rt分别为第一次反射、第二次反射和第三次反射的反射系数矩阵。折射系数矩阵和反射系数矩阵的非零元素满足以下关系式:
rfp=crp(if,zf);rfs=crs(if,zf);rep=crp(ie,ze);res=crs(ie,ze);
rtp=crp(it,zt);rts=crs(it,zt);
其中,(if,zf)、(ie,ze)和(it,zt)分别表示激光束在第一次反射、第二次反射和第三次反射时对应的入射角和折射角,if=arccos(A⊙N1),zf=arcsin(υsin if),ie=arccos(Bm⊙N2),ze=arcsin(υsin ie),it=arccos(Cm⊙N3),zt=arcsin(υsin it),式中,υ=nr/nf,nf为激光角反射器膜层的折射率。函数ctp和cts分别表示p波和s波的透射系数的计算公式,函数crp和crs分别表示p波和s波的反射系数的计算公式,它们满足以下关系式:
其中,表示函数的输入量。
(6.4.4)计算由入射激光偏振态所导致的出射激光束p波分量和s波分量的附加复振幅APout和ASout的表达形式:
式中,EPin和ESin分别为入射光场p波分量和s波分量的复振幅。
步骤6.5,计算出射激光束p波分量和s波分量的附加复振幅分布PEm和SEm,m=1、2、3…、6,满足以下关系式:
PEm=|PEm|exp(iφpm)=|APout|exp[iτ(ωrm+ωfm+ωam)+iωpm];
SEm=|SEm|exp(iφsm)=|ASout|exp[iτ(ωrm+ωfm+ωam)+iωsm];
式中,|PEm|和Фpm分别表示PEm的模和相位,|SEm|和Фsm分别表示SEm的模和相位,i表示虚数单位。τ=2π/λr,λr为入射激光束的波长,ωpm和ωsm分别为APout和ASout的幅角。
在上述的一种激光角反射器远场衍射强度的仿真方法,计算接收屏位置处的激光角反射器远场衍射强度分布,包括以下子步骤:
步骤7.1,将接收屏细分为μx×μy个离散接收点,每个接收点在接收屏上的横坐标xrcd和纵坐标yrcd分别满足:
xrcd=-Rw+2(c-1)Rw/μx,c=1、2、3…、μx;
yrcd=-Rw+2(d-1)Rw/μy,d=1、2、3…、μy;
其中,Rw为接收屏的半径大小。
步骤7.2,计算出射激光束的p波分量和s波分量的复振幅分布Pmab和Smab,满足以下关系式:
Pmab=exp{iτ[(dRmx+A1)xfab+(dRmy+A2)yfab]}·PEm;
Smab=exp{iτ[(dRmx+A1)xfab+(dRmy+A2)yfab]}·SEm;
步骤7.3,计算出射激光束p波分量和s波分量的远场衍射强度值IPcd和IScd,满足以下形式:
式中,d0为激光角反射器与接收屏之间的距离,i表示虚数单位,xgcd和ygcd可以表示为:
步骤7.4,计算接收屏位置处的激光角反射器远场衍射强度分布Icd,满足关系式:Icd=IPcd+IScd。
本发明具有如下优点:以激光角反射器的实际几何参数、二面直角误差、反射面与底面的面形误差等为仿真输入,在考虑倾光束斜入射和激光束偏振态的情况下,利用光线的追迹方法计算得到激光角反射器实际的衍射区域范围和出射激光束的复振幅分布,并通过离散的数值模拟出激光角反射器的远场衍射强度分布。所采用的仿真方法不仅仅以光束入射条件、二面直角误差、面形误差或激光偏振态等部分因素作为分析对象,扩充和完善了激光角反射器远场衍射强度的仿真范畴,特别适用于远距离条件下激光角反射器的使用性能指标的综合评估。
附图说明
图1(a)是激光角反射器的几何结构示意图。
图1(b)是圆形切割方式的激光角反射器结构。
图2是激光角反射器远场衍射仿真计算过程中所涉及的几种坐标框架。
图3是激光角反射器远场衍射强度数值仿真计算的流程。
图4是实施例激光角反射器的衍射区域分布。
图5是实施例出射激光束p波分量和s波分量的附加复振幅的模分布。
图6(a)是实施例出射激光束p波分量附加复振幅的相位分布。
图6(b)是实施例出射激光束s波分量附加复振幅的相位分布。
图7是实施例激光角反射器远场衍射强度的等高线分布。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例:
一、首先介绍一下本发明的具体方法流程,主要包括以下步骤:
1.输入激光角反射器远场衍射强度数值仿真的初始参数,包括下述过程:
(1.1)输入激光角反射器的几何参数,包括口径和倒棱宽度;
(1.2)输入激光角反射器的加工误差,包括二面直角误差、反射面的面形误差和底面的面形误差;
(1.3)输入激光角反射器材料的折射率和反射面膜层的折射率,若未反射面镀膜,则膜层折射率置为1;
(1.4)输入入射激光束的入射条件,包括入射激光束的入射角和方位角,入射光场p波与s波分量的复振幅分布;
(1.5)输入入射激光束的波长和用于激光角反射器面形检测的激光波长;
(1.6)输入接收屏的大小及其与激光角反射器之间的距离;
(1.7)输入激光角反射器底面和接收屏细分的点数。
2.计算激光束在激光角反射器内部的光线方向,包括下述过程:
(2.1)计算入射激光束的入射矢量R和折射光矢量A:
A=[A1,A2,A3]T=refr(-MR,N0);
式中,和θ分别为入射激光束在激光角反射器底面坐标系的入射角和方位角,N0为激光角反射器底面的法线方向,M=[M1;M2;M3], refr函数为光线折射定律的矢量形式,R1、R2和R3分别为入射光矢量在激光角反射器底面坐标系下的分量,A1、A2和A3分别为折射光矢量在激光角反射器反射面坐标系下的分量。
(2.2)计算激光束在激光角反射器内部经过6种不同反射次序后所对应的光线方向矢量Bm、Cm、Dm和Em,m=1、2、3…、6:
Bm=γ1A,m=1、2;Bm=γ2A,m=3、4;Bm=γ3A,m=5、6;
Cm=γ2Bm,m=1、6;Cm=γ3Bm,m=2、4;Cm=γ1Bm,m=3、5;
Dm=-A;Em=R;
其中,
3.计算激光束与激光角反射器表面的交点坐标,包括下述过程:
(3.1)将入射激光束细分为91×91根激光线,每根激光线在激光角反射器底面上的横坐标xfab和纵坐标yfab分别满足:
xfab=-r+2(a-1)r/91,a=1、2、3…、91;
yfab=-r+2(b-1)r/91,b=1、2、3…、91;
其中,r为激光角反射器底面的半径。
(3.2)利用激光角反射器底面坐标系和反射面坐标系之间的转换关系,获得每根激光线在反射面坐标系下的坐标(xlab,ylab,zlab):
[xlab,ylab,zlab]T=M[xfab,yfab,0]T+Mr,
式中,L为激光角反射器的棱长,通常情况下,
(3.3)根据入射激光束细化光线的起始坐标以及激光束在激光角反射器内部的矢量方向,计算激光束与激光角反射器底面的交点坐标:
[xomab,yomab,zomab]T=g(βk,Cm)·[xlab,ylab,zlab]T;
[xtmab,ytmab,ztmab]T=g(βj,Bm)·[xomab,yomab,zomab]T;
[xhmab,yhmab,zhmab]T=g(βi,A)·[xtmab,ytmab,ztmab]T;
式中,(xomab,yomab,zomab)、(xtmab,ytmab,ztmab)和(xhmab,yhmab,zhmab)分别表示激光束与激光角反射器反射面的交点坐标,(xcmab,ycmab,zcmab)表示激光束与激光角反射器底面的交点坐标。βi、βj、βk和β4表示系数矩阵,下标i、j和k的数值与光线反射次序有关,它们取值为:
当m=1、2时,i=1;当m=3、4时,i=2;当m=5、6时,i=3;
当m=3、5时,j=1;当m=1、6时,j=2;当m=2、4时,j=3;
当m=4、6时,k=1;当m=2、5时,k=2;当m=1、3时,k=3;
系数矩阵β1=[1,0,0],β2=[0,1,0],β3=[0,0,1],β4=[1,1,1],函数g(β,J)满足关系式:
其中,β和J表示函数的输入量,β=[βx,βy,βz],βx、βy和βz分别为β的分量;J=[Jx,Jy,Jz],Jx、Jy和Jz分别为J的分量。
4.计算激光角反射器的衍射区域范围,包括下述过程:
(4.1)激光束与反射面的交点必须限制在激光角反射器反射面范围内,则符合该条件的激光线序列号集合(Ra,Rb)必须满足以下形式:
(Ra,Rb)={(a,b)|T1∩T2∩T3}
其中,符号“∩”表示交集运算,集合T1、T2和T3表达式为:
式中,l表示激光角反射器的倒棱宽度。
(4.2)入射激光线和激光束与底面交点必须限制在激光角反射器底面范围内,则符合该条件的激光线序列号集合(Fa,Fb)必须满足:
(4.3)符合反射区域和底面区域限制要求的激光线所占的底面范围,即为激光角反射器的衍射区域(Ha,Hb),其表达形式为:
(Ha,Hb)=(Ra,Rb)∩(Fa,Fb)。
在本步骤中,激光角反射器的衍射区域范围综合考虑了激光角反射器圆形切割和倒棱等因素的影响。
5.计算出射激光束的附加复振幅分布,包括下述过程:
(5.1)计算由激光角反射器反射面面形误差所导致的出射激光束附加光程ωrm(m=1、2、3...、6),它满足以下关系:
式中,ξi(i=1、2、3)表示第i个反射面的面形光圈值,nr为表示激光角反射器材料的折射率,λc为面形检测的激光波长,Pi(i=1、2、3)可以表示为:
其中,
(5.2)计算由激光角反射器底面面形误差所导致的出射激光束附加光程ωfm(m=1、2、3...、6),满足以下关系:
式中,ξf为激光角反射器底面的面形光圈值,α为入射激光束在底面的折射角,
(5.3)计算由二面直角误差所导致的出射激光束附加光程,包括以下子过程:
(5.3.1)计算存在二面直角误差时,经过激光角反射器反射面三次反射后的反射光束的附加矢量dFm,m=1、2、3…、6:
dF1=[2(A2δ1+A3δ2),2(A3δ3-A1δ1),-2(A2δ3+A1δ2)]T;
dF2=[2(A3δ2+A2δ1),-2)A3δ3+A1δ1),2(A2δ3-A1δ2)]T;
dF3=[2(A3δ2-A2δ1),2(A1δ1+A3δ3),-2(A1δ2+A2δ3)]T;
dF4=-dF2;dF5=-dF3;dF6=-dF1,
式中,δ1、δ2和δ3为激光角反射器三个反射面两两之间的二面直角误差。
(5.3.2)计算反射光束的附加矢量dFm经过激光角反射器底面折射后,出射激光束的附加矢量dRm:
dRm=[dRmx,dRmy,dRmz]T=refr(dFm,-N0),m=1、2、3…、6。
(5.3.3)计算由二面直角误差导致的出射激光束附加光程ωam(m=1、2、3...、6),满足以下关系:
(5.4)计算由入射激光偏振态所导致的出射激光束的附加复振幅分布,包括以下子过程:
(5.4.1)计算入射激光束经过一次折射和三次反射后,在激光角反射器底面和反射面上的p波和s波的矢量方向,满足以下关系式:
式中,(pr,sr)和(pfm,sfm)分别为折射光矢量A在底面和第一次反射面上的p波和s波的矢量方向,(pem,sem)和(ptm,stm)分别为反射光矢量Bm在第一次反射面和第二次反射面上的p波和s波的矢量方向,(pum,sum)和(pvm,svm)分别为反射光矢量Cm在第二次反射面和第三次反射面上的p波和s波的矢量方向,(pgm,sgm)和(phm,shm)分别为反射光矢量Dm在第三次反射面和底面上的p波和s波的矢量方向,下标m为正整数,m=1、2、3…、6,符号表示向量的叉乘运算,N4为激光出射时底面的法线方向,N4=-N0,N1、N2和N3为反射面的法线方向,它们满足以下形式:
N1=[-1,0,0]T,m=1、2;N1=[0,-1,0]T,m=3、4;N1=[0,0,-1]T,m=5、6;
N2=[-1,0,0]T,m=3、5;N2=[0,-1,0]T,m=1、6;N2=[0,0,-1]T,m=2、4;
N3=[-1,0,0]T,m=4、6;N3=[0,-1,0]T,m=2、5;N3=[0,0,-1]T,m=1、3。
(5.4.2)计算折射光和反射光在激光角反射器表面的转换矩阵Γt(t=1、2、3、4),满足以下关系式:
式中,符号“⊙”表示矢量的点乘运算。
(5.4.3)计算折射光和反射光的p波和s波分量的反射系数矩阵和折射系数矩阵,满足以下关系式:
式中,Tf和Tl分别为第一次折射和最后一次折射的折射系数矩阵,Rf、Re和Rt分别为第一次反射、第二次反射和第三次反射的反射系数矩阵。折射系数矩阵和反射系数矩阵的非零元素满足以下关系式:
rfp=crp(if,zf);rfs=crs(if,zf);rep=crp(ie,ze);res=crs(ie,ze);
rtp=crp(it,zt);rts=crs(it,zt);
其中,(if,zf)、(ie,ze)和(it,zt)分别表示激光束在第一次反射、第二次反射和第三次反射时对应的入射角和折射角,if=arccos(A⊙N1),zf=arcsin(υsin if),ie=arccos(Bm⊙N2),ze=arcsin(υsin ie),it=arccos(Cm⊙N3),zt=arcsin(υsin it),式中,υ=nr/nf,nf为激光角反射器膜层的折射率。函数ctp和cts分别表示p波和s波的透射系数的计算公式,函数crp和crs分别表示p波和s波的反射系数的计算公式,它们满足以下关系式:
其中,表示函数的输入量。
(5.4.4)计算由入射激光偏振态所导致的出射激光束p波分量和s波分量的附加复振幅APout和ASout的表达形式:
式中,EPin和ESin分别为入射光场p波分量和s波分量的复振幅。
(5.5)计算出射激光束p波分量和s波分量的附加复振幅分布PEm和SEm,m=1、2、3…、6,满足以下关系式:
PEm=|PEm|exp(iφpm)=|APout|exp[iτ(ωrm+ωfm+ωam)+iωpm];
SEm=|SEm|exp(iφsm)=|ASout|exp[iτ(ωrm+ωfm+ωam)+iωsm];
式中,|PEm|和Фpm分别表示PEm的模和相位,|SEm|和Фsm分别表示SEm的模和相位,i表示虚数单位。τ=2π/λr,λr为入射激光束的波长,ωpm和ωsm分别为APout和ASout的幅角。
在本步骤中,出射激光束的附加复振幅分布综合考虑了二面直角误差、面形误差和入射激光偏振态等多重因素的影响。
6.计算接收屏位置处的激光角反射器远场衍射强度分布,包括下述过程:
(6.1)将接收屏细分为101×101个离散接收点,每个接收点在接收屏上的横坐标xrcd和纵坐标yrcd分别满足:
xrcd=-Rw+2(c-1)Rw/101,c=1、2、3…、101;
yrcd=-Rw+2(d-1)Rw/μy,d=1、2、3…、101;
其中,Rw为接收屏的半径大小。
(6.2)计算出射激光束的p波分量和s波分量的复振幅分布Pmab和Smab,满足以下关系式:
Pmab=exp{iτ[(dRmx+A1)xfab+(dRmy+A2)yfab]}·PEm;
Smab=exp{iτ[(dRmx+A1)xfab+(dRmy+A2)yfab]}·SEm;
(6.3)计算出射激光束p波分量和s波分量的远场衍射强度值IPcd和IScd,满足以下形式:
式中,d0为激光角反射器与接收屏之间的距离,i表示虚数单位,xgcd和ygcd可以表示为:
(6.4)计算接收屏位置处的激光角反射器远场衍射强度分布Icd,满足关系式:Icd=IPcd+IScd。
二、下面是本发明方法的一个具体实施例:
1.输入激光角反射器远场衍射强度数值仿真的初始参数,包括输入激光角反射器的口径和倒棱宽度、二面直角误差、反射面的面形误差和底面的面形误差、激光角反射器材料的折射率和反射面膜层的折射率、入射激光束的入射角和方位角、入射光场的p波与s波分量的复振幅分布等,其参数名称、符号与数值见表1所示;
表1激光角反射器远场衍射强度数值仿真的初始参数
2.计算激光束在激光角反射器内部的光线方向,包括入射激光矢量R、折射光矢量A、三个反射光矢量Bm、Cm、Dm和出射激光矢量Em,m=1、2、3…、6;
3.计算激光束与激光角反射器表面的交点坐标,包括入射激光束和出射激光束在底面的交点坐标(xlab,ylab,zlab)和(xcmab,ycmab,zcmab),激光束与激光角反射器反射面的交点坐标(xomab,yomab,zomab)、(xtmab,ytmab,ztmab)和(xhmab,yhmab,zhmab),其中,m=1、2、3…、6,a=1、2、3…、91,b=1、2、3…、91;
4.基于反射面交点必须落入反射面区域和底面交点必须落入底面区域的限制条件,计算激光角反射器的衍射区域范围(Ha,Hb),其分布见图4所示;
5.计算出射激光束p波分量和s波分量的附加复振幅分布PEm和SEm,m=1、2、3…、6,PEm和SEm的模和相位分别见图5和图6所示;
6.计算接收屏位置处的激光角反射器远场衍射强度分布Icd,即101×101个离散点的强度值,其等高线分布规律见图7所示。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (6)
1.一种激光角反射器远场衍射强度的仿真方法,其特征在于,包括以下过程:首先根据激光角反射器远场衍射强度数值仿真的初始参数计算激光束在激光角反射器内部的光线方向;然后利用光线方向计算激光束与激光角反射器表面的交点坐标;其次根据各交点坐标的限制条件计算激光角反射器的衍射区域范围;再次利用激光角反射器加工误差和入射激光偏振态计算出射激光束的附加复振幅分布;最后根据衍射区域范围和出射激光束的附加复振幅分布的结果,计算接收屏位置处的激光角反射器远场衍射强度分布;
所述的激光角反射器远场衍射强度数值仿真的初始参数包括:
激光角反射器的几何参数,包括口径和倒棱宽度;
激光角反射器的加工误差,包括二面直角误差、反射面的面形误差和底面的面形误差;
激光角反射器材料的折射率和反射面膜层的折射率,若未反射面镀膜,则膜层折射率置为1;
入射激光束的入射条件,包括入射激光束的入射角和方位角,入射光场p波与s波分量的复振幅分布;
入射激光束的波长和用于激光角反射器面形检测的激光波长;
接收屏的大小及其与激光角反射器之间的距离;
激光角反射器底面和接收屏细分的点数。
2.根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法,其特征在于,所述计算激光束在激光角反射器内部的光线方向,包括以下子步骤:
步骤3.1,计算入射激光束的入射矢量R和折射光矢量A:
式中,和θ分别为入射激光束在激光角反射器底面坐标系的入射角和方位角,N0为激光角反射器底面的法线方向,M=[M1;M2;M3], refr函数为光线折射定律的矢量形式,R1、R2和R3分别为入射光矢量在激光角反射器底面坐标系下的分量,A1、A2和A3分别为折射光矢量在激光角反射器反射面坐标系下的分量;
步骤3.2,计算激光束在激光角反射器内部经过6种不同反射次序后所对应的光线方向矢量Bm、Cm、Dm和Em,m=1、2、3…、6:
Bm=γ1A,m=1、2;Bm=γ2A,m=3、4;Bm=γ3A,m=5、6;
Cm=γ2Bm,m=1、6;Cm=γ3Bm,m=2、4;Cm=γ1Bm,m=3、5;
Dm=-A;Em=R;
其中,
3.根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法,其特征在于,所述计算激光束与激光角反射器表面的交点坐标,包括以下子步骤:
步骤4.1,将入射激光束细分为nx×ny根激光线,每根激光线在激光角反射器底面上的横坐标xfab和纵坐标yfab分别满足:
xfab=-r+2(a-1)r/nx,a=1、2、3…、nx;
yfab=-r+2(b-1)r/ny,b=1、2、3…、ny;
其中,r为激光角反射器底面的半径;
步骤4.2,利用激光角反射器底面坐标系和反射面坐标系之间的转换关系,获得每根激光线在反射面坐标系下的坐标(xlab,ylab,zlab):
[xlab,ylab,zlab]T=M[xfab,yfab,0]T+Mr,
式中,L为激光角反射器的棱长,通常情况下,
步骤4.3,根据入射激光束细化光线的起始坐标以及激光束在激光角反射器内部的矢量方向,计算激光束与激光角反射器底面的交点坐标:
[xomab,yomab,zomab]T=g(βk,Cm)·[xlab,ylab,zlab]T;
[xtmab,ytmab,ztmab]T=g(βj,Bm)·[xomab,yomab,zomab]T;
[xhmab,yhmab,zhmab]T=g(βi,A)·[xtmab,ytmab,ztmab]T;
<mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>xc</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>yc</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>zc</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>{</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&beta;</mi>
<mn>4</mn>
</msub>
<mo>,</mo>
<mi>A</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>xh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>yh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>zh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>+</mo>
<mfrac>
<mi>L</mi>
<mrow>
<mi>A</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>&beta;</mi>
<mn>4</mn>
</msub>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>A</mi>
<mn>3</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>r</mi>
</msub>
<mo>}</mo>
</mrow>
式中,(xomab,yomab,zomab)、(xtmab,ytmab,ztmab)和(xhmab,yhmab,zhmab)分别表示激光束与激光角反射器反射面的交点坐标,(xcmab,ycmab,zcmab)表示激光束与激光角反射器底面的交点坐标;βi、βj、βk和β4表示系数矩阵,下标i、j和k的数值与光线反射次序有关,它们取值为:
当m=1、2时,i=1;当m=3、4时,i=2;当m=5、6时,i=3;
当m=3、5时,j=1;当m=1、6时,j=2;当m=2、4时,j=3;
当m=4、6时,k=1;当m=2、5时,k=2;当m=1、3时,k=3;
系数矩阵β1=[1,0,0],β2=[0,1,0],β3=[0,0,1],β4=[1,1,1],函数g(β,J)满足关系式:
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>&beta;</mi>
<mo>,</mo>
<mi>J</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>J</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>x</mi>
</msub>
<mo>+</mo>
<msub>
<mi>J</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>y</mi>
</msub>
<mo>+</mo>
<msub>
<mi>J</mi>
<mi>z</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>z</mi>
</msub>
</mrow>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>J</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>y</mi>
</msub>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>z</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>J</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>J</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>J</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>x</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>J</mi>
<mi>x</mi>
</msub>
<msub>
<mi>l</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>J</mi>
<mi>z</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>J</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>J</mi>
<mi>z</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>x</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>J</mi>
<mi>z</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>J</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>x</mi>
</msub>
<mo>+</mo>
<msub>
<mi>J</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&beta;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,β和J表示函数的输入量,β=[βx,βy,βz],βx、βy和βz分别为β的分量;J=×[Jx,Jy,Jz],Jx、Jy和Jz分别为J的分量。
4.根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法,其特征在于,所述计算激光角反射器的衍射区域范围,包括以下子步骤:
步骤5.1,激光束与反射面的交点必须限制在激光角反射器反射面范围内,则符合该条件的激光线序列号集合(Ra,Rb)必须满足以下形式:
(Ra,Rb)={(a,b)|T1∩T2∩T3}
其中,符号“∩”表示交集运算,集合T1、T2和T3表达式为:
式中,l表示激光角反射器的倒棱宽度;
步骤5.2,入射激光线和激光束与底面交点必须限制在激光角反射器底面范围内,则符合该条件的激光线序列号集合(Fa,Fb)必须满足:
<mrow>
<mo>(</mo>
<mi>F</mi>
<mi>a</mi>
<mo>,</mo>
<mi>F</mi>
<mi>b</mi>
<mo>)</mo>
<mo>=</mo>
<mo>{</mo>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>|</mo>
<mo>(</mo>
<msubsup>
<mi>xf</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>yf</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>&le;</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
<mo>&cap;</mo>
<mo>(</mo>
<msubsup>
<mi>xc</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>yc</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>&le;</mo>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
<mo>}</mo>
<mo>,</mo>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3...</mn>
<mo>,</mo>
<mn>6</mn>
<mo>;</mo>
</mrow>
步骤5.3,符合反射区域和底面区域限制要求的激光线所占的底面范围,即为激光角反射器的衍射区域(Ha,Hb),其表达形式为:
(Ha,Hb)=(Ra,Rb)∩(Fa,Fb)。
5.根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法,其特征在于,所述计算出射激光束的附加复振幅分布,包括以下子步骤:
步骤6.1,计算由激光角反射器反射面面形误差所导致的出射激光束附加光程ωrm(m=1、2、3...、6),它满足以下关系:
<mrow>
<msub>
<mi>&omega;r</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<mfrac>
<mrow>
<mn>3</mn>
<mrow>
<mo>|</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
</mrow>
<msub>
<mi>&xi;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>n</mi>
<mi>r</mi>
</msub>
<msub>
<mi>&lambda;</mi>
<mi>c</mi>
</msub>
</mrow>
<mrow>
<mn>4</mn>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
式中,ξi(i=1、2、3)表示第i个反射面的面形光圈值,nr为表示激光角反射器材料的折射率,λc为面形检测的激光波长,Pi(i=1、2、3)可以表示为:
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>yo</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>zo</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>yt</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>zt</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>3</mn>
<mo>,</mo>
<mn>5</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>yh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>zh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>4</mn>
<mo>,</mo>
<mn>6</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>xo</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>zo</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>xt</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>zt</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>6</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>xh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>yh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>2</mn>
<mo>,</mo>
<mn>5</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>xo</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>yo</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>5</mn>
<mo>,</mo>
<mn>6</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>xt</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>yt</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>2</mn>
<mo>,</mo>
<mn>4</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>xh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>yh</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中,
步骤6.2,计算由激光角反射器底面面形误差所导致的出射激光束附加光程ωfm(m=1、2、3...、6),满足以下关系:
式中,ξf为激光角反射器底面的面形光圈值,α为入射激光束在底面的折射角,
步骤6.3,计算由二面直角误差所导致的出射激光束附加光程,包括以下子过程:
(6.3.1)计算存在二面直角误差时,经过激光角反射器反射面三次反射后的反射光束的附加矢量dFm,m=1、2、3…、6:
dF1=[2(A2δ1+A3δ2),2(A3δ3-A1δ1),-2(A2δ3+A1δ2)]T;
dF2=[2(A3δ2+A2δ1),-2(A3δ3+A1δ1),2(A2δ3-A1δ2)]T;
dF3=[2(A3δ2-A2δ1),2(A1δ1+A3δ3),-2(A1δ2+A2δ3)]T;
dF4=-dF2;dF5=-dF3;dF6=-dF1,
式中,δ1、δ2和δ3为激光角反射器三个反射面两两之间的二面直角误差;
(6.3.2)计算反射光束的附加矢量dFm经过激光角反射器底面折射后,出射激光束的附加矢量dRm:
dRm=[dRmx,dRmy,dRmz]T=refr(dFm,-N0),m=1、2、3…、6;
(6.3.3)计算由二面直角误差导致的出射激光束附加光程ωam(m=1、2、3...、6),满足以下关系:
<mrow>
<msub>
<mi>&omega;a</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<mi>r</mi>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>dR</mi>
<mrow>
<mi>m</mi>
<mi>x</mi>
</mrow>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>dR</mi>
<mrow>
<mi>m</mi>
<mi>y</mi>
</mrow>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3...</mn>
<mo>,</mo>
<mn>6</mn>
<mo>;</mo>
</mrow>
步骤6.4,计算由入射激光偏振态所导致的出射激光束的附加复振幅分布,包括以下子过程:
(6.4.1)计算入射激光束经过一次折射和三次反射后,在激光角反射器底面和反射面上的p波和s波的矢量方向,满足以下关系式:
<mrow>
<mi>s</mi>
<mi>r</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>A</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<mi>A</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
<mi>p</mi>
<mi>r</mi>
<mo>=</mo>
<mi>A</mi>
<mo>&CircleTimes;</mo>
<mi>s</mi>
<mi>r</mi>
<mo>;</mo>
<msub>
<mi>sf</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>A</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<mi>A</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>pf</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mi>A</mi>
<mo>&CircleTimes;</mo>
<msub>
<mi>sf</mi>
<mi>m</mi>
</msub>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>se</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>B</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>B</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>pe</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>se</mi>
<mi>m</mi>
</msub>
<mo>;</mo>
<msub>
<mi>st</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>B</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>B</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>pt</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>st</mi>
<mi>m</mi>
</msub>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>su</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>C</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>pu</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>su</mi>
<mi>m</mi>
</msub>
<mo>;</mo>
<msub>
<mi>sv</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>C</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>3</mn>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>pv</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>sv</mi>
<mi>m</mi>
</msub>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>sg</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>3</mn>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>pg</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msub>
<mi>D</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>sg</mi>
<mi>m</mi>
</msub>
<mo>;</mo>
<msub>
<mi>sh</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>4</mn>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>D</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>N</mi>
<mn>4</mn>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>ph</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msub>
<mi>D</mi>
<mi>m</mi>
</msub>
<mo>&CircleTimes;</mo>
<msub>
<mi>sh</mi>
<mi>m</mi>
</msub>
<mo>;</mo>
</mrow>
式中,(pr,sr)和(pfm,sfm)分别为折射光矢量A在底面和第一次反射面上的p波和s波的矢量方向,(pem,sem)和(ptm,stm)分别为反射光矢量Bm在第一次反射面和第二次反射面上的p波和s波的矢量方向,(pum,sum)和(pvm,svm)分别为反射光矢量Cm在第二次反射面和第三次反射面上的p波和s波的矢量方向,(pgm,sgm)和(phm,shm)分别为反射光矢量Dm在第三次反射面和底面上的p波和s波的矢量方向,下标m为正整数,m=1、2、3…、6,符号表示向量的叉乘运算,N4为激光出射时底面的法线方向,N4=-N0,N1、N2和N3为反射面的法线方向,它们满足以下形式:
N1=[-1,0,0]T,m=1、2;N1=[0,-1,0]T,m=3、4;N1=[0,0,-1]T,m=5、6;
N2=[-1,0,0]T,m=3、5;N2=[0,-1,0]T,m=1、6;N2=[0,0,-1]T,m=2、4;
N3=[-1,0,0]T,m=4、6;N3=[0,-1,0]T,m=2、5;N3=[0,0,-1]T,m=1、3;
(6.4.2)计算折射光和反射光在激光角反射器表面的转换矩阵Γt(t=1、2、3、4),满足以下关系式:
式中,符号“⊙”表示矢量的点乘运算;
(6.4.3)计算折射光和反射光的p波和s波分量的反射系数矩阵和折射系数矩阵,满足以下关系式:
<mrow>
<mi>T</mi>
<mi>f</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>tf</mi>
<mi>p</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>tf</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<mi>T</mi>
<mi>t</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>tl</mi>
<mi>p</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>tl</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
<mrow>
<mi>R</mi>
<mi>f</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>rf</mi>
<mi>p</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>rf</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<mi>Re</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>re</mi>
<mi>p</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>re</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<mi>R</mi>
<mi>t</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>rt</mi>
<mi>p</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>rt</mi>
<mi>s</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
式中,Tf和Tl分别为第一次折射和最后一次折射的折射系数矩阵,Rf、Re和Rt分别为第一次反射、第二次反射和第三次反射的反射系数矩阵;折射系数矩阵和反射系数矩阵的非零元素满足以下关系式:
rfp=crp(if,zf);rfs=crs(if,zf);rep=crp(ie,ze);res=crs(ie,ze);
rtp=crp(it,zt);rts=crs(it,zt);
其中,(if,zf)、(ie,ze)和(it,zt)分别表示激光束在第一次反射、第二次反射和第三次反射时对应的入射角和折射角,if=arccos(A⊙N1),zf=arcsin(υsinif),ie=arccos(Bm⊙N2),ze=arcsin(υsinie),it=arccos(Cm⊙N3),zt=arcsin(υsinit),式中,υ=nr/nf,nf为激光角反射器膜层的折射率;函数ctp和cts分别表示p波和s波的透射系数的计算公式,函数crp和crs分别表示p波和s波的反射系数的计算公式,它们满足以下关系式:
其中,表示函数的输入量;
(6.4.4)计算由入射激光偏振态所导致的出射激光束p波分量和s波分量的附加复振幅APout和ASout的表达形式:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>AP</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>AS</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msub>
<mi>Tl&Gamma;</mi>
<mn>4</mn>
</msub>
<msub>
<mi>Rt&Gamma;</mi>
<mn>3</mn>
</msub>
<msub>
<mi>Re&Gamma;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>Rf&Gamma;</mi>
<mn>1</mn>
</msub>
<mi>T</mi>
<mi>f</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>EP</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>EP</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
式中,EPin和ESin分别为入射光场p波分量和s波分量的复振幅;
步骤6.5,计算出射激光束p波分量和s波分量的附加复振幅分布PEm和SEm,m=1、2、3…、6,满足以下关系式:
PEm=|PEm|exp(iφpm)=|APout|exp[iτ(ωrm+ωfm+ωam)+iωpm];
SEm=|SEm|exp(iφsm)=|ASout|exp[iτ(ωrm+ωfm+ωam)+iωsm];
式中,|PEm|和Фpm分别表示PEm的模和相位,|SEm|和Фsm分别表示SEm的模和相位,i表示虚数单位;τ=2π/λr,λr为入射激光束的波长,ωpm和ωsm分别为APout和ASout的幅角。
6.根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法,其特征在于,所述计算接收屏位置处的激光角反射器远场衍射强度分布,包括以下子步骤:
步骤7.1,将接收屏细分为μx×μy个离散接收点,每个接收点在接收屏上的横坐标xrcd和纵坐标yrcd分别满足:
xrcd=-Rw+2(c-1)Rw/μx,c=1、2、3…、μx;
yrcd=-Rw+2(d-1)Rw/μy,d=1、2、3…、μy;
其中,Rw为接收屏的半径大小;
步骤7.2,计算出射激光束的p波分量和s波分量的复振幅分布Pmab和Smab,满足以下关系式:
Pmab=exp{iτ[(dRmx+A1)xfab+(dRmy+A2)yfab]}·PEm;
Smab=exp{iτ[(dRmx+A1)xfab+(dRmy+A2)yfab]}SEm;
步骤7.3,计算出射激光束p波分量和s波分量的远场衍射强度值IPcd和IScd,满足以下形式:
式中,d0为激光角反射器与接收屏之间的距离,i表示虚数单位,xgcd和ygcd可以表示为:
步骤7.4,计算接收屏位置处的激光角反射器远场衍射强度分布Icd,满足关系式:Icd=IPcd+IScd。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610545255.5A CN105929531B (zh) | 2016-07-12 | 2016-07-12 | 一种激光角反射器远场衍射强度的数值仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610545255.5A CN105929531B (zh) | 2016-07-12 | 2016-07-12 | 一种激光角反射器远场衍射强度的数值仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105929531A CN105929531A (zh) | 2016-09-07 |
CN105929531B true CN105929531B (zh) | 2018-04-20 |
Family
ID=56827753
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610545255.5A Expired - Fee Related CN105929531B (zh) | 2016-07-12 | 2016-07-12 | 一种激光角反射器远场衍射强度的数值仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105929531B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107480360B (zh) * | 2017-08-03 | 2020-04-24 | 中北大学 | 光束细分和相界面漫反射的激光烧熔的数值计算方法 |
CN110208777B (zh) * | 2019-06-11 | 2021-07-13 | 中国科学院电子学研究所 | 精确的角反射器几何误差测量方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5379115A (en) * | 1991-01-28 | 1995-01-03 | Excel Precision | Differential interferometer |
CN103207422A (zh) * | 2013-04-02 | 2013-07-17 | 中国科学院光电研究院 | 光学角反射靶球及其制造方法 |
-
2016
- 2016-07-12 CN CN201610545255.5A patent/CN105929531B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5379115A (en) * | 1991-01-28 | 1995-01-03 | Excel Precision | Differential interferometer |
CN103207422A (zh) * | 2013-04-02 | 2013-07-17 | 中国科学院光电研究院 | 光学角反射靶球及其制造方法 |
Non-Patent Citations (2)
Title |
---|
倒棱工艺对激光角反射器衍射特性的影响;石岩,周辉,文翔;《实验室研究与探索》;20130831;第32卷(第8期);128-132、136页 * |
卫星激光角反射器的远场衍射光强研究;钟声远,徐广平,吴键;《激光与红外》;20090228;第39卷(第2期);41-46页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105929531A (zh) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108983331B (zh) | 一种半透半逆反膜片和空中显示装置 | |
CN106773489B (zh) | 一种光学投影装置及深度相机 | |
CN108957911A (zh) | 散斑结构光投影模组及3d深度相机 | |
CN206020665U (zh) | 复合光束收发的差分吸收激光雷达系统 | |
CN105629449B (zh) | 一种新型菲涅耳光学天线发射系统 | |
CN105929531B (zh) | 一种激光角反射器远场衍射强度的数值仿真方法 | |
CN109579728A (zh) | 基于全息波导的散斑结构光投影模组 | |
CN107767420A (zh) | 一种水下立体视觉系统的标定方法 | |
CN109212773A (zh) | 紧凑式散斑投影模组及深度相机 | |
CN110400519A (zh) | 背光模组及其制作方法、显示装置 | |
CN105547649A (zh) | 一种短波红外镜头杂散辐射的检测方法 | |
Gong et al. | 3-D ice shape measurements using mid-infrared laser scanning | |
CN105511012B (zh) | 一种用于实现偏振准直面光源的导光结构 | |
CN105739101A (zh) | 匀光结构及匀光系统 | |
CN201285473Y (zh) | 一种基于散射的消相干匀场装置 | |
CN114296245A (zh) | 一种拉曼光束整形装置 | |
CN114428375A (zh) | 衍射光波导及ar眼镜 | |
CN101685181A (zh) | 一种基于散射的消相干匀场装置 | |
CN105572867B (zh) | 利用双楔形板和标准角锥棱镜实现反射光束偏转的方法 | |
CN103226241A (zh) | 一种基于能量守恒定律的光学系统分析设计方法 | |
CN111273261A (zh) | 一种基于离轴入射的同轴发射与接收激光雷达 | |
CN102628713B (zh) | 基于数字微镜器件的曲率波前传感器 | |
JP2017003766A (ja) | 光学素子及び表示装置 | |
US10018577B2 (en) | Methods and systems for imaging bulk motional velocities in plasmas | |
CN108693653A (zh) | 射频全息透射/光反射式共口径复合平面波场模拟装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180420 Termination date: 20200712 |
|
CF01 | Termination of patent right due to non-payment of annual fee |