CN103885031B - 基于搜索优化滤波的运动卫星干扰源定位方法 - Google Patents

Abstract

本发明公开了基于搜索优化滤波的运动卫星干扰源定位方法，属于卫星通信中的干扰源定位领域。针对地面定位系统利用双星转发信号对运动卫星干扰源定位的问题，提出充分利用观测信息，通过对时差线上位置参考点的搜索优化，得到运动干扰源位置的预估计值，再使用滤波算法对预估计值进行处理，实现对运动卫星干扰源的定位。该发明是对目前静止卫星干扰源定位技术的发展，为实现对运动卫星干扰源的定位提供了一种可行的思路与方法。

Description

基于搜索优化滤波的运动卫星干扰源定位方法

技术领域

本发明属于卫星通信中的干扰源定位领域，具体说是基于搜索优化滤波的运动卫星干扰源定位方法。

背景技术

当前，卫星通信已广泛应用于电视转播、通信广播、导航定位、资源探测、军事侦察等领域。随着卫星通信业务的快速发展，卫星所处的电磁环境也日益恶化，射频干扰事件频繁发生，已成为卫星通信所要面临的常见问题，并有愈演愈烈之势，特别是当今越来越依赖空间信息基础设施的信息时代以及美国提出“进入空间，利用空间与控制空间”战略的激烈竞争年代，卫星通信的安全性已成为卫星通信运营商及其用户最为关切的焦点问题。为保障卫星通信系统的信息安全，查找、识别干扰源就显得尤为重要，准确地对射频干扰信号进行定位查找与识别，并采取针对性措施对干扰卫星通信业务的行为进行惩罚，是当前卫星通信系统面临的新需求。

卫星干扰源定位技术一般采用双星定位方法。当干扰源对准被干扰卫星实施干扰时，干扰源的发射天线将天线主瓣对准被干扰卫星（称为主星），同时其副瓣不可避免的指向了与被干扰卫星临近的卫星（称为邻星）。基于这种情况，就可以利用两部接收天线对主星和邻星转发的干扰信号一起检测，从而实现对干扰源的定位。

国外对卫星干扰源定位的研究较早，最早开始于1982年，由英国国防评估研究局（DERA）利用时间差和频率差在相关频段内进行干扰源定位的研究，并在1993年对干扰源定位涉及到的关键技术进行了详细的研究，之后DERA联合英国MARLIN国际通讯有限公司研制推出了多卫星地面定位系统（SatID）。同时，美国的相关研究部门也对这一技术进行了研究，并就其关键技术在美国申请了多项专利。1999年美国的Interferometrics公司研制成功了用于商业化的干扰源定位系统（TLSModel2000）。SatID系统和TLSModel2000系统均解决了查找地面干扰源的问题，受到了全球各大卫星公司、政府和军事部门的关注和重视。国内对卫星干扰源定位的研究起步较晚，多所高校和研究所对卫星干扰源定位技术展开了研究。

目前的卫星干扰源定位技术主要针对地面静止干扰源，而随着卫星通信移动平台的普遍使用，也可能出现以飞机、舰船、汽车等移动通信平台为发射载体的干扰源，此种干扰源的运动特性将给卫星干扰源定位带来极大的挑战，为此很有必要对运动卫星干扰源定位问题进行研究。

发明内容

本发明的目的在于避免上述背景技术中的不足，提供一种能对运动卫星干扰源定位的基于搜索优化滤波的运动卫星干扰源定位方法。

本发明的目的是这样实现的：

一种基于搜索优化滤波的运动卫星干扰源定位方法，其特点是包括下述步骤：

a.利用当前时刻运动卫星干扰源信号到达主星和到达邻星的时间差，并辅以地球方程，得到与当前时刻运动卫星干扰源位置相对应的等时差线；其中，主星是受运动卫星干扰源信号干扰的地球同步轨道卫星，邻星是与主星邻近的地球同步轨道卫星；

