专利名称：一种基于实时卫星信号模拟器同步增减可见卫星装置的制作方法
技术领域：
本实用新型涉及一种卫星信号模拟器，尤其涉及一种基于实时卫星信号模拟器同 步增减可见卫星装置。它属于卫星定位导航技术领域。
背景技术：
在卫星导航定位系统中，导航接收机安装在载体上，其接收的信号含有载体动态 引起的多普勒频移，同时信号经过空间传播，不可避免地受到各种误差源的影响，信号在接 收时刻的状态已不同于发射时刻状态，这种差别与载体位置、动态特性、测量时间、应用环 境等有关。信号模拟器作为卫星导航系统和各种接收设备研制的关键环节和测试工具，能 够在卫星不在轨的条件下为导航系统中的地面系统研制、建设和测试提供一个可靠、稳定、 准确和易用的仿真和测试环境。卫星信号模拟器能够灵活模拟各类测试环境，检验载体武 器和地面系统的功能和性能指标是否满足研制要求、信息流程及处理算法是否合理、各模 块之间接口关系是否正确，从而为系统验收及实弹实验提供依据，同时对整个系统的能力 进行评估。在卫星信号实时模拟过程中，实时可见星判断和增减可见星操作是一项关键技 术难点。载体高动态环境下的可见星判断通过计算卫星与接收机之间的仰角和方位角实 现。通过星历中轨道参数和卫星初始参数计算得到卫星的位置，再根据接收机的运动轨迹 计算出二者的相对位置，利用卫星与接收机之间的仰角计算进行可见星判断。信号模拟器的可见星计算原理是以载体接收机天线位置为原点，建立相对坐标 系，计算卫星在该坐标系中的位置，从而计算出卫星相对于接收机天线的仰角，进而判断卫 星是否可见。这部分的步骤如下1.以导航接收机天线位置为原点，XoY平面与地球相切，Y轴指向正北，X轴在XoY 平面内与Y轴垂直指向东，Z轴与XoY平面垂直指向地球球体外侧。计算出卫星在该坐标 系中的位置(X，Y，z)。2.计算出卫星相对于接收机天线的仰角EL。
ζEL = arctan(-p= ) · 180 / 苁
^X2+Y23.由仰角判断该颗卫星是否可见。设立可见星仰角门限ELmin，如果EL彡ELmin，则该颗卫星可见，否则该颗卫星不可 见。由于当卫星仰角十分接近0°时，由电离层引起的测距误差会急剧增大，所以一般设置 ELmin在5 °左右。当模拟器进行真实卫星信号的长时间模拟时，需要根据载体运动状态实时计算卫 星与接收机的相对位置，并进行可见星判断。在仿真过程中必然会出现卫星增加或减少的 现象，在此情况下，只有保证新出现卫星或即将消失卫星信号模拟通道与其他正常信号模 拟通道时序完全同步且位置计算精确，才能使得整个模拟器系统工作正常，否则将出现信 号模拟误差和信号阶跃现象，会极大地影响模拟器的模拟精度和工作连续性。
4[0011]基于以上特点，本专利发明了一种具有实时信号模拟功能的同步增减可见卫星装 置及其构建的方法，实现高动态环境下卫星信号的长时间正确模拟。
实用新型内容1、目的本实用新型的目的在于提供一种基于实时卫星信号模拟器同步增减可见 卫星装置，以解决现有技术中的不足。2、技术方案本实用新型一种基于实时卫星信号模拟器同步增减可见卫星装置， 它主要可以分为实时可见星判断，模拟信息交互及通道同步开关状态控制三个部分。其中， 实时可见星判断部分主要实现模拟器工作时间内的可见卫星计算通过初始星历参数计算 卫星位置，根据在轨卫星和载体接收机的运动状态实时计算卫星和接收机的相对位置，计 算卫星仰角和方位角进行可见星判断，计算每个模拟通道中载波NCO(数字控制振荡器) 和码NCO的频率和相位控制字，实现模拟器运行过程中的实时加减星逻辑。模拟信息交互 部分主要实现仿真数据在数字信号处理器即DSP和现场可编程门阵列即FPGA�？橹涞� 交互通过总线接口驱动逻辑将载波NCO和码NCO的频率控制字和相位控制字同步传递到 FPGA的NCO计算模块中，并在系统仿真过程中，通过中断信号对FPGA计算�？榻卸� NCO控制字的实时更新；通道同步开关状态控制部分主要实现信号模拟器各通道之间的信 号同步和通道开关操作通过FPGA内部的同步脉冲信号，将DSP�？榉⑺偷募扑憧刂谱趾� 初始信息进行同步处理，确保各个信号调制通道的时序一致，针对模拟过程中不可避免的 出现某颗卫星中途出现或消失(即仰角低于可见星判断门限)的现象，在DSP和FPGA模 块之间建立各个模拟通道的同步开关逻辑，通过同步状态机的处理实现任意模拟通道的开 启、关闭和信息注入，并确保该通道与其他通道的时序同步。本实用新型采用以DSP+FPGA�？槲诵牡娜砑尴叩缣逑导芄梗ü凳笨杉� 星计算配合同步控制逻辑，使卫星信号模拟器能够在独立运行过程中实时改变模拟通道状 态，实现任意通道同步加星和去星操作，达到了卫星信号模拟器长时间独立工作的效果，具 有信号模拟精度高，长时间模拟无阶跃误差、算法复杂度低等优势，对于针对高动态环境下 长时间卫星信号模拟测试具有很高的创新型和实用意义。本实用新型一种基于实时卫星信号模拟器同步增减可见卫星装置，主要包括上位 机通信�？�、信息处理�？椤⑿藕糯砟？楹妥芟咝畔⒔换ツ？椋浣峁雇既缤�1所示。它们之间的连接逻辑关系是上位机通信�？橥ü齊S232串行接口与信息处理模 块连接，信息处理模块通过数字信号处理器即DSP和现场可编程门阵列即FPGA芯片之间的 总线接口逻辑与总线信息交互模块连接，信号处理�？橥ü鼶SP和FPGA芯片之间的总线接 口逻辑与总线信息交互�？榱�。所述上位机通信模块是在PC机上构建的系统显控界面软件，它的功能是进行模 拟器系统初始化操作，设置模拟初始参数和控制信息(包括模拟用户轨迹设定、星历文件 选取、初始用户位置和时刻设置等)，并通过串口逻辑将上述信息下发到信息处理�？椋� 时在工作过程中从信息处理�？橹薪邮漳Ｄ庑畔�(包括可见卫星信息，各通道模拟状态信 息等)，在显控界面中进行显示。该上位机通信�？椴捎孟钟屑际跏迪�。所述信息处理�？樵贒SP芯片中实现，它的功能是接收上位机通信�？榉⑺偷牟� 数和指令进行系统初始化和数据信息处理。信息处理�？橹饕晌佬俏恢糜朐靥褰邮栈�
5置实时计算子�？�、同步增减可见卫星实时判断子�？楹屯ǖ滥Ｄ庑畔⒓扑阕幽？楣钩�。所述卫星位置与载体接收机位置实时计算子�？椋涔δ苁歉菸佬切抢问� 载体运动轨迹参数实时计算所模拟卫星的位置和载体接收机的位置。该子�？榈慕峁褂尚� 历解析结构体、卫星位置计算函数和接收机轨迹计算函数构成，其算法原理采用现有技术 实现。所述同步增减可见卫星实时判断子�？榈慕峁褂晌佬茄鼋羌扑愫⒖杉佬橇� 表更新函数和通道状态设置函数构成。卫星仰角计算函数根据卫星位置与载体接收机位置 实时计算子�？槭涑龅奈佬俏恢煤徒邮栈恢檬菁扑忝靠盼佬窍喽杂谠靥褰邮栈难� 角信息，实时进行可见星座判断，将判决结果传递至可见星列表更新函数；可见星列表更新 函数读取此卫星在该计算周期和上一计算周期内的可见性判决结果后可得到该颗卫星属 于何种状态(卫星状态包括由不可见到可见、持续可见、由可见到不可见及持续不可见)， 以此状态更新可见卫星列表，并将此状态传递至通道状态设置函数；通道状态设置函数根 据卫星状态进行模拟通道的状态设置(模拟通道状态包括开启、持续模拟、关闭和持续关 闭)。所述通道模拟信息计算子�？榈墓δ苁歉萃皆黾蹩杉佬鞘凳迸卸献幽？� 建立的卫星列表信息，计算所有可见卫星的伪码/载波控制字信息、传输延迟、多普勒频移 及导航电文信息，将上述模拟信息通过总线信息交互�？榇涞叫藕糯砟？椋苯� 分模拟信息回传上位机通信�？榻邢允�。该子�？榈慕峁褂晌本嗉扑愫�、电离层/对 流层延迟计算函数、伪码/载波频率字计算函数和导航电文组帧函数构成，起算法原理采 用现有技术实现。所述信号处理�？樵贔PGA芯片中实现，它的功能是实现各模拟通道的状态控制 和信号处理，完成多通道卫星信号调制。信号处理�？橹饕赏ǖ雷刺娇刂谱幽？楹� 信号调制子�？楣钩�。所述通道状态同步控制子�？榈慕峁褂赏ǖ揽乜刂谱刺钩桑渥刺坡� 辑为停止状态通道保持关闭，相关数据寄存器复位(保持该状态)运行状态通道持续运行，相关数据寄存器实时更新(保持该状态)启动状态有新的可见卫星出现，该通道即将模拟该颗卫星信号(转移至运行状 态)关闭状态该通道目前模拟的卫星即将消失，通道即将关闭(转移至停止状态)通道状态同步控制子�？榈墓δ苁窍煊π畔⒋砟？榉⑺偷耐ǖ揽刂菩畔⒔� 模拟通道同步开关控制逻辑，实现可见卫星情况改变时模拟通道的准确变化，同时通过总 线信息交互�？槎酝ǖ滥Ｄ庑畔⒔刑崛『突捍妫瓿捎谛畔⒋砟？橹涞氖萁换ィ凰鲂藕诺髦仆ǖ赖墓δ苁峭ü隢CO和载波NCO设计实现动态伪码和载波的生 成，完成导航电文数据的扩频调制，生成模拟器发射信号。该子�？榈慕峁褂晌甭敕⑸髂� 块、码NCO�？�、载波NCO�？楹偷髦颇？楣钩桑渌惴ㄔ聿捎孟钟屑际跏迪�。所述总线信息交互�？橛勺芟咝畔⒔换ツ？橛蒁SP读数逻辑、DSP写数逻辑和中 断函数构成。通过DSP读数逻辑实现数据信息由信号处理�？榈叫畔⒋砟？榈拇莺突� 存，通过DSP写数逻辑实现数据信息由信息处理�？榈叫藕糯砟？榈拇莺突捍妫ü卸虾低持卸闲藕藕屯铰龀澹繁？杉佬乔榭龊湍Ｄ庑畔⒌氖凳备�。本实用新型一种基于实时卫星信号模拟器同步增减可见卫星装置的构建法，该方 法具体步骤如下步骤一系统初始化，DSP通过上位机通信�？榻邮粘跏寄Ｄ獠问�、星历参数和载 体接收机初始信息等。其中模拟卫星的星历参数通过读取星历文件获得，接收机初始信息 通过上位机设置，包括接收机初始位置、初始速度和初始加速度。根据星历信息计算每一颗 卫星的初始位置，对卫星和接收机的初始位置信息进行存储。根据GPS星历提供的轨道参数计算卫星位置的步骤如下所示，其中t是DSP系统 每次计算位置的时刻(1)根据星历数据参考时刻t。e计算归一化时间tktk = t-toe若 tk > 302400，tk = tk_604800若 tk < -302400，tk = tk+604800(2)计算卫星运行的平均角速度n0no = ^Igm / Ai = VGM /(Vl")3利用电文中平均角速度的修正量δ n，求得卫星运行的平均角速度为η = η0+Δη(3)计算观测时刻的卫星平近点角MkMk = M。+ntk式中，Mtl为电文中给出的相应卫星在参考时刻t。e平近点角。(4)计算偏近点角Ek和相对论效应误差REk = Mk+e sin Ek为加快收敛速度，采用微分迭代法解算Ek。经过8次迭代即可获得IEw-EiI < ε =1. OXlO-15的精度级别。
P _ J7 Mk-Ek+e·sinEk^ - K + !_e.cos^ E0 = Mk
>R = f -e-4A-smEK其中f 是 WGS 84 规定的校正参数，f = -4. 442809305 X 1(Γ10。(5)计算真近点角vk
r一Vl-e2 Sin^i.vk = arctan-k-
CosEk -e(6)计算升交距角ΦκΦκ = vk+ωω为卫星电文给出的近地点角距。(7)计算摄动改正项K、δ ” δ i 'Su=Cuccos(20k) + Cussm(20k)I Sr = Crc cos(20,) + Crx sin(20,) ^=C,ccos(20,) + C sin(20J
7[0058] 式中，5u> 5r> Si分别为因地球非球形和日月引力因素而引起的升交距角 摄动量，卫星失径r的摄动量和轨道倾角i的摄动量。
变化量
(8)计算经过摄动改正的升交距角uk、卫星矢径rk和轨道倾角ik rk = A(l -ecosEk) + 久.
(9)计算升交距角Φκ、经过摄动改正的升交距角uk、卫星矢径rk和轨道倾角ik的
,=Φ,+2.0·Φ,· (Cm . cos(2.0-OJ-Cuc. sin(2.0 ·Φ,)) . A-n-e- ην^
(10)计算卫星在轨道平面的位置
(11)计算卫星在轨道平面的速度
(12)计算观测时刻的升交点经度Qk
式中，为地球自转速率，Ω Pi^ttre均是卫星电文参数t
(13)计算卫星在WGS84坐标系中的位置
(14)计算卫星在WGS 84坐标系中的速度
步骤二 遍历计算每颗卫星的仰角和方位角信息，并由此判断该卫星是否可见，建 立可见卫星列表和空闲通道列表。 可见卫星计算原理是以载体接收机天线位置为原点，建立相对坐标系，计算卫星 在该坐标系中的位置，从而计算出卫星相对于接收机天线的仰角，进而判断卫星是否可见。
这部分的计算步骤如下1.以导航接收机天线位置为原点，XoY平面与地球相切，Y轴指向正北，X轴在XoY 平面内与Y轴垂直指向东，Z轴与XoY平面垂直指向地球球体外侧。计算出卫星在该坐标 系中的位置(X，Y，Z)。2.计算出卫星相对于接收机天线的仰角EL。 3.由仰角判断该颗卫星是否可见。设立可见星仰角门限ELmin，如果EL彡ELmin，则该颗卫星可见，否则该颗卫星不可 见。由于当卫星仰角十分接近0°时，由电离层引起的测距误差会急剧增大，所以一般设置 ELmin在5 °左右。步骤三信息处理�？橥ü芟咝畔⒔换ツ？橄蛐藕糯砟？榉⑺屯ǖ雷刺� 息，完成对信号处理�？橹懈鞲瞿Ｄ馔ǖ赖某跏蓟�。每个通道有4种状态，“00”表示通道关 闭，“01”表示通道正常工作，“10”表示通道即将启动，“11”表示通道即将关闭。初始状态 下，将所有分配了可见星的通道的对应状态置为“10”，其余的空闲通道状态置为“00”。初 始模拟信息包括初始载波相位、初始码片相位和导航电文数据等信息。步骤四在模拟器系统工作过程中，根据卫星号实时计算每一颗星的位置(具体 方法见步骤一)。根据载体运动模型实时计算载体接收机位置。根据二者相对位置判断该 颗星是否可见(具体方法见步骤二)。当可见卫星出现增加或减少状况时，同步改变相应通 道状态，并更新可见卫星列表信息。步骤五以中断脉冲信号作为触发，实时计算可见卫星的码NCO控制字和载波NCO 控制字，同时计算每个通道的初始码相位、初始载波相位和导航电文等通道数据信息。1.针对每一颗可见卫星，分别计算该卫星的传输环境延迟，主要包括调用对流层 延迟计Δτκ。ρ和电离层延迟Δ_计算。2.计算伪距 式中，ρ C0DE为对应的码伪距，P f为对应的载波伪距，XSV、XSV、ZSV为卫星的位置，Xu、 yu、Zu为接收机的位置，c为光速，ΔΜ为星钟误差，A1-为电离层误差，Atot为对流层误 差。3.计算相邻采样时刻的伪距变化量 式中，Δ ■、Δ CAEEIEE为码和载波对应的伪距变化量，P C0DE_0LD> P caeeiee_OLD为缓存的 上一个采样时刻码和载波对应的伪距。4.计算伪距变化量对应的偏移频率字 式中，Δ Ncode, Δ Ncaeeiee为计算出的码和载波的偏移控制字，fEF为GPS射频频率，fs 为采样频率，K为NCO深度，最后的xlOOO/m是为了对应ms计算周期，m为每一个计算周期 的ms数。5.将基准频率字和偏移频率字相加得到发送给FPGA的频率字。NCODE = Ncodo+ Δ NCODENcaeeiee 一 Ncaeeiee_b+ Δ NCAEEIEE式中，NroDE—Β和NfflmEK—Β分别为伪码和载波的基准频率。6.导航电文计算电文在模拟器中为初始计算一次，运行中每一主帧周期长度计算一次。初始计算 和运行时计算方法一样，都是根据收到的每颗星的星历参数和模拟时刻，对每一颗可见星 每一个子帧的电文逐个生成。步骤六响应中断脉冲信号，将步骤四一五中计算的信息数据(主要包括通道状 态，NCO控制字，导航电文)通过总线信息交互�？榉⑺椭列藕糯砟？椤ＤＤ馄髡Ｔ诵� 时步骤四和步骤五分别按照设计的周期循环计算。步骤四中所计算卫星号从1号卫星开始 并在本次计算结束后将卫星号加1，如果当前卫星号达到最大卫星号，则将卫星号置1重新 遍历计算。步骤五中所计算的频率控制字根据可见卫星列表状态，每次计算周期对所有状 态为“01”和“ 10”的通道进行计算。步骤一-六为模拟器信息处理部分，其流程原理如图 2所示。步骤七通过复位信号对信号处理�？橹械男藕拍Ｄ馔ǖ澜卸潦拇嫫鞯母次� 和数据锁存。复位逻辑完成后产生中断脉冲信号，并驱动总线接口�？椋畔⒋砟？榉� 送来的各种控制信息和参数进行锁存。具体方法为在FPGA(现场可编程门阵列)中判断总 线接口逻辑中的读数信号，当写片选信号为低时，在DSP(数字信号处理器)写使能信号上 升沿处进行触发，通过地址总线进行数据地址译码，通过数据总线将DSP发送来的各种控 制信息和参数进行锁存，DSP写数逻辑如图5所示。步骤八信号处理�？樵谥卸闲藕诺目刂葡拢凳倍寥⊥ǖ雷刺畔ⅲ礁赂� 个模拟通道的通道状态，同时将各通道的初始模拟信息置入相应寄存器。步骤九在通道同步脉冲的触发下，将通道状态控制字和需要回传信息处理�？� 的实时数据进行提取后赋与总线信息交互�？橹邢嘤Φ氖菟婕拇嫫�。具体方法为在 FPGA中判断总线接口逻辑中的写数信号，利用地址总线的地址译码将需要向信息处理�？� 发送的各种实时信息和参数所存在DSP读数据总线锁存寄存器中，当读片选信号和DSP读 使能信号同为低电平出进行触发，将锁存器中数据通过数据总线发送，将实时改变的通道 状态及仿真过程中产生的有效数据回传至信息处理模块。DSP读数逻辑如图6所示。步骤十通过同步通道开关控制状态机逻辑实现信号处理�？橹械耐ǖ雷刺� 改变，根据步骤八中描述的通道控制逻辑，将步骤七中锁存的通道状态控制信息置入同步 通道开关控制状态机，实现通道状态同步改变。将步骤七中锁存的模拟数据同步置入FPGA 相应通道。对于每个工作通道，在每次中断信号来临时重复执行步骤五-十，对通道内的模 拟参数进行同步更新，根据接收的控制字信息实时生成具有高动态多普勒特性的载波和伪 码，并完成导航电文数据调制，实现高动态信号模拟器的实时工作。步骤七_十为模拟器信 号处理部分，其流程原理如图3所示。
10[0109]所属步骤四中，实时可见星位置计算和可见星列表的处理流程详见附图2，属于本 实用新型创新保护点。其具体实现方法如下模拟器工作过程中，根据导航系统的卫星编 号Is为间隔循环遍历计算每一颗星的位置，同时根据载体运动模型实时计算载体接收机 位置，根据二者相对位置判断该颗星是否可见。如果该颗卫星可见，则查看可见卫星表，若 可见卫星表中已经有该颗卫星存在，则维持该颗卫星对应的通道状态为“01”不变，若可见 卫星表中没有该颗卫星，表示该卫星是即将可见的，则查看空闲通道列表，若有空闲通道， 则将该通道状态标记为“10”，并在空闲通道表中将该通道删除，若没有空闲通道，则维持现 有状态不变。如果该颗卫星不可见，则查看可见卫星表，若可见卫星表中没有该颗卫星，则 维持现有状态不变，若可见卫星表中存在该颗卫星，则表示该颗卫星即将不可见，将对应通 道状态标记为“11”，同时在可见卫星列表中删除该颗卫星。所属步骤八中，同步通道状态控制逻辑详见附图四，属于本实用新型创新保护点。 其具体实现方法如下信号处理模块在中断信号的控制下读取通道状态信息，并以本地同 步脉冲为触发更新各个模拟通道的通道状态。对于即将开始模拟的通道，则将对应的初始 模拟信息置入相应寄存器组。每个通道有4种状态，其状态转移逻辑为状态“00”:通道保持关闭，相关数据寄存器复位(保持状态1)状态“01”通道持续运行，相关数据寄存器实时更新(保持状态2)状态“10” 有新的可见卫星出现，该通道即将模拟该颗卫星信号(瞬时状态1)状态“11”该通道目前模拟的卫星即将消失，通道即将关闭(瞬时状态2)当模拟通道的状态为“00”时，表示该通道处于关闭状态，不对其内部的逻辑进行 处理，该通道不产生调制信号；当模拟通道的状态为“01”时，表示该通道处于工作状态，其 内部处理逻辑在每个系统中断时刻更新相关数据寄存器信息，持续产生调制信号；当模拟 通道的状态为“10”时，表示该通道即将开始工作，进入通道开启状态的跳转；当该模拟通 道的状态为“11”时，表示该通道即将停止工作，进入通道关闭状态的跳转。上述通道状态控制逻辑能够保证在模拟器出现可见卫星的增加或减少情况时，即 将开启或关闭的通道和正在运行的各通道严格保持同步，从而使得模拟器模拟的可见星星 座变化与真实情况一致，不会产生影响接收机的定位解算处理的误差。本实用新型一种基于实时卫星信号模拟器同步增减可见卫星装置及其构建方法 的优点在于(1)实时可见星判断，根据载体接收机的高动态运动轨�：臀佬窃硕旒Ｊ凳奔� 算卫星仰角和方位角，实时更新各通道模拟的卫星信号，保证了信号模拟器长时间的连续 正常工作；(2)高效率的同步信息交互逻辑可使信号模拟器的硬件逻辑资源充分利用，实现 对高动态载体轨�：投嗥绽掌狄频母呔饶Ｄ猓岣呶本喾直媛示龋�(3)采用通道同步开关状态控制，能够实现通道加减星操作和信息交互的时序同 步，保证在长时间信号模拟过程中的通道一致性和模拟误差消除。
图1所示为本实用新型实时卫星信号模拟器同步增减可见卫星装置的结构图图2所示为本实用新型信息处理�？楣ぷ髁鞒掏糩0123]图3所示为本实用新型信号处理�？楣ぷ髁鞒掏纪�4所示为同步通道开关控制状态机结构图图5所示为异步总线接口 DSP读逻辑示意图图6所示为异步总线接口 DSP写逻辑示意图图7所示为模拟器原始数据和商业接收机的定位数据比对结果示意图(纬度)图8所示为模拟器原始数据和商业接收机的定位数据比对结果示意图(径度)图9所示为模拟器原始数据和商业接收机的定位数据比对结果示意图(高度)
具体实施方式
下面将结合附图和实例对本实用新型作进一步的详细说明。 本实用新型应用实例基于的硬件平台选用Xi 1 inx公司生产的V4SX55FPGA芯片和 TI公司生产的C6713DSP芯片，系统工作时钟68MHz，模拟12通道GPS信号，动态环境为高 动态。本实用新型一种基于实时卫星信号模拟器同步增减可见卫星装置主要包括上位 机通信�？椋畔⒋砟？椋藕糯砟？楹妥芟咝畔⒔换ツ？椋湓斫峁雇既绺酵�1所
7J\ ο所述上位机通信�？槭窃赑C机上构建的系统显控界面软件，它的功能是进行模 拟器系统初始化操作，设置模拟初始参数和控制信息(包括模拟用户轨迹设定、星历文件 选取、初始用户位置和时刻设置等)，并通过串口逻辑将上述信息下发到信息处理�？椋� 时在工作过程中从信息处理�？橹薪邮漳Ｄ庑畔�(包括可见卫星信息，各通道模拟状态信 息等)，在显控界面中进行显示。所述信息处理�？樵贒SP C6713芯片中实现，它的功能是接收上位机通信�？榉� 送的参数和指令进行系统初始化和数据信息处理。信息处理�？橹饕晌佬俏恢糜朐靥褰� 收机位置实时计算子�？椤⑼皆黾蹩杉佬鞘凳迸卸献幽？楹屯ǖ滥Ｄ庑畔⒓扑阕幽？� 构成。所述卫星位置与载体接收机位置实时计算子�？椋涔δ苁歉菸佬切抢问� 载体运动轨迹参数实时计算所模拟卫星的位置和载体接收机的位置。该子模块的结构由星 历解析结构体、卫星位置计算函数和接收机轨迹计算函数构成，其算法原理采用现有技术 实现。所述同步增减可见卫星实时判断子�？榈慕峁褂晌佬茄鼋羌扑愫⒖杉佬橇� 表更新函数和通道状态设置函数构成。卫星仰角计算函数根据卫星位置与载体接收机位置 实时计算子�？槭涑龅奈佬俏恢煤徒邮栈恢檬菁扑忝靠盼佬窍喽杂谠靥褰邮栈难� 角信息，实时进行可见星座判断，将判决结果传递至可见星列表更新函数；可见星列表更新 函数读取此卫星在该计算周期和上一计算周期内的可见性判决结果后可得到该颗卫星属 于何种状态(卫星状态包括由不可见到可见、持续可见、由可见到不可见及持续不可见)， 以此状态更新可见卫星列表，并将此状态传递至通道状态设置函数；通道状态设置函数根 据卫星状态进行模拟通道的状态设置(模拟通道状态包括开启、持续模拟、关闭和持续关 闭)。所述通道模拟信息计算子�？榈墓δ苁歉萃皆黾蹩杉佬鞘凳迸卸献幽？�
12建立的卫星列表信息，计算所有可见卫星的伪码/载波控制字信息、传输延迟、多普勒频移 及导航电文信息，将上述模拟信息通过总线信息交互�？榇涞叫藕糯砟？椋苯� 分模拟信息回传上位机通信�？榻邢允�。该子�？榈慕峁褂晌本嗉扑愫�、电离层/对 流层延迟计算函数、伪码/载波频率字计算函数和GPS Ll频点导航电文组帧函数构成，起 算法原理采用现有技术实现。所述信号处理�？樵贔PGA芯片中实现，它的功能是实现12个同步GPS模拟通道 的状态控制和信号处理，完成多通道卫星信号调制。信号处理�？橹饕赏ǖ雷刺娇� 制子模块和信号调制子�？楣钩�。所述通道状态同步控制子�？榈慕峁褂赏ǖ揽乜刂谱刺钩桑渥刺坡� 辑为停止状态通道保持关闭，相关数据寄存器复位(保持该状态)运行状态通道持续运行，相关数据寄存器实时更新(保持该状态)启动状态有新的可见卫星出现，该通道即将模拟该颗卫星信号(转移至运行状 态)关闭状态该通道目前模拟的卫星即将消失，通道即将关闭(转移至停止状态)通道状态同步控制子�？榈墓δ苁窍煊π畔⒋砟？榉⑺偷耐ǖ揽刂菩畔⒔� 模拟通道同步开关控制逻辑，实现可见卫星情况改变时模拟通道的准确变化，同时通过总 线信息交互�？槎酝ǖ滥Ｄ庑畔⒔刑崛『突捍妫瓿捎谛畔⒋砟？橹涞氖萁换ィ凰鲂藕诺髦仆ǖ赖墓δ苁峭ü隢CO和载波NCO设计实现动态C/A码和数字载 波生成，完成GPS Ll频点导航电文数据的扩频调制，生成模拟器发射信号。该子模块的结 构由伪码发生器�？�、码NCO模块、载波NCO�？楹偷髦颇？楣钩桑渌惴ㄔ聿捎孟钟屑� 术实现。所述总线信息交互�？橛勺芟咝畔⒔换ツ？橛蒁SP读数逻辑、DSP写数逻辑和中 断函数构成。通过DSP读数逻辑实现数据信息由信号处理�？榈叫畔⒋砟？榈拇莺突� 存，通过DSP写数逻辑实现数据信息由信息处理�？榈叫藕糯砟？榈拇莺突捍妫ü� 中断函数产生系统中断信号和同步脉冲，确�？杉佬乔榭龊湍Ｄ庑畔⒌氖凳备�。本实用新型一种基于实时卫星信号模拟器同步增减可见卫星构建方法，该方法具 体步骤如下步骤一系统初始化，信息处理�？橥ü谥卸闲藕沤邮丈衔换ㄐ欧⑺偷哪� 拟时间和所有导航卫星的星历参数，并计算每一颗卫星的位置。DSP程序同时接收串口发送 的接收机位置信息，并进行存储。步骤二 由卫星号顺序计算每颗卫星相对于接收机天线的位置，由此计算出卫星 的仰角，并判断该颗卫星是否可见。如果该颗卫星可见，则写入可见卫星表中，给该颗卫星 分配一个通道，读取步骤一中存储结果，进行下一颗星的可见判断；否则将该颗卫星标记为 暂时不可见，则直接进入下一颗星的可见判断。如果可见卫星的数量已经达到模拟器支持 的最大通道数，则停止可见卫星的判断。如果遍历完毕所有卫星后可见卫星数量仍未达到 模拟器支持的最大通道数，则将空闲的通道写入空闲通道表。步骤三信息处理�？橥ü芟咝畔⒔换ツ？橄蛐藕糯砟？榉⑺屯ǖ雷刺统� 始模拟信息，通道状态为一个32bit的整数，由最低位开始，每2bit表示一个通道的不同状态，最多可表示16通道的状态。其中，用00表示通道关闭，01表示通道正常工作，10表示 通道即将启动，11表示通道即将关闭。初始状态下，将所有分配了可见星的通道的对应状态 置为10，其余的空闲通道状态置为00。初始模拟信息包括初始载波相位、初始码片相位和 导航电文数据等信息。步骤四模拟器工作过程中，根据导航系统的卫星编号遍历计算每一颗星的位置， 从编号1号星开始遍历到32号星，以Is为间隔进行循环计算。同时根据载体运动模型实 时计算载体接收机位置，根据二者相对位置判断该颗星是否可见。如果该颗卫星可见，则 查看可见卫星表，若可见卫星表中已经有该颗卫星存在，则维持该颗卫星对应的通道状态 为“01”不变，若可见卫星表中没有该颗卫星，表示该卫星是即将可见的，则查看空闲通道列 表，若有空闲通道，则将该通道状态标记为“10”，并在空闲通道表中将该通道删除，若没有 空闲通道，则维持现有状态不变。如果该颗卫星不可见，则查看可见卫星表，若可见卫星表 中没有该颗卫星，则维持现有状态不变，若可见卫星表中存在该颗卫星，则表示该颗卫星即 将不可见，将对应通道状态标记为“11”，同时在可见卫星列表中删除该颗卫星。步骤五以中断脉冲信号作为触发，实时计算可见卫星的码NCO控制字和载波NCO 控制字，同时计算每个通道的初始码相位、初始载波相位和导航电文等通道数据信息。具体 实例中采用周期为0. 4ms的中断脉冲信号，对于每一颗可见卫星，包括之前存在的和新出 现的，都计算NCO频率字，对于新出现的卫星，将提前进行2帧导航电文数据的计算。1.针对每一颗可见卫星，分别计算该卫星的传输环境延迟，主要包括调用对流层 延迟计Δτκ。ρ和电离层延迟Δ_计算。2.计算伪距 式中，ρ C0DE为对应的码伪距，P f为对应的载波伪距，xST、xST、zST为卫星的位置，Xu、 YpZu为接收机的位置，c为光速，ΔΜ为星钟误差，A1-为电离层误差，Atot为对流层误 差。3.计算相邻采样时刻的伪距变化量 式中，δ ■、δ CAEEIEE为码和载波对应的伪距变化量，p c0de_0ld> p caeeiee_0ld为缓存的 上一个采样时刻码和载波对应的伪距。4.计算伪距变化量对应的偏移频率字 式中，Δ Ncode, Δ Ncaeeiee为计算出的码和载波的偏移控制字，fEF为GPS射频频率，fs 为采样频率，K为NCO深度，最后的xlOOO/m是为了对应ms计算周期，m为每一个计算周期 的ms数。5.将基准频率字和偏移频率字相加得到发送给FPGA的频率字。Ncode = Ncodo+ Δ Ncode
14[0168]Ncaeeiee 一 Ncaeeiee_b+ Δ Ncaeeiee式中，NroDE—Β和Ν。ΑΚΚΙΕΚ—Β分别为伪码和载波的基准频率。6.导航电文计算电文在模拟器中为初始计算一次，运行中每30s计算一次。初始计算和运行时计 算方法一样，都是根据收到的每颗星的星历参数和模拟时刻，对每一颗可见星每一个子帧 的电文逐比特生成。步骤六通过总线信息交互�？橄煊χ卸闲藕牛葑芟呓涌诼呒街栉逯屑� 算的数据(主要包括NCO控制字，通道状态，导航电文)通过总线信息交互�？榉⑺椭列藕� 处理�？�。模拟器正常运行时步骤四和步骤五分别按照设计的周期循环计算。步骤四中所 计算卫星号从1号卫星开始并在本次计算结束后将卫星号加1，如果当前卫星号达到最大 卫星号，则将卫星号置1重新遍历计算。步骤五中所计算的频率控制字根据可见卫星列表 状态，每次计算周期对所有状态为“01”和“ 10”的通道进行计算。步骤一-六为模拟器信 息处理部分，其流程原理如图2所示。步骤七通过复位信号对信号处理模块中的GPS通道(1-12)进行读数寄存器的复 位。复位逻辑完成后产生中断脉冲信号，并驱动总线接口�？椋畔⒋砟？榉⑺屠吹母� 种控制信息和参数进行锁存。具体方法为在FPGA中判断总线接口逻辑中的读数信号，当写片选信号为低时，在 DSP写使能信号上升沿处进行触发，通过地址总线的地址译码，通过数据总线将DSP发送来 的各种控制信息和参数进行锁存，DSP写数逻辑如图5所示。其中DSP写入的数据包括通道码频率控制字CHx_C0DE_NC0_INC_F通道初始码相位控制字CHx_C0DE_NC0_INC_PHASE通道载波频率控制字CHx_CARRIER_NCO_INC_F通道载波相位控制字CHx_CARRIER_NCO_INC_PHASE通道待调制导航电文CHx_NAV_DATA通道初始参数(包括该通道模拟的卫星号，初始码片周期和导航电文计数延迟)CHx_sat_num通道状态控制字chan_sat_state_con步骤八信号处理�？樵谥卸闲藕诺目刂葡露寥⊥ǖ雷刺畔ⅲ⒁员镜赝铰� 冲为触发更新各个模拟通道的通道状态。对于即将开始模拟的通道，则将对应的初始模拟 信息置入相应寄存器组。通道实时通道状态用一个32bit的状态字表示，从低到高，每2bit 数据表示一个通道的状态信息。每个通道有4种状态，其状态转移逻辑为状态“00'状态“01'状态“10'状态“11'
通道保持关闭，相关数据寄存器复位(保持状态1) 通道持续运行，相关数据寄存器实时更新(保持状态2) 有新的可见卫星出现，该通道即将模拟该颗卫星信号(瞬时状态1) 该通道目前模拟的卫星即将消失，通道即将关闭(瞬时状态2) 当模拟通道的状态为“00”时，表示该通道处于关闭状态，不对其内部的逻辑进行 处理，该通道不产生调制信号；当模拟通道的状态为“01”时，表示该通道处于工作状态，其 内部处理逻辑在每个系统中断时刻更新相关数据寄存器信息，持续产生调制信号；当模拟 通道的状态为“ 10”时，表示该通道即将开始工作，进入通道开启状态的跳转；当该模拟通道的状态为“U”时，表示该通道即将停止工作，进入通道关闭状态的跳转。上述通道状态控制逻辑能够保证在模拟器出现可见卫星的增加或减少情况时，即 将开启或关闭的通道和正在运行的各通道严格保持同步，从而使得模拟器模拟的可见星星 座变化与真实情况一致，不会产生影响接收机的定位解算处理的误差。同步通道开关控制 状态机结构如图4所示。步骤九在通道同步脉冲的触发下，将通道状态控制字和需要回传信息处理�？� 的实时数据进行偏置调整后赋与相应的数据锁存寄存器；需要回传的实时数据包括通道实时载波NCO累加值CHx_carrier_phase通道实时码NCO累加值CHx_code_phase通道码片计数CHx_code_addr实时通道状态信息chan_sat_state判断总线接口逻辑中的写数信号，利用地址译码将需要向信息处理�？榉⑺偷母� 种实时信息和参数所存在DSP读数据总线锁存寄存器中，当读片选信号和DSP读使能信号 同为低电平出进行触发，将锁存器中数据通过地址总线发送，将实时改变的通道状态及仿 真过程中产生的有效数据回传至信息处理�？�。DSP读数逻辑如图6所示。步骤十通过同步通道开关控制状态机逻辑实现信号处理模块中的通道状态同步 改变，根据步骤八中描述的通道控制逻辑，将步骤七中锁存的通道状态控制信息置入同步 通道开关控制状态机，实现通道状态同步改变；将步骤七中锁存的模拟数据同步置入FPGA 相应通道。对于每个工作通道，在每次中断信号来临时重复执行步骤五-十，对通道内的模 拟参数进行同步更新，根据接收的控制字信息实时生成具有高动态多普勒特性的载波和伪 码，并完成导航电文数据调制，实现高动态信号模拟器的实时工作。步骤七_十为模拟器信 号处理部分，其流程原理如图3所示。结果表示本实用新型应用实例基于的硬件平台选用XilinxV4SX55FPGA芯片和TI公司生产 的C6713DSP芯片，系统工作时钟68MHz。模拟12通道GPS信号，动态环境为高动态。测试 载体为JAVAD JNS100商业接收机。随着载体运动轨迹的变化，载体和卫星之间的相对位置 也在不断变化，期间模拟器可见星座也随之改变。本实用新型中的构建方法使得在可见卫 星的增加或减少同步改变，模拟器的模拟数据不发生抖动和偏差。图7-图9是基于一组载 体绕地球高速运动轨迹的模拟器测试结果，分别截取了一段工作时间内模拟器模拟输出的 用户位置三维坐标(经纬高)和商业接收机定位解算输出的用户位置三维坐标。从图中可 以看出，在长时间模拟周期内，接收机输出的用户位置曲线较为平滑，其变化趋势与模拟器 模拟的原始运动轨迹整体一致，没有出现明显的误差和抖动。这表明模拟器的模拟数据整 体稳定，在可见卫星星座分布情况频繁改变时未出现模拟偏差，没有影响接收机的正确解 算，验证了基于该构建方法的模拟器长时间的正确工作。
权利要求一种基于实时卫星信号模拟器同步增减可见卫星装置，其特征在于它包括上位机通信�？�、信息处理�？椤⑿藕糯砟？楹妥芟咝畔⒔换ツ？椋簧衔换ㄐ拍？橥ü齊S232串行接口与信息处理�？榱樱畔⒋砟？橥ü中藕糯砥骷碊SP和现场可编程门阵列即FPGA芯片之间的总线接口逻辑与总线信息交互�？榱樱藕糯砟？橥ü鼶SP和FPGA芯片之间的总线接口逻辑与总线信息交互模块连接；所述上位机通信�？槭窃赑C机上构建的系统显控界面软件，它的功能是进行模拟器系统初始化操作，设置模拟初始参数和控制信息，包括模拟用户轨迹设定、星历文件选取、初始用户位置和时刻设置，并通过串口逻辑将上述信息下发到信息处理�？椋痹诠ぷ鞴讨写有畔⒋砟？橹薪邮漳Ｄ庑畔ⅲǹ杉佬切畔ⅲ魍ǖ滥Ｄ庾刺畔ⅲ谙钥亟缑嬷薪邢允荆凰鲂畔⒋砟？樵贒SP芯片中实现，它的功能是接收上位机通信�？榉⑺偷牟问椭噶罱邢低吵跏蓟褪菪畔⒋恚恍畔⒋砟？橛晌佬俏恢糜朐靥褰邮栈恢檬凳奔扑阕幽？�、同步增减可见卫星实时判断子�？楹屯ǖ滥Ｄ庑畔⒓扑阕幽？楣钩桑凰鑫佬俏恢糜朐靥褰邮栈恢檬凳奔扑阕幽？椋涔δ苁歉菸佬切抢问驮靥逶硕旒２问凳奔扑闼Ｄ馕佬堑奈恢煤驮靥褰邮栈奈恢茫桓米幽？榈慕峁褂尚抢馕鼋峁固�、卫星位置计算函数和接收机轨迹计算函数构成，其算法原理采用现有技术实现；所述同步增减可见卫星实时判断子模块的结构由卫星仰角计算函数、可见卫星列表更新函数和通道状态设置函数构成；卫星仰角计算函数根据卫星位置与载体接收机位置实时计算子�？槭涑龅奈佬俏恢煤徒邮栈恢檬菁扑忝靠盼佬窍喽杂谠靥褰邮栈难鼋切畔ⅲ凳苯锌杉亲卸希芯鼋峁葜量杉橇斜砀潞豢杉橇斜砀潞寥〈宋佬窃诟眉扑阒芷诤蜕弦患扑阒芷谀诘目杉耘芯鼋峁罂傻玫礁每盼佬鞘粲诤沃肿刺源俗刺驴杉佬橇斜恚⒔俗刺葜镣ǖ雷刺柚煤煌ǖ雷刺柚煤菸佬亲刺心Ｄ馔ǖ赖淖刺柚茫媚Ｄ馔ǖ雷刺ǹ�、持续模拟、关闭和持续关闭；所述通道模拟信息计算子�？榈墓δ苁歉萃皆黾蹩杉佬鞘凳迸卸献幽？榻⒌奈佬橇斜硇畔ⅲ扑闼锌杉佬堑奈甭�/载波控制字信息、传输延迟、多普勒频移及导航电文信息，将上述模拟信息通过总线信息交互�？榇涞叫藕糯砟？椋苯糠帜Ｄ庑畔⒒卮衔换ㄐ拍？榻邢允荆桓米幽？榈慕峁褂晌本嗉扑愫�、电离层/对流层延迟计算函数、伪码/载波频率字计算函数和导航电文组帧函数构成，起算法原理采用现有技术实现；所述信号处理�？樵贔PGA芯片中实现，它的功能是实现各模拟通道的状态控制和信号处理，完成多通道卫星信号调制；该信号处理�？橛赏ǖ雷刺娇刂谱幽？楹托藕诺髦谱幽？楣钩桑凰鐾ǖ雷刺娇刂谱幽？榈慕峁褂赏ǖ揽乜刂谱刺钩桑渥刺坡呒Ｖ棺刺ǖ辣３止乇眨喙厥菁拇嫫鞲次唬３指米刺辉诵凶刺ǖ莱中诵校喙厥菁拇嫫魇凳备拢３指米刺黄舳刺行碌目杉佬浅鱿郑猛ǖ兰唇Ｄ飧每盼佬切藕牛浦猎诵凶刺还乇兆刺猛ǖ滥壳澳Ｄ獾奈佬羌唇Вǖ兰唇乇眨浦镣Ｖ棺刺桓猛ǖ雷刺娇刂谱幽？榈墓δ苁窍煊π畔⒋砟？榉⑺偷耐ǖ揽刂菩畔⒔心Ｄ馔ǖ劳娇乜刂坡呒迪挚杉佬乔榭龈谋涫蹦Ｄ馔ǖ赖淖既繁浠蓖ü芟咝畔⒔换ツ？槎酝ǖ滥Ｄ庑畔⒔刑崛『突捍妫瓿捎谛畔⒋砟？橹涞氖萁换ィ凰鲂藕诺髦仆ǖ赖墓δ苁峭ü隢CO和载波NCO设计实现动态伪码和载波的生成，完成导航电文数据的扩频调制，生成模拟器发射信号；该子模块的结构由伪码发生器模块、码NCO�？�、载波NCO�？楹偷髦颇？楣钩桑渌惴ㄔ聿捎孟钟屑际跏迪郑凰鲎芟咝畔⒔换ツ？橛勺芟咝畔⒔换ツ？橛蒁SP读数逻辑、DSP写数逻辑和中断函数构成；通过DSP读数逻辑实现数据信息由信号处理模块到信息处理�？榈拇莺突捍妫ü鼶SP写数逻辑实现数据信息由信息处理�？榈叫藕糯砟？榈拇莺突捍妫ü卸虾低持卸闲藕藕屯铰龀澹繁？杉佬乔榭龊湍Ｄ庑畔⒌氖凳备�。
专利摘要一种基于实时卫星信号模拟器同步增减可见卫星装置，它包括上位机通信�？�、信息处理�？�、信号处理�？楹妥芟咝畔⒔换ツ？椋桓蒙衔换ㄐ拍？橥ü齊S232串行接口与信息处理�？榱樱眯畔⒋砟？橥ü中藕糯砥骷碊SP和现场可编程门阵列即FPGA芯片之间的总线接口逻辑与总线信息交互�？榱樱眯藕糯砟？橥ü鼶SP和FPGA芯片之间的总线接口逻辑与总线信息交互�？榱樱凰錾衔换ㄐ拍？槭窃赑C机上构建的系统显控界面软件，它的功能是进行模拟器系统初始化操作，设置模拟初始参数和控制信息；所述信息处理�？樵贒SP芯片中实现；所述信号处理模块在FPGA芯片中实现；所述总线信息交互�？橛勺芟咝畔⒔换ツ？橛蒁SP读数逻辑、DSP写数逻辑和中断函数构成。
文档编号G01S19/40GK201673260SQ20102018090
公开日2010年12月15日 申请日期2010年4月29日 优先权日2010年4月29日
发明者常青, 齐巍 申请人:北京航空航天大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
