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()