1.一种基于正弦调频SFM信号线性模型的伪码调相正弦调频复合信号伪码序列估计方法,其特征在于,包括以下步骤:对SFM信号进行第一类贝塞尔函数级数展开,将复正弦调频信号s(t)=A exp{j[ω0t+mfsin(ωmt)]}展开为以Bessel函数为系数的指数级数,取Bessel函数的非负阶数,得到SFM信号的单边线性信号将SFM信号的单边线性信号s′(t)复合伪码调相信号,建立基于SFM信号线性模型的伪码调相复合信号;采用平滑伪Wigner-Ville分布SPWVD对该伪码调相复合信号进行时频分析,寻找伪码相位跳变点处出现的尖锐负脉冲;采用基于奇异值分解SVD计算奇异值Sx=UΛxVH,H保留对角矩阵Λx第一个最大奇异值,其它奇异值置零得到Λx′,再恢复计算Sx′=UΛx′V 得到奇异值Sx′;求出Sx′沿频率轴的最小值切面图,在最小值切面图中搜索负尖峰对应的时刻,完成伪码调相复合信号PN码序列的估计,其中,K为第一类Bessel函数的最高阶数,A为常幅度,Jk(mf)为第k个载波分量幅度,mf为调制指数,j为虚单位,ω0为载波角频率,ωm为调制角频率,t为时间变量,Sx表示含有噪声的伪码复合信号时频分布,Sx′表示SVD去噪增强处理后的伪码复合信号时频分布,H表示共轭转置,U、V为左右奇异矩阵,Λx为奇异值对角矩阵。
2.根据权利要求1所述的估计方法,其特征在于,基于SFM线性模型的伪码调相复合信号是将伪码的频率搬移到频率f=f0+kfm上,其中,f0为载波分量,k为第一类Bessel函数的阶数,fm为调制频率。
3.根据权利要求1所述的估计方法,其特征在于,计算基于SFM线性模型的伪码调相复合信号的SPWVD变换得到信号的时频分布,对信号时频分布进行SVD分解后,保留奇异值对角矩阵的第一个最大奇异值,将其它奇异值置零,再由SVD恢复计算出经去噪增强后的时频分布,获取时频面沿频率轴的最小值切面,检测最小值切面的峰值,估计伪码调相复合信号的PN码序列。
4.根据权利要求3所述的估计方法,其特征在于,检测最小值切面的峰值具体为,采用门限法,将门限设为切面图中最小值的一半,存在连续的SPWVD值小于门限值的三个时刻点,且满足中间点为局部极小值的位置为最小值切面的峰值。