Changeset 6025 in ntrip for trunk/BNC/src/PPP


Ignore:
Timestamp:
Aug 21, 2014, 10:42:42 AM (10 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/PPP/pppSatObs.cpp

    r6024 r6025  
    151151  }
    152152}
     153
     154//
     155////////////////////////////////////////////////////////////////////////////
     156void t_pppSatObs::lcCoeff(t_lc::type tLC,
     157                          map<t_frequency::type, double>& codeCoeff,
     158                          map<t_frequency::type, double>& phaseCoeff) const {
     159
     160  codeCoeff.clear();
     161  phaseCoeff.clear();
     162
     163  t_frequency::type fType1 = t_lc::toFreq(_prn.system(),t_lc::l1);
     164  t_frequency::type fType2 = t_lc::toFreq(_prn.system(),t_lc::l2);
     165
     166  double f1 = t_CST::freq(fType1, _channel);
     167  double f2 = t_CST::freq(fType2, _channel);
     168
     169  switch (tLC) {
     170  case t_lc::l1:
     171    phaseCoeff[fType1] = 1.0; 
     172    return;
     173  case t_lc::l2: 
     174    phaseCoeff[fType2] = 1.0; 
     175    return;
     176  case t_lc::lIF:
     177    phaseCoeff[fType1] =  f1 * f1 / (f1 * f1 - f2 * f2);
     178    phaseCoeff[fType2] = -f2 * f2 / (f1 * f1 - f2 * f2);
     179    return;
     180  case t_lc::MW: 
     181    phaseCoeff[fType1] =  f1 / (f1 - f2);
     182    phaseCoeff[fType2] = -f2 / (f1 - f2);
     183    codeCoeff[fType1]  = -f1 / (f1 + f2);
     184    codeCoeff[fType2]  = -f2 / (f1 + f2);
     185    return;
     186  case t_lc::CL: 
     187    phaseCoeff[fType1] =  0.5;
     188    codeCoeff[fType1]  =  0.5;
     189    return;
     190  case t_lc::c1: 
     191    codeCoeff[fType1] = 1.0; 
     192    return;
     193  case t_lc::c2: 
     194    codeCoeff[fType2] = 1.0; 
     195    return;
     196  case t_lc::cIF:
     197    codeCoeff[fType1] =  f1 * f1 / (f1 * f1 - f2 * f2);
     198    codeCoeff[fType2] = -f2 * f2 / (f1 * f1 - f2 * f2);
     199    return;
     200  case t_lc::dummy:
     201  case t_lc::maxLc:
     202    return;
     203  }
     204}
     205
    153206
    154207//
Note: See TracChangeset for help on using the changeset viewer.