非线性方程法确定低空探测机载激光雷达消光系数边界值

提出了机载激光雷达进行低空探测时确定消光系数边界值的非线性方程法。首先根据机载激光雷达方程推导出消光系数边界值与机载激光雷达回波信号之间的函数关系,构建一个以消光系数边界值为变量的非线性方程。然后采用改进Jarratt法求此非线性方程的数值解,得到消光系数边界值。使用地基激光雷达真实回波信号数据模拟得到的机载激光雷达回波信号,进行了实验验证,并与斜坡法进行了比较。结果表明,本方法确定边界值的相对误差比斜坡法减小5.9%,由于其具有六阶的收敛特性,迭代次数仅为3次。该方法不需辅助测量设备,精度较高,收敛速度快,迭代次数少,具有良好的应用前景。     

一种新的大气消光系数边界值确定方法

本文以激光雷达测量机场跑道斜程能见度为应用背景,提出了一种基于不动点原理的大气消光系数边界值确定方法,将确定消光系数边界值的问题转化为求解函数不动点。首先在边界层底部选取参考点,假设一个消光系数边界值,在Klett反演算法的基础上,构建消光系数均值与边界值之间的函数关系。然后依据函数不动点存在性和唯一性的条件判断不动点...     

多线光谱的参数分组约束最小二乘拟合

本文在非线性最小二乘拟合中提出了参数分组约束优化方法,并以水汽的双线光声光谱分析为例显示了该方法的高效可靠性。     

改进的牛顿法确定大气消光系数边界值

在用激光雷达方程反演大气消光系数时,大气消光系数边界值对反演精度影响较大,而在低层大气中该值较难确定。文中提出了一种基于改进牛顿法的大气消光系数边界值确定方法,其核心思想是,把确定大气消光系数边界值的问题转化为求非线性方程的数值解。首先,根据大气消光系数边界值与激光雷达回波信号功率以及大气光学厚度之间的关系,假设大气消光系数边界值为x,构建一个非线性方程。其次,采用改进的牛顿法求非线性方程的数值解,得到大气消光系数边界值。使用香港天文台装置在香港国际机场的多普勒激光雷达回波信号数据,对该方法的可行性和可靠性进行了验证。结果表明:利用该方法确定边界值,可以较为准确地反演出低层大气消光系数。该方法收敛速度快,迭代次数少,并且不需要计算导数值,极大地减少了运算量,具有较强的实际应用价值。     

基于Fernald和Klett方法确定气溶胶消光系数边界值

基于Fernald方法和Klett方法,推导出一个具有明确物理意义的确定气溶胶消光系数边界值的表达式,该式比目前用于确定边界值的Collis斜率法表达式增加了两项：空气分子消光系数项和后向散射项,这两项与Collis斜率法的值符号相反。空气分子消光系数项较小,但后向散射项为后向散射系数的倒数和导数的乘积,绝对值能达到Collis斜率法的75.2 %。分析表明,考虑了新增两项反演的大气气溶胶光学厚度（AOD）与实测更接近,所以增加这两项是合理的、必要的。不同标高反演2007年9月20日的AOD在0.20~0.25之间,变化范围较小;反演的AOD方差为0.003,相对较小。说明新方法对标高的依赖较小且较稳定。分析424个时次晴空资料的反演结果可知,反演比实测大7.4%,反演与实测的相关系数为93.2%,相对误差和绝对误差分别为10.9%和0.03,反演的AOD方差为0.02,AOD小于0.45（占到资料总数的91.7%）时,反演结果较好。     

成像激光雷达大气消光系数的反演方法

针对成像激光雷达的应用背景,利用Klett迭代方法求解水平路程激光大气传播的消光系数,并采用Monte-Carlo方法对典型应用背景下的接收信号进行了模拟,结果表明：采用该方法求解的消光系数的理论误差小于6％;在考虑大气湍流起伏的情况下,消光系数的起伏变化明显,但是与理论模型的变化趋势吻合。     

成像激光雷达大气消光系数的反演方法

针对成像激光雷达的应用背景,利用Klett 迭代方法求解水平路程激光大气传播的消 光系数,并采用Monte2Carlo 方法对典型应用背景下的接收信号进行了模拟,结果表明:采用该方法求解的消光系数的理论误差小于6 %;在考虑大气湍流起伏的情况下,消光系数的起伏变化明显,但是与理论模型的变化趋势吻合。     

递推最小二乘法在射表拟合中的应用

本文用以（距离、俯仰角）为自变量的有理分式函数对射表拟合，采用递推最小二乘法求取拟合函数的参数。以得到满足精度要求的射表拟合函数。     

基于不动点原理的大气气溶胶消光系数边界值确定方法

提出了一种基于不动点原理的大气气溶胶消光系数边界值确定方法,其核心思想是将确定大气气溶胶消光系数边界值的问题转化为求解函数不动点。首先建立大气消光系数边界值与大气光学厚度和激光雷达回波信号之间的函数关系;其次依据函数不动点存在性和唯一性的条件估计不动点的存在,通过不动点迭代求得大气消光系数边界值,并由此值来确定大气气溶胶消光系数边界值。将本方法应用于实际激光雷达回波信号的反演中,得到低层大气气溶胶消光系数垂直廓线,并与在对流层顶使用洁净层法确定边界值所得的结果进行了对比。结果表明,利用本方法确定边界值,可以较为准确地反演出低层大气气溶胶消光系数。本方法可以预先估计不动点的存在区间、合理选取迭代初始值,具有收敛速度快、迭代次数少的优点,实际应用价值较强。     

