Changeset 6025 in ntrip for trunk/BNC/src/PPP/pppSatObs.cpp

Timestamp:

Aug 21, 2014, 10:42:42 AM (10 years ago)

Author:

mervart

Message:

 

File:

• trunk/BNC/src/PPP/pppSatObs.cpp (modified) (1 diff)

trunk/BNC/src/PPP/pppSatObs.cpp

@@ -151 +151 @@
   }
 }
+
+// 
+////////////////////////////////////////////////////////////////////////////
+void t_pppSatObs::lcCoeff(t_lc::type tLC, 
+                          map<t_frequency::type, double>& codeCoeff,
+                          map<t_frequency::type, double>& phaseCoeff) const {
+
+  codeCoeff.clear();
+  phaseCoeff.clear();
+
+  t_frequency::type fType1 = t_lc::toFreq(_prn.system(),t_lc::l1);
+  t_frequency::type fType2 = t_lc::toFreq(_prn.system(),t_lc::l2);
+
+  double f1 = t_CST::freq(fType1, _channel);
+  double f2 = t_CST::freq(fType2, _channel);
+
+  switch (tLC) {
+  case t_lc::l1:
+    phaseCoeff[fType1] = 1.0;  
+    return;
+  case t_lc::l2:  
+    phaseCoeff[fType2] = 1.0;  
+    return;
+  case t_lc::lIF: 
+    phaseCoeff[fType1] =  f1 * f1 / (f1 * f1 - f2 * f2);
+    phaseCoeff[fType2] = -f2 * f2 / (f1 * f1 - f2 * f2);
+    return;
+  case t_lc::MW:  
+    phaseCoeff[fType1] =  f1 / (f1 - f2);
+    phaseCoeff[fType2] = -f2 / (f1 - f2);
+    codeCoeff[fType1]  = -f1 / (f1 + f2);
+    codeCoeff[fType2]  = -f2 / (f1 + f2);
+    return;
+  case t_lc::CL:  
+    phaseCoeff[fType1] =  0.5;
+    codeCoeff[fType1]  =  0.5;
+    return;
+  case t_lc::c1:  
+    codeCoeff[fType1] = 1.0;  
+    return;
+  case t_lc::c2:  
+    codeCoeff[fType2] = 1.0;  
+    return;
+  case t_lc::cIF: 
+    codeCoeff[fType1] =  f1 * f1 / (f1 * f1 - f2 * f2);
+    codeCoeff[fType2] = -f2 * f2 / (f1 * f1 - f2 * f2);
+    return;
+  case t_lc::dummy: 
+  case t_lc::maxLc: 
+    return;
+  }
+}
+
 
 //