- Timestamp:
- Oct 7, 2015, 3:56:55 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/ephemeris.cpp
r7481 r7487 1530 1530 double zz = yp*sini; 1531 1531 1532 Matrix R 1= BNC_PPP::t_astro::rotX(-5.0 / 180.0 * M_PI);1533 Matrix R 2= BNC_PPP::t_astro::rotZ(ll);1532 Matrix RX = BNC_PPP::t_astro::rotX(-5.0 / 180.0 * M_PI); 1533 Matrix RZ = BNC_PPP::t_astro::rotZ(ll); 1534 1534 1535 1535 ColumnVector X1(3); X1 << xx << yy << zz; 1536 ColumnVector X2 = R 2*R1*X1;1536 ColumnVector X2 = RZ*RX*X1; 1537 1537 1538 1538 xc[0] = X2(1); … … 1553 1553 1554 1554 ColumnVector V(3); 1555 1556 1555 V[0] = cosom *vx - cosi*sinom *vy // dX / dr 1557 1556 - xp*sinom*dotom - yp*cosi*cosom*dotom // dX / dOMEGA … … 1564 1563 V[2] = sini *vy + yp*cosi *doti; 1565 1564 1566 Matrix R 3(3,3);1565 Matrix RdotZ(3,3); 1567 1566 double C = cos(ll); 1568 1567 double S = sin(ll); … … 1571 1570 UU[1][0] = -C; UU[1][1] = -S; UU[1][2] = 0.0; 1572 1571 UU[2][0] = 0.0; UU[2][1] = 0.0; UU[2][2] = 0.0; 1573 R 3= omegaBDS * UU;1572 RdotZ = omegaBDS * UU; 1574 1573 1575 1574 ColumnVector VV(3); 1576 1577 VV = R2*R1*V + R3*R1*X1; 1575 VV = RZ*RX*V + RdotZ*RX*X1; 1578 1576 1579 1577 vv[0] = VV(1);
Note:
See TracChangeset
for help on using the changeset viewer.