发明名称 多路GPS信号P(Y)码相位差的测量系统
摘要 本发明公开了一种多路GPS信号的P(Y)码相位差的测量系统,该系统包括有直射信号捕获跟踪模块、阵列乘法器模块、高度测量模块。依据GPS的正交调制特性提取接收到的卫星信号的正交和同相信号,分离出C/A码和P(Y)码,然后分别对其进行等间距相同的延迟,再将其延迟后的信号和反射信号进行互相关运算,不同的延迟码片导致不同的非相干积分值,得到不同路GPS信号的P(Y)码相位差。该测量系统降低了硬件接收设备的复杂度,同时测量精度也得到了提高。
申请公布号 CN103017731B 申请公布日期 2014.08.06
申请号 CN201210480241.1 申请日期 2012.11.22
申请人 北京航空航天大学 发明人 李伟强;杨东凯;秦瑾;张波;张帅
分类号 G01C5/00(2006.01)I 主分类号 G01C5/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 李有浩
主权项 一种多路GPS信号P(Y)码相位差的测量系统,其特征在于:该测量系统包括有直射信号捕获跟踪模块(10)、阵列乘法器模块(20)、高度测量模块(30); 直射信号捕获跟踪模块(10)用于对直射数字中频信号s<sub>D‑IF</sub>(n)进行捕获、跟踪及导航电文的解调,得到载波频率控制字F<sub>CW</sub>,反射卫星高度角ε,反‑直延迟值N,直‑基正交s<sub>D‑BB‑P</sub>(n),直‑基同相s<sub>D‑BB‑C</sub>(n);其中,反射卫星高度角ε输出给高度测量模块(30);载波频率控制字F<sub>CW</sub>、直‑基正交s<sub>D‑BB‑P</sub>(n)、直‑基同相s<sub>D‑BB‑C</sub>(n)、反‑直延迟值N分别输出给阵列乘法器模块(20); 阵列乘法器模块(20)利用多路延迟器阵列对反射数字中频信号s<sub>R‑IF</sub>(n)、直‑基正交s<sub>D‑BB‑P</sub>(n)和直‑基同相s<sub>D‑BB‑C</sub>(n)进行处理,得到直‑基正交互相关功率值X<sub>P</sub>(τ<sub>n</sub>)和直‑基同相互相关功率值X<sub>C</sub>(τ<sub>n</sub>),并输出到高度测量模块(30); 高度测量模块(30)包括反‑基正交功率值曲线拟合模块(304)、反‑基同相功率值曲线拟合模块(305)、P(Y)码峰值检测模块(302)、C/A码峰值检测模块(303)、高度计算模块(301); 反‑基正交功率值曲线拟合模块(304)将多个反‑基正交功率值,即反‑基正交初始功率值X<sub>P(0)</sub>、反‑基正交第一功率值X<sub>P(1)</sub>、反‑基正交第二功率值X<sub>P(2)</sub>、……、反‑基正交第M功率值X<sub>P(M)</sub>按照散射信号相关功率模型&lt;|YS(Δ,f)|<sup>2</sup>&gt;进行曲线拟合,输出P(Y)码的一维功率值‑时延函数X<sub>P(</sub><sub>τ</sub><sub>)</sub>; <img file="FDA0000511958410000011.GIF" wi="1600" he="169" />D<sup>2</sup>(r)表示天线增益的平方值,决定了天线覆盖区; Λ<sup>2</sup>[Δ(r)]表示反‑基正交功率值或者反‑基同相功率值,决定了等延迟区; S[Δf(r)]表示等多普勒区,其中Δf(r)=f<sub>c</sub>‑f<sub>0</sub>,即多普勒频差; <img file="FDA0000511958410000012.GIF" wi="141" he="70" />表示发射机到散射点的距离的平方;R<sub>r</sub><sup>2</sup>(r)表示接收机到散射点的距离的平方; σ<sub>0</sub>(r)d<sup>2</sup>r表示归一化双基散射截面,散射截面的大小决定了被称为照射区的海面区域;d<sup>2</sup>表示二维积分; r表示距离向量; A<sup>2</sup>表示幅度因子的平方值; T<sup>2</sup>表示积分时间的平方值; W<sup>2</sup>表示加性高斯白噪声; 反‑基同相功率值曲线拟合模块(305)将多个反‑基同相功率值,即反‑基同相初始功率值X<sub>C(0)</sub>、反‑基同相第一功率值X<sub>C(1)</sub>、反‑基同相第二功率值X<sub>C(2)</sub>……、反‑基同相第M功率值X<sub>C(M)</sub>按照散射信号相关功率模型&lt;|YS(Δ,f)|<sup>2</sup>&gt;进行曲线拟合,输出C/A码的一维功率值‑时延函数X<sub>C(</sub><sub>τ</sub><sub>)</sub>; P(Y)码峰值检测模块(302)对P(Y)码的一维功率值‑时延函数X<sub>P(</sub><sub>τ</sub><sub>)</sub>进行采用峰值跟踪法,拟合出峰值点,计算出峰值点对应的时延值,再将其转化为反射信号和直射信号之间的距离差ρ<sub>p</sub>; C/A码峰值检测模块(303)对C/A码的一维功率值‑时延函数X<sub>C(</sub><sub>τ</sub><sub>)</sub>进行采用峰值跟踪法,拟合出峰值点,计算出峰值点对应的时延值,再将其转化为反射信号和直射信号之间的距离差ρ<sub>c</sub>; 高度计算模块(301)对多路P(Y)码相位差计算直射信号和反射信号的距离差ρ<sub>p</sub>、多路C/A码相位差计算直射信号和反射信号的距离差ρ<sub>c</sub>、反射卫星高度角ε进行处理,根据海面反射的空间几何关系,信号处理平台与反射点的相对高程的计算公式为<img file="FDA0000511958410000021.GIF" wi="270" he="120" />其中h表示信号处理平台与反射点的相对高程,则<img file="FDA0000511958410000022.GIF" wi="294" he="127" />h<sub>p</sub>表示多路GPS信号P(Y)码的相位差的高度测量结果,<img file="FDA0000511958410000023.GIF" wi="290" he="120" />h<sub>c</sub>表示多路GPS信号C/A码的相位差的高度测量结果,两者的测高精度可进行对比。
地址 100191 北京市海淀区学院路37号