AdditionalEquations
public class PartialDerivativesEquations extends Object implements AdditionalEquations
additional equations
computing the partial derivatives
of the state (orbit) with respect to initial state and force models parameters.
This set of equations are automatically added to a numerical propagator
in order to compute partial derivatives of the orbit along with the orbit itself. This is
useful for example in orbit determination applications.
The partial derivatives with respect to initial state can be either dimension 6 (orbit only) or 7 (orbit and mass).
The partial derivatives with respect to force models parameters has a dimension
equal to the number of selected parameters. Parameters selection is implemented at
force models
level. Users must retrieve a parameter driver
using ForceModel.getParameterDriver(String)
and then
select it by calling setSelected(true)
.
If several force models provide different drivers
for the
same parameter name, selecting any of these drivers has the side effect of
selecting all the drivers for this shared parameter. In this case, the partial
derivatives will be the sum of the partial derivatives contributed by the
corresponding force models. This case typically arises for central attraction
coefficient, which has an influence on Newtonian attraction
, gravity field
, and relativity
.
Constructor | Description |
---|---|
PartialDerivativesEquations(String name,
NumericalPropagator propagator) |
Simple constructor.
|
Modifier and Type | Method | Description |
---|---|---|
double[] |
computeDerivatives(SpacecraftState s,
double[] pDot) |
Compute the derivatives related to the additional state parameters.
|
JacobiansMapper |
getMapper() |
Get a mapper between two-dimensional Jacobians and one-dimensional additional state.
|
String |
getName() |
Get the name of the additional state.
|
ParameterDriversList |
getSelectedParameters() |
Get the selected parameters, in Jacobian matrix column order.
|
SpacecraftState |
setInitialJacobians(SpacecraftState s0) |
Set the initial value of the Jacobian with respect to state and parameter.
|
SpacecraftState |
setInitialJacobians(SpacecraftState s1,
double[][] dY1dY0,
double[][] dY1dP) |
Set the initial value of the Jacobian with respect to state and parameter.
|
SpacecraftState |
setInitialJacobians(SpacecraftState s0,
int stateDimension) |
Deprecated.
as of 9.0, replaced by
setInitialJacobians(SpacecraftState) |
init
public PartialDerivativesEquations(String name, NumericalPropagator propagator) throws OrekitException
Upon construction, this set of equations is automatically added to
the propagator by calling its AbstractIntegratedPropagator.addAdditionalEquations(AdditionalEquations)
method. So
there is no need to call this method explicitly for these equations.
name
- name of the partial derivatives equationspropagator
- the propagator that will handle the orbit propagationOrekitException
- if a set of equations with the same name is already presentpublic String getName()
getName
in interface AdditionalEquations
public ParameterDriversList getSelectedParameters() throws OrekitException
The force models parameters for which partial derivatives are desired,
must have been selected
before this method is called, so the proper list is returned.
OrekitException
- if an existing driver for a
parameter throws one when its value is reset using the value
from another driver managing the same parameterpublic SpacecraftState setInitialJacobians(SpacecraftState s0) throws OrekitException
This method is equivalent to call setInitialJacobians(SpacecraftState,
double[][], double[][])
with dYdY0 set to the identity matrix and dYdP set
to a zero matrix.
The force models parameters for which partial derivatives are desired,
must have been selected
before this method is called, so proper matrices dimensions are used.
s0
- initial stateOrekitException
- if the partial equation has not been registered in
the propagator or if matrices dimensions are incorrectgetSelectedParameters()
@Deprecated public SpacecraftState setInitialJacobians(SpacecraftState s0, int stateDimension) throws OrekitException
setInitialJacobians(SpacecraftState)
This method is equivalent to call setInitialJacobians(SpacecraftState,
double[][], double[][])
with dYdY0 set to the identity matrix and dYdP set
to a zero matrix.
The force models parameters for which partial derivatives are desired,
must have been selected
before this method is called, so proper matrices dimensions are used.
s0
- initial statestateDimension
- state dimension, must be either 6 for orbit only or 7 for orbit and massOrekitException
- if the partial equation has not been registered in
the propagator or if matrices dimensions are incorrectgetSelectedParameters()
public SpacecraftState setInitialJacobians(SpacecraftState s1, double[][] dY1dY0, double[][] dY1dP) throws OrekitException
The returned state must be added to the propagator (it is not done automatically, as the user may need to add more states to it).
The force models parameters for which partial derivatives are desired,
must have been selected
before this method is called, and the dY1dP
matrix dimension must
be consistent with the selection.
s1
- current statedY1dY0
- Jacobian of current state at time t₁ with respect
to state at some previous time t₀ (must be 6x6)dY1dP
- Jacobian of current state at time t₁ with respect
to parameters (may be null if no parameters are selected)OrekitException
- if the partial equation has not been registered in
the propagator or if matrices dimensions are incorrectgetSelectedParameters()
public JacobiansMapper getMapper() throws OrekitException
OrekitException
- if the initial Jacobians have not been initialized yetsetInitialJacobians(SpacecraftState, int)
,
setInitialJacobians(SpacecraftState, double[][], double[][])
public double[] computeDerivatives(SpacecraftState s, double[] pDot) throws OrekitException
When this method is called, the spacecraft state contains the main
state (orbit, attitude and mass), all the states provided through
the additional
state providers
registered to the propagator, and the additional state
integrated using this equation. It does not contains any other
states to be integrated alongside during the same propagation.
computeDerivatives
in interface AdditionalEquations
s
- current state information: date, kinematics, attitude, and
additional statepDot
- placeholder where the derivatives of the additional parameters
should be putOrekitException
- if some specific error occursCopyright © 2002-2018 CS Systèmes d'information. All rights reserved.