CN103576166B - 基于单颗静止轨道卫星对干扰源进行定位的方法及系统 - Google Patents

基于单颗静止轨道卫星对干扰源进行定位的方法及系统 Download PDF

Info

Publication number
CN103576166B
CN103576166B CN201210277476.0A CN201210277476A CN103576166B CN 103576166 B CN103576166 B CN 103576166B CN 201210277476 A CN201210277476 A CN 201210277476A CN 103576166 B CN103576166 B CN 103576166B
Authority
CN
China
Prior art keywords
signal
satellite
doppler shift
frequency
source
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201210277476.0A
Other languages
English (en)
Other versions
CN103576166A (zh
Inventor
刘岩
李景春
张小飞
李安平
夏楠
唱亮
魏梅英
刘海洋
闫肃
赖幸君
崔晓曼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NATIONAL RADIO MONITORING CENTER
Original Assignee
NATIONAL RADIO MONITORING CENTER
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NATIONAL RADIO MONITORING CENTER filed Critical NATIONAL RADIO MONITORING CENTER
Priority to CN201210277476.0A priority Critical patent/CN103576166B/zh
Publication of CN103576166A publication Critical patent/CN103576166A/zh
Application granted granted Critical
Publication of CN103576166B publication Critical patent/CN103576166B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/015Arrangements for jamming, spoofing or other methods of denial of service of such systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position

Abstract

本发明实施例提供一种基于单颗静止轨道卫星对干扰源进行定位的方法及系统,所述的系统包括:信号采集装置,用于采集静止轨道卫星上的干扰信号;参考源信号发射装置,用于向所述的静止轨道卫星发送参考源信号,并接收所述的静止轨道卫星根据所述的参考源信号发送的参考源返回信号;信号处理装置,用于对所述的干扰信号以及所述的参考源返回信号的实际多普勒频移进行误差排除;干扰源定位装置,用于根据所述的干扰信号以及所述的参考源返回信号确定干扰源的位置。仅通过一颗静止轨道卫星即能实现对卫星干扰源的准确定位,大大地提升了卫星对干扰的查找能力,降低了干扰查找的诸多束缚条件。

Description

基于单颗静止轨道卫星对干扰源进行定位的方法及系统
技术领域
[0001] 本发明关于无线电技术,特别是关于静止卫星通信领域,具体的讲是一种基于单 颗静止轨道卫星对干扰源进行定位的方法及系统。
背景技术
[0002] 随着无线通信技术的发展,人们对于频谱资源的需求愈来愈大,"频谱枯竭"现象 日趋明显。静止卫星通信领域也不可避免地出现了由于频谱紧张等因素而引起的无线电干 扰事件,通信卫星越来越多地受到地面发射源的干扰,这些干扰给正常的卫星广播和通信 业务的正常开展带来了巨大的威胁,甚至可能产生无法估量的社会影响。
[0003] 为了应对这些射频干扰,通过受干扰卫星迅速而有效地确定地面干扰源位置是非 常必要的。当前,主要采用图1所示的双星定位系统对卫星干扰源进行定位,主要测定访问 地球同步人造通信卫星的上行链路发射机的位置。这是通过测量上行链路信号通过两个不 同的地球同步人造卫星而造成的到达接收站的时间差(TD0A)和到达频率差(FD0A)来实现 的。利用测量的TD0A、FD0A和准确的卫星轨道信息,计算上行发射机的位置。一个卫星上 行链路发射机访问一个地球同步人造卫星(主星),上行链路天线旁瓣的一些功率会辐射到 同步的邻近卫星。如果邻近卫星(邻星)有上行接收机,这个接收机频率和极化与发射的上 行信号一致,未知信号将在邻星相应的下行链路出现。
[0004] 上述的双星定位技术,是基于两颗静止卫星,通过计算干扰信号通过这两颗静止 卫星在接收站点所产生的时间差和位置差计算出干扰源的地理位置。但是,随着轨道位置 越来越密集、不同频段的大范围使用以及转发器的利用率逐年升高等因素,双星定位技术 也愈来愈受局限,即是否有符合条件的邻星来实现双星定位技术已经成为能否成功精确定 位的主要瓶颈之一。由于卫星轨位、频段、极化、发射天线口径、覆盖等参数的影响,很多时 候基于双星定位技术是不能够进行定位,或者不能满足前面提到双星定位条件,因此需要 一种全新的定位技术来解决上述干扰定位问题。
发明内容
[0005] 本发明提出一种基于单颗静止轨道卫星对干扰源进行定位的方法及系统,通过对 干扰信号进行采集,并通过一系列的后台处理,最终即可得到可靠的干扰源位置,用于解决 当前日趋增多的卫星干扰,对维护空中电波秩序,优化卫星频率资源,维护社会稳定有不可 替代的作用。
[0006] 本发明的目的之一是,提供了一种基于单颗静止轨道卫星对干扰源进行定位的方 法,所述的方法包括:采集静止轨道卫星上的干扰信号;向所述的静止轨道卫星发送参考 源信号;接收所述的静止轨道卫星根据所述的参考源信号发送的参考源返回信号;对所述 的干扰信号以及所述的参考源返回信号的实际多普勒频移进行误差排除;根据所述的干扰 信号以及所述的参考源返回信号确定干扰源的位置。
[0007] 本发明的目的之一是,提供一种基于单颗静止轨道卫星对干扰源进行定位的系 统,所述的系统包括:信号采集装置,用于采集静止轨道卫星上的干扰信号;参考源信号发 射装置,用于向所述的静止轨道卫星发送参考源信号,并接收所述的静止轨道卫星根据所 述的参考源信号发送的参考源返回信号;信号处理装置,用于对所述的干扰信号以及所述 的参考源返回信号的实际多普勒频移进行误差排除;干扰源定位装置,用于根据所述的干 扰信号以及所述的参考源返回信号确定干扰源的位置。
[0008] 本发明的有益效果在于,仅通过一颗静止轨道卫星即能实现对卫星干扰源的准确 定位,大大地提升了卫星对干扰的查找能力,降低了干扰查找的诸多束缚条件,弥补了当前 主流定位技术的不足,填补了我国在卫星定位系统的空白,极大的促进了无线电管理水平 的提尚。
附图说明
[0009] 为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本 发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可 以根据这些附图获得其他的附图。
[0010] 图1为现有双星定位系统的原理图;
[0011] 图2为本发明实施例提供的一种基于单颗静止轨道卫星对干扰源进行定位的系 统的原理图;
[0012] 图3为本发明实施例中的天线系统的组成原理图;
[0013] 图4为本发明实施例中的射频接收单元组成原理图;
[0014] 图5为频率测量精度验证平台的结构框图;
[0015] 图6为信号频率随时间的变化图;
[0016] 图7为信号的频率差随时间变化图;
[0017] 图8为步骤S7的具体流程图;
[0018] 图9为X、Y、Z三维地固坐标系统;
[0019] 图10为单星定位中地固坐标系中地球与卫星的模型;
[0020] 图11为zhongxinglO卫星一周内的运动轨迹图;
[0021] 图12为zhongxinglO卫星一天X、Y、Z方向下速度随时间的变化图;
[0022] 图13为基于四阶差分法计算出的多普勒频移随时间变化图;
[0023] 图14为基于四阶差分法计算出的多普勒频移变化率随时间变化图;
[0024] 图15为星上转发器频率随时间的变化曲线图;
[0025] 图16为小波去噪前后频率变化曲线图;
[0026] 图17为拟合值与真实值的误差图;
[0027] 图18为实际多普勒频移的数据筛选图;
[0028]图19为本发明实施例提供的一种基于单颗静止轨道卫星对干扰源进行定位的方 法的流程图;
[0029] 图20为上行多普勒频移理论与实测波形图;
[0030] 图21为经过滤波后的普勒频移与理论多普勒频移图;
[0031] 图22为截取2100点至2300点的数据图;
[0032] 图23为北京站经炜度定位结果图;
[0033] 图24为误差随噪声标准差变化曲线图;
[0034] 图25为成都站多普勒频率滤波结果图;
[0035] 图26为成都站经炜度定位结果图;
[0036] 图27为深圳站多普勒频率滤波结果图;
[0037] 图28为深圳站经炜度定位结果图。
具体实施方式
[0038] 下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完 整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于 本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例,都属于本发明保护的范围。
[0039] 本发明提供了一种基于单颗静止轨道卫星对干扰源进行定位的系统,如图2所 示,未知发射源向单颗静止轨道卫星发送干扰信号,所述的定位系统对应的定位方法具体 包括:
[0040]S1:采集静止轨道卫星上的干扰信号。
[0041] 在对卫星的干扰源进行定位时,首先需要保证能够接收到该卫星收到的干扰信 号。基于卫星通信的特点及不同频段的需要,在具体的实施方式中可通过不同口径及不同 馈源的抛物面天线来实现。例如,在该系统的一种实施方式中,只对C/Ku频段的信号进行 分析并定位,信号采集装置采用7. 3米C/Ku的双频段天线系统。该天线系统主要由天线主 副反射面、C/Ku双频段馈源(波纹喇叭和微波网络)、天线结构、方位一俯仰型全转台座架、 方位一俯仰驱动及其同步装置、天线控制器、射频接收单元(LNA、功分器、波导开关、C&Ku 跟踪下变频器)、跟踪接收机、安全保护装置等组成。天线系统组成原理框图如图3所示,射 频接收单元组成原理框图如图4所示。
[0042] 7. 3米C/Ku双频段天线系统的主要功能是接收来自地球静止轨道卫星的无线电 波信号,并将接收到的信号高质量的传输至信号处理装置。由于该定位系统中仅需要接收 C、Ku两个频段的信号,且对信号的信噪比、不同极化的信号、不同轨道位置的卫星上的干 扰、邻星泄露信号、系统噪声等方面有较高要求,因此7. 3米C/Ku双频段天线系统的主要技 术指标如下:
[0043] >工作频段:
[0044]C频段:3. 4GHz~4. 2GHz
[0045]Ku频段:10. 7GHz~12. 75GHz
[0046] >天线增益:
[0047]C频段:彡 47. 3+201g(f/3. 9)dBif单位:GHz
[0048]Ku频段:彡 57. 4+201g(f/2. 5)dBif单位:GHz
[0049] >极化方式:左、右旋圆极化/双线极化可选,电动线/圆极化转换,水平和垂直线 极化同时接收,电动线极化调整
[0050] >线极化调整范围:±50°
[0051] >天线噪声温度:(在晴天、微风10°仰角条件下)
[0052]C频段:Ta彡 45。K
[0053]Ku频段:Ta彡 50。K
[0054] >第一旁瓣:彡一 14dB
[0055] 宽角旁瓣:(90%峰值满足如下包络线):
[0056] 32-251g0 (dBi) 1° 彡 9 彡 48。
[0057] -10 (dBi) 0 >48。
[0058](式中0为偏离主瓣波束的角度)
[0059] >圆极化轴比:< 1 • 09
[0060] >线极化交叉极化隔离度:彡30dB(轴向)
[0061] >驻波(VSWR):彡 1. 3
[0062] >端口隔离:
[0063] 左、右圆极化端口:彡20dB
[0064] 水平、垂直线极化端口:彡35dB
[0065]C/Ku频段端 口:多35dB
[0066] >馈源插入损耗:
[0067] C频段:<0.3dB
[0068] Ku频段:彡0.45dB
[0069]C波段和Ku波段下变频器从天线接收单路射频输入,并且下变频到L波段用于匹 配处理速度。下变频器工作于扩展的RF输入范围。然而,每个基带变频器每次只能处理 500MHz频谱。因此,下变频器把RF频谱分成n个500MHz段中的一段,并把这一段下变频到 L波段。表1是有重叠的分段。
[0070]表1
Figure CN103576166BD00101
[0072] 射频切后单元是将切换后的L波段信号(包括Ku或C波段)作为输入,根据定位算 法的需要,将不同的输入端软件切换至不同的输出端。在本发明中,要求至少4*4路切换开 关。
[0073] 本发明的其他实施方式中,可采用了两路高速宽带实时频谱仪进行信号采集。天 线是射频前端,是用来接收射频信号的,而为了进一步对信号进行处理,需要用诸如高速宽 带实时频谱仪的采集装置对信号进行采集。为了满足足够的频率分辨率及定位算法的需 要,在信号采集过程中设置了流程化的采样。即步骤S1具体包括:
[0074] S101 :设置采样中心频率。由于地面信号经过移动的转发器后,信号频率会发生漂 移,因此,要设置合适的中心频率保证能够对信号频率偏移的情况下不丢失频率数据。
[0075]S102:设置采样带宽。根据信号可能的多普勒频移大小,预估计一个带宽,取一定 的冗余,保证在此带宽下,不丢失频率数据。本发明中,至少需要保证采样带宽为预测的多 普勒频移峰值的2. 5倍。
[0076]S103:设置采样时长。由于在计算多普勒频移的时候,需要保证多普勒频移的精度 及分辨率,因此过短的采样时长不能保证分辨率,过长的时长易引起时间延迟并导致时间 不同步。并且,由于测频算法中,需要根据所设置的采样时间,综合多种因素,选择不同的滤 波方法,保证频率的精度和分辨率。
[0077]S104:设置采样数据存储方式。由于对于不同的信号,其采样速率和采样带宽都不 一样,因此根据定位算法的需要,需要设置多长时间存储一个数据文件并设置要采样多少 个文件参与定位。本具体的实施方式中,根据中星10的经炜度计算,在Ku频段的单载波, 信号的多普勒频移大约为300Hz-400Hz之间,因此设置了采样带宽为lKHz,采样时长为1分 钟为一个单位数据文件。
[0078] S105 :采集静止轨道卫星上的干扰信号;
[0079]S106 :对所述的干扰信号进行频率估计;
[0080] S107 :对所述的干扰信号进行筛选。
[0081]图19为本发明实施例提供的一种基于单颗静止轨道卫星对干扰源进行定位的方 法的流程图,由图19可知,基于单颗静止轨道卫星对干扰源进行定位的方法还包括:
[0082]S2:向所述的静止轨道卫星发送参考源信号;
[0083]S3:接收所述的静止轨道卫星根据所述的参考源信号发送的参考源返回信号;
[0084]S4:确定所述的干扰信号以及所述的参考源返回信号对应的理论多普勒频移以及 实际多普勒频移;
[0085]S5:对所述的干扰信号以及所述的参考源返回信号对应的实际多普勒频移进行误 差排除;
[0086]S6:根据所述的干扰信号以及所述的参考源返回信号确定干扰源的位置。除了上 述步骤,在本发明的其他实施方式中,本发明还包括S7:确定所述的静止轨道卫星的频率 漂移。图8为步骤S7的具体流程图,由图8可知,该步骤具体包括:
[0087]S71:向所述的静止轨道卫星发送标准频率信号;
[0088]S72:接收所述的静止轨道卫星根据所述的标准频率信号发送的标准频率返回信 号;
[0089]S73:根据所述的标准频率返回信号确定其对应的频率。该步骤与计算干扰信号的 频率采用相同的方法,即通过采集到的标准频率返回信号,采用四阶差分法计算出精确的 频率值。
[0090]S74:确定所述的标准频率返回信号的多普勒频移。通过步骤S73计算出精确的频 率值后,恢复出有多普勒和卫星本振漂移产生的频率随时间的变化值。
[0091]S75:根据所述的标准频率返回信号对应的频率以及多普勒频移确定所述的静止 轨道卫星的频率漂移。在具体的实施方式中,该步骤可通过将标准频率信号与标准频率返 回信号相减,即可得到卫星的频率偏差。
[0092] 标准频率是提供一个精确的频率,为其他信号提供参考,如接收的干扰信号。标准 频率是用来衡量其他频率的尺度。
[0093] 步骤S2具体包括:
[0094]S21:对所述的静止轨道卫星进行星历校准;
[0095]S22:对所述的静止轨道卫星进行星历数据转换;
[0096] S23:向所述的静止轨道卫星发送参考源信号。
[0097] 步骤S4中确定所述的干扰信号以及所述的参考源返回信号对应的理论多普勒频 移以及实际多普勒频移。下面介绍多普勒频移的确定原理。❖地球模型的建立
[0098] 众所周知,一个物体的位置和它的运动情况只能相对一个参考系来确定。适用惯 性定律(即牛顿第一定律)的参考系称为惯性参考系。
[0099] 在天体力学中,研宄天体运动,参考系必须是惯性的。从绝对意义上来说,任意地 心坐标系统都不是惯性的。因为地球绕太阳公转,而太阳又在做空间加速运动。但是,为了 研宄人造卫星运动,卫星实际上是被认为是地球携带着,地球的这一运动将不含卫星运动 所涉及的力学问题。因此,总是采用地心系统。然而,有为了表达卫星的空间位置,又要用 到站心坐标系。
[0100] 定义一个空间坐标系,应包含三个要素:坐标原点、参考平面(XY平面)和参考平 面上的主方向(X轴方向)。
[0101] 在该实施例中,采用了如图9所示的地固坐标系统,即认为地球是相对静止的。
[0102] 在卫星通信中,静止轨道卫星的轨位一般都是以经度来表示。本实施例中,以中星 10为实际干扰定位的卫星载体。中星10官方轨位位于东经110. 5度,根据此轨道位置,可 以根据地心坐标系转换成地固坐标系下的经炜度。设定坐标系的零经线为格林威治子午 线。
[0103] 设椭球长半轴为a,短半轴为b,变率倒数为1/f,那么
Figure CN103576166BD00121
[0105] 式中:v为维度识处的卯酉圈曲率半径,v=or/(l-e2sin2史)a5史和A分别为坐标点 的维度和经度,
[0106] h为相对椭球面的高度,
[0107]e为椭球第一偏心率,e2=(a2-b2)/a2=2f_f2.
[0108] 该实施例中,中星10的轨道经度为110. 5°,可以预测出该卫星在地固坐标系中 的位置为:
[0109]X=-14766144.0234
[0110]Y=39493846. 18746944
[0111]Z=0.00000
[0112] 单位:米
[0113] 接收站点的地理位置为:
[0114]经度:116. 25311
[0115]维度:39. 660222
[0116] 经过计算,接收站点在此坐标系下的坐标为:
[0117]X=-2172438.491219755
[0118]Y=4403331.893923063
[0119]Z=4070634. 2545
[0120] 根据上述模型,可以得出如图10所示的直观单星定位图形:
[0121] ❖卫星位置在坐标系中的表示
[0122] 建立上述坐标系后,需要根据卫星的实时位置及速度信息,计算出不同时刻的信 号经过卫星转发器后所产生的多普勒频移。因此,如果计算出卫星的实时坐标是决定定位 问题的基础。本实施例中采取星历(即卫星的实时信息)转换的方式获取精确星历数据。根 据参考网络定期提供的六根轨道根数,将其转换成地固坐标系下的六参数:位置X、Y、Z坐 标及速度Vx、Vy、Vz。可以计算出中星10 (E110. 5度)卫星的运动轨迹及速度变化轨迹。如 图12和图13所示。此处提及的参考网络和因特网不是一个概念,在需要的情况下,还将根 据现有的参考源网络(参考源分布)(即全国各地都有的参考源发射用于计算卫星的精确轨 道)。
[0123] 从图12可以看出,卫星在其轨迹附近做类似螺旋运动,运动幅度在数十公里以 内。卫星在其轨迹附近运动时,其速度在不同方向(X、Y、Z三个方向)是变化的,速度幅度 约为2-3米/秒。zhongxinglO卫星一周内的运动轨迹如图11所示。
[0124] ❖理论多普勒频移及变化率的计算
[0125] 在地固坐标系下,建立了地球和卫星的状态方程后,并根据卫星的实时轨道位置 参数,根据多普勒频移计算公式能够分别计算出上行和下行产生的多普勒频移:
Figure CN103576166BD00131
[0128]dopp1er_total=dopp1er_fu+dopp1er_fd
[0129] 其中dopp1er_fu和dopp1er_fd分别是上行多普勒频移和下行多普勒频移, doppler_total为卫星移动产生的总的多普勒频移。C为光速,大小为2. 99997*108,fu为上 行信号的中心频率,fd经卫星转发器转后下行信号的中心频率,;卫星的速度矢量,I为地 球张指向卫星的单位方向矢量。根据每一分钟卫星的星历数据,每一分钟可以计算出一个 多普勒频移。
[0130] ❖实际多普勒频移及变化率的测量
[0131] 在单星定位初期,采用单载波信号进行仿真,若为调制信号,则可以将其转化为单 载波信号后进行估计。目前单载波信号的估计算法较多且较成熟。本实施例中采用了时域 中四阶差分法。
[0132] 时域法的基本思想是对信号瞬时相位进行差分,来得到频率的估计值。由于一阶 差分会产生较大的误差,一般多采用多阶差分,利用其平滑滤波效应来降低差分噪声。比较 典型的是最优多阶差分滤波器,这里采用4阶差分滤波器,其表达式为:
[0133]y(k) =0. 0464A:+〇, 0846A2+〇. 1084A3+〇. 1148A4
[0134]其中Ak= x(n+k)-x(n_k),{x(n)}为原始信号相位序列。
[0135] 该方法相对计算量小,实现简单,但对各种噪声的影响非常敏感,将会降低信号的 估计方差。图13和图14为实际测量并通过四阶差分算法计算出的多普勒频移及变化率随 时间变化图。
[0136] 因此,步骤S4具体包括S41:根据所述的干扰信号以及所述的参考源返回信号确 定实际多普勒频移;
[0137]S42:根据所述的静止轨道卫星的轨道位置参数确定所述的静止轨道卫星的理论 多普勒频移;
[0138]S43:将所述的实际多普勒频移以及所述的理论多普勒频移进行同步。
[0139] 步骤S5具体包括:
[0140]S51:排除所述的静止轨道卫星的频率漂移引起的误差;
[0141]S52:排除所述的静止轨道卫星由于温度变化引起的误差;
[0142]S53:排除所述的静止轨道卫星由于星历引起的误差;
[0143]S54:排除由所述的干扰信号以及所述的参考源返回信号确定实际多普勒频移时 产生的误差;
[0144]S55:对所述的实际多普勒频移进行筛选。
[0145] 实际多普勒频移的误差主要包括:静止轨道卫星漂移引起的误差;星历准确性引 起的误差;测频算法引起的频率误差;接收端引起的频率误差。下面逐一介绍对各种误差 的排除方式。
[0146] a.静止轨道卫星漂移引起的误差,主要包括如下两种:卫星本身漂移引起的误差; 卫星随温度周期变化引起的误差。
[0147] 1.卫星本身漂移引起的误差
[0148] 本发明中,对卫星移动引起的多普勒频移的量级及大小要求非常高,在精确计算 有卫星运动引起的多普频差外,还需要剔除有本振本身引起的误差。可以通过两信道同步 的方式去除一部分。即在接收干扰源的同时,另一路信道同时采集由自参考源发射的标准 频率校准信号,并对该两路信号进行严格时间同步,计算下行接收站信号的频率对时间的 变化并同时计算出该时段的理论多普勒频移(仅由卫星运动引起的频率偏移),两个量做差 运算,即可得卫星本振引起的频率漂移。星上转发器频率随时间的变化曲线如图15所示。
[0149] 本振随温度周期变化引起的误差
[0150] 在卫星收发系统中,同步卫星和地面上的收发机之间的相对运动引起多普勒效 应,使得通信信号频率产生多普勒频移。多普勒频移信号中包含收发机地理坐标、运动速度 等重要信息。因此对该信号捕获、跟踪以及有效地分析会对无线电监测、定位起到至关重要 的作用。但是在宽温度范围估计多普勒频移的过程中,发射机以及接收机本振频率的稳定 度是测量误差的影响因素之一。由于晶体振荡器本身的不稳定性,特别是在温度变化时频 率漂移较大,制约了估计结果的准确度。这里以上行多普勒信号为例,表示为
[0151] fup(t)=fd(t) +A f (t)
[0152] 式中fd(t)表示卫星相对于地面发射端的多普勒频移理论值,A f(t)表示由温度 引起的频移偏差。可以看出如何消除温度引起的频率漂移是实现多普勒频移精确跟踪的关 键。通过大量实测与仿真实验,发现对于一个固定的发射端,国内多个监测站所获得的多普 勒频移的理论波形相差十分微小。因此在该实施例中将北京监测站的多普勒频移理论波形 作为参考,将其与实际接收波形做差,即
Figure CN103576166BD00151
[0154]然后对频偏Af(t)进行拟合得到&_(/)在采样多项式拟合时,不要选择过高的 阶数,通常选择4-5阶为宜。最后用接收波形减去拟合的结果
Figure CN103576166BD00152
[0156] 得到校准后的波形。此处校准的波形是用算法来剔除奇异点,并将频率值做均值 处理,使处理后的数据平坦。
[0157]b.星历准确性引起的误差
[0158] 为了精确计算多普勒频移,星历的准确性是基础。本发明中,一方面采用定期从网 站对星历进行更新,另一方面,采用测轨系统定期对卫星进行测控,及时更新星历。
[0159]c•测频算法引起的频率误差
[0160] 采用四阶最优差分滤波的方法作为频率估计方法。由于多普勒频率随时间较慢的 连续变化,且估计结果会存在随机噪声,因此必须对估计序列进行处理。
[0161] 考虑到多普勒频率最大变化率为1Hz/分钟,而定位所需精度至少为KT^KTHz的 量级。为最大程度避免估计不确定度的影响,将Is内的估计结果进行平均,由于在Is内的 最大Doppler频率变化为1(T2量级,可以认为其对定位结果不构成主要影响。
[0162] 尽管采用Is内估计结果进行平均的处理方法可以很大程度缓解估计随机噪声的 影响,但是处理后的结果仍存在一定分量的随机噪声。为便于准确获得多普勒变化率等参 数,需要采用进一步的方法进行处理来得到平滑的频率变化曲线。
[0163] 本实施例中采用基于多项式拟合的方法来进行,采用10阶的多项式进行拟合,可 以得到较正确的平滑曲线如图16和图17所示,由图16可以看出,实测值基本分布在拟合 值的两边,可以较好的描绘出真实值的变化趋势。拟合值可以较好的反映真实值的变化趋 势。在Is的数据下并经拟合处理后,其估计精度至少应可以保持在104抱数量级。
[0164]d.接收端引起的频率误差评估与消除
[0165] 任何硬件,由于热噪声、晶振的稳定度等影响,其频率响应特性都不可能恒定。因 此为了避免各个环节所引起的频率漂移引起的定位误差,在接收端采用了标准频率校准接 收机及采集设备的方式,保证接收端接收到的信号误差在1(T3的精度。
[0166] 本实施例中,核心内容之一就是精确估计某一信号经过单颗同步地球轨道卫星并 转发到地球站时产生的多普勒频移。由于地球运动速率仅为数米每秒,产生的多普勒频移 量级(Ku频段时)为150Hz左右。因此,在经度方向上,定位精度要在十公里以内,需要测频 精度为〇. 1Hz的量级。在维度方向上,定位精度要在十公里以内,需要测频精度为0.01Hz 的量级。因此,对于标准C频段(下行3. 4GHz-4. 2GHz、上行5. 625GHz-6. 425GHz)和标准Ku 频段(下行12. 25GHz-12. 75GH、下行14. 475GH-14. 975GHz),要保证维度下的0. 01Hz的经 炜度量级,则需要参考源信号发射装置本身带来的频率误差小于〇. 01Hz。基于此,在具体的 实施方式中,还可包括标准频率模块,其功能是为参考源信号发射装置提供一个频率误差 小于0. 01Hz标准频率参考。
[0167] 对于Ku频段,要满足0. 01Hz的频率精度,则需要引入标准参考频率的精度量级 为1.0E-12,对于C频段,要满足0. 1Hz的频率精度,则需要引入标准参考频率的精度为 5E-12。本发明中,引入了铷钟作为定位的标准频率参考。为了保证进一步的定位精度,可 以引入铯钟作为信号参考。
[0168] 某款铯钟(5071A)的技术参数如下:
[0169] 性能参数
[0170] 精度:±1.0E_12
[0171] 预热时间(典型值):30minutes
[0172] 重复性:±2.0E_13
[0173] 设置参数
[0174] 范围:±1.0E_9
[0175] 分辨率:1.0E-15
[0176] 稳定性
[0177] 不同时间(s)的稳定性
[0178] 1 < 1.2E-11
[0179] 10 < 8.5E-12
[0180] 100 < 2.7E-12
[0181] 1,000 < 8.5E-13
[0182] 10,000 < 2.7E-13
[0183] 100,000 < 8.5E-14
[0184] 长期 〈 5.0E-14
[0185] 相位噪声
[0186] Offset(Hz) 5MHz输出
[0187] 1 < -95dBc
[0188] 10 < -130dBc
[0189] 100 < -145dBc
[0190] 1,000 < -155dBc
[0191] 10,000 < -155dBc
[0192] 100,000 < -160dBc
[0193] 为了保证能够达到定位精度的要求,可通过下述的验证系统对频率测量精度进行 验证,验证平台的示意图如图5所示。发射频率为一个4GHz的单载波信号,该信号发生器采 用外部高精度参考源进行频率校准,并采用两台某型号的高速采集设备对信号进行采集, 采样速率为lOkbs,每一分钟存储一个文件,采样时间为2个小时。经过实际多普勒频移及 变化率的计算算法,计算出的两路信号的频率随时间变化如图6所示,两路信号的频率变 化的差值如下图7所示。
[0194] 从图6可以看出,在对干扰信号采集时,需要保证信号采集设备、标准频率参考等 需要高精度频率的设备有足够的时间预热。预热时间要保证半个小时以上。由于两台设备 固有本振的差别,计算出的两路信号其频率变化规律严格一致。
[0195] 图7给出的是对两路信号计算的频率相减,以得到信号频率随时间的变化大小。 由图可以看出,两路信号的频率差变化在0. 01Hz的量级,这一方面说明目前的采集设备具 有一定的采集精度,另一方面说明频率估计算法有较高的精度保证。
[0196] 步骤S6具体包括:
[0197]S61:对所述的实际多普勒频移进行滤波处理;
[0198]S62:从所述的滤波后的实际多普勒频移中提取有效数据;
[0199]S63:根据所述的有效数据确定与所述的干扰信息对应的干扰源的位置。
[0200] 首先对实际多普勒频移进行数据筛选。由于采集的下行信号的频率变化有以下几 个部分组成:上行频率漂移、实际多普勒频移、卫星本振频移、接收端设备频率漂移等。在这 四部分构成中,由于上行频率漂移是不可控的,其他几部分的频率变化是可以根据一定的 技术和算法进行消除。因此,为了保证从接收端的信号中精确得到卫星运动引起的多普勒 频移,且受上行频率漂移影响最小,则需要引入信号筛选机制和方法,剔除上行信号频率变 化较大的部分,并形成成型的算法。本发明中,提出了滑动均值滤波的方法,并根据多普勒 频移的峰值来生产相应的滤波参照函数,并采用最大似然估计的方法,求出多普勒频移最 好的一段,并与理论多普勒频移同步,并最终实现定位算法,进一步提升定位精度。实际多 普勒频移的数据筛选图如图18所示。
[0201] 基于粒子滤波的定位算法是近年来提出的一种适于非线性非高斯场合的滤波方 法。由于多普勒频移是一个随时间变化的量,但是由于除了固有多普勒频移变化之外,还有 很多因素都会引起频率随时间变化,需要用一种适合于该变化的滤波方法对频率变化进行 滤波。粒子滤波的核心思想就是利用大量随机样本及其相应权值毕竟所要求的后验概率分 布。从贝叶斯理论来看,跟踪问题就是在给定观测数据z1:k(下标l:k表示1到k时所有数 据)的情况下,递归计算当前的系统状态XK。粒子滤波就是要利用序贯蒙特卡洛的方法毕 竟状态后验分布P(XK|z1:k),即通过从后验概率分布采集带有权重的样本集,用 该样本集表不后验分布:
Figure CN103576166BD00171
[0203] 式中S(•)为狄拉克函数;N为所用粒子数目。
[0204] 本发明中,根据卫星直角坐标系和卫星各项参数,并利用粒子滤波的算法建立状 态方程和粒子,如下:
[0205]a.多普勒频移计算公式
[0206]k时刻卫星直角坐标xs=[xs,k,ys,k,zs,k]T,卫星速度u=[vx,k,vy,k,vz,k]T,接收端直 角坐标[Xnk,ynk,Zpk]T,目标经炜度为(民口),转换成直角坐标为xt= [xt,yt,zt]T, 其中
Figure CN103576166BD00172
[0208],且v=6378000为地球半径,h=35786000为轨位离地面高度。
[0209]设ds -t=[xs,k-xt,ys,k-yt,zs,k_zt]T,并对其归一化,有ds -t=ds -1/| |ds -1| |2。
[0210] 设ds -r=[xs,k-xr,ys,k-yr,zs,k_zr]T,并对其归一化,有ds -r=ds -r/| |ds -r| |2。
[0211] 计算上行多普勒频移,有
Figure CN103576166BD00181
其中心是卫星上行中心频率,C(l= 3e8
Figure CN103576166BD00182
表示光速。计算下行多普勒频移,有 其中fd(l是卫星下行中心频 J 率。贝1J总多普勒频移可表示为ft()tal=fup+fdOTn。
[0212] b.基于粒子滤波算法的单星定位
[0213] (1)运用粒子滤波算法对多普勒频移信号进行滤波;
[0214] (2)对经过滤波后的信号进行有效数据提取,并将其再次经过粒子滤波得到目标 经炜度的估计。
[0215] 具体算法如下:
[0216] 利用粒子滤波算法对单星定位问题进行建模,构建状态空间方程及观测方程。
[0217] 状态空间方程:
[0218] 由于干扰源位置固定不变,则有:
Figure CN103576166BD00183
[0220] 其中,<和<是均值为0,方差分别为<和< 的高斯白噪声。观测方程:
Figure CN103576166BD00184
[0222] 其中ftotal,k表示k时刻地面接收端实测多普勒频移,A(死爲)表示多普勒频移理 论值,其可表示为与干扰源经炜度有关的函数,%表示噪声,其中包括频率漂移及环境噪 声。
[0223] 产生粒子0
Figure CN103576166BD00185
,其中i= 1,2, . . .,N,N表示粒子 数量。
[0224] 更新粒子状态:
Figure CN103576166BD00186
[0226] 计算得到重要性权重
Figure CN103576166BD00187
[0228] 并归一化,对粒子进行重采样。将表示经炜度的粒子0和#带入函数h(□)得到 •.朽并对其进行加权平均,得到对多普勒频移的估计
[0229]
Figure CN103576166BD00191
[0230]接下来,需要对滤波器输出信号进行数据截取,选择多普勒频移变化率最大 的一段数据。对截取后的信号再次进行粒子滤波,这一次的目标是对干扰源的经炜度€和 <进行估计,同样需要产生粒子
Figure CN103576166BD00192
更新粒子状态、计算 重要性权重及粒子重采样,最后得到干扰源经炜度估计为:
Figure CN103576166BD00193
[0233] 本发明还提供了一种基于单颗静止轨道卫星对干扰源进行定位的系统,如图2所 示,未知发射源向单颗静止轨道卫星发送干扰信号,所述的定位系统具体包括:
[0234] 信号采集装置,用于采集静止轨道卫星上的干扰信号;
[0235] 参考源信号发射装置,用于向所述的静止轨道卫星发送参考源信号,并接收所述 的静止轨道卫星根据所述的参考源信号发送的参考源返回信号;
[0236] 多普勒频移确定装置,用于确定所述的干扰信号以及所述的参考源返回信号对应 的理论多普勒频移以及实际多普勒频移;
[0237] 信号处理装置,用于对所述的干扰信号以及所述的参考源返回信号的实际多普勒 频移进行误差排除;
[0238] 干扰源定位装置,用于根据所述的干扰信号以及所述的参考源返回信号确定干扰 源的位置。
[0239] 所述的参考源信号发射装置包括:
[0240] 标准频率信号发射模块,用于向所述的静止轨道卫星发送标准频率信号;
[0241] 标准频率信号接收模块,用于接收所述的静止轨道卫星根据所述的标准频率信号 发送的标准频率返回信号;
[0242] 所述的信号处理装置包括:
[0243] 第一确定模块,用于根据所述的标准频率返回信号确定其对应的频率;
[0244] 第二确定模块,用于确定所述的标准频率返回信号的多普勒频移;
[0245] 第三确定模块,用于根据所述的标准频率返回信号对应的频率以及多普勒频移确 定所述的静止轨道卫星的频率漂移。
[0246] 所述信号采集装置具体包括:
[0247] 第一设置模块,用于设置采样中心频率以及采样带宽;
[0248] 第二设置模块,用于设置采样时长以及采样数据存储方式;
[0249] 采集模块,用于采集静止轨道卫星上的干扰信号;
[0250] 所述的信号处理装置还包括:
[0251] 频率估计模块,用于对所述的干扰信号进行频率估计;
[0252] 筛选模块,用于对所述的干扰信号进行筛选。
[0253] 所述的系统还包括:
[0254] 星历校准装置,用于对所述的静止轨道卫星进行星历校准;
[0255] 星历转换装置,用于对所述的静止轨道卫星进行星历数据转换。
[0256] 所述的多普勒频移确定装置具体包括:
[0257] 实际多普勒频移确定模块,用于根据所述的干扰信号以及所述的参考源返回信号 确定实际多普勒频移;
[0258] 理论多普勒频移确定模块,用于根据所述的静止轨道卫星的轨道位置参数确定所 述的静止轨道卫星的理论多普勒频移;
[0259]同步模块,用于将所述的实际多普勒频移以及所述的理论多普勒频移进行同步。
[0260] 所述的信号处理装置包括:
[0261] 第一误差排除模块,用于排除所述的静止轨道卫星的频率漂移引起的误差;
[0262] 第二误差排除模块,用于排除所述的静止轨道卫星由于温度变化引起的误差;
[0263] 第三误差排除模块,用于排除所述的静止轨道卫星由于星历引起的误差;
[0264] 第四误差排除模块,用于排除由所述的干扰信号以及所述的参考源返回信号确定 实际多普勒频移时产生的误差;
[0265] 筛选模块,用于对所述的实际多普勒频移进行筛选。
[0266] 干扰源定位装置具体包括:
[0267] 滤波模块,用于对所述的实际多普勒频移进行滤波处理;
[0268] 提取模块,用于从所述的滤波后的实际多普勒频移中提取有效数据;定位模块,用 于根据所述的有效数据确定与所述的干扰信息对应的干扰源的位置。
[0269] 以下结合具体的实施例,详细介绍本发明的技术方案。在北京、成都、深圳三个 城市发射干扰源,并经过zhongxinglO(E110.5° )转发,在北京(经度116.25311,维度 39. 660222)对干扰源进行采集并发射参考源信号,最终定位,定位结果可靠,精度大约在 100公里左右。结合卫星上行台站分布及地面干扰快速查找的能力,能够有效对干扰进行确 定并锁定。
[0270] 1)北京站定位结果统计
[0271] 经过信号优化及筛选算法后,从多普勒频移信号中筛选2100点-2300点作为实际 多普勒频移进行定位。图20为上行多普勒频移理论与实测波形图;图21为经过滤波后的 普勒频移与理论多普勒频移图;图22为截取2100点至2300点的数据图;图23为北京站经 炜度定位结果图;图24为误差随噪声标准差变化曲线图。
[0272] 表1北京站定位误差
Figure CN103576166BD00211
[0274] 2)成都站数据结果
[0275] 经过信号优化及筛选算法后,从多普勒频移信号中筛选2200点-2400点作为实际 多普勒频移进行定位。图25为成都站多普勒频率滤波结果图;图26为成都站经炜度定位 结果图。
[0276] 表2成都站定位误差
Figure CN103576166BD00212
[0278] 3)深圳站数据结果
[0279] 经过信号优化及筛选算法后,从多普勒频移信号中筛选2200点-2400点作为实际 多普勒频移进行定位。图27为深圳站多普勒频率滤波结果图,图28为深圳站经炜度定位 结果图。
[0280] 表3深圳站定位误差
Figure CN103576166BD00221
[0282] 综上所述,本发明的有益成果是:提供了一种基于单颗静止轨道卫星对干扰源进 行定位的方法及系统,通过仅仅一颗静止轨道卫星即能实现卫星干扰源定位,大大地提升 了卫星干扰的查找能力,降低了干扰查找的诸多束缚条件,弥补了当前主流定位技术的不 足,填补了我国在卫星定位系统的空白,极大的促进了无线电管理水平的提高。
[0283] 本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例 的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员, 依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内 容不应理解为对本发明的限制。

Claims (15)

1. 一种基于单颗静止轨道卫星对干扰源进行定位的方法,其特征是,所述的方法包 括: 采集静止轨道卫星上的干扰信号; 向所述的静止轨道卫星发送参考源信号; 接收所述的静止轨道卫星根据所述的参考源信号发送的参考源返回信号; 确定所述的干扰信号以及所述的参考源返回信号对应的理论多普勒频移以及实际多 普勒频移; 对所述的干扰信号以及所述的参考源返回信号对应的实际多普勒频移进行误差排除; 以及 根据所述的干扰信号以及所述的参考源返回信号确定干扰源的位置; 其中,基于粒子滤波算法的单星定位包括: (1) 运用粒子滤波算法对多普勒频移信号进行滤波; (2) 对经过滤波后的信号进行有效数据提取,并将其再次经过粒子滤波得到目标经炜 度的估计; 具体算法如下: 利用粒子滤波算法对单星定位问题进行建模,构建状态空间方程及观测方程, 状态空间方程: 由于干扰源位置固定不变,则有:
Figure CN103576166BC00021
其中,< 和 < 是均值为〇,方差分别为4和 < 的高斯白噪声, 观测方程: .ftotal 'k=咐,%)+";. 其中表示k时刻地面接收端实测多普勒频移,,约)表示多普勒频移理论值, 其可表示为与干扰源经炜度有关的函数,Uk表示噪声,其中包括频率漂移及环境噪声; 产生粒子€ ~RUmax)和< ~"(Umax),其中i= 1,2, . . .,N,N表示粒子数量; 更新粒子状态: e'k =GU+4- ,(pX+ni 计算得到重要性权重:
Figure CN103576166BC00022
并归一化,对粒子进行重采样,将表示经炜度的粒子贫和S带入函数h( •)得到 并对其进行加权平均,得到对多普勒频移的估计
Figure CN103576166BC00023
接下来,需要对滤波器输出信号进行数据截取,选择多普勒频移变化率最大的一 段数据,对截取后的信号再次进行粒子滤波,这一次的目标是对干扰源的经炜度穴和#进 行估计,同样需要产生粒子€ ~VWminU和W,更新粒子状态、计算重要 性权重及粒子重采样,最后得到干扰源经炜度估计为:
Figure CN103576166BC00031
2. 根据权利要求1所述的方法,其特征是,所述方法还包括确定所述的静止轨道卫星 的频率漂移。
3. 根据权利要求2所述的方法,其特征是,确定所述的静止轨道卫星的频率漂移具体 包括: 向所述的静止轨道卫星发送标准频率信号; 接收所述的静止轨道卫星根据所述的标准频率信号发送的标准频率返回信号; 根据所述的标准频率返回信号确定其对应的频率; 确定所述的标准频率返回信号的多普勒频移; 根据所述的标准频率返回信号对应的频率以及多普勒频移确定所述的静止轨道卫星 的频率漂移。
4. 根据权利要求1所述的方法,其特征是,所述采集静止轨道卫星上的干扰信号具体 包括: 设置采样中心频率以及采样带宽; 设置采样时长以及采样数据存储方式; 采集静止轨道卫星上的干扰信号; 对所述的干扰信号进行频率估计; 对所述的干扰信号进行筛选。
5. 根据权利要求1所述的方法,其特征是,向所述的静止轨道卫星发送参考源信号具 体包括: 对所述的静止轨道卫星进行星历校准; 对所述的静止轨道卫星进行星历数据转换; 向所述的静止轨道卫星发送参考源信号。
6. 根据权利要求4所述的方法,其特征是,确定所述的干扰信号以及所述的参考源返 回信号对应的理论多普勒频移以及实际多普勒频移具体包括: 根据所述的干扰信号以及所述的参考源返回信号确定实际多普勒频移; 根据所述的静止轨道卫星的轨道位置参数确定所述的静止轨道卫星的理论多普勒频 移; 将所述的实际多普勒频移以及所述的理论多普勒频移进行同步。
7. 根据权利要求1所述的方法,其特征是,对所述的干扰信号以及所述的参考源返回 信号的实际多普勒频移进行误差排除具体包括: 排除所述的静止轨道卫星的频率漂移引起的误差; 排除所述的静止轨道卫星由于温度变化引起的误差; 排除所述的静止轨道卫星由于星历引起的误差; 排除由所述的干扰信号以及所述的参考源返回信号确定实际多普勒频移时产生的误 差; 对所述的实际多普勒频移进行筛选。
8. 根据权利要求1所述的方法,其特征是,根据所述的干扰信号以及所述的参考源返 回信号确定干扰源的位置具体包括: 对所述的实际多普勒频移进行滤波处理; 从所述的滤波后的实际多普勒频移中提取有效数据; 根据所述的有效数据确定与所述的干扰信息对应的干扰源的位置。
9. 一种基于单颗静止轨道卫星对干扰源进行定位的系统,其特征是,所述的系统包 括: 信号采集装置,用于采集静止轨道卫星上的干扰信号; 参考源信号发射装置,用于向所述的静止轨道卫星发送参考源信号,并接收所述的静 止轨道卫星根据所述的参考源信号发送的参考源返回信号; 多普勒频移确定装置,用于确定所述的干扰信号以及所述的参考源返回信号对应的理 论多普勒频移以及实际多普勒频移; 信号处理装置,用于对所述的干扰信号以及所述的参考源返回信号的实际多普勒频移 进行误差排除;以及 干扰源定位装置,用于根据所述的干扰信号以及所述的参考源返回信号确定干扰源的 位置; 其中,基于粒子滤波算法的单星定位包括: (1) 运用粒子滤波算法对多普勒频移信号进行滤波; (2) 对经过滤波后的信号进行有效数据提取,并将其再次经过粒子滤波得到目标经炜 度的估计; 具体算法如下: 利用粒子滤波算法对单星定位问题进行建模,构建状态空间方程及观测方程, 状态空间方程: 由于干扰源位置固定不变,则有:
Figure CN103576166BC00041
其中,< 和 < 是均值为〇,方差分别为W和 < 的高斯白噪声, 观测方程: ft〇tal,k=h{0^(P!,) +Uk 其中表示k时刻地面接收端实测多普勒频移,MA#/.I表示多普勒频移理论值, 其可表示为与干扰源经炜度有关的函数,Uk表示噪声,其中包括频率漂移及环境噪声; 产生粒子€,其中i=1,2,. . .,N,N表示粒子数量; 更新粒子状态: 0K+noka+nl 计算得到重要性权重:
Figure CN103576166BC00051
并归一化,对粒子进行重采样,将表示经炜度的粒子贫和4带入函数h(〇得到 并对其进行加权平均,得到对多普勒频移的估计
Figure CN103576166BC00052
接下来,需要对滤波器输出信号无进行数据截取,选择多普勒频移变化率最大的一 段数据,对截取后的信号再次进行粒子滤波,这一次的目标是对干扰源的经炜度贫和#进 行估计,同样需要产生粒子《和#~叭Pmin,化,更新粒子状态、计算重要 性权重及粒子重采样,最后得到干扰源经炜度估计为:
Figure CN103576166BC00053
10. 根据权利要求9所述的系统,其特征是,所述的参考源信号发射装置包括: 标准频率信号发射模块,用于向所述的静止轨道卫星发送标准频率信号; 标准频率信号接收模块,用于接收所述的静止轨道卫星根据所述的标准频率信号发送 的标准频率返回信号; 所述的信号处理装置包括: 第一确定模块,用于根据所述的标准频率返回信号确定其对应的频率; 第二确定模块,用于确定所述的标准频率返回信号的多普勒频移; 第三确定模块,用于根据所述的标准频率返回信号对应的频率以及多普勒频移确定所 述的静止轨道卫星的频率漂移。
11. 根据权利要求9所述的系统,其特征是,所述信号采集装置具体包括: 第一设置模块,用于设置采样中心频率以及采样带宽; 第二设置模块,用于设置采样时长以及采样数据存储方式; 采集模块,用于采集静止轨道卫星上的干扰信号; 所述的信号处理装置还包括: 频率估计模块,用于对所述的干扰信号进行频率估计; 第一筛选模块,用于对所述的干扰信号进行筛选。
12. 根据权利要求9所述的系统,其特征是,所述的系统还包括: 星历校准装置,用于对所述的静止轨道卫星进行星历校准; 星历转换装置,用于对所述的静止轨道卫星进行星历数据转换。
13. 根据权利要求10所述的系统,其特征是,所述的多普勒频移确定装置具体包括: 实际多普勒频移确定模块,用于根据所述的干扰信号以及所述的参考源返回信号确定 实际多普勒频移; 理论多普勒频移确定模块,用于根据所述的静止轨道卫星的轨道位置参数确定所述的 静止轨道卫星的理论多普勒频移; 同步模块,用于将所述的实际多普勒频移以及所述的理论多普勒频移进行同步。
14. 根据权利要求9所述的系统,其特征是,所述的信号处理装置包括: 第一误差排除模块,用于排除所述的静止轨道卫星的频率漂移引起的误差; 第二误差排除模块,用于排除所述的静止轨道卫星由于温度变化引起的误差; 第三误差排除模块,用于排除所述的静止轨道卫星由于星历引起的误差; 第四误差排除模块,用于排除由所述的干扰信号以及所述的参考源返回信号确定实际 多普勒频移时产生的误差; 第二筛选模块,用于对所述的实际多普勒频移进行筛选。
15. 根据权利要求9所述的系统,其特征是,干扰源定位装置具体包括: 滤波模块,用于对所述的实际多普勒频移进行滤波处理; 提取模块,用于从所述的滤波后的实际多普勒频移中提取有效数据; 定位模块,用于根据所述的有效数据确定与所述的干扰信息对应的干扰源的位置。
CN201210277476.0A 2012-08-06 2012-08-06 基于单颗静止轨道卫星对干扰源进行定位的方法及系统 Active CN103576166B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210277476.0A CN103576166B (zh) 2012-08-06 2012-08-06 基于单颗静止轨道卫星对干扰源进行定位的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210277476.0A CN103576166B (zh) 2012-08-06 2012-08-06 基于单颗静止轨道卫星对干扰源进行定位的方法及系统

Publications (2)

Publication Number Publication Date
CN103576166A CN103576166A (zh) 2014-02-12
CN103576166B true CN103576166B (zh) 2015-08-19

Family

ID=50048337

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210277476.0A Active CN103576166B (zh) 2012-08-06 2012-08-06 基于单颗静止轨道卫星对干扰源进行定位的方法及系统

Country Status (1)

Country Link
CN (1) CN103576166B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103885031B (zh) * 2014-04-08 2016-03-23 中国电子科技集团公司第五十四研究所 基于搜索优化滤波的运动卫星干扰源定位方法
CN104297723A (zh) * 2014-11-06 2015-01-21 中国电子科技集团公司第五十四研究所 基于搜索优化和非线性滤波的运动卫星干扰源定位方法
CN104635203B (zh) * 2015-02-12 2017-04-19 国家无线电监测中心 一种基于粒子滤波算法的无线电干扰源测向定位方法
CN106842256B (zh) * 2016-12-12 2019-03-12 西安空间无线电技术研究所 一种利用gnss单星信号的导航定位方法
CN107947849B (zh) * 2017-11-21 2019-10-29 中国电子科技集团公司第五十四研究所 一种多信关站同步方法
CN109946725A (zh) * 2017-12-20 2019-06-28 慧众行知科技(北京)有限公司 一种卫星定位方法及系统
CN110531315B (zh) * 2019-07-15 2020-08-21 南京航空航天大学 一种基于信号强度变化率的卫星干扰源直接定位方法及定位装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5088676B2 (ja) * 2007-07-31 2012-12-05 セイコーエプソン株式会社 時刻修正装置、時刻修正装置付き計時装置及び時刻修正方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
利用互模糊函数实现卫星干扰源定位;孙正波 等;《电波科学学报》;20041031;第19卷(第5期);第525-529页 *
卫星通信电磁干扰定位分析与干扰查找;李安平 等;《电波科学学报》;20101031;第25卷;第113-116页 *
通信卫星干扰定位的BRF神经网络方法;童新海 等;《通信学报》;20001130;第21卷(第11期);第35-40页 *

Also Published As

Publication number Publication date
CN103576166A (zh) 2014-02-12

Similar Documents

Publication Publication Date Title
CN103576166B (zh) 基于单颗静止轨道卫星对干扰源进行定位的方法及系统
US8456353B2 (en) Method and system for determining clock corrections
US10148345B2 (en) Accurate timing distribution by high-frequency radio
Ivanov et al. Chirp ionosonde and its application in the ionospheric research
CN102215558B (zh) 由通信广播卫星信号辅助的地面移动通信网定位方法
CN104270190A (zh) 基于电离层资料的同步自适应短波通信选频方法
Brown A physical meteor-burst propagation model and some significant results for communication system design
Kalaivaanan et al. Evaluation of ka-band rain attenuation for satellite communication in tropical regions through a measurement of multiple antenna sizes
CN103220774A (zh) 一种借助蜂窝网进行精密授时的方法及系统
CN103744082B (zh) 基于dmb信号的被动雷达水汽探测方法
Morton et al. Ionospheric effects, monitoring, and mitigation techniques
EP1955446A2 (en) Hf radio network frequency management
Hatsuda et al. Ku-band long distance site-diversity (SD) characteristics using new measuring system
Storek et al. Measurements of phase fluctuations for reliable MIMO space communications
Giesbrecht et al. An empirical study of the probability density function of HF noise
CN105425262B (zh) 一种实现卫星转发导航系统载波相位精密测量的方法
Yao et al. The NTSC VLBI System and its application in UT1 measurement
Storek et al. Interferometer for measurements of the MIMO satellite channel at Ku-band
Li et al. Over-the-horizon time and frequency synchronization for maneuverable radar system based on troposcatter
Luini et al. Joint effects of clouds and rain on Ka-band earth observation data downlink systems
Boithias et al. Propagation due to tropospheric inhomogeneities
CN111866754A (zh) 一种无线广播授时信息处理方法
Sun et al. Analysis of two way satellite time and frequency transfer system
Wilk-Jakubowski Sun transits in geo satellite systems in the aspect of radio waves propagation
Xie et al. Selection of Communication Frequency Points for Shortwave Radio Network Based on VOACAP

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
C14 Grant of patent or utility model