Changeset 2934 in ntrip for trunk/BNC/combination


Ignore:
Timestamp:
Jan 29, 2011, 8:08:19 PM (14 years ago)
Author:
mervart
Message:
 
Location:
trunk/BNC/combination
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/combination/bnccomb.cpp

    r2933 r2934  
    4646double cmbParam::partial(const QString& acIn, t_corr* corr) {
    4747 
    48   if      (type == offset) {
     48  if      (type == AC_offset) {
     49    if (AC == acIn) {
     50      return 1.0;
     51    }
     52  }
     53  else if (type == Sat_offset) {
    4954    if (AC == acIn && prn == corr->prn) {
    5055      return 1.0;
     
    99104      it.next();
    100105      cmbAC* AC = it.value();
    101       _params.push_back(new cmbParam(cmbParam::offset, ++nextPar, AC->name,prn));
     106      if (iGps == 1) {
     107        _params.push_back(new cmbParam(cmbParam::AC_offset, ++nextPar, AC->name, ""));
     108      }
     109      _params.push_back(new cmbParam(cmbParam::Sat_offset, ++nextPar, AC->name, prn));
    102110    }
    103111  }
     
    107115  _QQ = 0.0;
    108116
    109   _sigClk0 = 100.0;
    110   _sigOff0 = 100.0;
     117  _sigACOff  = 100.0;
     118  _sigSatOff = 100.0;
     119  _sigClk    = 100.0;
    111120
    112121  for (int iPar = 1; iPar <= _params.size(); iPar++) {
    113122    cmbParam* pp = _params[iPar-1];
    114     if      (pp->type == cmbParam::clk) {
    115       _QQ(iPar,iPar) = _sigClk0 * _sigClk0;
    116     }
    117     else if (pp->type == cmbParam::offset) {
    118       _QQ(iPar,iPar) = _sigOff0 * _sigOff0;
     123    if      (pp->type == cmbParam::AC_offset) {
     124      _QQ(iPar,iPar) = _sigACOff * _sigACOff;
     125    }
     126    else if (pp->type == cmbParam::Sat_offset) {
     127      _QQ(iPar,iPar) = _sigSatOff * _sigSatOff;
     128    }
     129    else if (pp->type == cmbParam::clk) {
     130      _QQ(iPar,iPar) = _sigClk * _sigClk;
    119131    }
    120132  }
     
    314326  for (int iPar = 1; iPar <= _params.size(); iPar++) {
    315327    cmbParam* pp = _params[iPar-1];
    316     if      (pp->type == cmbParam::clk) {
     328    if  (pp->type == cmbParam::AC_offset || pp->type == cmbParam::clk) {
    317329      pp->xx = 0.0;
    318330      for (int jj = 1; jj <= _params.size(); jj++) {
    319331        _QQ(iPar, jj) = 0.0;
    320332      }
    321      _QQ(iPar,iPar) = _sigClk0 * _sigClk0;
     333    }
     334    if      (pp->type == cmbParam::AC_offset) {
     335     _QQ(iPar,iPar) = _sigACOff * _sigACOff;
     336    }
     337    else if (pp->type == cmbParam::clk) {
     338     _QQ(iPar,iPar) = _sigClk * _sigClk;
    322339    }
    323340  }
  • trunk/BNC/combination/bnccomb.h

    r2933 r2934  
    1010class cmbParam {
    1111 public:
    12   enum parType {offset, clk};
     12  enum parType {AC_offset, Sat_offset, clk};
    1313  cmbParam(parType typeIn, int indexIn,
    1414           const QString& acIn, const QString& prnIn);
     
    7979  cmbCaster*            _caster;
    8080  QVector<cmbParam*>    _params;
    81   double                _sigClk0;
    82   double                _sigOff0;
     81  double                _sigACOff;
     82  double                _sigSatOff;
     83  double                _sigClk;
    8384  SymmetricMatrix       _QQ;
    8485};
Note: See TracChangeset for help on using the changeset viewer.