b.在步骤a中得到的等时差线上取若干对经纬度值，称为参考点，并在设定范围内对运动卫星干扰源的航向值和速度值进行搜索，由参考点、航向值和速度值组成参考矢量矩阵，参考矢量矩阵的每一个元素即为包含了运动卫星干扰源的经纬度、航向及速度的矢量；

c.利用参考矢量矩阵，构造与参考矢量矩阵对应的观测量的预测矩阵，并将预测矩阵中的元素与实际观测值进行比较，优化得到最接近实际观测量的预测矩阵元素，再以该元素为线索，倒推出其所对应的参考矢量矩阵中的位置参考点，即为运动卫星干扰源的位置预估计值；

d.对步骤c中得到的预估计值进行滤波处理，修正因搜索优化不准而产生较大误差的估计值，得到当前时刻的运动卫星干扰源位置的估计值；

e.更新当前时刻，重复步骤a～d，得到更新后时刻的运动卫星干扰源位置的估计值；

完成对运动卫星干扰源的定位。

其中，所述的运动卫星干扰源是指对卫星实施干扰的干扰源是运动的。

其中，所述的地面定位系统包括接收两颗卫星转发信号的地面接收站和一个参考辐射源发射站。

其中，步骤c中所述的观测量包括运动卫星干扰源信号到达主星、邻星的到达时间差和到达频率差以及地面站接收到的运动卫星干扰源信号经主星转发后的信号频率。

其中，所述运动卫星干扰源信号到达主星、邻星的到达时间差TDOA为：

$\mathrm{TDOA}=\frac{1}{\hat{c}}\left[(l_2-l_1)\right]=\mathrm{DTO}\left(\mathrm{unk}\right)-\mathrm{DTO}\left(\mathrm{ref}\right)+\frac{1}{\hat{c}}(l_2^r-l_1^r)---\left(1\right)$

其中，表示电磁波传播速度，l₁和l₂分别表示运动卫星干扰源到达主星、邻星的距离，和分别表示已知参考辐射源到达主星、邻星的距离，DTO(unk)和DTO(ref)分别表示运动卫星干扰源信号和参考辐射源信号经双星转发到达接收站的差分时差，用公式表示为：

$\mathrm{DTO}\left(\mathrm{unk}\right)=\frac{1}{\hat{c}}[(l_2+l_{m2})-(l_1+l_{m1})]=\frac{1}{\hat{c}}(l_2-l_1)+\frac{1}{\hat{c}}(l_{m2}-l_{m1})$

$\mathrm{DTO}\left(\mathrm{ref}\right)=\frac{1}{\hat{c}}[(l_2^r+l_{m2})-(l_1^r+l_{m1})]=\frac{1}{\hat{c}}(l_2^r-l_1^r)+\frac{1}{\hat{c}}(l_{m2}-l_{m1})$

其中，l_m1和l_m2分别为主星、邻星到达地面接收站的距离。

其中，所述的运动卫星干扰源信号到达主星、邻星的到达频率差FDOA为：

$\mathrm{FDOA}=\mathrm{DFO}\left(\mathrm{unk}\right)-\mathrm{DFO}\left(\mathrm{ref}\right)+(\frac{f_c}{\hat{c}}-\frac{f_0}{\hat{c}})(v_2\left({\stackrel{\→}{r}}_m\right)-v_1\left({\stackrel{\→}{r}}_m\right))-\frac{f_0}{\hat{c}}{v}_{21}\left({\stackrel{\→}{r}}_0\right)---\left(2\right)$

其中，f_c和f₀分别为干扰源和已知参考源的信号发射频率；为主星速度矢量在地面接收站径向上的分量，为邻星速度矢量在地面接收站径向上的分量，为两颗卫星速度矢量在参考源径向上速度分量的差值，其中和分别为主星、邻星的位置矢量和速度矢量，和分别为地面接收站的两个接收天线和参考源的位置矢量；DFO(unk)和DFO(ref)分别表示干扰源信号和参考辐射源信号经双星转发到达接收站的差分频差，用公式表示为：

