基于FRFT的低信噪比LFM信号参数快速估计算法.pdf

（）指挥控制与仿真 引用格式：东锦鹏，陈世文，杨锦程，等基于 的低信噪比 信号参数快速估计算法指挥控制与仿真，（）：，（）：基于 的低信噪比 信号参数快速估计算法东锦鹏，陈世文，杨锦程，韩 啸（信息工程大学，河南 郑州；部队，河南 济源）摘 要：基于分数阶傅里叶变换（，）对线性调频（，）信号参数进行估计，问题关键是确定 最佳阶数，根据误差迭代思想提出新的参数估计算法，该算法利用归一化带宽和旋转角的转化关系，由估计误差推算角度差值，有效降低了运算量，不需要调频斜率正负的先验信息，改进的对数搜索算法可以进一步提高参数估计结果的稳定性和可靠性。仿真结果表明，信噪比在 以上时该方法在高效率的前提下仍具有良好的参数估计性能，平均估计误差在 以内，估计结果接近 下限，满足工程实时处理需求。关键词：低信噪比；分数阶傅里叶变换；线性调频信号；参数估计中图分类号：文献标志码：，（，；，）：（）（），：（）；（）；（）；收稿日期：修回日期：作者简介：东锦鹏（），男，硕士研究生，研究方向为目标信息获取与处理。陈世文（），男，博导，教授。线性调频（，）信号作为基本低截获概率雷达信号，在工程中广泛应用，如何在低信噪比战场环境下对其实时处理，完成参数估计有着重要的研究前景和现实意义。于 年首次将分数阶傅里叶变换（，）应用于信号分析领域，不同于二次型时频分析工具，是一种线性变换，用分数域中的单一变量来表示信号的时频信息且没有交叉项的干扰，可以充分利用 信号的稀疏性特点，在低信噪比的条件下仍保持较高的检测概率和参数估计准确率，已有学者证明其估计结果为最小方差估计，推导了 下限（）公式。传统的 参数估计方法为二维搜索，运算速度较慢，不满足实时处理需求。为提高运算速度，大量学者对其离散算法优化进行了研究，同样在参数估计算法上，针对求取最佳变换阶数这一关键问题，文献根据 和时频分布关系，通过旋转角度等参数可计算得出最佳阶数，但信噪比较低时估计性能急剧下降。文献在此基础上引入了四阶原点矩作为目标优化函数，改进了传统的定步长搜索算法，文献将算法进行改进，应用于多分量信号，均取得了良好的估计效果。文献分别以大津法、最小范数、信息熵为目标优化函数，提高了估计性能但在噪声能量较大时三者表现欠佳，文献基于 原理，提出了改进型算法 及，同样存在运算量较大的缺点。以上快速算法大都需要在运算前判断未知信号调频斜率的正负，具有一定的局限性，信号实时处理能力和参数估计结果的稳 东锦鹏，等：基于 的低信噪比 信号参数快速估计算法第 卷定性仍需进一步改进。针对低信噪比环境下 信号参数估计结果准确度不高、现有算法运算量大的问题，本文提出快速算法，按照粗精两步估计思路。本算法包括根据误差迭代思想的快速估计和对数搜索的精确估计两部分，快速估计部分通过 归一化带宽、旋转角差值、旋转角估计值三者关系，根据误差自适应调整迭代步长，快速求得旋转角估计值；精确计算部分提出改良的对数搜索算法，通过 次运算，可将搜索步长减小至 ，进一步提高了参数估计精度和稳定性。信号分数阶频谱特性分析 信号分数阶频谱特性 信号一般表示为（）（）（）（）其中，为 信号的初始频率，为调频斜率，（）为高斯白噪声。文献分析了 信号时频分布与 的联系，如图 所示，运算可以理解为将信号时频面坐标轴绕原点逆时针转动某一特定角度，频谱即为转动后形成的分数域上的时频投影。图 投影示意图 在实际工程应用中需要进行归一化步骤，详细理论过程可参考文献，在此给出归一化量纲归一化因子 的表达式为 （）其中，分别为 信号的时宽和采样频率。归一化后的横（纵）坐标总长度为（）当旋转角度 接近最优旋转角时，归一化频率的带宽 趋近于，即 变换近似于冲激信号。无论是传统算法还是改进的高效算法，关键步骤均是通过搜索（估计）取得最优旋转角，之后的最优变换阶数 由下式求得（）根据变换阶数 以及最大峰值的坐标，从归一化变换关系估计得出信号的调频斜率和载频。功率谱平滑在估计 信号归一化频率带宽时，噪声对估计精度影响较大，直接影响参数估计结果的可靠性。为降低噪声对估计结果的影响，对信号频谱进行光滑处理，采用以下平滑滤波公式可以充分利用原功率谱特征信息（）（，）()（）其中，为信号长度，为平滑窗长度，（）为信号原频谱，（）为平滑功率谱，为取余运算。根据文献 的论证说明，频谱平滑窗长度取附近的整数值，归一化频谱则选取 附近的整数值。四阶原点矩本文选择分数阶频谱的四阶原点矩作为搜索目标优化函数。信号（）的分数阶频谱四阶原点矩定义为（）（）（）最佳旋转角度 下的四阶原点矩为（）（）（）（）其中，为信号幅值，为信号调制周期。若旋转角度与最佳值有差值，其四阶原点矩为（）（）（）其中（）为常数，显然，越小，四阶原点矩越大，即最佳旋转角度 对应阶数的 频谱取得最大四阶原点矩。四阶原点矩作为搜索目标优化函数相比于频谱幅度特征，在接近最佳旋转角度 时，特征变化更加平滑明显，能有效提升低信噪比环境下的最值判断准确度。基于 的 信号参数估计为提高 信号参数估计的速度和稳定性，本算法流程分为粗细估计两部分。基于误差迭代的快速算法基本原理：通过 与时频分布图像旋转对应关系，图 中的归一化频率的带宽 可表示为第 期指挥控制与仿真 （）（）其中，为归一化后的横（纵）坐标总长度，为最优阶数对应的旋转角，为当前旋转角。式（）给出了当前旋转角与最优旋转角的差值即估计误差 与当前 频谱带宽 之间的关系，根据 变换的旋转可加性，可以通过加（减）差值来快速计算最优旋转角。图 归一化带宽与旋转角差值的关系 特别的，当 时，实质上为傅里叶变换，对应带宽，归一化调频斜率为，角度差值为 （）（）为避免差值选取过大，消除功率谱平滑导致估计带宽偏大的不利影响，提高迭代精度，取 （）其中，为接近于 的常数系数。调频斜率的正负决定了最佳旋转角及 阶数的区间，调频斜率为负值则旋转角为，阶数为；为正值则旋转角在 ，阶数为。当前有的快速算法将信号分成两段，做 次，对比中心频率来确定调频斜率的正负，增加了运算量，计算过程较为繁琐。本文对此流程进行了改进，在求得当前旋转角与最优旋转角的差值后，分别将角度 和 代入式（）对应阶数的 变换，考虑到更接近最佳旋转角度的 变换四阶原点矩更大，取两者中较大四阶原点矩的角度为旋转角的较精确计算值 ，即 （），（）（）此时不需要调频斜率正负的先验信息或进行 运算，可自动调整取值以更接近最佳旋转角度，具体过程见图。利用上述 运算作用于 信号的基本机理，通过功率谱平滑计算归一化带宽，得出旋转角差值，更新旋转角计算值，再计算旋转角估计值下的 频谱归一化带宽，进而求得更精确的旋转角差值，根据 变换的旋转可加性，逐次迭代，直到旋转角估计值的四阶原点矩达到最大，在对计算速度要求较高的情况下此部分可作为快速估计算法单独使用。算法中迭代步长 为自适应取得，递减速度极快，轮即可取得旋转角最优值。图 逐步接近最优旋转角的 频谱 快速估计步骤及流程图以下为快速估计算法步骤，流程如图 所示。图 快速估计流程图 步骤 通过 得到信号频谱。步骤 确定信号带宽的粗计算值。根据式（）做功率谱平滑，降低噪声影响后求得归一化频率下的信号带宽的粗估计值。步骤 求得初始旋转差值。根据式（）可得 （）其中，为归一化后的横（纵）坐标总长度。代入式（）得（）东锦鹏，等：基于 的低信噪比 信号参数快速估计算法第 卷步骤 确定旋转角的较精确估计值。由于调频斜率的正负未知，有 （）比较四阶原点矩大小，获得旋转角的较精确估计值为 （）（）（）为旋转角 下的频谱四阶原点矩。步骤 更新归一化频率下的信号带宽。上步求得的角度下的 频谱做功率谱平滑，求得归一化频率下的信号带宽的粗估计值。步骤重复步骤 。由于逐步接近最佳旋转角度，归一化频率带宽 趋近于零，旋转差值 随之减小，直到旋转角估计值 下的四阶原点矩达到最大，此时对应最优旋转角。最优旋转角由式（）迭代求得。（）（）（）改进的对数搜索算法根据快速估计取得的旋转角估计值，在选定区间内以使用时间复杂度较低的对数搜索算法快速搜索，进一步提高参数准确度和稳定性，适合在低信噪比或对误差要求更严格的情况下使用。改进的对数搜索算法原理为：取旋转角度区间，初始步长 取较大数值以确保在初始区间内存在最优解。计算区间两端角度对应的 变换及其四阶原点矩，若两端的四阶原点矩小于中心角度的四阶原点矩，则保持搜索中心角度不变，新步长缩减至原步长的一半；若左侧四阶原点矩较大则将左侧端点设为搜索中心角度，右侧四阶原点矩较大则将右侧端点设为搜索中心角度，新步长缩减至原步长的一半。（），（），（）（）为第 次迭代后更精确的旋转角度。每次搜索计算 次，步长缩减为原来的一半，经过 轮搜索，搜索步长可缩小 倍，即为初始步长的 。精确计算步骤及流程为进一步接近最佳旋转角度，算法的精确估计步骤如下，基于改良对数搜索算法的精确估计流程如图 所示。步骤 确定初始步长。步骤 区间内搜索。利用快速估计得到的较精确的旋转角度，取区间，分别计算区间端点的四阶原点矩，比较大小后更新旋转角估图 精确估计流程图 计值。步骤 更新搜索步长。判断此时步长是否小于目标步长，不满足则更新步长，反之则结束搜索。步骤 估计 信号参数。取得最优旋转角度估计值 和最佳变换阶数 后，设 变换结果为（），信号的调频斜率、带宽及载频由式（）求得。（）（）（）其中，为归一化量纲归一化因子。算法运算量及稳定性分析本文中 变换采用文献提出的数值计算算法，设信号的采样点数为，则进行次变换需要的计算量为（），进行 次四阶原点矩计算需要的运算量为（）。第一部分作为快速算法，最多进行 轮迭代，每轮迭代进行 次 变换和 次四阶原点矩计算，运算量为（）。第二部分通过对数搜索，进一步提高参数准确度。初始步长 取 ，目标步长取 ，根据前文分析，缩小至需进行 轮搜索，共进行了 次 变换和四阶原点矩计算，运算量为（），粗精两步估计总运算量为（）。传统二维搜索算法同样达到 步长精度，需进行 次 运算，运算量为（），可见，本文算法在低信噪比环境下有较高参数估计准确率的同第 期指挥控制与仿真 时，运算量得到显著降低。仿真实验及分析实验 低信噪比环境的有效性验证。设定 为，选取的 信号为（）（），其中载频 ，调频斜率 ，采样频率 ，参数 取 。实验中，经功率谱平滑滤波，改进算法的快速估计部分的归一化长度变化如图 所示，在迭代过程中，信号能量更加集中，经过 轮迭代计算，最佳阶数取得估计值 ，调频斜率估计结果相对误差 ，载频估计相对误差 ，估计误差较低，说明此信噪比下，本文算法的快速估计部分能够对 信号参数进行有效估计。图 迭代后的归一化频谱 图 四种算法的仿真实验对比 算法的第二部分精确计算环节可以进一步提高参数估计结果的准确性、可靠性，精确估计部分对数搜索算法的初始步长取 ，目标步长取 ，计算结果见图。经进一步计算，将 最佳阶数的估计值由 进一步提升至 ，载频误差和调频斜率误差进一步减小至 及 。图 估计结果 实验 算法性能评估。设定 范围为 ，每个信噪比进行 次 实验，以文献改进高效 算法与文献黄金分割搜索算 法 作 为 对 比，选 取 的 信 号 为（）（），其中载频 ，调频斜率 ，采样频率 。本文算法参数 取 ，搜索初始步长 取 。图 为四种算法的载频和调频斜率估计误差均值以及均方误差 的对比图。实验结果表明，与改进高效 算法及黄金分割算法相比，本文算法的 东锦鹏，等：基于 的低信噪比 信号参数快速估计算法第 卷快速估计部分，既有速度优势也有精度优势，且在增加精确搜索部分后，在 低信噪比环境下仍能较为准确地估计出 信号参数，算法提升效果更加明显。从仿真结果的误差方差分布与 结合来看，在低信噪比下本文算法的估计结果仍然较为稳定，更贴近 信号的，算法鲁棒性好，结果可靠性高。实验 运算量评估。为验证低信噪比下本文算法的高效性，将本文算法与文献的改进高效 算法及传统算法进行分析比较，黄金分割搜索算法在低信噪比环境下已失效，不再列入对比。设定信噪比为，选取的 信号及算法参数不变，目标搜索步长取 、，进行 次 实验，取调频斜率相对误差均值 和载频相对误差均值，实验结果见图，运算量统计见表。从参数估计精度角度分析，信噪比为 时，本文算法快速部分基于迭代仍可进行准确计算，包含精确计算步骤后，也取得了预期的更准确的结果，进一步提高了参数值的准确性和稳定性，而改进高效 算法在快速估计归一化带宽时误差较大，难以搜索取得较准确的最佳变换阶数。从算法运算量角度分析，相比改进高效 算法，本文所提算法的快速估计部分以及对数搜索精确估计部分不需要进行 运算来分析调频斜率正负，较为高效便捷，另外，在阶次误差要求更严苛，目标搜索步长更小时，本文的改进算法在精确度、运算量和可靠性等方面具有更突出的优势。图 不同目标步长下的估计结果 表 不同目标步长下的运算量对比 序号目标搜索步长 本文快速估计部分 次 次 传统搜索算法 次 次 次 改进高效 次 次 次 次 次 次 本文快速估计对数搜索 次 次 次 次 次 次 结束语本文提出了一种基于 的 信号参数估计算法，运用粗精两步计算思想，快速计算部分基于误差迭代，利用了 与时频变换的旋转关系，通过迭代快速求得信号参数，不需要调频斜率正负的先验信息，稳定性高准确度好，适合于对参数估计时效性要求高的场合；精确计算部分利用快速估计结果，结合改进的对数搜索算法，仅通过 次运算，可将搜索步长减小至 ，进一步提高参数估计精度。实验证明，本文算法在信噪比 环境下仍能保持良好的参数估计性能，算法能够满足对效率和精度不同要求的场景，相比 经典二维搜索算法及其他快速算法，准确性和可靠性优势明显，可用于信号实时处理分析的工程实践。参考文献：，（）：，（）：，第 期指挥控制与仿真 ，（）：，（）：，（）：，（）：平先军，陶然，周思永，等 一种新的分数阶傅立叶变换快速算法 电子学报，（）：，（）：，（）：黄响，唐世阳，张林让，等 一种基于高效 的 信号检测与参数估计快速算法 电子与信息学报，（）：，（）：，（）：刘利民，李豪欣，李琦，等 基于分数阶傅里叶变换的低信噪比线性调频信号参数快速估计算法 电子与信息学报，（）：，（）：，（）：，（）：，（）：，（）：，（）：赵兴浩，邓兵，陶然 分数阶傅里叶变换数值计算中的量纲归一化 北京理工大学学报，（）：，（）：齐林，陶然，周思永，等 基于分数阶 变换的多分量 信号的检测和参数估计 中国科学 辑：技术科学，（）：，（）：（责任编辑：胡前进）