船载激光雷达观测青岛海域气溶胶光学性质实验研究

研制了一套能连续探测海洋大气气溶胶光学性质垂直分布的船载小型米散射激光雷达系统,用该系统于2015年5月和9月在青岛海域开展了海洋大气气溶胶观测实验.以同船搭载的Microtops II型手持式太阳光度计同步测量得到的气溶胶光学厚度作对比数据,比较了两种仪器的气溶胶光学厚度测量结果,得到观测海域柱平均激光雷达比为32.4sr,标准偏差4.6sr.由Fernald方法得到了实验海域夏季15 km以下的气溶胶消光系数廓线.连续观测实验表明研制的船载米散射激光雷达结构紧凑,易于在船载平台操作,能对海洋大气气溶胶的光学特性进行连续有效的高时空分辨率廓线测量.     

船载激光雷达测量水体光学参数的仿真模拟研究

激光雷达能够高效获取海洋光学特性的垂直剖面信息,是海洋光学探测的重要手段之一。利用蒙特卡罗仿真方法,基于Gordon(1982)的机载激光雷达测量水体光学参数模型,研究了船载激光雷达在水中的传输过程和水中光场分布。特别研究考虑了接收视场角和望远镜半径等参数的影响,建立了适用于船载海洋激光雷达的模拟系统。在激光雷达的传输等效为太阳光传输后,该模拟系统与常用的HydroLight的模拟进行了比对印证并获得了一致的结果。在此基础上,通过模拟得到的激光雷达回波信号分析了不同激光雷达测量模式及典型水体条件下激光雷达消光系数α和海水光学参数之间的关系。船载激光雷达结果表明,在窄接收视场角情况下,激光雷达消光系数α趋向于水体光束衰减系数c;在宽接收视场角情况下,α趋近于水体的向下辐照度漫射衰减系数K_d。相比机载观测,船载观测的α趋近K_d的速度变缓。在垂直分层水体中,激光雷达在下层水体中测量的α值会向上层水体的α值偏移。该结果为研究海洋激光雷达测量参数与海洋光学参数之间的关系提供了进一步的认知。     

气溶胶后向散射消光对数比对消光系数反演的影响研究

使用Mie散射激光雷达在阴、雾和晴三种天气状况下采集了大气回波信号,从所测数据反演气溶胶消光系数的过程中发现,气溶胶后向散射消光对数比对反演结果存在较大影响.在阴天天气下,当后向散射消光对数比从0.7变化到1.0时,气溶胶消光系数反演结果相差近5倍;对晴天天气,反演结果相差近3倍;而在有雾天气,其反演结果变化不大.通过参考大气能见度因子对消光系数进行修正,并从而获得对应的后向散射消光对数比,给出了一种根据天气状况确定气溶胶后向散射消光对数比的新方法.实际反演计算表明,用此方法可获得与实际更为接近的气溶胶消光系数.     

云的多次散射对激光雷达测量结果影响的研究

为了更加准确地获得大气消光系数和后向散射系数,利用半解析Monte Carlo方法对云的大气多次散射激光雷达回波信号进行了模拟计算。分析了激光雷达接收视场角及光学厚度对多次散射回波信号的影响及水云和卷云多次散射因子与光学厚度的关系。激光雷达测量数据的对比结果表明,多次散射对卷云和水云消光系数影响较明显,而对后向散射系数的影响可以忽略。     

微脉冲激光雷达测量大气水平能见度

为了改进乃至革新大气能见度的探测技术,探索新型、快速、准确测量大气能见度的方法,采用自行研制的能见度微脉冲激光雷达,进行了水平大气能见度测量实验,并同已经商品化的能见度仪NQ-1进行了同步对比测量。实验结果表明,在水平能见度较差时,微脉冲激光雷达和能见度仪具有很高的相关性,相关系数达到了0.88;在水平能见度较好的时候,相关系数保持在0.73以上。说明能见度微脉冲激光雷达系统可以应用于大气能见度的探测。     

Computationally efficient DOA estimation for coprime linear array: a successive signal subspace fitting algorithm

ABSTRACT

In this paper, direction of arrival (DOA) estimation of multiple signals with coprime array is investigated and signal subspace fitting (SSF) method is linked to the coprime array, which achieves a better DOA estimation performance than the traditional uniform array. While the SSF method requires expensive computational cost in the case of multiple signals due to the multidimensional global angular searching, we propose a successive SSF (S-SSF) algorithm from a computationally efficient perspective. In the proposed algorithm, we employ rotational invariance and coprime property to obtain the initial estimates. Then, via a successive scheme, we transform the traditional multidimensional global angular searching problem into one-dimensional partial angular searching one. Consequently, the computational complexity has been significantly reduced. Specifically, the proposed S-SSF algorithm can obtain almost the same DOA estimation performance as SSF but with remarkably lower complexity. Finally, Cramer-Rao Bound (CRB) is provided and numerical simulations demonstrate the effectiveness of the proposed algorithm.