$\begin{array}{c}\mathrm{DFO}\left(\mathrm{unk}\right)=-\frac{f_c}{\hat{c}}({\stackrel{\→}{v}}_{s2}-\stackrel{\→}{v})\·\frac{{\stackrel{\→}{r}}_{s2}-\stackrel{\→}{r}}{\left|\right|{\stackrel{\→}{r}}_{s2}-\stackrel{\→}{r}\left|\right|}-\frac{f_c-f_{T2}}{\hat{c}}\·{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}\left|\right|}\\ +\frac{f_c}{\hat{c}}({\stackrel{\→}{v}}_{s1}-\stackrel{\→}{v})\·\frac{{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}}{\left|\right|{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}\left|\right|}+\frac{f_c-f_{T1}}{\hat{c}}\·{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}\left|\right|}+(f_{T1}-f_{T2})\end{array}$

$\begin{array}{c}\mathrm{DFO}\left(\mathrm{ref}\right)=-\frac{f_0}{\hat{c}}{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_0}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_0\left|\right|}-\frac{f_0-f_{T2}}{\hat{c}}\·{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}\left|\right|}\\ +\frac{f_0}{\hat{c}}{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_0}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_0\left|\right|}+\frac{f_0-f_{T1}}{\hat{c}}\·{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}\left|\right|}+(f_{T1}-f_{T2})\end{array}$

其中，和分别为运动卫星干扰源的位置矢量和速度矢量，f_T1和f_T2分别为主星和邻星的信号转发频率。

其中，所述的地面站接收到的运动卫星干扰源信号经主星转发后的信号频率为：

$\tilde{f}=f_c+f_d+\mathrm{\Δcf}=f_c-\frac{f_c}{\hat{c}}({\stackrel{\→}{v}}_{s1}-\stackrel{\→}{v})\·\frac{{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}}{\left|\right|{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}\left|\right|}+\mathrm{\Δcf}---\left(3\right)$

其中，f_c为干扰源的信号发射频率，表示电磁波传播速度，f_d为运动卫星干扰源运动引起的多普勒，即Δcf为运动卫星干扰源信号频率测量误差，和分别为主星的位置矢量和速度矢量，和分别为运动卫星干扰源的位置矢量和速度矢量。

其中，所述步骤c中的观测量的预测矩阵是由参考矢量矩阵中的每一个元素，通过公式(1)～(3)计算得到的值来构成的。

其中，所述步骤c中的优化得到最接近实际观测量的预测矩阵元素，具体为：

利用参考矢量矩阵中的每一个元素，按照公式(1)～(3)计算与之对应的观测量的预测矩阵构造函数并求将该优化问题描述如下：

其中，|·|表示求绝对值，为实际观测量， $\hat{Z}\left(k\right)={[{\mathrm{\Λ}}_{\mathrm{TDOA}}\left(k\right),{\mathrm{\Λ}}_{\mathrm{FDOA}}\left(k\right),{\mathrm{\Λ}}_{\mathrm{Doppler}}\left(k\right)]}_{(n,m,l)}$ 为观测量的预测矩阵，N为参考点χ_n的个数，M为航向值的搜索个数，L为速度值V_l的搜索个数；所对应的即为运动卫星干扰源位置的预估计值。

其中，步骤d中所述的滤波处理，是采用卡尔曼滤波算法对运动卫星干扰源位置的预估计值进行实时处理的。

本发明与现有技术相比的有益效果为：

现有的卫星干扰源定位技术可以有效解决对地面静止干扰源的定位问题，但对运动卫星干扰源的定位问题却无能为力。本发明提出了一种基于搜索优化滤波的运动卫星干扰源定位方法，为解决对运动卫星干扰源定位问题提供了一种可行的思路与方法。

附图说明

图1是本发明方法的流程图；

图2是本发明方法所借助的地面定位系统的定位示意图；

图3是本发明方法的搜索示意图。

