Changeset 4270 in ntrip


Ignore:
Timestamp:
Jun 22, 2012, 12:40:24 PM (12 years ago)
Author:
mervart
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/rinex/reqcanalyze.cpp

    r4269 r4270  
    156156    QString          prn     = it.key();
    157157    const t_satStat& satStat = it.value();
    158 
    159     if (satStat.MP1.size() > 1) {
    160       double mean = 0.0;
    161       for (int ii = 0; ii < satStat.MP1.size(); ii++) {
    162         mean += satStat.MP1[ii];
     158 
     159    for (int im = 1; im <= 2; im++) {
     160      const QVector<double>& MP = (im == 1) ? satStat.MP1 : satStat.MP2;
     161      if (MP.size() > 1) {
     162        double mean = 0.0;
     163        for (int ii = 0; ii < MP.size(); ii++) {
     164          mean += MP[ii];
     165        }
     166        mean /= MP.size();
     167        double stddev = 0.0;
     168        for (int ii = 0; ii < MP.size(); ii++) {
     169          double diff = MP[ii] - mean;
     170          stddev += diff * diff;
     171        }
     172        double multipath = sqrt(stddev / (MP.size()-1));
     173     
     174        *_log << "MP" << im << " " << prn << " " << multipath << endl;
    163175      }
    164       mean /= satStat.MP1.size();
    165       double stddev = 0.0;
    166       for (int ii = 0; ii < satStat.MP1.size(); ii++) {
    167         double diff = satStat.MP1[ii] - mean;
    168         stddev += diff * diff;
    169       }
    170       double MP1 = sqrt(stddev / (satStat.MP1.size()-1));
    171 
    172       *_log << "MP1 " << prn << " " << MP1 << endl;
    173176    }
    174177
Note: See TracChangeset for help on using the changeset viewer.