def search(x, chan, doppler, ca_code_phase):
n = int(fs * 0.004)
w = nco.nco(-(562500 * chan + doppler) / fs, 0, n)
m_metric, m_k = 0, 0
for k in range(1000):
q = 0
cp = 5110 * k + 10 * ca_code_phase
for block in range(20):
incr = 5110000.0 / fs
c = p.code(0, cp, incr, n)
xp = x[n * block:n * (block + 1)] * c * w
q = q + np.absolute(np.sum(xp))
cp += n * incr
print('%f %f' % (k, q))
if q > m_metric:
m_metric = q
m_k = k
return m_metric, m_k
def search(x,chan,doppler,ca_code_phase):
n = int(fs*0.004)
w = nco.nco(-(437500*chan+doppler)/fs,0,n)
m_metric,m_k = 0,0
for k in range(1000):
q = 0
cp = 5110*k + 10*ca_code_phase
for block in range(20):
incr = 5110000.0/fs
c = p.code(0,cp,incr,n)
xp = x[n*block:n*(block+1)]*c*w
q = q + np.absolute(np.sum(xp))
cp += n*incr
print('%f %f'%(k,q))
if q>m_metric:
m_metric = q
m_k = k
return m_metric,m_k