专利名称：一种跳频信号时差估计方法
技术领域：
本发明属于信号处理范畴，涉及时差估计方法，尤其涉及跳频信号的高精度、快 速时差估计方法。
技术背景高精度时差估计在雷达、声纳、卫星和通信等系统的电磁波辐射源定位中具有重要 的应用价值。由于不同位置处接收的信号之间的时差与信号包络之间的互相关函数的峰 值位置之间存在一一对应的关系，因此可以利用该关系设计时差估计技术，通过对不同 位置处接收的信号进行互相关处理，实现信号时差估计。在基于这种原理的现有时差估 计技术中，典型的是广义互相关(General Cross Correlation)时差估计技术。由信号 处理和参数估计的基本理论可知，这类时差估计的方差的下界与信号带宽的三次方成反 比，与信噪比和信号的时间长度成反比。因此，针对窄带信号的时差估计问题，需要较 大的信噪比和较长的信号时间长度，才能获得较高精度的窄带信号时差估计。当信号是跳频信号时，接收的信号在整个信号观测时间内按跳频周期随机的改变跳 频频偏。由于具有相同跳频频偏的信号仅占整个信号观测时间的一部分，使得对应该跳 频频偏的信号在整个信号观测时间内的平均信息速率远小于每个跳频周期内的信息速 率，即对应每个跳频频偏的信号的基带带宽实际上变得更窄了。此时如果仍利用不同位 置处接收的信号之间的时差与信号包络之间的互相关函数的峰值位置之间存在的一一对 应关系设计跳频信号时差估计方法，时差估计的方差的下界将难以满足实际应用的需要。若要避免对应每个跳频频偏的信号的基带带宽变得更窄，则只能利用每个跳频周期 的接收信号估计信号之间的时差。由于一个跳频周期的时间是有限的，将信号观测时间 限制为一个跳频周期也将使时差估计的方差的下界难以满足实际应用的需要。可见，如能利用较短时间内在两个不同位置处接收的跳频信号得到高精度的时差估计，将在保持高精度跳频信号时差估计性能的同时，有效的提高跳频信号时差估计的速 度，从而对跳频辐射信号源定位具有重要意义。 发明内容本发明的目的是为跳频辐射信弓源的精确、快速定位提供一种高精度、快速的时差估 计方法。本发明的目的是这样达到的一种跳频信号时差估计方法，其特征在于利用两个不同位置处接收的跳频信号之 间的时差与这两个不同位置处接收的跳频信号载波之间的相位差存在的线性对应关系， 通过测量这两个不同位置处接收的跳频信号载波之间的多个跳频周期内的相位差，估计 这两个不同位置处接收的跳频信号之间的时差。其具体步骤是根据一个跳频周期(r)的载波频率和信号带宽，对两个不同位置处接收的跳频信号进行带通滤波；确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变换； 利用Hilbert变换结果分别确定两个不同位置处接收的跳频信号在一个跳频周期内 的相位；确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差； 确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值； 依次确定两个不同位置处接收的跳频信号在K (大于1的整数)个跳频周期内的相 位差的中值序列；由K个跳频周期内相位差的中值序列对应的线性方程组的最小二乘解确定跳频信号 的时差估计。所述根据一个跳频周期内的载波频率和带宽，对两个不同位置处接收的跳频信号X, (z+;tr)和x2 (r+奵)分别进行带通滤波，结果为-<formula>formula see original document page 6</formula>和<formula>formula see original document page 7</formula>其中，-7V2^fS772， a(m)为根据一个跳频周期内的载波频率和带宽设计的带通滤波 器的系数，w = l，2，.-.,M。为-所述确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变换分别<formula>formula see original document page 7</formula>和<formula>formula see original document page 7</formula>其中，乂 = 7^是单位虚数￡所述利用Hilbert变换结果分别确定两个不同位置处接收的跳频信号在一个跳频周 期内的相位分别为 <formula>formula see original document page 7</formula>和 p2<formula>formula see original document page 7</formula>其中，angle()表示取复数的相角。 所述确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差为W+々r) =+ at)o+at)其中，-772S"772。 所述确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值为-.A = median k(,+")] 其中'median[^+w)]表示^(f+^)在一"2s""2范围内的中值。-772^772所述依次确定两个不同位置处接收的跳频信号在^个跳频周期的相位差的中值序列为<formula>formula see original document page 7</formula>所述由尺个跳频周期内相位差的中值序列对应的线性方程组:2:/2 124 1的解2;r乂 1 2#2 1力为第A个跳频周期内相对确定跳频信号的时差估计为f，其中，矩阵力=2tt/k 1—于载频的跳频频偏，()'和()—'分别表示矩阵的转置和矩阵的逆。本发明的积极效果为跳频辐射信号源的精确、快速定位提供了一种高精度、快速 的时差估计方法。满足了在雷达、声纳、卫星和通信等系统的电磁波辐射源定位中高精 度时差估计的需要，为电磁波监测技术的发展和应用拓展了新的发展空间。根据本发明的方法进行的实验表明，将本发明提出的跳频信号时差估计方法应用于每秒250跳，驻 留时间3. 5毫秒，跳频频偏依次为-25kHz， 225 kHz， 0kHz， -250 kHz， 250 kHz, 775 kHz 的跳频信号，在两个不同位置处接收的跳频信号时差为493纳秒，基带信噪比15dB的情 况下，估计的时差为534纳秒，时间为6个跳频周期，即24毫秒。
图1为本发明一种跳频信号时差估计方法的具体步骤流程图； 图2为本发明一种跳频信号时差估计方法的相位差模型拟合图。
具体实施方式
附图给出了本发明的一个实施例。以下结合附图对本发明的具体实施方式
进行说明。术语与符号定义AW和A(/) 两个不同位置处接收的跳频信号r 跳频周期r 跳频信号到达两个不同位置处的时差K 跳频周期数(大于l的整数)Hilbert变换 印)的Hilbert变换是XO-^) +丄义. 源信号的载频/和；T /时刻和Z-r时刻相对于载频/的跳频频偏%和"2 两个不同位置处接收的跳频信号的幅度变化3和《 两个不同位置处接收的跳频信号的初相v,(/)和vj/) 加性白噪声angle() 表示复数的相角median[州+^)] 表示州+w)在-:r / 2 a ^ r / 2范围内的中间值-772^772以两个不同位置处按收的跳频信号:c^)和A(0之间的时差估计为例子。假设跳频 信号到达两个不同位置处的时差为r ，在f时刻两个不同位置处接收的跳频信号x,(f)和 x2(0分别为乇0)-a^0)cos(2;r(/c+乂V + e,) + v!(y) (1)和x2(f) = a^G —r)cos(2;r(y; + /_r)(, —"+《)+v2(0 (2) 其中，力)为信号包络，X为源信号的载频，/和L为？时刻和卜r时刻相对于载频义的 跳频频偏，a,和A为两个不同位置处接收的跳频信号的幅度变化，0,和《为两个不同位 置处接收的跳频信号的初相，K(/)和v2W为加性白噪声。在f时刻两个不同位置处接收的跳频信号的相位差为州)=[2;r(义,+ L )(f — r) +《]一 [2;r(X ③ =一2tt(rr + 2;r(/;_r — ,)f +《一《一 2;r/cr9记跳频周期为r，由于时差r一般都是远小于跳频周期r，所以在一个跳频周期内的大部 分时刻两路接收信号的跳频频偏都满足于是，式(3)又可写成<formula>formula see original document page 10</formula> (4) 与得到(3)式的道理相同，在/ + ^时刻两个不同位置处接收的跳频信号的相位差为 <formula>formula see original document page 10</formula>(5) 其中，A-0，l,2,…,《-l， K是接收信号期间跳频信号的跳频次数或跳频周期数。由(3)-(5)式可见，两个不同位置处接收的跳频信号之间的时差与这两个不同位置处 接收的跳频信号载波之间的相位差之间存在线性对应关系，该线性对应关系受相对于载 频的跳频频偏的影响。若已知或预先估计相对于载频的跳频频偏，则可利用两个不同位 置处接收的跳频信号之间的时差与这两个不同位置处接收的跳频信号载波之间的相位差 存在的线性对应关系，通过测量这两个不同位置处接收的跳频信号载波之间的多个跳频 周期内的相位差，估计这两个不同位置处接收的跳频信号之间的时差。图1示出根据本发明的设计跳频信号时差估计方法的流程图。流程开始于步骤101。 在步骤102，根据一个跳频周期内的载波频率和带宽，对两个不同位置处接收的跳频信 号;c, (/ + W)和x2 + 分别进行带通滤波，结果为<formula>formula see original document page 10</formula> (6)和<formula>formula see original document page 10</formula> (7)其中，-772S/^r/2, "(w)为根据一个跳频周期内的载波频率和带宽设计的带通滤波 器的系数，》i = l,2，...,A/。在步骤1031,确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变 换分别为<formula>formula see original document page 10</formula> (8)和<formula>formula see original document page 11</formula>丌Jm其中，y-V^I是单位虚数。(9)和在步骤1032，确定两个不同位置处接收的跳频信号在一个跳频周期内的相位分别为<formula>formula see original document page 11</formula> 其中，angle()表示取复数的相角。在步骤1041，确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差为-0++at) (12)其中，-772<"r/2。在步骤1042，确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中 值为A = median [萍+^)] (13) 其中'median[州+w)]表示-("w)在-"2^""2范围内的中值。在步骤105，依次确定两个不同位置处接收的跳频信号在《个跳频周期的相位差的中值序列为<formula>formula see original document page 11</formula>(14)在步骤106，由《个跳频周期内相位差的中值序列对应的线性方程组:<formula>formula see original document page 12</formula>(15)的解:<formula>formula see original document page 12</formula>(16)确定跳频信号的时差估计为r"，其中，矩阵<formula>formula see original document page 12</formula>于载频的跳频频偏，()'和()—'分别表示矩阵的转置和矩阵的逆。 根据本发明的设计跳频信号时差估计方法的流程结束于步骤107 定义相位差中值序列拟合误差为y;为第A个跳频周期内相对<formula>formula see original document page 12</formula>(17)将本发明提出的跳频信号时差估计方法应用于每秒250跳，驻留时间3.5毫秒，跳 频频偏依次为-25kHz, 225 kHz， 0 kHz， -250 kHz, 250 kHz， 775 kHz的跳频信号，在 两个不同位置处接收的跳频信号时差为493纳秒，夂=6，带内信噪比15dB的情况下， 相位差中值序列拟合误差如图2所示。根据本发明的跳频信号时差估计方法估计的时差 为534纳秒，时间为6个跳频周期，即24毫秒。虽然已经参考附图对本发明的跳频信号时差估计方法以举例方式进行了描述，但是 本发明不限于上述这些细节，本申请含盖权利要求范围之内的各种变型或改变。工业应用性将本发明提出的跳频信号时差估计方法应用于跳频辐射信号源时差定位系统，满足 跳频辐射信号源时差定位系统对快速时差估计、高精度时差估计性能的要求。
权利要求
1、一种跳频信号时差估计方法，其特征在于利用两个不同位置处接收的跳频信号之间的时差与这两个不同位置处接收的跳频信号载波之间的相位差存在的线性对应关系，通过测量这两个不同位置处接收的跳频信号载波之间的多个跳频周期内的相位差，估计这两个不同位置处接收的跳频信号之间的时差；其具体步骤是根据一个跳频周期的载波频率和信号带宽，对两个不同位置处接收的跳频信号进行带通滤波；确定两个不同位置处接收的跳频信号在一个跳频周期内的Hilbert变换；利用Hilbert变换结果分别确定两个不同位置处接收的跳频信号在一个跳频周期内的相位；确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差；确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值；依次确定两个不同位置处接收的跳频信号在K(大于1的整数)个跳频周期内的相位差的中值序列；由K个跳频周期内相位差的中值序列对应的线性方程组的最小二乘解确定跳频信号的时差估计。
2、如权利要求1所述的跳频信号时差估计方法，其特征在于-所述根据一个跳频周期内的载波频率和带宽，对两个不同位置处接收的跳频信号分别 进行带通滤波是指对跳频信号x, + AT1)和x2+ W)进行带通滤波，其中<formula>formula see original document page 2</formula> a(w)为根据一个跳频周期内的载波频率和带宽设计的带通滤波器的系数，w = l,2,...，M。
3、如权利要求1所述的跳频信号时差估计方法，其特征在于所述确定两个不 同位置处接收的跳频信号在一个跳频周期内的Hilbert变换分别为其滤波结果为:<formula>formula see original document page 2</formula>和<formula>formula see original document page 2</formula>其中，y'-VTT是单位虛数。
4、 如权利要求1所述的跳频信号时差估计方法，其特征在于所述利用 Hilbert变换结果分别确定两个不同位置处接收的跳频信号在一个跳频周期内的 相位分别为a G <formula>formula see original document page 2</formula>其中，angle()表示取复数的相角。
5、 如权利要求1所述的跳频信号时差估计方法，其特征在于所述确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差为<formula>formula see original document page 2</formula>
6、 如权利要求1所述的跳频信号时差估计方法，其特征在于所述确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差的中值为<formula>formula see original document page 2</formula>其中，mediant"w)]表示W+^)在—7V2《"r/2范围内的中值。
7、 如权利要求1所述的跳频信号时差估计方法，其特征在于所述依次确定两个不同位置处接收的跳频信号在/:个跳频周期的相位差的中值序列为<formula>formula see original document page 2</formula>
8、如权利要求i所述的跳频信号时差估计方法，其特征在于所述由/:个跳频周期内相位差的中值序列对应的线性方程组 <formula>formula see original document page 2</formula>的解: <formula>formula see original document page 4</formula>最后确定跳频信号的时差估计为f，其中，矩阵j-<formula>formula see original document page 4</formula>，A为第A个跳频周期内相对于载频的跳频频偏，()'和(广分别表示矩阵的转置和矩阵的逆。
全文摘要
跳频信号时差估计方法属信号处理范畴。利用两个不同位置处接收的跳频信号间的时差与这两个位置接收的跳频信号载波之间的相位差的线性对应关系，通过测量这两个不同位置处接收的跳频信号载波之间的多个跳频周期内的相位差，估计这两个不同位置处接收的跳频信号之间的时差。包括对两个不同位置处接收的跳频信号进行带通滤波；在一个跳频周期内的Hilbert变换、计算相位；确定两个不同位置处接收的跳频信号在一个跳频周期内的相位差、相位差的中值；在K个跳频周期内的相位差中值序列；解出中值序列对应的线性方程组的最小二乘解确定跳频信号的时差估计等步骤。本方法可对跳频信号的高精度、快速的估计时差，广泛应用于跳频辐射信号源定位系统。
文档编号H04B1/713GK101247141SQ20071005070
公开日2008年8月20日 申请日期2007年12月5日 优先权日2007年12月5日
发明者群 万, 军 王, 琰 王, 罗来源, 魏合文 申请人:电子科技大学;西南电子电信技术研究所