具体实施方式

以下结合附图举例对本发明做进一步的说明。

基于搜索优化滤波的运动卫星干扰源定位方法的流程图如图1所示；所借助的地面定位系统如图2所示，其中，地面定位系统包括接收两颗卫星转发信号的地面接收站和一个已知参考辐射源发射站。地面定位系统所利用的卫星转发信号，是来自两颗临近的地球同步轨道卫星，其中受干扰的卫星称为主星，另一颗卫星称为邻星。

该发明方法的基本步骤如下：

a.利用当前时刻k的到达时间差（TDOA）观测量，并辅以地球方程，得到与当前时刻运动卫星干扰源位置相对应的等时差线；

b.参照图3，在等时差线上取若干个参考点χ（即潜在的运动卫星干扰源位置的经纬度），并在设定范围内对运动卫星干扰源的航向值和速度值V进行遍历搜索，由参考点、航向值和速度值组成参考矢量矩阵（矩阵的每一个元素即为包含了运动卫星干扰源的经纬度、航向及速度的矢量）。其中，航向的搜索范围可设为全向，即范围为0°～360°；速度值的搜索范围可由运动卫星干扰源运动特性的先验知识来设定；

参考矢量矩阵为：

其中，χ_N为从时差线上选取的N个参考点，为运动卫星干扰源的航向搜索值，V_L为运动卫星干扰源的速度大小搜索值，M为航向值的搜索个数，L为速度值的搜索个数；

c.利用参考矢量矩阵，构造与之对应的观测量的预测矩阵，并将其中的元素与实际观测值进行比较，优化得到最接近实际观测量的预测矩阵元素，再以该元素为线索，倒推出其所对应的参考点，即为运动卫星干扰源位置的预估计值。具体实现步骤如下：

c1.观测量包括到TDOA、到达频率差（FDOA）、地面站接收到的运动干扰源信号经主星转发后的信号频率，具体表述如下：

运动卫星干扰源信号到达主星、邻星的到达时间差TDOA为：

$\mathrm{TDOA}=\frac{1}{\hat{c}}\left[(l_2-l_1)\right]=\mathrm{DTO}\left(\mathrm{unk}\right)-\mathrm{DTO}\left(\mathrm{ref}\right)+\frac{1}{\hat{c}}(l_2^r-l_1^r)---\left(1\right)$

其中，表示电磁波传播速度，l₁和l₂分别表示干扰源到达主星、邻星的距离，和分别表示已知参考辐射源到达主星、邻星的距离，DTO(unk)和DTO(ref)分别表示干扰源信号和参考辐射源信号经双星转发到达接收站的差分时差，用公式表示为：

$\mathrm{DTO}\left(\mathrm{unk}\right)=\frac{1}{\hat{c}}[(l_2+l_{m2})-(l_1+l_{m1})]=\frac{1}{\hat{c}}(l_2-l_1)+\frac{1}{\hat{c}}(l_{m2}-l_{m1})$

$\mathrm{DTO}\left(\mathrm{ref}\right)=\frac{1}{\hat{c}}[(l_2^r+l_{m2})-(l_1^r+l_{m1})]=\frac{1}{\hat{c}}(l_2^r-l_1^r)+\frac{1}{\hat{c}}(l_{m2}-l_{m1})$

其中，l_m1和l_m2分别为主星、邻星到达地面接收站的距离。

运动卫星干扰源信号到达主星、邻星的到达时间差TDOA为：

$\mathrm{FDOA}=\mathrm{DFO}\left(\mathrm{unk}\right)-\mathrm{DFO}\left(\mathrm{ref}\right)+(\frac{f_c}{\hat{c}}-\frac{f_0}{\hat{c}})(v_2\left({\stackrel{\→}{r}}_m\right)-v_1\left({\stackrel{\→}{r}}_m\right))-\frac{f_0}{\hat{c}}{v}_{21}\left({\stackrel{\→}{r}}_0\right)---\left(2\right)$

