Changeset 4270 in ntrip
- Timestamp:
- Jun 22, 2012, 12:40:24 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/rinex/reqcanalyze.cpp
r4269 r4270 156 156 QString prn = it.key(); 157 157 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; 163 175 } 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;173 176 } 174 177
Note:
See TracChangeset
for help on using the changeset viewer.