一种激光角反射器远场衍射强度的数值仿真方法
【专利摘要】本发明涉及一种激光角反射器远场衍射强度数值仿真方法，属于光学设计与仿真领域。本发明以激光角反射器实际加工过程中的二面直角误差、面形误差和倒棱等误差参数、实际使用过程中的光束倾斜效应和激光偏振态效应等因素为输入条件，以激光角反射器的衍射区域、出射激光束的复振幅分布和简化的远场衍射强度模型为基础，通过入射激光束和接收屏的离散划分以及光线的追迹方法，仿真得到综合影响因素条件下激光角反射器远场衍射强度的数值分布。所采用的仿真方法充分考虑了影响激光角反射器远场强度分布的所有因素，扩充和完善了激光角反射器远场衍射强度的仿真范畴，特别适用于远距离条件下激光角反射器的使用性能指标的综合评估。
【专利说明】
一种激光角反射器远场衍射强度的数值仿真方法
技术领域
[0001] 本发明属于光学设计与仿真领域，具体涉及一种激光角反射器远场衍射强度数值 仿真方法，适用于激光测距合作目标的性能指标综合评价。
【背景技术】
[0002] 激光角反射器是由三个相互垂直的直角反射面和一个底面组成的实心四面体，其 几何结构见图1(a)所示。平面0BC、0AC和0AB为激光角反射器的直角反射面，分别以数字1、2 和3来表示，平面ABC为激光角反射器的底面。通常情况下，激光角反射器底面采用圆形切割 以便于实际装配，见图1 (b)所示。当入射激光束从激光角反射器底面入射时，将先后经过底 面的第一次折射、反射面的三次反射和底面的再次折射，以六种不同的反射次序：1-2-3、 1-3-2、2-1-3、2-3-1、3-1 -2、3-2-1，从底面出射。对于理想的激光角反射器而 言，考虑到其特殊的几何结构，出射激光束将与入射激光束的方向平行，即理想激光角反射 器具有定向反射特性，因此，通常其作为激光测距中的合作目标，用于增加测距仪接收机位 置处的激光强度。
[0003] 激光角反射器在使用过程中，其与激光测距机之间距离远大于激光角反射器的有 效口径，因此，出射激光束必然受到衍射效应的制约，即在激光测距仪接收机位置处接收到 的激光强度是激光角反射器远场衍射强度。由于激光角反射器在加工制造过程会出现二面 直角误差和面形误差，若偏振入射激光束倾斜入射至激光角反射器底面时，则激光角反射 器远场衍射强度的分布将受到激光角反射器几何参数、光束入射条件和接收区域有关。
[0004] 建立如图2所示的坐标框架，包括激光角反射器棱坐标系和底面坐标系、接收屏坐 标系，其中，激光角反射器棱坐标系的原点位于顶点处，其xl轴、yl轴和zl轴分别与激光角 反射器直角边重合;底面坐标系的原点位于底面中心，zf轴与底面的法线方向平行，yf轴是 yl轴在底面的投影，xf?轴与yf?轴和zf?轴构成右手螺旋法则;接收屏坐标系的原点是入射激 光束中心与接收屏的交点，zr轴与入射激光束方向平行，xr轴和yr轴位于接收屏范围内。因 此，在距离激光角反射器底面do位置处所接收到的远场衍射强度可以表示为：
[0006] 其中，P为入射激光束在激光角反射器底面上的入射角，A为入射激光束的波长，k =2JT八，S为激光角反射器的有效衍射区域，EP (xf，yf)和Es (xf，yf)分别为出射激光束在激 光角反射器底面上的P波和s波光场分量，r是激光角反射器底面上的点(xf，yf)与接收屏上 的点(xr，yr)之间的距离。考虑到激光角反射器远场衍射强度的计算模型不存在解析解，并 且其受到的影响因素较多，因此，构建完善的激光角反射器远场衍射强度的综合仿真方法 是十分复杂和困难的。
[0007] 近些年来，已经有一些关于激光角反射器远场衍射强度仿真的报道，如在垂直入 射并忽略入射激光束偏振态条件下的远场衍射强度分析(聂辉，等.角锥棱镜的远场衍射特 性.光学学报，23(12): 1470-1474，2003)、在忽略面形误差和入射激光束偏振效应的远场衍 射强度数值模拟(叶一东，等.角锥后向反射器的数值模拟研究.光学学报，23(4): 485-490， 2003)、在垂直入射并忽略激光角反射器加工误差的远场衍射强度仿真(Murphy T W，et al.Polarization and far-field diffraction patterns of total internal reflection corner cubes.Applied optics,52(2): 117-126,2013.)，这些仿真方法是在 仅考虑影响激光角反射器远场衍射强度分布的部分因素前提下，利用相应的假设近似来获 得激光角反射器远场衍射强度的分布规律。采用这些方法仅能仿真得到满足某些特定激光 角反射器的几何参数或测量条件的远场衍射强度的仿真，无法实现激光角反射器实际使用 过程中的远场衍射强度的真实仿真。
【发明内容】
[0008] 本发明主要是解决现有仿真方法所存在的问题:提出一种综合考虑激光角反射器 实际加工过程中的二面直角误差、面形误差和倒棱等误差参数，实际使用过程中的光束倾 斜效应和激光偏振态效应等因素，通过入射激光束和接收屏的离散划分以及光线的追迹方 法，仿真得到非理想条件下激光角反射器远场衍射强度的数值，以实现对激光角反射器实 际使用性能的精确评估。
[0009] 本发明的上述技术问题主要是通过下述技术方案得以解决的：
[0010] -种激光角反射器远场衍射强度的仿真方法，包括以下过程:首先根据激光角反 射器远场衍射强度数值仿真的初始参数计算激光束在激光角反射器内部的光线方向；然后 利用光线方向计算激光束与激光角反射器表面的交点坐标;其次根据各交点坐标的限制条 件计算激光角反射器的衍射区域范围；再次利用激光角反射器加工误差和入射激光偏振态 计算出射激光束的附加复振幅分布;最后根据衍射区域范围和出射激光束的附加复振幅分 布的结果，计算接收屏位置处的激光角反射器远场衍射强度分布。
[0011] 在上述的一种激光角反射器远场衍射强度的仿真方法，所述的激光角反射器远场 衍射强度数值仿真的初始参数包括：
[0012] 激光角反射器的几何参数，包括口径和倒棱宽度；
[0013] 激光角反射器的加工误差，包括二面直角误差、反射面的面形误差和底面的面形 误差；
[0014] 激光角反射器材料的折射率和反射面膜层的折射率，若未反射面镀膜，则膜层折 射率置为1;
[0015] 入射激光束的入射条件，包括入射激光束的入射角和方位角，入射光场p波与s波 分量的复振幅分布；
[0016] 入射激光束的波长和用于激光角反射器面形检测的激光波长；
[0017] 接收屏的大小及其与激光角反射器之间的距离；
[0018] 激光角反射器底面和接收屏细分的点数。
[0019] 在上述的一种激光角反射器远场衍射强度的仿真方法，计算激光束在激光角反射 器内部的光线方向，包括以下子步骤：
[0020] 步骤3.1，计算入射激光束的入射矢量R和折射光矢量A:
[0021 ] |' - | sin ^ cos 0, sin ^9 sin 0, cos(p f ; A = [Ai, A2, A3 ]T=refr (-MR, No)；
[0022] 式中，货和0分别为入射激光束在激光角反射器底面坐标系的入射角和方位角，No 为激光角反射器底面的法线方向，乂 =[-U-K-l]/V? ? M = [ M 1 ； M 2 ; M 3 ],
>refr函数为光线折 射定律的矢量形式，Ri、R2和R3分别为入射光矢量在激光角反射器底面坐标系下的分量，Ai、 A2和A3分别为折射光矢量在激光角反射器反射面坐标系下的分量。
[0023] 步骤3.2,计算激光束在激光角反射器内部经过6种不同反射次序后所对应的光线 方向矢量 Bm、Cm、Dm 和 Em，m=l、2、3…、6:
[0024] Bm= y iA，m=l、2;Bm= y 2A，m=3、4;Bm= y 3A，m = 5、6;
[0025] Cm= 丫 2Bm，m=l、6;Cm= 丫 3Bm，m = 2、4;Cm= y iBm，m = 3、5;
[0026] Dm=-A；Em=R； "-1： 0 on 「100] 「:i 0 0 -
[0027] 其中，A: 〇 1 0 ? f2 = 0 -1 0 ? f3 = 〇 1 〇 。 0 0: 1」 [_0 0: 1」 l0: 0 -1
[0028] 在上述的一种激光角反射器远场衍射强度的仿真方法，计算激光束与激光角反射 器表面的交点坐标，包括以下子步骤：
[0029] 步骤4.1，将入射激光束细分为nxXny根激光线，每根激光线在激光角反射器底面 上的横坐标Xfab和纵坐标yfab分别满足：
[0030] xfab = -r+2(a-l)r/nx，a = l、2、3...、nx;
[0031 ] yfab = -r+2(b-l)r/ny，b = l、2、3...、ny;
[0032] 其中，r为激光角反射器底面的半径。
[0033] 步骤4.2,利用激光角反射器底面坐标系和反射面坐标系之间的转换关系，获得每 根激光线在反射面坐标系下的坐标(xlab，yl ab，zlab):
[0034] [xlab,ylab,zlab]T = M[xfab,yfab,0]T+Mr,
[0035] 式中
山为激光角反射器的棱长，通常情况下，l = 。
[0036] 步骤4.3,根据入射激光束细化光线的起始坐标以及激光束在激光角反射器内部 的矢量方向，计算激光束与激光角反射器底面的交点坐标：
[0037] [ XOmab , y0mab , ZOmab ] -g (0k , Cm) * [xlab,ylab,zlab]；
[0038] [ Xtmab , y tmab , Ztmab ] -g(0j)Bm) * [ XOmab , y0mab , ZOmab ]；
[0039] [ xhmab , yhmab , zhmab ] -g(0i,A) * [xtmab , ytmab , Ztmab]；
[0041 ]式中，（XOmab，y〇mab，ZOmab )、（ Xtmab，y tmab，Z tmab )和（xhmab，yhmab，zhmab )分力丨」表激 束与激光角反射器反射面的交点坐标，（XCmab，yCmab，ZCmab)表示激光束与激光角反射器底面 的交点坐标。队、比、&和04表示系数矩阵，下标i、j和k的数值与光线反射次序有关，它们取 值为：
[0042] 当m = 1、2时，i = 1;当m = 3、4时，i = 2;当m = 5、6时，i = 3;
[0043] 当m = 3、5时，j = l;当m=l、6时，j = 2;当m = 2、4时，j = 3;
[0044] 当m = 4、6时，k = 1;当m = 2、5时，k = 2;当m = 1、3时，k = 3;
[0045] 系数矩阵 01 = [1，0,0]，02=[0，1，0]，&3=[0,0，1]，04=[1，1，1]，函数8(0，了)满足 关系式：
[0047]其中，时PJ表示函数的输入量，P H队]，瓜為和队分别为抑勺分量;J = [ Jx， Jy，Jz]，JX、Jy和Jz分别为J的分量。
[0048] 在上述的一种激光角反射器远场衍射强度的仿真方法，计算激光角反射器的衍射 区域范围，包括以下子步骤：
[0049] 步骤5.1，激光束与反射面的交点必须限制在激光角反射器反射面范围内，则符合 该条件的激光线序列号集合(Ra，Rb)必须满足以下形式：
[0050] (Ra,Rb) = {(a,b) |TinT2nT3}
[0051] 其中，符号" n "表示交集运算，集合Ti、t#pt3表达式为： wmuA>/&zomt(6>/ m=l,2
[0052] fx = < xomah > / & zomnh > / m = 3,4 ； x〇mab>!&yomab>! m = 5,6 m = 3,5 [yKab^-^^mub^1 m = 4,6
[0053] 7； = < ximah > / & ztmab > / /w = 1,6 ; T3=< xhmah > / & zhnMb >1 in = 2,5 , xtmah^! rn = 2A {xhmab>! &yhnwb> I w = l,3
[0054]式中，1表示激光角反射器的倒棱宽度。
[0055] 步骤5.2,入射激光线和激光束与底面交点必须限制在激光角反射器底面范围内， 则符合该条件的激光线序列号集合(Fa，Fb)必须满足：
[0056] (Fo,F/))-{(a,6)|(x/；t + j/；fr< r2)n(^^ m=\ s 2, 3,.., 6,
[0057] 步骤5.3,符合反射区域和底面区域限制要求的激光线所占的底面范围，即为激光 角反射器的衍射区域(Ha，Hb)，其表达形式为：
[0058] (Ha，Hb) = (Ra，Rb) n (Fa，Fb)。
[0059] 在上述的一种激光角反射器远场衍射强度的仿真方法，计算出射激光束的附加复 振幅分布，包括以下子步骤：
[0060] 步骤6.1，计算由激光角反射器反射面面形误差所导致的出射激光束附加光程co rm(m=l、2、3. . .、6)，它满足以下关系：
[0062]式中，l(i = l、2、3)表示第i个反射面的面形光圈值，nr为表示激光角反射器材料 的折射率，A。为面形检测的激光波长，Pi (i = 1、2、3)可以表示为：
[0066] 其中，A = ^"4。
[0067] 步骤6.2,计算由激光角反射器底面面形误差所导致的出射激光束附加光程《 fm (m=l、2、3. . .、6)，满足以下关系：
[0069]式中，If为激光角反射器底面的面形光圈值，a为入射激光束在底面的折射角， a = arcsin (siiu/?/// .) e
[0070]步骤6.3,计算由二面直角误差所导致的出射激光束附加光程，包括以下子过程： [0071 ] (6.3.1)计算存在二面直角误差时，经过激光角反射器反射面三次反射后的反射 光束的附加矢量dFm，m=l、2、3…、6:
[0072] dFi = [2(A25i+A352)，2(A353-Ai5i)，_2(A253+Ai52) ]T;
[0073] dF2 = [2(A352+A25i)，_2(A353+Ai5i)，2(A253~Ai52) ]T;
[0074] dF3 = [2(A352-A25i)，2(Ai5i+A353)，-2(Ai52+A253) ]T;
[0075] dF4=-dF2 ；dF5 = -dF3； dF6 = -dFi,
[0076]式中，和S3为激光角反射器三个反射面两两之间的二面直角误差。
[0077] (6.3.2)计算反射光束的附加矢量dFmg过激光角反射器底面折射后，出射激光束 的附加矢量dRm:
[0078] dRm= [dRmx，dRmy，dRmZ]T = refr(dFm，-No)，m= 1、2、3...、6〇
[0079] (6.3.3)计算由二面直角误差导致的出射激光束附加光程《 am(m= 1、2、3. . .、6)， 满足以下关系：
[0080] 〇}am =-^2rVdna^a'Riir co^.0 + dR,ll.s\n0^j, M=1 > 2> 3..,\ 6〇
[0081 ]步骤6.4，计算由入射激光偏振态所导致的出射激光束的附加复振幅分布，包括以 下子过程：
[0082] (6.4.1)计算入射激光束经过一次折射和三次反射后，在激光角反射器底面和反 射面上的P波和s波的矢量方向，满足以下关系式：
[0087] 式中，（pr，sr)和(pfm，sfm)分别为折射光矢量A在底面和第一次反射面上的p波和s 波的矢量方向，（pem，sem)和(ptm，stm)分别为反射光矢量Bm在第一次反射面和第二次反射面 上的P波和S波的矢量方向，（pUm，SUm)和(pVm，SVm)分别为反射光矢量Cm在第二次反射面和第 三次反射面上的P波和S波的矢量方向，化8^吨〇和化^，8^)分别为反射光矢量0?在第三次 反射面和底面上的P波和S波的矢量方向，下标m为正整数，m=l、2、3…、6,符号"?"表示向 量的叉乘运算，N4为激光出射时底面的法线方向，仏=_他，他、犯和犯为反射面的法线方向， 它们满足以下形式：
[0088] Ni=[-l，0,0]T，m=l、2;Ni = [0，-l，0]T，m=3、4;Ni=[0,0，-l]T，m = 5、6;
[0089] 他=卜1，0,0]T，m=3、5;N2 = [0 厂1，0]T，m=l、6;N2 = [0,0 厂l]T，m = 2、4;
[0090] N3=[-l，0,0]T，m=4、6;N3 = [0，-l，0]T，m=2、5;N3=[0,0，-l]T，m=l、3。
[0091] (6.4.2)计算折射光和反射光在激光角反射器表面的转换矩阵1^(〖=1、2、3、4)， 满足以下关系式： 「_91 r Jpr〇pfm sr〇Pf'?]. r \Pem(^Ptmsem〇Ptm]^
[prQsfm sr0sfm\ 2 [pem0stm sem0stm\
[0093 ] r = pUm opVm SUm〇pVi" . r = p8m °phm S8m °phm ? lPum〇svm sum〇svmj， 4 [pgm0shm % G)為,」'
[0094]式中，符号" 〇 "表示矢量的点乘运算。
[0095] (6.4.3)计算折射光和反射光的p波和s波分量的反射系数矩阵和折射系数矩阵， 满足以下关系式： 「r/； 0]
[0096] Tf= 1 I Tl= 1 ; L〇 c」 L():( rf 0 1 「re" 0 .0:
[0097] Jp' ； Re= p m = 1 5 _〇 ，:d _ L° L° rh_
[0098] 式中，Tf?和II分别为第一次折射和最后一次折射的折射系数矩阵，Rf、Re和Rt分别 为第一次反射、第二次反射和第三次反射的反射系数矩阵。折射系数矩阵和反射系数矩阵 的非零元素满足以下关系式：
[0099] fp=^ctp{jp,a)i tfs^cts{(p,a)% tlp^-dp{a.(p) t tl^-cts{cx,(p)%
[0100] rfP = crp(if ,zf) ；rfs = crs(if ,zf) ；reP = crp( ie, ze) ；res = crs(ie ,ze)；
[0101 ] rtP = crp(it,zt) ；rts = crs(it,zt)；
[0102]其中，（if，zf)、（ie，ze)和（it，zt)分别表示激光束在第一次反射、第二次反射和 第三次反射时对应的入射角和折射角，if = arccos(A 〇 Ni)，zf = arcsin(usin if)，ie = arccos(Bm 〇 N2) ,ze = arcsin(usin ie), it = arccos(Cm 〇 N3), zt = arcsin(usin it)，式 中，u = nr/nf，nf为激光角反射器膜层的折射率。函数ctp和cts分别表示p波和s波的透射系 数的计算公式，函数crp和crs分别表示p波和s波的反射系数的计算公式，它们满足以下关 系式：
[0105] 其中，(抖，识2)表示函数的输入量。
[0106] (6.4.4)计算由入射激光偏振态所导致的出射激光束p波分量和s波分量的附加复 振幅AP〇ut和AS〇ut的表达形式： 1 [EB ~
[0108] 式中，EPin和ESin分别为入射光场p波分量和s波分量的复振幅。
[0109] 步骤6.5，计算出射激光束p波分量和S波分量的附加复振幅分布PEjPSEm，m= 1、2、 3…、6,满足以下关系式：
[0110] PEm= | PEm| exp(i <}) Pm) = | APout | exp[ix( to rm+w fm+ o am)+i w Pm]；
[0111] SEm= | SEm| exp(i <}) Sm) = | AS〇ut | exp[ix( to rm+w fm+ o am)+i w sm]；
[0112] 式中，|PEm|和〇_分别表示PEm的模和相位，|SE m|和〇观分别表示SEm的模和相位，i 表示虚数单位。1 = 231/%，\为入射激光束的波长，〇^"1和《8111分别为4?。111；和43。1 11；的幅角。
[0113] 在上述的一种激光角反射器远场衍射强度的仿真方法，计算接收屏位置处的激光 角反射器远场衍射强度分布，包括以下子步骤：
[0114] 步骤7.1，将接收屏细分为fcXb个离散接收点，每个接收点在接收屏上的横坐标 xrcxi和纵坐标yrcxi分别满足：
[0115] xrcd = -Rw+2(c-l)Rw/lix，c = 1、2、3...、yx;
[0116] yrcd = -Rw+2(d-l)Rw/liy，d = 1、2、3...、liy;
[0117]其中，1为接收屏的半径大小。
[01 18]步骤7.2,计算出射激光束的P波分量和S波分量的复振幅分布Pmab和Smab，满足以下 关系式：
[0119] Pmab = exp{iT[(dRmx+Al)xfab+(dRmy+A2)yfab]} 1 PEm ；
[0120] Smab = exp{ix[(dRmx+Al)xfab+(dRmy+A2)yfab]} 1 SEm ；
[0121]步骤7.3，计算出射激光束p波分量和s波分量的远场衍射强度值IPcd和IScd，满足 以下形式：
[0124] 式中，do为激光角反射器与接收屏之间的距离，i表示虚数单位，xgcd和ygcd可以表 示为：
[0125] Xg.cil = cos^9cos0 '.xred - sin 0 yrctl + ditR^ :
[0126] 嫂= cospsin沒?尤匕 + cod ? .v/:" +
[0127] 步骤7.4,计算接收屏位置处的激光角反射器远场衍射强度分布Id，满足关系式： led 一 IPcd+I Sed 〇
[0128] 本发明具有如下优点：以激光角反射器的实际几何参数、二面直角误差、反射面与 底面的面形误差等为仿真输入，在考虑倾光束斜入射和激光束偏振态的情况下，利用光线 的追迹方法计算得到激光角反射器实际的衍射区域范围和出射激光束的复振幅分布，并通 过离散的数值模拟出激光角反射器的远场衍射强度分布。所采用的仿真方法不仅仅以光束 入射条件、二面直角误差、面形误差或激光偏振态等部分因素作为分析对象，扩充和完善了 激光角反射器远场衍射强度的仿真范畴，特别适用于远距离条件下激光角反射器的使用性 能指标的综合评估。
【附图说明】
[0129] 图1(a)是激光角反射器的几何结构示意图。
[0130] 图1(b)是圆形切割方式的激光角反射器结构。
[0131 ]图2是激光角反射器远场衍射仿真计算过程中所涉及的几种坐标框架。
[0132] 图3是激光角反射器远场衍射强度数值仿真计算的流程。
[0133] 图4是实施例激光角反射器的衍射区域分布。
[0134] 图5是实施例出射激光束p波分量和s波分量的附加复振幅的模分布。
[0135] 图6(a)是实施例出射激光束p波分量附加复振幅的相位分布。
[0136] 图6(b)是实施例出射激光束s波分量附加复振幅的相位分布。
[0137] 图7是实施例激光角反射器远场衍射强度的等高线分布。
【具体实施方式】
[0138] 下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。
[0139] 实施例：
[0140] -、首先介绍一下本发明的具体方法流程，主要包括以下步骤：
[0141] 1.输入激光角反射器远场衍射强度数值仿真的初始参数，包括下述过程：
[0142] (1.1)输入激光角反射器的几何参数，包括口径和倒棱宽度；
[0143] (1.2)输入激光角反射器的加工误差，包括二面直角误差、反射面的面形误差和底 面的面形误差；
[0144] (1.3)输入激光角反射器材料的折射率和反射面膜层的折射率，若未反射面镀膜， 则膜层折射率置为1;
[0145] (1.4)输入入射激光束的入射条件，包括入射激光束的入射角和方位角，入射光场 P波与S波分量的复振幅分布；
[0146] (1.5)输入入射激光束的波长和用于激光角反射器面形检测的激光波长；
[0147] (1.6)输入接收屏的大小及其与激光角反射器之间的距离；
[0148] (1.7)输入激光角反射器底面和接收屏细分的点数。
[0149] 2.计算激光束在激光角反射器内部的光线方向，包括下述过程：
[0150] (2.1)计算入射激光束的入射矢量R和折射光矢量A:
[0151 ] /? = [Sl，/?:，A']/=[sinpa)s0，sin(/?sinAcos^)] r;:A=[Ai，A2，A3]T=ref;r(-MR，N〇);
[0152] 式中，识和0分别为入射激光束在激光角反射器底面坐标系的入射角和方位角，No为 激光角反射器底面的法线方向，％ = [-1，-1-1丨/ #，M= [Mi ;M2 ;M3]，
，refr函数为光线折射定律的矢量形式，Ri、R2 和R3分别为入射光矢量在激光角反射器底面坐标系下的分量，Al、A2和A3分别为折射光矢量 在激光角反射器反射面坐标系下的分量。
[0153] (2.2)计算激光束在激光角反射器内部经过6种不同反射次序后所对应的光线方 向矢量 Bm、Cm、Dm 和 Em，m=l、2、3…、6:
[0154] Bm= y iA，m=l、2;Bm= y 2A，m=3、4;Bm= y 3A，m = 5、6;
[0155] Cm= y 2Bm，m=l、6;Cm= y 3Bm，m = 2、4;Cm= y lBm，m = 3、5;
[0156] Dm=-A；Em=R； "-1 0 0] 「10 0] 「1 0 .0 -
[0157] 其中，K = 〇 l () ;: 〇 -l 〇 ; r3 = 〇 l 〇 。 0 0 1」 [_0 0 1」 |_0 0 -1
[0158] 3.计算激光束与激光角反射器表面的交点坐标，包括下述过程：
[0159] (3.1)将入射激光束细分为91X91根激光线，每根激光线在激光角反射器底面上 的横坐标xfab和纵坐标yf ab分别满足：
[0160] xfab = -r+2(a_l )r/91，a=l、2、3...、91;
[0161 ] yfab = -r+2(b_l )r/91，b=l、2、3…、91;
[0162] 其中，r为激光角反射器底面的半径。
[0163] (3.2)利用激光角反射器底面坐标系和反射面坐标系之间的转换关系，获得每根 激光线在反射面坐标系下的坐标(11 ￡11)，71￡11)，21￡11 )):
[0164] [xlab,ylab,zlab]T = M[xfab,yfab,0]T+Mr,
[0165] 式中，
丄为激光角反射器的棱长，通常情况下，Z 。
[0166] (3.3)根据入射激光束细化光线的起始坐标以及激光束在激光角反射器内部的矢 量方向，计算激光束与激光角反射器底面的交点坐标：
[0167] [ XOmab，y0mab，ZOmab ] -g (0k，Cm) ? [ xlab，y lab，Z lab ];
[01 68] [ Xtmab , y tmab , Ztmab ] -g(0j)Bm) * [ XOmab , y0mab , ZOmab ]；
[0169] [ xhmab , yhmab , zhmab ] -g(0i,A) * [xtmab , ytmab , Ztmab]；
[01 71 ]式中，（XOmab，yOmab，ZOmab )、（ Xtmab，ytmab，Z tmab )和（xhmab，y hmab，zhmab )分力 U 表 激 束与激光角反射器反射面的交点坐标，（XCmab，yCmab，ZCmab)表示激光束与激光角反射器底面 的交点坐标。队、比、&和04表示系数矩阵，下标i、j和k的数值与光线反射次序有关，它们取 值为：
[0172]当m = 1、2时，i = 1;当m = 3、4时，i = 2;当m = 5、6时，i = 3;
[0173 ]当m = 3、5时，j = 1;当m = 1、6时，j = 2;当m = 2、4时，j = 3;
[0174] 当m = 4、6时，k = 1;当m = 2、5时，k = 2;当m = 1、3时，k = 3;
[0175] 系数矩阵 01 = [1，0,0]，02=[0，1，0]，&3=[0,0，1]，04=[1，1，1]，函数8(0，了)满足 关系式：
[0177] 其中，时PJ表示函数的输入量，0 H队]，瓜為和队分别为抑勺分量;J = [ Jx， Jy，Jz]，JX、Jy和Jz分别为J的分量。
[0178] 4.计算激光角反射器的衍射区域范围，包括下述过程：
[0179] (4.1)激光束与反射面的交点必须限制在激光角反射器反射面范围内，则符合该 条件的激光线序列号集合(Ra，Rb)必须满足以下形式：
[0180] (Ra,Rb) = {(a,b) |TinT2nT3}
[0181] 其中，符号" n"表示交集运算，集合、T#PT3表达式为： 3'^,,；, ^/&^!(!/) >/ m = \,2
[0182] - < xomab > / &zomab > l /? = 3,4 ? ^〇mub ^ 1 & > / m = 5,6 yimih>! &ztmab>! w = 3,5 \yhmab>lSLzhmab>l m = 4,6
[0183] T2 = < xtmab > l &ztmab >1 /? = 1,6 ; T3 = < xhmob > l &zhntab >1 m = 2,5 , m = 2,4 [xhmib>! &yhmab> l m = l,3
[0184] 式中，1表示激光角反射器的倒棱宽度。
[0185] (4.2)入射激光线和激光束与底面交点必须限制在激光角反射器底面范围内，则 符合该条件的激光线序列号集合(Fa，Fb)必须满足：
[0186] (Fa,Fh) = {(^/, h) \ (.v/；,； + yf；h < r ) R (.vt：；,.,, + rr,；,,, < r )}，m= 1 , 2> 3..., 6 〇
[0187] (4.3)符合反射区域和底面区域限制要求的激光线所占的底面范围，即为激光角 反射器的衍射区域(Ha，Hb)，其表达形式为：
[0188] (Ha，Hb) = (Ra，Rb) n (Fa，Fb)。
[0189] 在本步骤中，激光角反射器的衍射区域范围综合考虑了激光角反射器圆形切割和 倒棱等因素的影响。
[0190] 5.计算出射激光束的附加复振幅分布，包括下述过程：
[0191] (5.1)计算由激光角反射器反射面面形误差所导致的出射激光束附加光程c〇rm( m = 1、2、3. . .、6)，它满足以下关系：
[0193] 式中，l(i = l、2、3)表示第i个反射面的面形光圈值，nr为表示激光角反射器材料 的折射率，A。为面形检测的激光波长，Pi (i = 1、2、3)可以表示为： {y°mab _A) + (ZOmc,b ~ A) m -
[0194] Pl = <(ytmab-A)1 +(ztnmb- Af m - 3,5 ; {yKab - A)2 + (zhmab - A)2 m = 4,6 (?TO,- - A)2 + (功,--A)2 m = 3,4
[0195] P?=< (xt)mb~ A)2 ^(ztmab~hf m-1,6 ; 、(U)2 + (九6 - A)2 /" = 2,5 (x〇me,b-Ay+{y〇mab- Ay 削=5,6
[0196] = < (xtnwh - A)2 + (ytmab - A)2 m = 2A , {Anab ~AT + {y\ah ~ A)' m =
[0197] 其中,
[0198] (5.2)计算由激光角反射器底面面形误差所导致的出射激光束附加光程co fm(m = 1、2、3. . .、6)，丨两足以下关系：
[0200] 式中，If为激光角反射器底面的面形光圈值，a为入射激光束在底面的折射角， er= tircsin (siiip/"J。
[0201] (5.3)计算由二面直角误差所导致的出射激光束附加光程，包括以下子过程：
[0202] (5.3.1)计算存在二面直角误差时，经过激光角反射器反射面三次反射后的反射 光束的附加矢量dF m，m=l、2、3…、6:
[0203] dFi=[2(A28i+A382),2(A383-Ai8i),-2(A283+Ai82)] T；
[0204] dF2 = [2(A352+A25i) ,_2)A353+Ai5i) , 2(A253~Ai52) ]T;
[0205] dF3=[2(A3S2-A2S1)，2(A 1S1+A3S3)，-2(A1S 2+A2S3)]T;
[0206] dF4=-dF2 ；dF5 = -dF3； dF6 = -dFi,
[0207]式中，和S3为激光角反射器三个反射面两两之间的二面直角误差。
[0208] (5.3.2)计算反射光束的附加矢量dFmg过激光角反射器底面折射后，出射激光束 的附加矢量dRm:
[0209] dRm= [dRmx，dRmy，dRmZ]T = refr(dFm，-No)，m= 1、2、3…、6〇
[0210] (5.3.3)计算由二面直角误差导致的出射激光束附加光程c〇am(m=l、2、3. ..、6)， 满足以下关系：
[0211] =-?s/^rtana(c//?,',< cos沒十，爪=1、2、3…、6。
[0212] (5.4)计算由入射激光偏振态所导致的出射激光束的附加复振幅分布，包括以下 子过程：
[0213] (5.4.1)计算入射激光束经过一次折射和三次反射后，在激光角反射器底面和反 射面上的P波和s波的矢量方向，满足以下关系式：
[0218] 式中，（pr，sr)和(pfm，sfm)分别为折射光矢量A在底面和第一次反射面上的p波和s 波的矢量方向，（pem，sem)和(ptm，stm)分别为反射光矢量Bm在第一次反射面和第二次反射面 上的P波和S波的矢量方向，（pUm，SUm)和(pVm，SVm)分别为反射光矢量Cm在第二次反射面和第 三次反射面上的P波和S波的矢量方向，化8^吨〇和化^，8^)分别为反射光矢量0?在第三次 反射面和底面上的P波和S波的矢量方向，下标m为正整数，m=l、2、3…、6,符号"表示向 量的叉乘运算，N4为激光出射时底面的法线方向，仏=_他，他、犯和犯为反射面的法线方向， 它们满足以下形式：
[0219] Ni=[-l，0,0]T，m=l、2;Ni = [0，-l，0]T，m=3、4;Ni=[0,0，-l]T，m = 5、6;
[0220] 他=卜1，0,0]T，m=3、5;N2 = [0，-l，0]T，m=l、6;N2 = [0,0，-l]T，m = 2、4;
[0221 ] N3=[-l，0,0]T，m=4、6;N3 = [0，-l，0]T，m=2、5;N3=[0,0，-l]T，m=l、3。
[0222] (5.4.2)计算折射光和反射光在激光角反射器表面的转换矩阵1^(七=1、2、3、4)， 满足以下关系式： rnooo1 _ \pr〇pfm srOpfl JpemOptm semQptm |_0223 J x y - $ - .， iprQsfn, wOC」 -sem〇stm _ 「…厂 r 「风御^ s^ \PSm〇PK sgmQphm
[0224] r, = ； r4 = ., ； ' _PSm〇shm
[0225] 式中，符号" 0 "表示矢量的点乘运算。
[0226] (5.4.3)计算折射光和反射光的p波和s波分量的反射系数矩阵和折射系数矩阵， 满足以下关系式：
[。227] 叫。、上
[0228] Rf= p ; Re= p ; Rt = ^ ; ? [0 <」 L0 Lq rts.-
[0229] 式中，Tf?和II分别为第一次折射和最后一次折射的折射系数矩阵，Rf、Re和Rt分别 为第一次反射、第二次反射和第三次反射的反射系数矩阵。折射系数矩阵和反射系数矩阵 的非零元素满足以下关系式：
[0230] tfp = ctp((p,a); (f\ = cls(0,0;.); //,, = dp(a,(p) \ // = cis(a,(p);
[0231 ] rfP = crp(if ,zf) ；rfs = crs(if ,zf) ；reP = crp( ie, ze) ；res = crs(ie ,ze)；
[0232] rtP = crp(it,zt) ；rts = crs(it,zt)；
[0233] 其中，（if ,zf)、（ie,ze)和（it,zt)分别表示激光束在第一次反射、第二次反射和 第三次反射时对应的入射角和折射角，if = arccos(A 〇 Ni)，zf = arcsin(usin if)，ie = arccos(Bm 〇 N2) ,ze = arcsin(usin ie), it = arccos(Cm 〇 N3),zt = arcsin(usin it)，式 中，u = nr/nf，nf为激光角反射器膜层的折射率。函数ctp和cts分别表示p波和s波的透射系 数的计算公式，函数crp和crs分别表示p波和s波的反射系数的计算公式，它们满足以下关 系式：
[0236] 其中，(仍，妒2)表示函数的输入量。
[0237] (5.4.4)计算由入射激光偏振态所导致的出射激光束p波分量和s波分量的附加复 振幅AP〇ut和AS〇ut的表达形式： ~ 1 「￡P. _
[0238] Jui =Tir4Rtr,Rer2Rfr{rf pn , _ . om _ -…in _
[O239]式中，EPin和ESin分别为入射光场p波分量和S波分量的复振幅。
[0240] (5.5)计算出射激光束p波分量和s波分量的附加复振幅分布TOdPSEm，m= 1、2、 3…、6,满足以下关系式：
[0241 ] PEm= | PEm| exp(i <}) Pm) = | AP〇ut | exp[ix( to rm+w fm+ o am)+i w Pm]；
[0242] SEm= | SEm| exp(i sm) = | AS〇ut | exp[iT( 0 rm+w fm+ ?am)+i w sm];
[0243] 式中，|PEm|和〇_分别表示PEm的模和相位，|SE m|和〇观分别表示SEm的模和相位，i 表示虚数单位。1 = 231/%，\为入射激光束的波长，〇^"1和《8111分别为4?。111；和43。1 11；的幅角。
[0244] 在本步骤中，出射激光束的附加复振幅分布综合考虑了二面直角误差、面形误差 和入射激光偏振态等多重因素的影响。
[0245] 6.计算接收屏位置处的激光角反射器远场衍射强度分布，包括下述过程：
[0246] (6.1)将接收屏细分为101X101个离散接收点，每个接收点在接收屏上的横坐标 xrcxi和纵坐标yrcxi分别满足：
[0247] xrcd = -Rw+2(c-l)Rw/101，c = l、2、3…、101;
[0248] yrcd = -Rw+2(d-l)Rw/yy，d=l、2、3…、101;
[0249] 其中，1为接收屏的半径大小。
[0250] (6.2)计算出射激光束的P波分量和S波分量的复振幅分布Pmab和Smab，满足以下关 系式：
[0251] Pmab = exp{iT[(dRmx+Al)xfab+(dRmy+A2)yfab]} 1 PEm ；
[0252] Smab = exp{ix[(dRmx+Al)xfab+(dRmy+A2)yfab]} ? SEm；
[0253 ] (6.3)计算出射激光束p波分量和s波分量的远场衍射强度值IPcd和I Scd，满足以下 形式：
[0256]式中，do为激光角反射器与接收屏之间的距离，i表示虚数单位，xgcd和ygcd可以表 示为：
[0257] .vg (/ - coS(/?cos^-.vrf/ -sin^- v/;,. +dvRt ;
[0258] .观喊=c〇s</?sin0 cos0 ? .vr^ + 。
[0259] (6.4)计算接收屏位置处的激光角反射器远场衍射强度分布Id，满足关系式：led 一 IPcd+I Sed 〇
[0260] 二、下面是本发明方法的一个具体实施例：
[0261] 1.输入激光角反射器远场衍射强度数值仿真的初始参数，包括输入激光角反射器 的口径和倒棱宽度、二面直角误差、反射面的面形误差和底面的面形误差、激光角反射器材 料的折射率和反射面膜层的折射率、入射激光束的入射角和方位角、入射光场的P波与s波 分量的复振幅分布等，其参数名称、符号与数值见表1所示；
[0262] 表1激光角反射器远场衍射强度数值仿真的初始参数
[0264]
[0265] 2.计算激光束在激光角反射器内部的光线方向，包括入射激光矢量R、折射光矢量 A、三个反射光矢量Bm、Cm、DjP出射激光矢量Em，m=l、2、3…、6;
[0266] 3.计算激光束与激光角反射器表面的交点坐标，包括入射激光束和出射激光束在 底面的父点坐标(xlab，y lab，zlab)和(XCmab，yCmab，ZCmab)，激光束与激光角反射器反射面的父 点丰不(XOmab，y〇mab，ZOmab )、( Xtmab，y tmab，Z tmab )矛口 ( xhmab，yhmab，zhmab )，中，HI - l、2、3'..、6，￡l = 1、2、3…、91，b = l、2、3…、91;
[0267] 4.基于反射面交点必须落入反射面区域和底面交点必须落入底面区域的限制条 件，计算激光角反射器的衍射区域范围(Ha，Hb)，其分布见图4所示；
[0268] 5.计算出射激光束p波分量和s波分量的附加复振幅分布PEjPSEm，m= 1、2、3…、6， PE4PSEm的模和相位分别见图5和图6所示；
[0269] 6.计算接收屏位置处的激光角反射器远场衍射强度分布led，即101X101个离散点 的强度值，其等高线分布规律见图7所示。
[0270]本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领 域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替 代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
【主权项】
1. 一种激光角反射器远场衍射强度的仿真方法，其特征在于，包括以下过程:首先根据 激光角反射器远场衍射强度数值仿真的初始参数计算激光束在激光角反射器内部的光线 方向；然后利用光线方向计算激光束与激光角反射器表面的交点坐标;其次根据各交点坐 标的限制条件计算激光角反射器的衍射区域范围；再次利用激光角反射器加工误差和入射 激光偏振态计算出射激光束的附加复振幅分布;最后根据衍射区域范围和出射激光束的附 加复振幅分布的结果，计算接收屏位置处的激光角反射器远场衍射强度分布。2. 根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法，其特征在于，所 述的激光角反射器远场衍射强度数值仿真的初始参数包括： 激光角反射器的几何参数，包括口径和倒棱宽度； 激光角反射器的加工误差，包括二面直角误差、反射面的面形误差和底面的面形误差； 激光角反射器材料的折射率和反射面膜层的折射率，若未反射面镀膜，则膜层折射率 置为1; 入射激光束的入射条件，包括入射激光束的入射角和方位角，入射光场P波与s波分量 的复振幅分布； 入射激光束的波长和用于激光角反射器面形检测的激光波长； 接收屏的大小及其与激光角反射器之间的距离； 激光角反射器底面和接收屏细分的点数。3. 根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法，其特征在于，所 述计算激光束在激光角反射器内部的光线方向，包括以下子步骤： 步骤3.1，计算入射激光束的入射矢量R和折射光矢量A: 及== [sin供cos设,.sinpsini^cosy :;.A= [Ai，A2，A3]T=refr(-MR，N〇); 式中，P和Θ分别为入射激光束在激光角反射器底面坐标系的入射角和方位角，No为激光角反射器底面的法线方向，refr函数为光线折射定律的矢量形式，Ri、R2和 R3分别为入射光矢量在激光角反射器底面坐标系下的分量，Al、A2和A3分别为折射光矢量在 激光角反射器反射面坐标系下的分量； 步骤3.2,计算激光束在激光角反射器内部经过6种不同反射次序后所对应的光线方向 矢量Bm、Cm、Dm和Em，m= 1、2、3···、6: Bm= γ iA，m=l、2;Bm= γ 2A，m=3、4;Bm= γ 3A，m=5、6; Cm- Υ 2Bm?Π 1 -l、6;Cm- Y 3Bm?Π 1 -2、4;Cm- Y lBm?Π 1 -3、5; Dm = -A;Em=R；4. 根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法，其特征在于，所 述计算激光束与激光角反射器表面的交点坐标，包括以下子步骤： 步骤4.1，将入射激光束细分为nx X ny根激光线，每根激光线在激光角反射器底面上的 横坐标Xfab和纵坐标yfab分别满足： xfab = -r+2(a-l)r/nx，a= 1、2、3·"、nx; yfab = -r+2(b-l)r/ny，b = 1、2、3···、ny; 其中，r为激光角反射器底面的半径； 步骤4.2，利用激光角反射器底面坐标系和反射面坐标系之间的转换关系，获得每根激 光线在反射面坐标系下的坐标(11￡^，71￡11)，21￡11 )): [xlab,ylab,zlab]T = M[xfab,yfab,〇]T+Mr,L为激光角反射器的棱长，通常情况下步骤4.3，根据入射激光束细化光线的起始坐标以及激光束在激光角反射器内部的矢 量方向，计算激光束与激光角反射器底面的交点坐标： [XOmab , y Omab , Z Omab ] -g ( Pk , Cm ) * [xlab,ylab,zlab]; [X tmab , y tmab , Z tmab ] -g(Pj,Bm) * [ XOmab , y Omab , Z Omab ]; [xhmab , yhmab , zhmab ] -§(β?,Α) * [ X tmab , y tmab , Z tmab ]；式中，（XOmab，y〇mab，ZOmab )、（ Xtmab，y tmab，Z tmab )和（xhmab，yhmab，zhmab )分力U表激7?束与 激光角反射器反射面的交点坐标，（1(^1)，5^11131)，2(31 1131))表不激光束与激光角反射器底面的交 点坐标、比、&和β4表示系数矩阵，下标i、j和k的数值与光线反射次序有关，它们取值为： 当m=l、2时，i = l;当m=3、4时，i = 2;当m=5、6时，i = 3; 当m=3、5时，j = l;当m=l、6时，j = 2;当m=2、4时，j = 3; 当m=4、6时，k = 1;当m=2、5时，k = 2;当m= 1、3时，k = 3; 系数矩阵 01=[1，〇,〇]，02=[〇，1，〇]，&3=[〇,〇，1]，04=[1，1，1]，函数8(0，了)满足关系 式：其中，β和J表示函数的输入量，β= [βχ,βγ,βζ]，βχ、βγ和分别为β的分量;J= [jx, Jy, Jz]，Jx、Jy和Jz分别为J的分量。 5 .根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法，其特征在于，所 述计算激光角反射器的衍射区域范围，包括以下子步骤： 步骤5.1，激光束与反射面的交点必须限制在激光角反射器反射面范围内，则符合该条 件的激光线序列号集合(Ra，Rb)必须满足以下形式： (Ra,Rb) = {(a,b) |ΤιηΤ2ηΤ3} 其中，符号" η "表示交集运算，集合??、Τ#ΡΤ3表达式为：式中，1表示激光角反射器的倒棱宽度； 步骤5.2，入射激光线和激光束与底面交点必须限制在激光角反射器底面范围内，则符 合该条件的激光线序列号集合(Fa，Fb)必须满足：步骤5.3,符合反射区域和底面区域限制要求的激光线所占的底面范围，即为激光角反 射器的衍射区域(Ha,Hb)，其表达形式为： (Ha，Hb) = (Ra，Rb) n (Fa，Fb)。6.根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法，其特征在于，所 述计算出射激光束的附加复振幅分布，包括以下子步骤： 步骤6.1，计算由激光角反射器反射面面形误差所导致的出射激光束附加光程ω rm(m = 1、 2、3. . .、6)，它满足以下关系：式中，= 表示第i个反射面的面形光圈值，nr为表示激光角反射器材料的折 射率，λ。为面形检测的激光波长，h (i = 1、2、3)可以表示为：步骤6.2，计算由激光角反射器底面面形误差所导致的出射激光束附加光程ω fm(m= 1、 2、 3. . .、6)，?两足以下关系：式中，U为激光角反射器底面的面形光圈值，α为入射激光束在底面的折射角， a = arcsin (sin</)//；(.)； 步骤6.3,计算由二面直角误差所导致的出射激光束附加光程，包括以下子过程： (6.3.1) 计算存在二面直角误差时，经过激光角反射器反射面三次反射后的反射光束 的附加矢量 dFm，m=l、2、3"_、6: dFi= [2(Α2δι+Α3δ2)，2(Α3δ3~Αιδ!)，_2(Α2δ3+Αιδ2) ]τ; dF2= [2(Α3δ2+Α2δι)，_2(Α3δ3+Αιδ!)，2(Α2δ3~Αιδ2) ]τ; dF3= [2(Α3δ2~Α2δι)，2(Αιδ1ι+Α3δ3)，-2(Αιδ2+Α2δ3) ]τ; dF4 = ~dF2；dFs = ~dF3 ； dF6 = ~dFi, 式中，为激光角反射器三个反射面两两之间的二面直角误差； (6.3.2) 计算反射光束的附加矢量dFmg过激光角反射器底面折射后，出射激光束的附 加矢量dRm: dRm= [dRmx，dRmy，dRmZ]T=refr(dFm，-N〇)，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)分别为反射光矢量仏在第一次反射面和第二次反射面上的 P波和s波的矢量方向，（pum，sum)和(pVm，SVm)分别为反射光矢量C m在第二次反射面和第三次 反射面上的P波和s波的矢量方向，（pgm，sgm)和(ph m，shm)分别为反射光矢量Dm在第三次反射 面和底面上的P波和s波的矢量方向，下标m为正整数，111=1、2、3-、6，符号"@"表示向量的 叉乘运算，N4为激光出射时底面的法线方向，N4 = -No，Λ、犯和犯为反射面的法线方向，它们 满足以下形式： Ni=[-l，0,0]T，m=l、2;Ni=[0，-l，0]T，m=3、4;Ni=[0,0，-l] T，m=5、6; N2 = [-1，0,0]T，m=3、5;N2 = [0，-l，0]T，m=l、6;N2 = [0,0，-l]T，m=2、4; N3 = [-1，0,0]T，m=4、6;N3 = [0，-l，0]T，m=2、5;N3 = [0,0，-l]T，m=l、3; (6.4.2) 计算折射光和反射光在激光角反射器表面的转换矩阵Γ t(t = 1、2、3、4)，满足 以下关系式：式中，符号" 0 "表示矢量的点乘运算； (6.4.3) 计算折射光和反射光的p波和s波分量的反射系数矩阵和折射系数矩阵，满足 以下关系式：式中，Tf和T1分别为第一次折射和最后一次折射的折射系数矩阵，Rf、Re和Rt分别为第 一次反射、第二次反射和第三次反射的反射系数矩阵;折射系数矩阵和反射系数矩阵的非 零元素满足以下关系式： tfp = ο?ρ(φ,α) ; lj\ = cls{(p"a); //" = ν?ρ(α.φ) ; tis =cts{a,(p); 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)和分别表示激光束在第一次反射、第二次反射和第三 次反射时对应的入射角和折射角，if = arccos(A θ Ni)，zf = arcsin(usinif)，ie = arccos (Bm Θ N2) ,ze = arcsin(usinie), it = -arccos(Cm 〇 N3), zt = arcsin(usinit) ,?ζψ ,υ = ηΓ/ nf，nf为激光角反射器膜层的折射率;函数ctp和cts分别表示p波和s波的透射系数的计算公 式，函数crp和crs分别表示p波和s波的反射系数的计算公式，它们满足以下关系式：其中，:〇】，识2)表示函数的输入量； (6.4.4) 计算由入射激光偏振态所导致的出射激光束p波分量和s波分量的附加复振幅 AP0Ut和AS〇ut的表达形式：式中，EPin和ESin分别为入射光场P波分量和S波分量的复振幅； 步骤6.5，计算出射激光束p波分量和s波分量的附加复振幅分布PEjPSEm，m= 1、2、3···、 6,满足以下关系式： PEm- PEm I ΘΧρ ( ? Φ pm) - ΑΡout | ΘΧρ [ ? Τ ( Q !Γιη+ Q fm't' ^ Bm) + ? CO Pm]; SEm- SEm I ΘΧρ ( ? Φ sm) - ASout | ΘΧρ [ ? T ( Q Q fm't' ^ Bm) + ? CO Sm]； 式中，I PEm I和Φ[?分别表示PEm的模和相位，I SEm I和〇sm分别表示SEm的模和相位，i表示 虚数单位;τ = 2π/λι·，入射激光束的波长，ω pjP ω sm分别为APmjt和ASmjt的幅角。7.根据权利要求1所述的一种激光角反射器远场衍射强度的仿真方法，其特征在于，所 述计算接收屏位置处的激光角反射器远场衍射强度分布，包括以下子步骤： 步骤7.1，将接收屏细分为μχΧμγ个离散接收点，每个接收点在接收屏上的横坐标x^d和 纵坐标yrcxi分别满足： xr cd - _Rw+2 (c_l) Rw/l4x，c -1、2、3···、]4x; yr cd - _Rw+2 (d_ 1) Rw/l^y，d -1、2、3···、1?; 其中，Rw为接收屏的半径大小； 步骤7.2,计算出射激光束的P波分量和S波分量的复振幅分布Pmab和Smab，满足以下关系 式： Pmab = exp { i τ [( dRmx+Al ) xfab+( dRmy+A2 ) yfab ]}· PEm; Smab = -exp{iT[(dRmx+Al)xfab+(dRmy+A2)yfab]} · SEm; 步骤7.3，计算出射激光束p波分量和s波分量的远场衍射强度值IPcd和I Scd，满足以下形 式：式中，do为激光角反射器与接收屏之间的距离，i表示虚数单位，xgcd和ygcd可以表示为：步骤7.4,计算接收屏位置处的激光角反射器远场衍射强度分布led，满足关系式：Icd = IPcd+I Scd 〇
【文档编号】G02B27/00GK105929531SQ201610545255
【公开日】2016年9月7日
【申请日】2016年7月12日
【发明人】周辉, 李松, 郑国兴, 田昕
【申请人】武汉大学