Changeset 10010 in ntrip


Ignore:
Timestamp:
Mar 17, 2023, 2:27:30 PM (13 months ago)
Author:
stuerze
Message:

minor changes

Location:
trunk/BNC
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/BNC/scripts/pppPlot.pl

    r10005 r10010  
    198198                                                    title   => "Tropospheric Delay, RMS + -$str_rms_trp m",
    199199                                                    timefmt => '%s',
    200                                                     style   => "linespoints",
     200                                                    style   => "points",
    201201        );
    202202        $chartTRP->plot2d($dataTRP);
    … …  
    228228                                                     title   => "Receiver clock",
    229229                                                     timefmt => '%s',
    230                                                      style   => "linespoints",
     230                                                     style   => "points",
    231231            );
    232232            $chartCLK->plot2d($dataset);
    … …  
    261261                                                         title   => "Receiver Offset GLONASS",
    262262                                                         timefmt => '%s',
    263                                                          style   => "linespoints",
     263                                                         style   => "points",
    264264                );
    265265                $chartOFF_GLO->plot2d($dataset);    #system ("display $pngName&");
    … …  
    294294                                                         title   => "Receiver Offset Galileo",
    295295                                                         timefmt => '%s',
    296                                                          style   => "linespoints",
     296                                                         style   => "points",
    297297                );
    298298                $chartOFF_GAL->plot2d($dataset);    #system ("display $pngName&");
    … …  
    327327                                                         title   => "Receiver Offset Beidou",
    328328                                                         timefmt => '%s',
    329                                                          style   => "linespoints",
     329                                                         style   => "points",
    330330                );
    331331                $chartOFF_BDS->plot2d($dataset);    #system ("display $pngName&");
    … …  
    376376                    timefmt => '%s',
    377377
    378                     #style   => "linespoints",
    379                     style => "linespoints",
     378                    #style   => "points",
     379                    style => "points",
    380380                );
    381381                push ( @datasets, $dataset );
    … …  
    451451                        title => "$key_sat",
    452452                        timefmt => '%s',
    453                         style   => "linespoints",
    454 
    455                         #style   => "linespoints",
     453                        style   => "points",
     454
     455                        #style   => "points",
    456456                    );
    457457                    push ( @datasets_amb, $dataset_amb );
    … …  
    464464                        timefmt => '%s',
    465465
    466                         #style   => "linespoints",
    467                         style => "linespoints",
     466                        #style   => "points",
     467                        style => "points",
    468468                    );
    469469                    push ( @datasets_epo, $dataset_epo );
    … …  
    530530                        timefmt => '%s',
    531531
    532                         #style   => " linespoints ",
    533                         style => "linespoints",
     532                        #style   => " points ",
     533                        style => "points",
    534534                    );
    535535                    push ( @datasets, $dataset );
    … …  
    580580                                                    title   => "$key_sys",
    581581                                                    timefmt => '%s',
    582                                                     style   => "linespoints",
     582                                                    style   => "points",
    583583                      );
    584584                    push ( @datasets, $dataset );
    … …  
    643643                        title   => "$key_sat",
    644644                        timefmt => '%s',
    645 
    646                         #style   => "linespoints",
    647                         style => "linespoints",
     645                        style => "points",
    648646                    );
    649647                    push ( @datasets, $dataset );

  • trunk/BNC/src/PPP/pppFilter.cpp

    r10009 r10010  
    7272
    7373  const QMap<char, t_pppRefSat*> &refSatMap = epoch->refSatMap();
    74 
     74  const QList<char> &usedSystems = _parlist.usedSystems();
    7575  //--
    7676  // Set Parameters
    … …  
    7878    return failure;
    7979  }
    80 
     80#ifdef BNC_DEBUG_PPP
    8181  if (OPT->_obsModelType == OPT->DCMcodeBias ||
    82   OPT->_obsModelType == OPT->DCMphaseBias) {
    83 #ifdef BNC_DEBUG_PPP
    84   _parlist.printParams(_epoTime);
     82      OPT->_obsModelType == OPT->DCMphaseBias) {
     83    _parlist.printParams(_epoTime);
     84  }
    8585#endif
    86   }
     86
    8787  // Status Vector, Variance-Covariance Matrix
    8888  // -----------------------------------------
    8989  ColumnVector xFltOld = _xFlt;
    9090  SymmetricMatrix QFltOld = _QFlt;
    91   setStateVectorAndVarCovMatrix(xFltOld, QFltOld);
     91  bool setNeuNoiseToZero = false;
     92  for (int iSys = 0; iSys < usedSystems.size(); iSys++) {
     93    char sys = usedSystems[iSys];
     94    int num = 0;
     95    for (unsigned jj = 0; jj < allObs.size(); jj++) {
     96      if (allObs[jj]->prn().system() == sys) {
     97        ++num;
     98      }
     99    }
     100    if (num < OPT->_minObs) {
     101      setNeuNoiseToZero = true;
     102    }
     103  }
     104  setStateVectorAndVarCovMatrix(xFltOld, QFltOld, setNeuNoiseToZero);
    92105
    93106  // Pre-Process Satellite Systems separately
    … …  
    97110  OPT->_obsModelType == OPT->DCMphaseBias) {
    98111    preProcessing = true;
    99     const QList<char> &usedSystems = _parlist.usedSystems();
    100112    for (int iSys = 0; iSys < usedSystems.size(); iSys++) {
    101113      char sys = usedSystems[iSys];
    … …  
    137149  // ------------------------------------
    138150  preProcessing = false;
    139   const QList<char> &usedSystems = _parlist.usedSystems();
    140151  for (int iSys = 0; iSys < usedSystems.size(); iSys++) {
    141152    char sys = usedSystems[iSys];
    … …  
    594605//
    595606////////////////////////////////////////////////////////////////////////////
    596 void t_pppFilter::predictCovCrdPart(const SymmetricMatrix &QFltOld) {
     607void t_pppFilter::predictCovCrdPart(const SymmetricMatrix &QFltOld, bool setNeuNoiseToZero) {
    597608
    598609  const vector<t_pppParam*> &params = _parlist.params();
    … …  
    622633  } else {
    623634    double dt = _epoTime - _firstEpoTime;
    624     if (dt < OPT->_seedingTime) {
     635    if (dt < OPT->_seedingTime || setNeuNoiseToZero) {
    625636      _QFlt.SymSubMatrix(1, 3) = QFltOld.SymSubMatrix(1, 3);
    626637    } else {
    … …  
    633644////////////////////////////////////////////////////////////////////////////
    634645void t_pppFilter::setStateVectorAndVarCovMatrix(const ColumnVector &xFltOld,
    635     const SymmetricMatrix &QFltOld) {
     646    const SymmetricMatrix &QFltOld, bool setNeuNoiseToZero) {
    636647
    637648  const vector<t_pppParam*> &params = _parlist.params();
    … …  
    666677    }
    667678  }
    668   predictCovCrdPart(QFltOld);
     679  predictCovCrdPart(QFltOld, setNeuNoiseToZero);
    669680}
    670681

  • trunk/BNC/src/PPP/pppFilter.h

    r10008 r10010  
    206206              const QMap<char, t_pppRefSat*>& refSatMap);
    207207
    208   void setStateVectorAndVarCovMatrix(const ColumnVector& xFltOld, const SymmetricMatrix& QFltOld);
    209 
    210   void predictCovCrdPart(const SymmetricMatrix& QFltOld);
     208  void setStateVectorAndVarCovMatrix(const ColumnVector& xFltOld, const SymmetricMatrix& QFltOld,
     209                                     bool setNeuNoiseToZero);
     210
     211  void predictCovCrdPart(const SymmetricMatrix& QFltOld, bool setNeuNoiseToZero);
    211212
    212213  t_irc addNoiseToPar(t_pppParam::e_type parType, char sys);

  • trunk/BNC/src/PPP/pppSatObs.cpp

    r10009 r10010  
    3333using namespace std;
    3434
    35 const double   GLO_WEIGHT_FACTOR = 1.0;
    36 const double   BDS_WEIGHT_FACTOR = 1.0;
     35const double   GLO_WEIGHT_FACTOR = 2.0;
     36const double   BDS_WEIGHT_FACTOR = 2.0;
    3737
    3838// Constructor
    … …  
    355355  // De-Weight GLO+BDS
    356356  // -----------------
    357   if (_prn.system() == 'R') {
     357  if (_prn.system() == 'R' && t_lc::includesCode(tLC)) {
    358358    retVal *= GLO_WEIGHT_FACTOR;
    359359  }
    360   if (_prn.system() == 'C') {
     360  if (_prn.system() == 'C' && t_lc::includesCode(tLC)){
    361361    retVal *= BDS_WEIGHT_FACTOR;
    362362  }
    … …  
    398398
    399399  retVal = sqrt(retVal);
     400
     401  // De-Weight GLO+BDS
     402  // -----------------
     403  if (_prn.system() == 'R' && t_lc::includesCode(tLC)) {
     404    retVal *= GLO_WEIGHT_FACTOR;
     405  }
     406  if (_prn.system() == 'C' && t_lc::includesCode(tLC)){
     407    retVal *= BDS_WEIGHT_FACTOR;
     408  }
     409
    400410
    401411  return retVal;
Note: See TracChangeset for help on using the changeset viewer.