专利名称：多路gps信号p(y)码相位差的测量系统的制作方法
技术领域：
本发明涉及一种相位码的测量，更具体地说，是指一种利用导航卫星的直射信号 代替本地信号与反射信号进行互相关运算，利用互相关结果进行多路GPS信号P (Y)码相位 差的估计，进而将估计的结果应用在海面高度测量、海面风场反演、飞行器测姿、形变监测 等等应用方面。
背景技术：
近年来，基于全球导航卫星系统反射信号的微波遥感技术以其全天候、全天时、多 信号源、宽覆盖、高时空分辨率等应用优势，在遥感领域展现出广阔的应用前景。迅速发展 起来的GPS测量技术使传统大地测量学的面貌焕然一新，给测绘领域带来一场深刻的技术 革命。GPS测量技术是在海、陆、空任意测站点，通过天线连续接收GPS卫星的直射或者反射 信号，并计算测站点的三维坐标、海拔高度、运动载体的三维速度及时间等。
目前在该技术领域一般采用GPS接收机设备完成，在GPS接收机中由于C/A码周 期短，实现相关算法方便，很多研究人员已经通过测量多路GPS信号的C/A码相位差来进行 海面高度测量。在岸基实验方面，Martin-Neira等人在最早进行的大桥实验中对直射信号 和反射信号分别进行处理，在仅利用C/A码相位的情况下达到了 3m的测高精度Jreuhaft 等人在火山口湖实验中利用本地信号与直射和反射信号进行相关运算，利用非线性参数估 计的方法提取了直射和反射信号的载波相位差，并在Is时间内进行平均，获得了 2cm的测 高精度；Ruffini等人和Rivas等人分别在海港中平静海面的条件下利用码信号的载波相 位获得3.1cm和5cm的测闻精度。在机载实验方面，Lowe在实验中获得了 60cm的测闻精 度，并认为长时间的平均可使测高精度提高至5cm量级。
通过上述实验，可见利用C/A码来进行多路GPS信号相位差的测高精度并不高。针 对上述问题，本发明提出一种利用直射信号与反射信号互相关的海面测高方法，该方法利 用GPS信号结构和特点，对直射信号的同相(包含C/A码)分量和正交(包含P(Y)码)分量 进行正交分离，采用正交分离后的直射信号代替本地信号与反射信号分别进行相关运算， 利用互相关结果进行多路GPS信号的相位差估计，从而实现接收机所在位置距离反射面的 高度测量。该系统避免了反射信号处理中本地信号的产生，降低了接收设备的复杂度，并利 用具有更高扩频增益的P(Y)码信号提高了相位差的估计精度。发明内容
本发明的目的是提供一种多路GPS信号P (Y)码相位差的测量系统，该系统通过 在FPGA芯片上用Verilog HDL语言设计得到硬件的功效。在本发明中，采用时域伪码串行 载波串行的捕获跟踪方法对直射数字中频信号进行捕获跟踪；采用阵列乘法器进行多路直 射数字基频信号和反射数字基频信号的相关功率值的运算；采用跟踪峰值方法对多路不同 延时的直射数字基频信号和反射数字基频信号的相关功率值进行峰值检测处理；并根据海 面反射的空间几何关系对直射信号和反射信号之间的距离差和卫星高度角进行处理，得到信号处理平台至反射面的相对高程。
本发明的一种多路GPS信号P (Y)码相位差的测量系统，该测量系统包括有直射信号捕获跟踪�？�(10 )、阵列乘法器�？�(20 )、高度测量�？�(30 );
直射信号捕获跟踪�？�(10)用于对直射数字中频信号sD_IF(n)进行捕获、跟踪及导航电文的解调，得到载波频率控制字Few，反射卫星高度角ε，反-直延迟值N，直-基正交sD_BB_P(n)，直-基同相sD_BB_e(n);其中，反射卫星高度角ε输出给高度测量�？�(30);载波频率控制字Fcw、直-基正交sD_BB_P(n)、直-基同相sD_BB_c(n)、反-直延迟值N分别输出给阵列乘法器�？�(20)；
阵列乘法器模块(20)利用多路延迟器阵列对反射数字中频信号sK_IF(n)、直-基正交sD_BB_P(n)和直-基同相sD_BB_e(n)进行处理，得到直-基正交互相关功率值ΧΡ( τ n)和直-基同相互相关功率值XJ τ n)，并输出到高度测量�？�(30)；
高度测量�？�(30)包括反-基正交功率值曲线拟合�？�(304)、反-基同相功率值曲线拟合�？�(305 )、P (Y)码峰值检测模块(302)、C/A码峰值检测�？�(303)、高度计算�？�(301)；
反-基正交功率值曲线拟合�？�(304)将多个反-基正交功率值，即反-基正交初始功率值Xpft0、反-基正交第一功率值Xp(1)、反-基正交第二功率值Xp(2)、……、反-基正交第M功率值XP(M)按照散射信号相关功率模型<|YD(A，f) |2>进行曲线拟合，输出P (Y)码的一维功率值-时延函数ΧΡ( τ)。
反-基同相功率值曲线拟合�？�(305)将多个反-基同相功率值，即反-基同相初始功率值χ_)、反-基同相第一功率值xc(1)、反-基同相第二功率值xc(2)……、反-基同相第M功率值Xe(M)按照散射信号相关功率模型〈IYS (△，f) 12>进行曲线拟合，输出C码的一维功率值-时延函数\(τ)。
P (Y)码峰值检测模(302)块对P (Y)码的一维功率值-时延函数ΧΡ(τ)进行采用峰值跟踪法，拟合出峰值点 ，计算出峰值点对应的时延值，再将其转化为反射信号和直射信号之间的距离差P ρ。
C/A码峰值检测�？�(303)对C/A码的一维功率值-时延函数Xe(T)进行采用峰值跟踪法，拟合出峰值点，计算出峰值点对应的时延值，再将其转化为反射信号和直射信号之间的距离差P。。
高度计算�？�(301)对多路P (Y)码相位差计算直射信号和反射信号的距离差 Pp、多路码相位差计算直射信号和反射信号的距离差P。、反射卫星高度角ε进行处理，根据海面反射的空间几何关系，信号处理平台与反射点的相对高程的计算公式为IsinE其中h表示信号处理平台与反射点的相对高程，则\ =，hp表示多路GPS信号P⑴2sin^码的相位差的高度测量结果，<，h。表示多路GPS信号C/A码的相位差的高度测2smi*量结果，两者的测高精度可进行对比。
本发明多路GPS信号P (Y)码相位差的测量系统的优点在于
①利用P(Y)码相位差进行测高可以满足不同观测平台的需求，同时也满足高精度测高的需求。
②利用P(Y)码相位差进行测高对环境的要求较宽松，应用范围不受限制，具有遥感卫星信号测量的特点，即宽覆盖、高时空分辨率等。
③在信号处理中，利用接收到的直射信号代替本地信号和反射信号进行互相关运算的处理方式，避免了产生本地信号，降低了接收设备的硬件复杂度。
④利用多路GPS信号的C/A码相位差进行测高获得的测量结果，可以与利用多路 GPS信号的P (Y)码相位差进行测高获得的测量结果进行精度对比。对比可知P (Y)码的测量精度更高。
图1是本发明利用多路GPS信号的P(Y)码相位差、C/A码相位差进行互相关测高装置的结构示图。
图2是本发明阵列乘法器�？榻峁故就�。
图3是本发明高度测量模块结构示图。
具体实施方式
GPS信号上存在着C/A码和P (Y)码两种测距码，C/A码只调制在LI载波信号上， 周期为1ms，码率为1.023Mcps ;而？化)码同时调制在LI和L2载波信号上，周期为7天，码率为 10. 23Mcps。
参见图1所示，本发明是一种利用导航卫星的直射信号代替本地信号与反射信号进行互相关运算，利用不同的互相关功率值进行多路GPS信号的P(Y)码相位差估计，进而实现海面高度测量、海面风场反演、飞行器测姿、形变监测等等的测量系统，该多路GPS信号的P码相位差的测量系统包括有直射信号捕获跟踪�？�10、阵列乘法器�？�20、高度测量�？�30。该系统通过在GPS接收机的FPGA芯片上用Verilog HDL语言设计得到硬件的功效。
(一)直射信号捕获跟踪模块10
直射信号捕获跟踪模块10用于快速对接收到的GPS卫星的直射数字中频信号 sD_IF(η)信号进行捕获、跟踪及导航电文的解调。
在本发明中，直射信号捕获跟踪�？�10采用了在2006年7月，第27卷第4期，《遥测�？亍饭摹癎PS接收机基带信号处理�？榈腇PGA实现”中提到的算法，分别得到载波频率控制字Few，反射卫星高度角ε，反射信号相对于直射信号的延迟相位估计值Ν(简称为反-直延迟值N)，直射数字基频信号的正交分量(包含P (Y)码)sD_BB_P(n)(简称为直-基正交sD_BB_P(n))和直射数字基频信号的同相分量(包含C/A码)sD_BB_c(n)(简称为直-基同相 sD-BB-c(η));其中，反射卫星高度角ε输出给高度测量�？�30 ;载波频率控制字Few、直-基正交sD_BB_P(n)、直-基同相sD_BB_c(n)、反-直延迟值N分别输出给阵列乘法器�？�20。
sD_BB_P(n)和 sD_BB_c(n)中的(η)表示离散时间。
(二)阵列乘法器模块20
参加图2所示，阵列乘法器模块20利用多路延迟器阵列对接收到的反射数字中频信号sK_IF (η)、直-基正交sD_BB_P (η)和直-基同相sD_BB_。(η)进行处理，得到直-基正交互相关功率值ΧΡ(τη)和直-基同相互相关功率值Xe(Tn),并输出到高度测量�？�30。τη表示相关时间的取值范围，即O、Λ、......、ΜΛ。
所述直-基正交互相关功率值Xp ( τη)是指反射数字中频信号81；1(11)与直-基正交sD,-P(n)的不同时延间距信号的互相关功率值。
所述直-基同相互相关功率值\(τη)是指反射数字中频信号^1 (η)与直-基同相sD,-c(n)的不同时延间距信号的互相关功率值。
具体地，阵列乘法器模块20对反射数字中频信号sK_IF (η)、直-基正交sD_BB_P (η)和直-基同相sD_BB_c(n)的处理步骤如下
步骤201 :载波NCO根据载波频率控制字Few产生对应频率的正弦波和余弦波，产生的波形取决于ROM中存放的是正弦表还是余弦表；
步骤202 :正弦乘法器中将反射数字中频信号sK_IF(n)和正弦表输出的正弦信号进行相乘，输出反射数字基频信号的正交分量(包含P (Y)码)sK-P(n)(简称为反-基正交 sE-p (η))；
余弦乘法器将反射数字中频信号sK_IF(n)和余弦表输出的余弦信号进行相乘，输出反射数字基频信号的同相分量(包含C/A码)sE_c(n)(简称为反-基同相SK_e(n))。
步骤203 :延迟器P将直-基正交sD_BB_P(n)延迟N个时间单位，即信号将转换成延迟直-基正交sD_BB_P(n_NA ) ；N是直射信号捕获跟踪�？槭涑龅牟捎弥鄙浔栈贩瓷淇返姆绞蕉苑瓷湫藕畔喽杂谥鄙湫藕诺难映傧辔坏墓兰浦担魇茄映俚氖奔洳匠ぃ煌砜傻�
延迟器pi将延迟直-基正交sD_BB_P(n-NA )延迟I个时间步长Λ，即信号将转换成第一延迟直-基正交sD_BB_P(η-ΝΔ-Δ)；
延迟器ρ2将第一延迟直-基正交sD_BB_P(n-NA-A)延迟I个时间步长Λ，即信号将转换成第二延迟直-基正交sD_BB_P(η-ΝΔ-2Δ)；……；
延迟器ρΜ将第M-1延迟直-基正交sD_BB_P (η-Ν Λ-(M-1) Δ)延迟I个 时间步长Λ， 即信号将转换成第M延迟直-基正交sD_BB_P(η-ΝΔ-ΜΔ)；
延迟器c将直-基同相sD_BB_c(n)延迟N个时间单位，即信号将转换成延迟直_基同相sD_BB_Jn-NA) ;N是直射信号捕获跟踪�？槭涑龅牟捎弥鄙浔栈贩瓷淇返姆绞蕉苑瓷湫藕畔喽杂谥鄙湫藕诺难映傧辔坏墓兰浦担魇茄映俚氖奔洳匠ぃ煌砜傻�
延迟器Cl将延迟直-基同相sD_BB_c(n_NA)延迟I个时间步长Λ，即信号将转换成第一延迟直-基同相sD_BB_c(n-NA-A)；
延迟器c2将第一延迟直-基同相sD_BB_c(n-NA-A)延迟I个时间步长Λ，即信号将转换成第二延迟直-基同相sD_BB_c(η-ΝΔ-2Δ)；……；
延迟器CM将第M-1延迟直-基同相sD_BB_c (η-Ν Λ-(M-1) Δ)延迟I个时间步长Λ， 即信号将转换成第M延迟直-基同相sD_BB_c(η-ΝΔ-ΜΔ)；
步骤204 :乘法器p0将反-基正交sK_P (η)与延迟直_基正交sD_BB_P (η-Ν Δ )进行相乘，得到反-基正交初始互相关值SK_D_PW ;同理可得
乘法器pi将反-基正交sK_P (η)与第一延迟直-基正交sD_BB_P(n-NA-A)进行相乘，得到反-基正交第一互相关值 Sr-D-P(I);
乘法器P2将反-基正交sK_P (η)与第二延迟直_基正交sD_BB_P (η-Ν Δ-2 Δ )进行相乘，得到反-基正交第二互相关值 Sr-D-P (2);......;
乘法器PM将反-基正交sK_P (η)与第M延迟直-基正交sD_BB_P (η-Ν Δ-ΜΔ)进行相乘，得到反-基正交第M互相关值Sr-D-P(M);
乘法器CO将反-基同相sK_c(n)与延迟直-基同相sD_BB_c(η_ΝΛ )进行相乘，得到反-基同相初始互相关值 Sr-D-C(O);
乘法器Cl将反-基同相sK_c (η)与第一延迟直-基同相sD_BB_c(n-NA-A)进行相乘，得到反-基同相第一互相关值 Sr-D-C(I);
乘法器c2将反-基同相sK_c(n)与第二延迟直_基同相sD_BB_c (η_Ν Λ _2 Λ )进行相乘，得到反-基同相第二互相关值 Sr-D-C (2);......;
乘法器CM将反-基同相sK_c(n)与第M延迟直-基同相sD_BB_c(n_NA _ΜΛ )进行相乘，得到反-基同相第M互相关值Sk_d_c(M)；
步骤205 :积分累加单元PO将反-基正交初始互相关值SK_D_P(CI)在预检积分时间内进行累加，得到未截断反-基正交初始功率值Ypftl)；
积分累加单元pi将反-基正交第一互相关值SK_D_P(1)在预检积分时间内进行累加， 得到未截断反-基正交第一功率值YP(1)；
积分累加单元ρ2将反-基正交第二互相关值SK_D_P(2)在预检积分时间内进行累加， 得到未截断反-基正交第二功率值YP(2)；……；
积分累加单元ρΜ将反-基正交第M互相关值SK_D_P(M)在预检积分时间内进行累加， 得到未截断反-基正交第M功率值YP(M)；
积分累加单元CO将反-基同相初始互相关值SK_D_e((l)在预检积分时间内进行累加， 得到未截断反-基同相初始功率值Y_)；
积分累加单元Cl将反-基同相第一互相关值SK_D_e(1)在预检积分时间内进行累加， 得到未截断反-基同相第一功率值Yc(1)；
积分累加单元c2将反-基同相第二互相关值SK_D_e(2)在预检积分时间内进行累加， 得到未截断反-基同相第二功率值Yc(2)；……；
积分累加单元CM将反-基同相第M互相关值SK_D__在预检积分时间内进行累加， 得到未截断反-基同相第M功率值Ye(M)；
步骤206 :截断器pO将未截断反-基正交初始功率值YP(Q)进行12比特的截断处理，得到反-基正交初始功率值 Xp(O);
截断器pi将未截断反-基正交第一功率值YP(1)进行12比特的截断处理，得到反-基正交第一功率值Xp(1)；
截断器p2将未截断反-基正交第二功率值Yp⑵进行12比特的截断处理，得到反-基正交第二功率值Xp(2)；……；
截断器ρΜ将未截断反-基正交第M功率值YP(M)进行12比特的截断处理，得到反-基正交第M功率值XP(M)；
截断器CO将未截断反-基同相初始功率值Y_)进行12比特的截断处理，得到反-基同相初始功率值χ_)；
截断器Cl将未截断反-基同相第一功率值Ye(1)进行12比特的截断处理，得到反-基同相第一功率值xc(1)；
截断器c2将未截断反-基同相第二功率值Y。⑵进行12比特的截断处理，得到反-基同相第二功率值xc(2)；……；
截断器cM将未截断反-基同相第M功率值YC(M)进行12比特的截断处理，得到反-基同相第M功率值XC(M)。
在本发明中，采用阵列形式排布乘法器、延迟器、积分累加器、截断器进行多路直射数字基频信号和反射数字基频信号的相关功率值的运算，实现了在信号处理中，利用接收到的直射信号代替本地信号和反射信号进行互相关运算的处理方式，避免了产生本地信号，降低了接收设备的硬件复杂度。
(三)高度测量�？�30
参见图3所示，高度测量�？�30包括有反-基正交功率值曲线拟合�？�304、 反-基同相功率值曲线拟合模块305、P码峰值检测�？�302、C码峰值检测�？�303、高度计算�？�301。
在本发明中，记录保存下多个截断器输出的反-基正交功率值和反-基同相功率值，应用散射信号相关功率模型〈IYS ( △，f) 12>得到二维功率值-时延-多普勒函数。而采用的散射信号相关功率模型〈I YS( Λ，f) |2>为时间延迟步长Λ和多普勒频率f两变量的函数。/,|Λf.D2(r)xk2 ΓΔ( ) χ|^Γδ/(γ)12
(|Κ(Δ,/)| ) = ^2χΓ2χ||-4^(r)xi ,'(r)-^{r)d2r^W2·,
D2(r)表示天线增益的平方值，决定了天线覆盖区；
Λ2[Δ (r)]表示反-基正交功率值或者反-基同相功率值，决定了等延迟区；
S[Af(r)]表示等多普勒区，其中Af(r) = 即多普勒频差；
杧(r)表示发射机到散射点的距离的平方；
Rr2 (r)表示接收机到散射点的距离的平方；
OtlOOd2r表示归一化双基散射截面，散射截面的大小决定了被称为照射区的海面区域；d2表示二维积分；
r表示距离向量；
A2表示幅度因子的平方值；
T2表示积分时间的平方值；
W2表不加性闻斯白噪声。
反-基正交功率值曲线拟合�？�304将多个反-基正交功率值，即反-基正交初始功率值Xpft0、反-基正交第一功率值Xp(1)、反-基正交第二功率值Xp(2)、……、反-基正交第M功率值XP(M)按照散射信号相关功率模型〈IYS (△，f) 12>进行曲线拟合，输出P码的一维功率值-时延函数ΧΡ(τ)。
反-基同相功率值曲线拟合�？�305将多个反-基同相功率值，即反-基同相初始功率值χ_)、反-基同相第一功率值xc(1)、反-基同相第二功率值xc(2)……、反-基同相第 M功率值Xe(M)按照散射信号相关功率模型〈IYS (△，f) 12>进行曲线拟合，输出C码的一维功率值-时延函数X&)。
P (Y)码峰值检测�？�302对P (Y)码的一维功率值-时延函数ΧΡ(τ)进行采用峰值跟踪法，拟合出峰值点，计算出峰值点对应的时延值，再将其转化为反射信号和直射信号之间的距离差P ρ。
C/A码峰值检测�？�303对C/A码的一维功率值-时延函数\认τ)进行采用峰值跟踪法，拟合出峰值点，计算出峰值点对应的时延值，再将其转化为反射信号和直射信号之间的距离差P。。
高度计算模块301对多路P (Y)码相位差计算直射信号和反射信号的距离差P p (P (Y)码峰值检测模块输出的)、多路C/A码相位差计算直射信号和反射信号的距离差P。 (C码峰值检测�？槭涑龅�)、反射卫星高度角ε (直射信号捕获跟踪模块输出的)进行处理，根据海面反射的空间几何关系，信号处理平台与反射点的相对高程的计算公式为*=^7，其中h表示信号处理平台与反射点的相对高程，则有\K表示多路GPS信号P(Y)码的相位差的高度测量结果A，h。表示多路GPS信号C/A码的相位差的高度测量结果，两者的测高精度可进行对比。
在本发明中，高度计算�？�301计算出的反射信号和直射信号之间的相位差可进行多种应用，如海面高度测量、海面风场反演、飞行器测姿、形变监测等等，本发明中以海面高度测量为例。
本发明涉及的利用多路GPS信号的P (Y)码相位差进行测高的测量系统，实质是利用导航卫星的直射信号代替本地信号与反射信号进行互相关运算，利用GPS信号结构和特点，对直射信号的同相(包含C/A码)分量和正交(包含P(Y)码)分量进行正交分离，采用正交分离后的直射信号代替本地信号与反射信号分别进行相关运算，利用互相关结果进行多路GPS信号的相位差估计，从而实现接收机所在位置距离反射面的高度测量。该系统避免了反射信号处理中本地信号的产生，降低了接收设备的复杂度，并利用具有更高扩频增益的P(Y)码信号提高了相位差的估计精度。随着导航卫星系统的快速发展，空间中将有多个导航卫星系统并存，导航卫星信号资源日益丰富，利用导航卫星的直射信号和反射信号进行互 相关运算来计算出多路GPS信号的相位差的技术实施可行性与有效性将越来越强。
权利要求
1.一种多路GPS信号P (Y)码相位差的测量系统，其特征在于该测量系统包括有直射信号捕获跟踪模块(10 )、阵列乘法器�？�(20 )、高度测量�？�(30 ); 直射信号捕获跟踪�？�(10)用于对直射数字中频信号sD_IF(n)进行捕获、跟踪及导航电文的解调，得到载波频率控制字F ，反射卫星高度角ε，反-直延迟值N，直-基正交sD_BB_P(n)，直-基同相sD_BB_e(n);其中，反射卫星高度角ε输出给高度测量�？�(30);载波频率控制字Fcw、直-基正交sD_BB_P(n)、直-基同相sD_BB_c(n)、反-直延迟值N分别输出给阵列乘法器模块(20)； 阵列乘法器�？�(20)利用多路延迟器阵列对反射数字中频信号sK_IF(n)、直-基正交sD-BB-P(η)和直-基同相sD_BB_。(η)进行处理，得到直-基正交互相关功率值Xp( τ η)和直-基同相互相关功率值XJ τη)，并输出到高度测量�？�(30)； 高度测量模块(30)包括反-基正交功率值曲线拟合模块(304)、反-基同相功率值曲线拟合�？�(305)、Ρ (Y)码峰值检测模块(302)、C/A码峰值检测模块(303)、高度计算�？�(301)； 反-基正交功率值曲线拟合�？�(304)将多个反-基正交功率值，即反-基正交初始功率值Xpftl)、反-基正交第一功率值XP(1)、反-基正交第二功率值ΧΡ(2)、……、反-基正交第M功率值ΧΡ(Μ)按照散射信号相关功率模型<|YS(A，f) |2>进行曲线拟合，输出P (Y)码的一维功率值-时延函数ΧΡ(τ)。
反-基同相功率值曲线拟合�？�(305)将多个反-基同相功率值，即反-基同相初始功率值Χ_)、反-基同相第一功率值Xc(1)、反-基同相第二功率值Xc(2)……、反-基同相第M功率值Xew按照散射信号相关功率模型〈IYS (△，f) 12>进行曲线拟合，输出C码的一维功率值-时延函数X&)。
P (Y)码峰值检测模(302)块对P (Y)码的一维功率值-时延函数ΧΡ(τ)进行采用峰值跟踪法，拟合出峰值点，计算出峰值点对应的时延值，再将其转化为反射信号和直射信号之间的距离差P ρ。
C/A码峰值检测�？�(303)对C/A码的一维功率值-时延函数进行采用峰值跟踪法，拟合出峰值点，计算出峰值点对应的时延值，再将其转化为反射信号和直射信号之间的距离差P co 高度计算�？�(301)对多路P (Y)码相位差计算直射信号和反射信号的距离差Pp、多路码相位差计算直射信号和反射信号的距离差P。、反射卫星高度角ε进行处理，根据海面反射的空间几何关系，信号处理平台与反射点的相对高程的计算公式为
2.根据权利要求1所述的多路GPS信号P(Y)码相位差的测量系统，其特征在于阵列乘法器�？�(20)对反射数字中频信号sK_IF(n)、直-基正交sD_BB_P(n)和直-基同相sD_BB_c(n)的处理步骤如下步骤201 :载波NCO根据载波频率控制字Few产生对应频率的正弦波和余弦波，产生的波形取决于ROM中存放的是正弦表还是余弦表； 步骤202 :正弦乘法器中将反射数字中频信号sK_IF(n)和正弦表输出的正弦信号进行相乘，输出反射数字基频信号的正交分量，即反-基正交sK_P(n)； 余弦乘法器将反射数字中频信号sK_IF(n)和余弦表输出的余弦信号进行相乘，输出反射数字基频信号的同相分量，即反-基同相sK_c(n)。
步骤203 :延迟器P将直-基正交sD_BB_P(n)延迟N个时间单位，即信号将转换成延迟直-基正交sD_BB_P(n-NA) ;N是直射信号捕获跟踪�？槭涑龅牟捎弥鄙浔栈贩瓷淇返姆绞蕉苑瓷湫藕畔喽杂谥鄙湫藕诺难映傧辔坏墓兰浦担魇茄映俚氖奔洳匠ぃ煌砜傻醚映倨鱌l将延迟直-基正交sD_BB_P(n_NA)延迟I个时间步长Λ，即信号将转换成第一延迟直-基正交sD_BB_P(η-ΝΔ-Δ)； 延迟器ρ2将第一延迟直-基正交sD_BB_P(n-NA-A)延迟I个时间步长Λ，即信号将转换成第二延迟直-基正交s請_ρ(η-ΝΔ-2Δ)；……； 延迟器PM将第M-1延 迟直-基正交81),_1)(11州八-(1-1)八)延迟I个时间步长Λ，ΒΡ信号将转换成第M延迟直-基正交sD_BB_P(η-ΝΔ-ΜΔ)； 延迟器c将直-基同相sD_BB_c(n)延迟N个时间单位，即信号将转换成延迟直-基同相sD_BB_c(n-NA) #是直射信号捕获跟踪�？槭涑龅牟捎弥鄙浔栈贩瓷淇返姆绞蕉苑瓷湫藕畔喽杂谥鄙湫藕诺难映傧辔坏墓兰浦担茄映俚氖奔洳匠ぃ煌砜傻� 延迟器Cl将延迟直-基同相sD-BB-c(n-NA)延迟I个时间步长Λ，即信号将转换成第一延迟直-基同相s請_c(η-ΝΔ-Δ)； 延迟器c2将第一延迟直-基同相sD_BB_c(n-NA-A)延迟I个时间步长Λ，即信号将转换成第二延迟直-基同相sD_BB_c(η-ΝΔ-2Δ)；……； 延迟器CM将第M-1延迟直-基同相sD_BB_c(n-NA-(M_l) Λ)延迟I个时间步长Λ，即信号将转换成第M延迟直-基同相sD_BB_c(η-ΝΔ-ΜΔ)； 步骤204 :乘法器p0将反-基正交sK_P(n)与延迟直-基正交sD_BB_P(n-NA )进行相乘，得到反-基正交初始互相关值SK_D_P(CI);同理可得 乘法器Pl将反-基正交sK_P (η)与第一延迟直-基正交sD_BB_P (η-Ν Δ-Δ)进行相乘，得到反-基正交第一互相关值 Sr-D-P(I); 乘法器Ρ2将反-基正交sK_P(n)与第二延迟直-基正交sD_BB_P(n-NA -2 Δ )进行相乘，得到反-基正交第二互相关值 Sr-D-P (2);......; 乘法器PM将反-基正交sK_P (η)与第M延迟直-基正交sD_BB_P (η-Ν Δ-ΜΔ)进行相乘，得到反-基正交第M互相关值Sr-D-P(M); 乘法器CO将反-基同相sK_c (η)与延迟直-基同相sD_BB_c(n-NA)进行相乘，得到反-基同相初始互相关值 Sr-D-C(O); 乘法器Cl将反-基同相sK-c (η)与第一延迟直-基同相sD_BB_c (η-Ν Δ-Δ)进行相乘，得到反-基同相第一互相关值 Sr-D-C(I); 乘法器c2将反-基同相sK_c(n)与第二延迟直-基同相sD_BB_c(n-NA-2A)进行相乘，得到反-基同相第二互相关值 Sr-D-C (2);......; 乘法器CM将反-基同相sK_c (η)与第M延迟直-基同相sD_BB_c (η-Ν Δ-ΜΔ)进行相乘，得到反-基同相第M互相关值Sk_d_c(M)； 步骤205 :积分累加单元PO将反-基正交初始互相关值SK_D_P(CI)在预检积分时间内进行累加，得到未截断反-基正交初始功率值YP(CI)； 积分累加单元Pl将反-基正交第一互相关值SK_D_P(1)在预检积分时间内进行累加，得到未截断反-基正交第一功率值YP(1)； 积分累加单元Ρ2将反-基正交第二互相关值SK_D_P⑵在预检积分时间内进行累加，得到未截断反-基正交第二功率值YP(2)；……； 积分累加单元PM将反-基正交第M互相关值SK_D_P(M)在预检积分时间内进行累加，得到未截断反-基正交第M功率值YP(M)； 积分累加单元CO将反-基同相初始互相关值SK_D__)在预检积分时间内进行累加，得到未截断反-基同相初始功率值; 积分累加单元Cl将反-基同相第一互相关值Sk_m(1)在预检积分时间内进行累加，得到未截断反-基同相第一功率值Yc(1); 积分累加单元c2将反-基同相第二互相关值SK_D_。⑵在预检积分时间内进行累加，得到未截断反-基同相第二功率值Yc(2)；……； 积分累加单元CM将反-基同相第M互相关值Sk_m(m)在预检积分时间内进行累加，得到未截断反-基同相第M功率值Ye(M)； 步骤206 :截断器pO将未截断反-基正交初始功率值Ypftl)进行12比特的截断处理，得到反-基正交初始功率值Xpftl); 截断器Pl将未截断反-基正交第一功率值YP(1)进行12比特的截断处理，得到反-基正交第一功率值XP(1); 截断器p2将未截断反-基正交第二功率值Yp⑵进行12比特的截断处理，得到反-基正交第二功率值Xp⑵;……； 截断器PM将未截断反-基正交第M功率值YP(M)进行12比特的截断处理，得到反-基正交第M功率值XP(M)； 截断器CO将未截断反-基同相初始功率值Y_)进行12比特的截断处理，得到反-基同相初始功率值Xc(O); 截断器Cl将未截断反-基同相第一功率值Yc(1)进行12比特的截断处理，得到反-基同相第一功率值Xca)； 截断器c2将未截断反-基同相第二功率值Y。⑵进行12比特的 截断处理，得到反-基同相第二功率值Xc(2)；……； 截断器CM将未截断反-基同相第M功率值Y_)进行12比特的截断处理，得到反-基同相第M功率值Xc(Μ)。
3.根据权利要求1所述的多路GPS信号P (Y)码相位差的测量系统，其特征在于该系统通过在FPGA芯片上用Verilog HDL语言设计得到硬件的功效。
全文摘要
本发明公开了一种多路GPS信号的P(Y)码相位差的测量系统，该系统包括有直射信号捕获跟踪模块、阵列乘法器�？�、高度测量�？�。依据GPS的正交调制特性提取接收到的卫星信号的正交和同相信号，分离出C/A码和P(Y)码，然后分别对其进行等间距相同的延迟，再将其延迟后的信号和反射信号进行互相关运算，不同的延迟码片导致不同的非相干积分值，得到不同路GPS信号的P(Y)码相位差。该测量系统降低了硬件接收设备的复杂度，同时测量精度也得到了提高。
文档编号G01C5/00GK103017731SQ20121048024
公开日2013年4月3日 申请日期2012年11月22日 优先权日2012年11月22日
发明者李伟强, 杨东凯, 秦瑾, 张波, 张帅 申请人:北京航空航天大学
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
