HTML
-
通常理想目标表面材料的光散射特性复杂,理论上由粗糙面散射理论可获知实验数值,但实际上很难提供计算所需的粗糙面元微观特性如相关长度、统计分布等。因此,通常使用双向反射分布函数BRDF来描述目标面元的光散射特征。BRDF所描述的光散射特征分布和粗糙面元的微观特性、材料介电常数、入射光波长等因素紧密联系。通常粗糙面元会被看作是由大量微面元构成,且所有微面元的反射都遵循菲涅耳效应[10-11]。
如图 2所示,对于小面元dS,入射光源方向为(θi, φi),探测器的观测方向为(θr, φr)。其中θ是天顶角,φ是方位角,z是代表粗糙面元的法线。BRDF定义为入射到粗糙表面的光辐射照度与光辐射亮度相联系,在图中表示为dLr(θi,φi,θr,φr)与沿(θi,φi)方向入射到被照射面元的亮度dEi(θi, φi)之比,即:
式中, fr(θi, φi, θr, φr)可以简化为fr(θi, θr, φ),φ=φr-φi为小面元坐标系内入射角与反射角的差值。
-
根据粗糙面元的BRDF样片材料,可以获知样片表面的光散射特征。但由于直接获取样片BRDF很困难,实验次数有限,不能多次重复实验获取数据。一般情况下是采用所得到的部分BRDF实验值拟合五参量统计模型,五参量模型基于粗糙面元的目标光散射计算获得,最终给出了BRDF的五参量统计模型[12]:
式中, kb,kd,kr,a,b为待定参量。kb为镜反射系数,kd为漫反射系数,kr为与微面元表面法线分布函数相关的参量,a和b是与表面菲涅耳函数相关的系数。根据能量守恒定律,kd+kb≤1,θi, θr, φr分别为入射角、散射角和相对方位角,γ为微观平面上本地坐标系的入射角。(2)式前半部分表示的是双向反射分布函数的相干散射,其中粗糙表面上偏振二向分布函数用kr2×cosα/[1+(kr2-1)cosα]标记,α为微面元平面法线与粗糙面元平面法线的夹角[13-14],菲涅耳反射系数用exp[b(1-cosγ)]来表示,遮蔽系数用G(θi, θr, φr)来表示。
遗传算法的收敛速度极快,它是一种全局优化概率算法,可以很快找到最佳解。而且各态历经性使其能够有效进行全局优化搜索,减少计算时间,保证算法有效性基础上得到参量值。因此可以通过遗传算法来计算五参量模型的参量值[15]。选取样片1和2的实测BRDF数据来进行优化建模,波长和五参量如表 1所示。
wavelength/μm kb kd kr a b error/% 0.905 1.436 0.109 26.36 0.411 -12.85 3.07 0.86 3.201 0.047 0.681 0.595 -98.04 4.71 Table 1. Five parameters data table of sample piece
选取不同入射波长的BRDF实测数据来进行优化建模并与实测数据相比较[16],入射波长分别为0.86μm和0.905μm。如图 3所示,实线为实测数据,虚线为五参量模型计算结果。可知对于不同材料的BRDF,目标回波的形状基本相似,幅值不同。而且随着入射角的增大,峰值功率逐渐增大。由于脉冲激光引信的探索过程是高速变化的,其入射角的波动范围值较广。因此一般检验时应选择在较大入射角的情况下进行[17-18],以便回波信号可以在能量较少时也能够成功实现检测。模拟曲线能很好地符合实验数据,证明了建模结果的正确性。
3.1. 双向反射分布函数
3.2. 五参量模型
-
针对引信工作中的实际环境,如在云烟环境中,结合脉冲激光发射系统信号处理单元对目标回波的判别,提出了基于回波包络上升速率的区域联合判别方式,并通过实际的实验测试得出数据。测试实验中用表面反射率接近实体目标的白板代替,小型目标表面反射率约为0.33,大中型目标表面反射率约为0.9,具体以不同尺寸白板来区分各型号的探测目标,云烟环境用烟雾机产生的烟雾模拟。由于研究的引信探测视场呈环形,则每一路形成的探测视场角为90°。因此,设置3个测试角度,分别为中心点0°、左侧45°和右侧45°。
如图 5所示,通过背景板可以看到,激光光束照射的视场角以对称形式左右可分别到达45°,测试得到的光束光斑分布均匀。系统设计的小型目标作用距离为3m,中型目标作用距离为5m,大型目标作用距离为7m。在没有干扰下, 不同场景目标测试记录数据如表 2和表 3所示。
targettype range number of measurement echo signal amplitude/mV left 45° center 0° right 45° small object 3m No.1 240 460 250 No.2 260 400 270 No.3 220 280 310 medium target 5m No.1 210 300 130 No.2 150 300 110 No.3 460 600 280 large target 7m No.1 130 200 160 No.2 110 110 110 No.3 280 350 250 Table 2. Objective test data record sheet under ideal environment
target type range number of measurement echo signal amplitude/mV left 45° center 0° right 45° small object 3m No.1 250 360 210 No.2 260 280 250 No.3 210 220 210 medium target 5m No.1 310 260 230 No.2 150 210 180 No.3 160 230 240 large target 7m No.1 200 190 180 No.2 130 130 130 No.3 130 140 130 Table 3. Objective test data record sheet under cloud and smoke environment
将模拟目标白板在不同作用距离点上下移动,找到反射点,通过示波器显示反射回波,读取回波幅值。利用数据记录仪记录回波包络图如图 6a所示。
从记录结果可以看到, 在中心位置即发射激光直射到目标时, 发射回波的幅值最大,位于探测视场两侧的光束能量略有下降。由于作用距离变大,反射回波幅值会降低,这符合激光传输能量的变化趋势。图中矩形框内包络的上升速率呈递增趋势,到达峰值后稳定持续一段时间,包络走势满足对目标识别的判定准则。
使用烟雾机在装置探测视场内间隔喷出烟雾来模拟实际环境,测试结果如表 3所示。
在实际环境中,系统正常工作,比较表 2数据,回波幅值大小有所降低,偏差约为65.3%。这是因为激光脉冲抵抗干扰而损失能量,但在判定准则的应用下依然实现对目标的判定。利用数据记录仪记录此测试环境下的回波包络如图 6b所示。
与图 6a相比,系统探测到的回波包络图出现明显的差异,在探测光路中脉冲激光由于实际环境如云烟的影响,激光反射率较低,因此, 系统接收到此部分反射回波后作出的反应会出现小幅的的包络响应。当激光穿过云烟到达目标表面时,反射回波到达接收系统后会出现幅度较大的反应目标的包络图形。从包络图中可以明显看出, 目标形成的包络在前沿上升速率较快,且达到峰值后稳定持续一段时间。对比图中方框1和方框2可知,干扰形成的包络在达到某一峰值后会迅速下降,不会形成持续稳定得包络现象,因此, 系统在探测过程中依据判定准则较好地实现了抗干扰。