技术领域本发明属于循环平稳信号处理技术领域，涉及到时延和多普勒频移联合估计的方法，特别涉及一种使用广义Sigmoid变换循环模糊函数的时延与多普勒频移联合估计的方法。
背景技术：
大多数通信、雷达信号都具有循环平稳特性，利用这一特性可以有效去除同频带干扰信号的影响。循环模糊函数就是应用了循环平稳理论的一个经典方法，用于同频干扰下的时延和多普勒频移联合估计具有很好的效果。但是该方法是基于噪声服从高斯分布的假设下的，当噪声中包含脉冲特性时算法性能将严重退化或失效。基于分数低阶统计量理论的分数低阶循环模糊函数，可以同时在脉冲噪声和同频干扰下工作，但也有其固有的缺点：分数低阶的阶数需要噪声的先验知识，否则阶数的不当选择会影响方法的估计效果，另外，分数低阶统计量在脉冲性较强时，对脉冲噪声抑制能力不足，从而导致估计误差增大。故本发明提出一种可调Sigmoid变换循环相关，在可调Sigmoid循环相关的基础上，应用广义Sigmoid变换循环模糊函数在脉冲噪声和同频干扰共存条件下进行时延和多普勒频移的联合估计，并且可以通过调节可调Sigmoid变换循环相关的倾斜参数，调整对不同信号的衰减程度，以提高方法的适用范围。
技术实现要素：
针对现有技术的不足，本发明提供一种利用广义Sigmoid变换循环模糊函数进行时延和多普勒频移的联合估计方法。该方法采用可调Sigmoid变换可以通过倾斜参数的调整有针对性地对噪声中的脉冲性进行抑制，利用同载频不同波特率的信号具有不同循环频率的特性，可以通过循环频率有效区分待测信号和同频干扰信号。故本方法可以在脉冲噪声和同频干扰同时存在时对时延和多普勒频移进行联合估计。本发明的技术方案为：一种基于广义Sigmoid变换循环模糊函数的时延与多普勒频移联合估计的方法，主要包括以下步骤：第一步，计算广义Sigmoid变换循环模糊函数1.1)利用可调Sigmoid循环自相关公式计算参考信号的可调Sigmoid循环自相关；1.2)利用可调Sigmoid循环自相关公式计算待测信号的可调Sigmoid循环自相关；1.3)利用可调Sigmoid循环互相关公式计算参考信号与待测信号之间的可调Sigmoid循环互相关。1.4)利用可调Sigmoid循环互相关公式计算待测信号与参考信号之间的可调Sigmoid循环互相关。1.5)根据广义Sigmoid变换循环模糊函数公式计算获得广义Sigmoid变换循环模糊函数。第二步，参数估计2.1计算1.5)结果中广义Sigmoid变换循环模糊函数的绝对值；2.2搜索2.1)结果所对应的时延值和多普勒频移值，即为所求的估计值。本方法能够在脉冲噪声和同频干扰同时存在的条件下对时延和多普勒频移正确估计，且抗脉冲噪声能力较强，满足实际需求。附图说明图1是本发明的算法流程图；图2是本发明中广义Sigmoid变换循环模糊函数以BPSK信号为例的三维图；图3(a)是本发明中广义Sigmoid变换循环模糊函数在频率等于真实频移时的时间轴截面；图3(b)是本发明中广义Sigmoid变换循环模糊函数在时延等于真实时延时的频率轴截面。具体实施方式为使本发明实施例的目的、技术方案及其优点更加清楚，下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚完整的描述，整体算法流程图如图1所示：第一步，计算广义Sigmoid变换循环模糊函数；1.1利用可调Sigmoid变换循环自相关公式计算参考信号x(t)在时间区间[-T/2,T/2]上的可调Sigmoid变换循环自相关如公式(1)所示：Rx,TSϵ(τ)=ΔlimT→∞1T∫-T/2T/2Sigmoid[x(t+τ/2)]Sigmoid*[x(t-τ/2)]e-j2πϵtdt---(1)]]>采集数据为有限长度时的可调Sigmoid变换循环自相关可由公式(2)估计得到；R^x,TSϵ(τ)=1T∫-T/2T/2Sigmoid[x(t+τ/2)]Sigmoid*[x(t-τ/2)]e-j2πϵtdt---(2)]]>其中，ε表示循环频率，由s(t)的循环频率确定；τ表示时延；x(t)为接收到的参考信号，x(t)＝s(t)+w1(t)+sI(t)；w1(t)为加性噪声，设加性噪声为服从位置参数a＝0、对称参数β＝0的Alpha稳定分布噪声，若加性噪声为高斯分布或无干扰信号，不影响后续的估计步骤及估计结果；s(t)为具有循环平稳特性的发射信号；sI(t)为干扰信号，sI(t)和s(t)的载频相同，波特率不同；加性噪声项w1(t)和干扰项sI(t)为脉冲噪声与同频干扰并存下的信号模型；Sigmoid[x(t)]表示如公式(3)所示的对x(t)进行可调Sigmoid变换；Sigmoid[x(t)]=21+exp[-λx(t)]-1---(3)]]>1.2)与1.1)类似，利用可调Sigmoid变换循环自相关公式计算待测信号y(t)在时间区间[-T/2,T/2]上的可调Sigmoid变换循环自相关如公式(4)所示：Ry,TSϵ(τ)=ΔlimT→∞1T∫-T/2T/2Sigmoid[y(t+τ/2)]Sigmoid*[y(t-τ/2)]e-j2πϵtdt---(4)]]>采集数据为有限长度时的可调Sigmoid变换循环自相关可由公式(5)估计得到；R^y,TSϵ(τ)=1T∫-T/2T/2Sigmoid[y(t+τ/2)]Sigmoid*[y(t-τ/2)]e-j2πϵtdt---(5)]]>其中，为待测信号；D为待估计的时间延迟；fd为待估计的多普勒频移；w2(t)为加性噪声；1.3利用可调Sigmoid循环互相关公式计算参考信号与待测信号之间的可调Sigmoid循环互相关Rxy,TSϵ(τ)=ΔlimT→∞1T∫-T/2T/2Sigmoid[y(t+τ/2)]Sigmoid*[x(t-τ/2)]e-j2πϵtdt---(6)]]>采集数据为有限长度时的估计公式为：R^xy,TSϵ(τ)=1T∫-T/2T/2Sigmoid[y(t+τ/2)]Sigmoid*[x(t-τ/2)]e-j2πϵtdt---(7)]]>1.4利用可调Sigmoid循环互相关公式计算待测信号与参考信号之间的可调Sigmoid循环互相关Ryx,TSϵ(τ)=ΔlimT→∞1T∫-T/2T/2Sigmoid[x(t+τ/2)]Sigmoid*[y(t-τ/2)]e-j2πϵtdt---(8)]]>采集数据为有限长度时的估计公式为R^yx,TSϵ(τ)=1T∫-T/2T/2Sigmoid[x(t+τ/2)]Sigmoid*[y(t-τ/2)]e-j2πϵtdt---(9)]]>1.5根据广义Sigmoid循环模糊函数公式，计算得到广义Sigmoid循环模糊函数u和f分别代表时延和多普勒频移；Cyx,GSϵ(u,f)=Δ∫Rx,Sigmoidϵ(τ+u)[Ryx,Sigmoidϵ-f(τ)]*e-jπfτdτ+∫Ry,Sigmoidϵ(τ+u)[Ryx,Sigmoidϵ-f(τ)]*e-jπfτdτ+∫Rxy,Sigmoidϵ(τ+u)[Rx,Sigmoidϵ-f(τ)]*e-jπfτdτ+∫Rxy,Sigmoidϵ(τ+u)[Ry,Sigmoidϵ-f(τ)]*e-jπfτdτ---(10)]]>第二步，参数估计2.1)计算公式(10)结果中广义Sigmoid循环模糊函数的绝对值；2.2)搜索广义Sigmoid循环模糊函数的绝对值最大值的一组时延值u和多普勒频移f，为待测信号的估计值如公式(11)所示。(D^,f^d)=argmax|Cyx,GSϵ(u,f)|---(11)]]>