Changeset 2084 in ntrip for trunk/BNC/bncmodel.cpp
- Timestamp:
- Dec 6, 2009, 4:57:58 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/bncmodel.cpp
r2083 r2084 56 56 const double sig_crd_p = 100.0; 57 57 const double sig_clk_0 = 1000.0; 58 const double sig_trp_0 = 0.01; 59 const double sig_trp_p = 1e-6; 58 60 const double sig_amb_0 = 100.0; 59 61 const double sig_P3 = 1.0; 60 const double sig_L3 = 0.01;62 const double sig_L3 = 0.01; 61 63 62 64 // Constructor … … 92 94 return 1.0; 93 95 } 96 else if (type == TROPO) { 97 return 1.0 / sin(satData->eleSat); 98 } 94 99 else if (type == AMB_L3) { 95 100 if (prnIn == prn) { … … 106 111 //////////////////////////////////////////////////////////////////////////// 107 112 bncModel::bncModel() { 108 _xcBanc.ReSize(4); _xcBanc = 0.0; 113 114 bncSettings settings; 115 116 _static = false; 117 if ( Qt::CheckState(settings.value("pppStatic").toInt()) == Qt::Checked) { 118 _static = true; 119 } 120 121 _usePhase = false; 122 if ( Qt::CheckState(settings.value("pppUsePhase").toInt()) == Qt::Checked) { 123 _usePhase = true; 124 } 125 126 _estTropo = false; 127 if ( Qt::CheckState(settings.value("pppEstTropo").toInt()) == Qt::Checked) { 128 _estTropo = true; 129 } 130 131 _xcBanc.ReSize(4); _xcBanc = 0.0; 132 _ellBanc.ReSize(3); _ellBanc = 0.0; 133 109 134 _params.push_back(new bncParam(bncParam::CRD_X, 1, "")); 110 135 _params.push_back(new bncParam(bncParam::CRD_Y, 2, "")); 111 136 _params.push_back(new bncParam(bncParam::CRD_Z, 3, "")); 112 137 _params.push_back(new bncParam(bncParam::RECCLK, 4, "")); 113 _ellBanc.ReSize(3); 138 if (_estTropo) { 139 _params.push_back(new bncParam(bncParam::TROPO, 5, "")); 140 } 114 141 115 142 unsigned nPar = _params.size(); 143 116 144 _QQ.ReSize(nPar); 145 _xx.ReSize(nPar); 117 146 _QQ = 0.0; 147 _xx = 0.0; 118 148 119 149 _QQ(1,1) = sig_crd_0 * sig_crd_0; … … 121 151 _QQ(3,3) = sig_crd_0 * sig_crd_0; 122 152 _QQ(4,4) = sig_clk_0 * sig_clk_0; 123 124 _xx.ReSize(nPar); 125 _xx = 0.0; 126 127 bncSettings settings; 128 129 _static = false; 130 if ( Qt::CheckState(settings.value("pppStatic").toInt()) == Qt::Checked) { 131 _static = true; 132 } 133 134 _usePhase = false; 135 if ( Qt::CheckState(settings.value("pppUsePhase").toInt()) == Qt::Checked) { 136 _usePhase = true; 153 if (_estTropo) { 154 _QQ(5,5) = sig_trp_0 * sig_trp_0; 137 155 } 138 156 } … … 219 237 satData->rho = (satData->xx - xRec).norm_Frobenius(); 220 238 221 double tropDelay = delay_saast(satData->eleSat); 239 double tropDelay = delay_saast(satData->eleSat) + 240 trp() / sin(satData->eleSat); 222 241 223 242 return satData->rho + clk() - satData->clk + tropDelay; … … 371 390 _QQ(4,4) = sig_clk_0 * sig_clk_0; 372 391 392 // Tropospheric Delay 393 // ------------------ 394 if (_estTropo) { 395 _params[4]->x0 += _params[4]->xx; 396 _QQ(5,5) += sig_trp_p * sig_trp_p; 397 } 398 373 399 // Ambiguities 374 400 // -----------
Note:
See TracChangeset
for help on using the changeset viewer.