其中，f_c和f₀分别为干扰源和已知参考源的信号发射频率。为主星速度矢量在地面接收站径向上的分量，为邻星速度矢量在地面接收站径向上的分量，为两颗卫星速度矢量在参考源径向上速度分量的差值，其中和分别为主星、邻星的位置矢量和速度矢量，和分别为地面接收站的两个接收天线和参考源的位置矢量。DFO(unk)和DFO(ref)分别表示干扰源信号和参考辐射源信号经双星转发到达接收站的差分频差，用公式表示为：

$\begin{array}{c}\mathrm{DFO}\left(\mathrm{unk}\right)=-\frac{f_c}{\hat{c}}({\stackrel{\→}{v}}_{s2}-\stackrel{\→}{v})\·\frac{{\stackrel{\→}{r}}_{s2}-\stackrel{\→}{r}}{\left|\right|{\stackrel{\→}{r}}_{s2}-\stackrel{\→}{r}\left|\right|}-\frac{f_c-f_{T2}}{\hat{c}}\·{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}\left|\right|}\\ +\frac{f_c}{\hat{c}}({\stackrel{\→}{v}}_{s1}-\stackrel{\→}{v})\·\frac{{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}}{\left|\right|{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}\left|\right|}+\frac{f_c-f_{T1}}{\hat{c}}\·{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}\left|\right|}+(f_{T1}-f_{T2})\end{array}$

$\begin{array}{c}\mathrm{DFO}\left(\mathrm{ref}\right)=-\frac{f_0}{\hat{c}}{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_0}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_0\left|\right|}-\frac{f_0-f_{T2}}{\hat{c}}\·{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}\left|\right|}\\ +\frac{f_0}{\hat{c}}{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_0}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_0\left|\right|}+\frac{f_0-f_{T1}}{\hat{c}}\·{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}\left|\right|}+(f_{T1}-f_{T2})\end{array}$

其中，和分别为运动卫星干扰源的位置矢量和速度矢量，f_T1和f_T2分别为主星和邻星的信号转发频率。

地面接收到的运动干扰源信号经主星转发后的信号频率为：

$\tilde{f}=f_c+f_d+\mathrm{\Δcf}=f_c-\frac{f_c}{\hat{c}}({\stackrel{\→}{v}}_{s1}-\stackrel{\→}{v})\·\frac{{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}}{\left|\right|{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}\left|\right|}+\mathrm{\Δcf}$

其中，f_c为干扰源的信号发射频率，表示电磁波传播速度，f_d为运动卫星干扰源运动引起的多普勒，即Δcf为运动卫星干扰源信号频率测量误差；

c2.利用参考矢量矩阵中的每一个元素，按照公式(1)～(3)计算与之对应的观测量的预测矩阵构造函数并求将该优化问题描述如下：

其中，|·|表示求绝对值，为实际观测量， $\hat{Z}\left(k\right)={[{\mathrm{\Λ}}_{\mathrm{TDOA}}\left(k\right),{\mathrm{\Λ}}_{\mathrm{FDOA}}\left(k\right),{\mathrm{\Λ}}_{\mathrm{Doppler}}\left(k\right)]}_{(n,m,l)}$ 为预测量的观测值矩阵；所对应的即为运动卫星干扰源位置的预估计值；

d.对步骤c中得到的运动卫星干扰源位置的预估计值进行卡尔曼滤波实时处理，修正因搜索优化不准而产生较大误差的估计值，输出运动干扰源运动卫星干扰源位置的估计值；

e.更新当前时刻，重复步骤a～d，得到更新后时刻的运动卫星干扰源位置的估计值；

完成对运动卫星干扰源的定位。

Claims (10)

1.基于搜索优化滤波的运动卫星干扰源定位方法，是通过借助地面定位系统来实现的，其特征在于包括以下步骤：

