1.一种空中目标光子计数激光雷达探测仿真的方法,其特征在于,所述方法包括如下步骤:步骤S1、计算空中目标所处的大气透过率,以确定所述空中目标的背景噪声;步骤S2、对所述空中目标进行小面元拆分,得到小面元集;步骤S3、对所述小面元集中每个小面元进行姿态变换;步骤S4、根据姿态变换结果,获取所述空中目标不同时刻的激光脉冲波形;步骤S5、根据所述空中目标的背景噪声以及所述空中目标不同时刻的激光脉冲波形,计算不同时刻的信号光子数和噪声光子数;步骤S6、根据所述不同时刻的信号光子数和噪声光子数,对所述光子计数激光雷达中的单光子探测器的回波光子接收过程进行仿真,得到单光子探测器阵列中各个单元的光子回波概率分布;在所述步骤S6中,所述仿真的具体过程包括:步骤S61、设置产生光子事件的时间t=1;步骤S62、根据t时刻的信号光子数和噪声光子数,产生一个t时刻的泊松分布的随机数;步骤S63、根据所述随机数,判断光子事件是否产生,如是,则将所述光子事件对应的时间进行存储,并令t=t+t dead ,进入步骤S64;如否,则令t=t+1,返回步骤S62;步骤S64、判断t是否大于脉冲周期,如是,则进入步骤S65;如否,则令t=t+1,返回步骤S62;其中,t dead 为所述单光子探测器的死时间;步骤S65、判断t是否大于总探测时间,如是,则输出存储的所述光子事件对应的时间,以获取所述单光子探测器阵列各个单元的光子回波概率分布,结束;如否,则t=t+1,返回步骤S62。
2.根据权利要求1所述的方法,其特征在于,在步骤S1中,所述大气透过率为: ; ;其中, T 为大气透过率;λ为激光波长; λ 0 为常数; V 为大气能见度,单位为 km ;X为激光传输距离,单位为 km 。
3.根据权利要求2所述的方法,其特征在于,在步骤S3中,按照如下公式进行姿态变换:S ( x , y , z )= RS 0 ( x 0 , y 0 , z 0 )+ M ( x m , y m , z m ); ;其中, S ( x , y , z )为姿态变换后的各个小面元坐标; S 0 ( x 0 , y 0 , z 0 )为姿态变换前的各个小面元坐标; M ( x m , y m , z m )为空中目标的中心坐标; R 为姿态变换矩阵; θ pitch 、 θ roll 和 θ yaw 分别为空中目标的俯仰角、翻滚角和偏航角。
4.根据权利要求3所述的方法,其特征在于,在所述步骤S4中,所述获取所述空中目标不同时刻的激光脉冲波形的具体过程包括:步骤S41、根据姿态变换结果,计算所述各个小面元到所述光子计数激光雷达中的激光器发射的激光中心轴的距离,以选取距离小于激光光斑半径的小面元并放入到新的小面元集中;步骤S42、对所述激光器发射的激光进行分束;并从分束后得到的子激光中选取出与所述新的小面元集中小面元相交的子激光作为新子激光;步骤S43、获取所述新的小面元集各个新的小面元与各个新子激光的交点,构成交点集;步骤S44、从所述交点集中抽取出与同一个新子激光的交点,并选择距离激光器最近的交点所对应的小面元在不同时刻的波形作为对应新子激光不同时刻的脉冲波形;步骤S45、对所有的新子激光不同时刻的脉冲波形进行求和,得到所述空中目标不同时刻的激光脉冲波形。
5.根据权利要求4所述的方法,其特征在于,所述方法还包括:步骤S7、根据所述单光子探测器阵列各个单元的光子回波概率分布,计算所述单光子探测器阵列各个单元获得的目标距离值,以获得所述空中目标的三维光子点云分布。
6.一种空中目标光子计数激光雷达探测仿真的系统,其特征在于,所述系统包括:第一计算模块,用于计算空中目标所处的大气透过率,以确定所述空中目标的背景噪声;拆分模块,用于对所述空中目标进行小面元拆分,得到小面元集;姿态变换模块,用于对所述小面元集中每个小面元进行姿态变换;获取模块,用于根据姿态变换结果,获取所述空中目标不同时刻的激光脉冲波形;第二计算模块,用于根据所述空中目标的背景噪声以及所述空中目标不同时刻的激光脉冲波形,计算不同时刻的信号光子数和噪声光子数;仿真模块,用于根据所述不同时刻的信号光子数和噪声光子数,对所述光子计数激光雷达中的单光子探测器的回波光子接收过程进行仿真,得到所述单光子探测器阵列各个单元的光子回波概率分布;所述仿真模块包括:设置子模块,用于设置产生光子事件的时间t=1;产生子模块,用于根据t时刻的信号光子数和噪声光子数,产生一个t时刻的泊松分布的随机数;第一判断子模块,用于根据所述随机数,判断光子事件是否产生,如是,则将所述光子事件对应的时间进行存储,并令t=t+t dead 后传输给第二判断子模块;如否,则令t=t+1后传输给产生子模块;第二判断子模块,用于判断t是否大于脉冲周期,如是,则将t传输给第三判断子模块;如否,则令t=t+1后传输给产生子模块;其中,t dead 为所述单光子探测器的死时间;第三判断子模块,用于判断t是否大于总探测时间,如是,则输出存储的所述光子事件对应的时间,以获取所述单光子探测器阵列各个单元的光子回波概率分布,结束;如否,则t=t+1后传输给产生子模块。
7.根据权利要求6所述的系统,其特征在于,所述获取模块包括:计算子模块,用于根据姿态变换结果,计算所述各个小面元到所述光子计数激光雷达中的激光器发射的激光中心轴的距离,以选取距离小于激光光斑半径的小面元并放入到新的小面元集中;选取子模块,用于对所述激光器发射的激光进行分束;并从分束后得到的子激光中选取出与所述新的小面元集中小面元相交的子激光作为新子激光;构成子模块,用于获取所述新的小面元集各个新的小面元与各个新子激光的交点,构成交点集;抽取子模块,用于从所述交点集中抽取出与同一个新子激光的交点,并选择距离激光器最近的交点所对应的小面元在不同时刻的波形作为对应新子激光不同时刻的脉冲波形;求和子模块,用于对所有的新子激光不同时刻的脉冲波形进行求和,得到所述空中目标不同时刻的激光脉冲波形。
8.根据权利要求7所述的系统,其特征在于,所述系统还包括:第三计算模块,用于根据所述单光子探测器阵列各个单元的光子回波概率分布,计算所述单光子探测器阵列各个单元获得的目标距离值,以获得所述空中目标的三维光子点云分布。