Changeset 4269 in ntrip

Timestamp:

Jun 22, 2012, 12:31:03 PM (12 years ago)

Author:

mervart

Message:

 

File:

• trunk/BNC/rinex/reqcanalyze.cpp (modified) (1 diff)

trunk/BNC/rinex/reqcanalyze.cpp

@@ -156 +156 @@
     QString          prn     = it.key();
     const t_satStat& satStat = it.value();
-    if (satStat.MP1.size()) {
+
+    if (satStat.MP1.size() > 1) {
+      double mean = 0.0;
       for (int ii = 0; ii < satStat.MP1.size(); ii++) {
-        *_log << "MP1 " << prn << " " << satStat.MP1[ii] << endl;
+        mean += satStat.MP1[ii];
       }
-    }
-    if (satStat.MP2.size()) {
-      for (int ii = 0; ii < satStat.MP2.size(); ii++) {
-        *_log << "MP2 " << prn << " " << satStat.MP2[ii] << endl;
+      mean /= satStat.MP1.size();
+      double stddev = 0.0;
+      for (int ii = 0; ii < satStat.MP1.size(); ii++) {
+        double diff = satStat.MP1[ii] - mean;
+        stddev += diff * diff;
       }
-    }
+      double MP1 = sqrt(stddev / (satStat.MP1.size()-1));
+
+      *_log << "MP1 " << prn << " " << MP1 << endl;
+    }
+
   }