a.利用当前时刻运动卫星干扰源信号到达主星和到达邻星的时间差，并辅以地球方程，得到与当前时刻运动卫星干扰源位置相对应的等时差线；其中，主星是受运动卫星干扰源信号干扰的地球同步轨道卫星，邻星是与主星邻近的地球同步轨道卫星；

b.在步骤a中得到的等时差线上取若干对经纬度值，称为参考点，并在设定范围内对运动卫星干扰源的航向值和速度值进行搜索，由参考点、航向值和速度值组成参考矢量矩阵，参考矢量矩阵的每一个元素即为包含了运动卫星干扰源的经纬度、航向及速度的矢量；

c.利用参考矢量矩阵，构造与参考矢量矩阵对应的观测量的预测矩阵，并将预测矩阵中的元素与实际观测量进行比较，优化得到最接近实际观测量的预测矩阵元素，再以该元素为线索，倒推出其所对应的参考矢量矩阵中的位置参考点，即为运动卫星干扰源的位置预估计值；

d.对步骤c中得到的预估计值进行滤波处理，修正因搜索优化不准而产生较大误差的预估计值，得到当前时刻的运动卫星干扰源的位置估计值；

e.更新当前时刻，重复步骤a～d，得到更新后时刻的运动卫星干扰源位置估计值；

完成对运动卫星干扰源的定位。

2.根据权利要求1所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：所述的运动卫星干扰源是指对卫星实施干扰的干扰源是运动的。

3.根据权利要求1所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：所述的地面定位系统包括接收两颗卫星转发信号的地面接收站和一个参考辐射源发射站。

4.根据权利要求1所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：步骤c中所述的观测量包括主星、邻星相对于运动卫星干扰源的到达时间差和到达频率差以及地面站接收到的运动卫星干扰源信号经主星转发后的信号频率。

5.根据权利要求4所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：所述的主星、邻星相对于运动卫星干扰源的到达时间差TDOA为：

$TDOA=\frac{1}{c}\[\left(l_2-l_1\right)\]=DTO\left(unk\right)-DTO\left(ref\right)+\frac{1}{c}\left(l_2^r-l_1^r\right)---\left(1\right)$

其中，c表示电磁波传播速度，l₁和l₂分别表示运动卫星干扰源到达主星、邻星的距离，和分别表示已知参考辐射源到达主星、邻星的距离，DTO(unk)和DTO(ref)分别表示运动卫星干扰源信号和参考辐射源信号经双星转发到达接收站的差分时差，用公式表示为：

$DTO\left(unk\right)=\frac{1}{c}\[(l_2+l_{m2})-(l_1+l_{m1})\]=\frac{1}{c}(l_2-l_1)+\frac{1}{c}(l_{m2}-l_{m1})$

$DTO\left(ref\right)=\frac{1}{c}\[(l_2^r+l_{m2})-(l_1^r+l_{m1})\]=\frac{1}{c}(l_2^r-l_1^r)+\frac{1}{c}(l_{m2}-l_{m1})$

其中，l_m1和l_m2分别为主星、邻星到达地面接收站的距离。

6.根据权利要求5所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：所述的主星、邻星相对于运动卫星干扰源的到达频率差FDOA为：

$FDOA=DFO\left(unk\right)-DFO\left(ref\right)+(\frac{f_c}{c}-\frac{f_0}{c})\left(v_2\right({\stackrel{\→}{r}}_m)-v_1({\stackrel{\→}{r}}_m\left)\right)-\frac{f_0}{c}{v}_{21}\left({\stackrel{\→}{r}}_0\right)---\left(2\right)$

其中，f_c和f₀分别为干扰源和已知参考源的信号发射频率；为主星速度矢量在地面接收站径向上的分量，为邻星速度矢量在地面接收站径向上的分量，为两颗卫星速度矢量在参考源径向上速度分量的差值，其中和分别为主星、邻星的位置矢量和速度矢量，和分别为地面接收站的两个接收天线和参考源的位置矢量；DFO(unk)和DFO(ref)分别表示干扰源信号和参考辐射源信号经双星转发到达接收站的差分频差，用公式表示为：

