找到校验不通过的原因:
将相对码还原到数据码设置的实现出错
从书本(《全球导航卫星系统原理》,P172)中可以看到,GLONASS将数据码转变为相对码的公式为:
c_i = b_i xor c_i - 1 (i=2,3,4,…,85,而c_i=0)
据此反推将相对码还原到数据码的公式:
b_i = c_i xor c_i-1
这里b_1同样为0
实际实现过程中,先设置i初值为2,将c_2与c_1异或得到b_2,然后i++,逐步得到i=2,3,4,…,85的b_i。然而,由于相对码和数据码使用的是同一个数组,导致迭代过程中b_i的值覆盖了c_i的值,产生错误,得到的数据码与真实数据码不一致,导致其不能通过校验。改变实现方法,先设i的初值为85,然后向前逐步迭代得到i=85,84,83,…,2的b_i。
通过python仿真验证了上述过程,随机生成长度为85的0/1数组,使用向后迭代方法,得到的数组与原始不一致,而向前迭代则完全一致,重复10000次该过程,验证上述结果。
附上相应代码
import numpy as np
def xor(a, b):
if a == b:
return 0
else:
return 1
cnt = 0
N = 1000
for i0 in range(N):
array0 = np.random.randint(0,2,85)
array0[0] = 0
array1 = np.zeros(85)
for i in range(1, 85): # 变换得到相对码 c_i = b_i xor c_i-1
array1[i] = xor(array1[i-1], array0[i])
array2 = array1.copy()
array3 = array1.copy()
for i in range(1, 85): # 尝试反变换回去
array2[i] = xor(array2[i], array2[i-1]) # 向后迭代
array3[85-i] = xor(array3[85-i], array3[84-i]) #向前迭代
if False in (array0 == array2).tolist(): # 与原始数组比较,需要验证array3时相应修改即可
print("i=%d" % i0)
print(array0)
print(array2)
else:
cnt += 1
print("总比较成功次数:%d" % cnt)