Changeset 10168 in ntrip


Ignore:
Timestamp:
Aug 22, 2023, 6:39:48 PM (10 months ago)
Author:
stuerze
Message:
 
Location:
trunk/BNC/src/PPP
Files:
2 edited

Legend:

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

    r10167 r10168  
    268268      }
    269269      if (par) {
    270 //        if (par->ambResetCandidate()) {
     270        if (par->ambResetCandidate()) {
    271271          resetAmb(par->prn(), obsVector, maxOutlierLC, &QSav, &xSav);
    272 //        }
    273 //        else {
    274 //          par->setAmbResetCandidate();
    275 //          obs->setOutlier();
    276 //        }
     272        }
     273        else {
     274          par->setAmbResetCandidate();
     275          obs->setOutlier();
     276        }
     277        addNoiseToPar(t_pppParam::ion, par->prn(), 0.001);
    277278      }
    278279      else {
     
    391392      }
    392393      t_lc::type tLC = par->tLC();
    393 //      if (tLC != lc) {continue;}
     394      if (tLC != lc) {continue;}
    394395      LOG << string(_epoTime) << " RESET " << par->toString() << endl;
    395396      delete par; par = new t_pppParam(t_pppParam::amb, prn, tLC, &obsVector);
     
    413414      }
    414415      _x0[ind] = par->x0();
     416      irc = success;
     417    }
     418  }
     419
     420  return irc;
     421}
     422
     423// Add noise to individual parameter
     424////////////////////////////////////////////////////////////////////////////
     425t_irc t_pppFilter::addNoiseToPar(t_pppParam::e_type parType, t_prn prn, double noise) {
     426  t_irc irc = failure;
     427  vector<t_pppParam*> &params = _parlist->params();
     428  for (unsigned iPar = 0; iPar < params.size(); iPar++) {
     429    t_pppParam *par = params[iPar];
     430    if (par->type() == parType && par->prn() == prn) {
     431      int ind = par->indexNew();
     432      LOG << string(_epoTime) << " ADD NOISE TO " << par->toString() << endl;
     433      par->setIndex(ind);
     434      _QFlt(ind + 1, ind + 1) = noise * noise;
    415435      irc = success;
    416436    }
  • trunk/BNC/src/PPP/pppFilter.h

    r10034 r10168  
    8383                 SymmetricMatrix* QSav = 0, ColumnVector* xSav = 0);
    8484
     85  t_irc addNoiseToPar(t_pppParam::e_type parType, t_prn prn, double noise);
     86
    8587  void cmpDOP(const std::vector<t_pppSatObs*>& obsVector);
    8688
Note: See TracChangeset for help on using the changeset viewer.