def compareRecievers():
channelDict = Scenario.readNMEAFile()
satelliteDict = Scenario.readSimulatorMotionFile()
startTime = 406800 #second of week
CNO = 48
PLLBW = 2
FLLBW = 10
for SVNum in [1,3,6,11,16,18,19,22,31]:
T,SignalIn,stateHistory = \
GenSignal.createScenarioFromSpirentFile(SVNum,CNO)
coherentIntegrationTime = 4
numOutputSamples = T.size/coherentIntegrationTime
print('Running Reciever')
Z,DcoFreqArray,LockStateValues,PhaseError = \
runReceiver(T,SignalIn,18,FLLBW)
PolarisFreq = DcoFreqArray[500:1000]
plt.plot(PolarisFreq,'r')
plt.ylabel('Doppler Shift (Hz)')
plt.xlabel('Time (s)')
plt.ylim(-4,4)
plt.title('Bn=18Hz')
plt.savefig(str(SVNum)+'18Polaris.eps',format='eps', dpi=1000)
plt.close()
Z,DcoFreqArray,LockStateValues,PhaseError = \
runReceiver(T,SignalIn,32,FLLBW)
PolarisFreq = DcoFreqArray[500:1000]
plt.plot(PolarisFreq,'r')
plt.ylim(-4,4)
plt.ylabel('Doppler Shift (Hz)')
plt.xlabel('Time (s)')
plt.title('Bn=32Hz')
plt.savefig(str(SVNum)+'32Polaris.eps',format='eps', dpi=1000)
plt.close()
