1   /* Copyright 2002-2024 Thales Alenia Space
2    * Licensed to CS GROUP (CS) under one or more
3    * contributor license agreements.  See the NOTICE file distributed with
4    * this work for additional information regarding copyright ownership.
5    * CS licenses this file to You under the Apache License, Version 2.0
6    * (the "License"); you may not use this file except in compliance with
7    * the License.  You may obtain a copy of the License at
8    *
9    *   http://www.apache.org/licenses/LICENSE-2.0
10   *
11   * Unless required by applicable law or agreed to in writing, software
12   * distributed under the License is distributed on an "AS IS" BASIS,
13   * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14   * See the License for the specific language governing permissions and
15   * limitations under the License.
16   */
17  package org.orekit.estimation.measurements.modifiers;
18  
19  import org.orekit.estimation.measurements.EstimatedMeasurementBase;
20  import org.orekit.estimation.measurements.gnss.OneWayGNSSRangeRate;
21  import org.orekit.propagation.SpacecraftState;
22  import org.orekit.utils.PVCoordinates;
23  
24  /** Class modifying theoretical range-rate measurement with relativistic frequency deviation.
25   * <p>
26   * Relativistic clock correction is caused by the motion of the satellite as well as
27   * the change in the gravitational potential
28   * </p>
29   * @author Luc Maisonobe
30   * @since 12.1
31   *
32   * @see "Teunissen, Peter, and Oliver Montenbruck, eds. Springer handbook of global navigation
33   * satellite systems. Chapter 19.2. Springer, 2017."
34   */
35  public class RelativisticClockOneWayGNSSRangeRateModifier
36      extends AbstractRelativisticClockOnBoardRangeRateModifier<OneWayGNSSRangeRate> {
37  
38      /** Simple constructor.
39       * @param gm gravitational constant for main body in signal path vicinity.
40       */
41      public RelativisticClockOneWayGNSSRangeRateModifier(final double gm) {
42          super(gm);
43      }
44  
45      /** {@inheritDoc} */
46      @Override
47      public void modifyWithoutDerivatives(final EstimatedMeasurementBase<OneWayGNSSRangeRate> estimated) {
48  
49          // local satellite
50          final SpacecraftState local = estimated.getStates()[0];
51  
52          // compute semi major axis of remote (transmitter) satellite
53          final PVCoordinates remote  = estimated.getParticipants()[0];
54          final double        rRemote = remote.getPosition().getNorm();
55          final double        vRemote = remote.getVelocity().getNorm();
56          final double        aRemote = 1.0 / (2 / rRemote - vRemote * vRemote / getGm());
57  
58          modifyWithoutDerivatives(estimated,
59                                   local.getA(), local.getPosition().getNorm(),
60                                   aRemote, rRemote);
61  
62      }
63  
64  }