[发明专利]一种任意时刻P码发生方法

权利要求书

1.一种任意时刻P码发生方法，其特征在于，所述任意时刻P码发生方法包括如下步骤：

步骤A，建立移位寄存器状态表：在P码发生器中，将X1A、X1B、X2A和X2B四路移位寄存器各自一个周期内相位对应的所有状态值存储为一个状态表；该状态表以寄存器相位为索引，存储对应的寄存器状态值；

步骤B，设置P码发生器参数：设置卫星星号，及X1A、X1B、X2A和X2B四路移位寄存器的初始状态；本步骤通过对时间辅助信息的解算得到移位寄存器的相位，由该相位作为索引，在预存储的寄存器状态表中读取对应状态值，该步骤分为以下5个子步骤：

步骤B-1，由于P码发生周期为一周即七天，一周内P码的相位均与时间一一对应，P码发生器读取本地时间信息，转换成时间信息控制字，其中包含相对于GPS周起始的相对时间t及产生的码元总数N；

步骤B-2，X1、X2序列的相位计算：读取时间信息控制字，计算X1、X2序列的相位；

步骤B-3，移位寄存器周期数及相位计算：根据X1、X2序列的相位，对X1A、X1B、X2A和X2B四路移位寄存器的周期数与相位进行计算；

步骤B-4，读取移位寄存器状态：利用四路移位寄存器的相位，索引本地存储的寄存器状态表，读取当前相位对应的各移位寄存器的状态值；

步骤B-5，设置移位寄存器状态：利用读取的各移位寄存器状态值，设置当前移位寄存器的状态；

步骤C，产生X1、X2序列：发出控制信号使四路移位寄存器X1A、X1B及X2A、X2B同时产生相应的码序列，将X1A、X1B移位寄存器输出序列异或相加得到X1序列，X2A、X2B移位寄存器输出序列异或相加得到X2序列；

步骤D，P码序列产生：加载卫星号，作为X2序列的延迟长度，将X2序列延迟得到X2i序列，X1序列和X2i序列进行模二加运算，产生当前时刻的P码。

2.如权利要求1所述的一种任意时刻P码发生方法，其特征在于：所述步骤A中，四路移位寄存器X1A、X1B、X2A和X2B，其特征多项式分别如下：

X1A:1+x⁶+x⁸+x¹¹+x¹²

X1B:1+x¹+x²+x⁵+x⁸+x⁹+x¹⁰+x¹¹+x¹²

X2A:1+x¹+x³+x⁴+x⁵+x⁷+x⁸+x⁹+x¹⁰+x¹¹+x¹²

X2B:1+x²+x³+x⁴+x⁸+x⁹+x¹²。

3.如权利要求1所述的一种任意时刻P码发生方法，其特征在于：所述步骤B-1中，所述时间信息控制字计算方法如下：

t=d*24*60*60+h*60*60+m*60+s+i/1000

N＝t*10.23*10⁶

其中，d为本周经历天数，h代表小时，m代表分钟，s代表秒，i表示毫秒。

4.如权利要求1所述的一种任意时刻P码发生方法，其特征在于：所述步骤B-2中，所述X1、X2序列的相位计算，其计算公式如下：

N_X1=N%15345000

N_X2=N%15345037

其中，N_X1、N_X1为X1、X2序列的当前相位，%为取余数运算。

5.如权利要求4所述的一种任意时刻P码发生方法，其特征在于：所述步骤B-3中，移位寄存器相位计算步骤如下：

步骤B-3-1，移位寄存器周期数计算，其公式为：

X1A_counter=int(N_X1/4092)

X1B_counter=int(N_X1/4093)

X2A_counter=int(N_X2/4092)

X2B_counter=int(N_X2/4093)

X1A_counter、X1B_counter、X2A_counter、X2B_counter分别为四路移位寄存器当前的周期数，int(*)为取整运算；

步骤B-3-2，移位寄存器相位计算，其公式为：

X1A_shift=N_X1%4092

X1B_shift=N_X1%4093X1B_counter<37490X1B_counter≥3749]]>

X2A_shift=N_X2%4092X2A_counter<37500X2A_counter≥3750]]>

X2B_shift=N_X2%4093X2B_counter<37490X2B_counter≥3749]]>

X1A_shift，X1B_shift,X2A_shift,X2B_shift分别为四路移位寄存器当前的相位值。