$\begin{array}{c}DFO\left(unk\right)=-\frac{f_c}{c}\left({\stackrel{\→}{v}}_{s2}-\stackrel{\→}{v}\right)\·\frac{{\stackrel{\→}{r}}_{s2}-\stackrel{\→}{r}}{\left|\right|{\stackrel{\→}{r}}_{s2}-\stackrel{\→}{r}\left|\right|}-\frac{f_c-f_{T2}}{c}\·{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}\left|\right|}\\ +\frac{f_c}{c}\left({\stackrel{\→}{v}}_{s1}-\stackrel{\→}{v}\right)\·\frac{{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}}{\left|\right|{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}\left|\right|}+\frac{f_c-f_{T1}}{c}\·{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}\left|\right|}+\left(f_{T1}-f_{T2}\right)\end{array}$

$\begin{array}{c}DFO\left(ref\right)=-\frac{f_0}{c}{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_0}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_0\left|\right|}-\frac{f_0-f_{T2}}{c}\·{\stackrel{\→}{v}}_{s2}\·\frac{{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}}{\left|\right|{\stackrel{\→}{r}}_{s2}-{\stackrel{\→}{r}}_{m2}\left|\right|}\\ +\frac{f_0}{c}{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_0}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_0\left|\right|}+\frac{f_0-f_{T1}}{c}\·{\stackrel{\→}{v}}_{s1}\·\frac{{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}}{\left|\right|{\stackrel{\→}{r}}_{s1}-{\stackrel{\→}{r}}_{m1}\left|\right|}+\left(f_{T1}-f_{T2}\right)\end{array}$

其中，f_T1和f_T2分别为主星和邻星的信号转发频率。

7.根据权利要求6所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：所述的地面站接收到的运动卫星干扰源信号经主星转发后的信号频率为：

$\tilde{f}=f_c+f_d+\mathrm{\Δ}cf=f_c-\frac{f_c}{c}({\stackrel{\→}{v}}_{s1}-\stackrel{\→}{v})\·\left|\right|{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}\left|\right|+\mathrm{\Δ}cf---\left(3\right)$

其中，Δcf为运动干扰源信号频率测量误差，f_c为干扰源的信号发射频率，c表示电磁波传播速度，和分别为主星的位置矢量和速度矢量，f_d为运动卫星干扰源运动引起的多普勒，即 $f_d=-\frac{f_c}{c}\left({\stackrel{\→}{v}}_{s1}-\stackrel{\→}{v}\right)\·\left|\right|{\stackrel{\→}{r}}_{s1}-\stackrel{\→}{r}\left|\right|.$

8.根据权利要求7所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：所述步骤c中的观测量的预测矩阵是由参考矢量矩阵中的每一个元素，通过公式(1)～(3)计算得到的值来构成的。

9.根据权利要求8所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：所述步骤c中的优化得到最接近实际观测量的预测矩阵元素，具体为：

利用参考矢量矩阵中的每一个元素，按照公式(1)～(3)计算与之对应的观测量的预测矩阵构造函数并求ming(χ_n,θ_m，V_l)，将该优化问题描述如下：

$\left\{\begin{array}{c}\min g\left({\mathrm{\χ}}_n,{\mathrm{\θ}}_m,V_l\right)\\ \begin{array}{cc}s.t.& n\∈1,2,...,N\\ & m\∈1,2,...,M\\ & l\∈1,2,...,L\end{array}\end{array}\right.---\left(4\right)$

其中，|·|表示求绝对值，为实际观测量矩阵，为预测观测值矩阵，ming(χ_n,θ_m，V_l)所对应的(χ_n,θ_m，V_l)即为运动卫星干扰源位置的预估计值。

10.根据权利要求1所述的基于搜索优化滤波的运动卫星干扰源定位方法，其特征在于：步骤d中所述的滤波处理，是采用卡尔曼滤波算法对预估计值进行实时处理。