Cislunar Space Beginner's GuideCislunar Space Beginner's Guide
  • Satellite Simulation
Cislunar Glossary
Resources & Tools
Space News
AI Q&A
Forum
Home
Gitee
GitHub
  • 简体中文
  • English
  • Satellite Simulation
Cislunar Glossary
Resources & Tools
Space News
AI Q&A
Forum
Home
Gitee
GitHub
  • 简体中文
  • English
  • Site map

    • Home (overview)
    • What is cislunar space
    • Spacecraft trajectories
    • Directions & labs
    • Glossary · terms & definitions
    • Data & code
    • Space industry archive
  • Cislunar glossary (terms & definitions)

    • Cislunar Space Glossary
    • Fundamentals

      • Absolute Range
      • Aerodynamic Coefficient
      • Aerodynamic Moment
      • Aerospace Vehicle
      • Allan Deviation (ADEV)
      • Ballistic Coefficient
      • Bi-Elliptic Transfer
      • Body Frame
      • Celestial Coordinate System
      • Celestial Sphere
      • Characteristic Velocity
      • Coverage Angle
      • Dual One-Way Ranging (DOWR)
      • Earth Ellipsoid
      • Earth Oblateness Perturbation
      • Earth-Centered Earth-Fixed Frame (ECEF)
      • Einstein Equivalence Principle (EEP)
      • Energy Parameter
      • Earth Observation (EO)
      • Finite Thrust Maneuver
      • Free-Flight Phase
      • Free-Flight Trajectory
      • Frozen Orbit
      • Gaussian Perturbation Equations
      • Geocentric Inertial Frame
      • GPS Time
      • Gravitational Potential
      • Gravitational Redshift
      • Gravity Turn
      • Gravity vs Gravitation
      • High Altitude Airship (HAA)
      • Hit Equation
      • Hohmann Transfer
      • Inertial Navigation System
      • Instantaneous Balance Assumption
      • In-Situ Resource Utilization (ISRU)
      • Julian Date
      • Kepler's Equation
      • Korea Multi-Purpose Satellite (KOMPSAT)
      • Lagrangian Perturbation Equations
      • Launch Azimuth
      • Launch Window
      • Lift-to-Drag Ratio
      • Load Factor
      • Longitudinal and Lateral Motion
      • Lunar Lander
      • Minimum Energy Trajectory
      • Near-space
      • Newton's Iteration Method
      • Nuri (KSLV-II)
      • Nutation
      • Optimal Velocity Inclination
      • Orbit Capture
      • Orbit Insertion Conditions
      • Orbital Elements
      • Orbital Equation
      • Orbital Maneuver
      • Orbital Phase
      • Orbital Transfer Vehicle
      • Passive Hydrogen Maser (PHM)
      • Perturbation Motion
      • Phasing Orbit
      • Pitch Program Angle
      • Powered Phase
      • Precession
      • Center of Pressure
      • Range Error Coefficient
      • Reentry Corridor
      • Reentry Phase
      • Repeat Ground Track Orbit
      • Reusable Launch Vehicle
      • Synthetic Aperture Radar (SAR)
      • Satellite Ring
      • Sequential Quadratic Programming
      • Skip Reentry
      • Solar Exposure Factor
      • Specific Angular Momentum
      • Specific Impulse
      • Stagnation Heat Flux
      • Standard Atmosphere
      • Stratospheric Airship
      • Subsatellite Track
      • Sun-Synchronous Orbit
      • Thrust-to-Weight Ratio
      • Thrust
      • Total Angle of Attack
      • Trajectory Equation
      • Trajectory Optimization
      • Trim Angle of Attack
      • True Anomaly
      • Tsiolkovsky Rocket Equation
      • Powered Phase Turning Process
      • Two-Body Problem
      • Coordinated Universal Time
      • Variation of Parameters
      • Velocity Frame
      • Velocity Inclination Angle
      • Vis-Viva Equation
      • Very Low Earth Orbit (VLEO)
      • Walker Constellation
      • Zero-Angle-of-Attack Reentry
    • Dynamics & math

      • A* Search Algorithm (A* Search)
      • A2PPO (Attention-Augmented Proximal Policy Optimization)
      • Action-Angle Variables
      • Backstepping Sliding Mode Control
      • Backward Stability Set
      • Bang-bang Control (Bang-bang Control)
      • Barycentric Synodic Coordinate System
      • Batch Deployment (Batch Deployment)
      • Bicircular Four-Body Problem
      • Birkhoff-Gustavson Normal Form
      • Buoyancy-weight Imbalance
      • Capture Set
      • Central Manifold
      • Chaos Effect
      • Clohessy-Wiltshire (CW) Equation
      • Co-state Normalization (Co-state Normalization)
      • Co-state Variables
      • Coasting Arc (Coasting Arc)
      • Continuation Method (Parameter Continuation)
      • Continuation
      • Cooperative Agent (CA)
      • CR3BP with Low-Thrust (CR3BP-LT)
      • Circular Restricted Three-Body Problem (CR3BP)
      • Curriculum Learning
      • Deep Deterministic Policy Gradient (DDPG)
      • Deep Reinforcement Learning
      • Detection Graph
      • Differential Correction
      • Differential Evolution (DE) Algorithm
      • Differential Games (Differential Games)
      • Direct Collocation
      • Dynamic Programming (Dynamic Programming)
      • Dynamic Target Method
      • Ephemeris Model
      • Equinoctial Orbital Elements (Equinoctial Orbital Elements)
      • Earth Restricted Three-Body Problem (ERTBP)
      • Fuel-optimal Control
      • Fuzzy Backstepping Control
      • Generalized Advantage Estimation (GAE)
      • Gaussian Process Regression
      • Geocentric Rotating Coordinate System (GRC)
      • Hamiltonian
      • Hybrid Cluster Particle Swarm Optimization (HCPSO)
      • Heteroclinic Orbit Transfer (Heteroclinic Orbit Transfer)
      • Hill Three-Body Problem
      • Homotopy Method (Homotopy Method)
      • Improved Baseline Control-Point Method (Improved Baseline Control-Point Method)
      • Impulsive Maneuver
      • Initial Value Optimization
      • Invariant Manifold (Invariant Manifold)
      • J2000 Geocentric Equatorial Coordinate System (J2000 Geocentric Equatorial Coordinate System)
      • Jacobi Constant (Jacobi Integral)
      • K-Means Clustering (K-Means Clustering)
      • K-Medoids Clustering (K-Medoids Clustering)
      • KD-Tree (KD-Tree)
      • Libration Point (Equilibrium Point)
      • Libration Point Spacecraft Body Coordinate System (Libration Point Spacecraft Body Coordinate System)
      • Libration Point Spacecraft Orbital Coordinate System (Libration Point Spacecraft Orbital Coordinate System)
      • Lindstedt-Poincare Method (Lindstedt-Poincare Method)
      • L2-centered Rotating Coordinate System (L2-centered Rotating Coordinate System, LRC)
      • LSTM Neural Network
      • Low-Thrust Transfer MDP Formulation
      • Mass Discontinuity (Mass Discontinuity)
      • Multi-Objective Monte Carlo Tree Search (MO-MCTS)
      • Modal Analysis
      • Monodromy Matrix
      • Monte Carlo Tree Search
      • Newton-Euler Equations
      • NSGA II (Non-dominated Sorting Genetic Algorithm II)
      • Pareto Optimality
      • Particle Swarm Optimization
      • Patch Point (Splicing Point)
      • Patched Method
      • Poincaré Map
      • Poincaré Section
      • Pontryagin's Maximum Principle
      • Pseudo-Arclength Continuation
      • Spacecraft Pursuit-Evasion Game
      • Q-Law Control Law
      • Quasi-Bicircular Problem (QBCP)
      • Quasi-Bicircular Four-Body Problem
      • Reachable Set
      • Reduced-Order Dynamic Equations
      • Regional Station-keeping Control
      • Regularization
      • Reinforcement Learning Enhanced Particle Swarm Optimization (RLEPSO)
      • Saddle-Point Strategy
      • Seven-node Model
      • Shooting Method
      • Six-DOF Motion Equations
      • Sliding Mode Control
      • Solar Radiation Pressure (SRP)
      • Stability Index
      • Stability Set
      • State-Dependent Traveling Salesman Problem (SDTSP)
      • State Transition Matrix (STM)
      • Static Lift
      • Strobe Map
      • Switching Function
      • Targeting Method
      • Thermo-mechanical Coupling Model
      • Thermodynamic Model
      • Two-Point Boundary Value Problem (TPBVP)
      • Trim Condition
      • Two-Dominant Invariant Manifold Method
      • Two-Level Differential Correction Method
      • Two-node Model
      • Variational Mode Decomposition
      • Zero-Effort Miss (ZEM)
      • Zero-Velocity Surface
    • Mission orbits

      • Apolune
      • Axial Orbit
      • Ballistic Capture Orbit
      • Butterfly Orbit
      • Cycler Trajectory
      • Distant Prograde Orbit (DPO)
      • DRO Constellation
      • Distant Retrograde Orbit (DRO)
      • Earth-Moon L1/L2 Halo Orbit (EML1/EML2 Halo)
      • Free-Return Trajectory
      • Full Lunar Surface Coverage Orbit
      • Halo Orbit
      • Heteroclinic Connection
      • Horseshoe Orbit
      • Hub-and-Spoke
      • Lissajous Orbit
      • Long Period Orbit
      • Low Prograde Orbit (LoPO)
      • Low-Energy Transfer Orbit
      • Low-Thrust Transfer Orbit
      • Lyapunov Orbit
      • Multi-Revolution Halo Orbit
      • Near-Rectilinear Halo Orbit (NRHO)
      • Orbit Identification
      • Orbit Keeping (Station-Keeping)
      • Parking Orbit
      • Perilune
      • Polynomial Constraint Station-Keeping
      • Primary Impulse Orbit Transfer
      • Prograde
      • Quasi-Periodic Orbit
      • Resonance Orbit
      • Retrograde
      • Short Period Orbit
      • Transfer Orbit
      • Triangular Libration Points
      • Vertical Orbit
    • Navigation & systems

      • Altitude Regulation
      • Autonomous Navigation
      • Cislunar Spatiotemporal Reference
      • Earth-Moon Hybrid Navigation
      • Extended Kalman Filter (EKF)
      • GPS Aided GEO Augmented Navigation (GAGAN)
      • Earth GNSS Weak Signal Navigation
      • Inter-Satellite Link Navigation
      • Indian Regional Navigation Satellite System (IRNSS)
      • LEO Navigation Augmentation
      • LiAISON Navigation
      • LunaNet (Lunar Network)
      • Lunar Navigation Constellation
      • Moonlight Initiative
      • Observability
      • Positioning, Navigation, and Timing (PNT)
      • Sun-Earth-Moon Autonomous Navigation
      • Tiandu-1
      • Trajectory Planning
      • X-ray Pulsar Navigation
    • Astronomy & observation

      • Astrometry
      • Background Star Elimination
      • Cislunar Moving Objects
      • Continuous Coverage (CP)
      • Earth Albedo
      • Ephemeris Correlation
      • Hot Pixel
      • Illumination Constraint
      • Image Registration
      • Image Stacking
      • Infrared Radiation
      • Lunar Glare Zone
      • Pointing Constraint
      • Quasi-zero Wind Layer
      • Segmentation Map
      • Shift-and-Add (SAA)
      • Sidereal Tracking
      • Signal-to-Noise Ratio (SNR)
      • Solar Radiation
      • Source Extraction
      • Synthetic Tracking
      • Zonal Wind
    • Military space doctrine

      • Anti-Satellite Test (ASAT)
      • Cislunar Space Situational Awareness
      • Civil-Military Integration
      • Competitive Endurance
      • Component Field Commands
      • Commander, Space Forces (COMSPACEFOR)
      • Counterspace Operations
      • Directed Energy Weapon (DEW)
      • Distributed Architecture
      • DOTMLPF-P Framework
      • Force Design
      • Force Development
      • Force Employment
      • Force Generation
      • Golden Dome
      • Kinetic Weapon
      • Mission Command
      • Mission Delta (MD)
      • Operational Test and Training Infrastructure (OTTI)
      • Persistent Detection Corridor (PDC)
      • Resilience Map
      • Resilient/Disaggregated Architecture
      • Space Domain Awareness (SDA)
      • Space Mission Task Force (SMTF)
      • Space Superiority
      • Space Force Generation Process (SPAFORGEN)
      • System Delta (SYD)
    • Organizations

      • Anduril Industries
      • Booz Allen Hamilton
      • Danuri Lunar Orbiter
      • General Dynamics Mission Systems
      • GITAI USA
      • Indian Space Research Organisation
      • Korea Aerospace Administration
      • Lockheed Martin
      • Northrop Grumman
      • Quindar
      • Raytheon Missiles & Defense
      • Sci-Tec
      • SpaceX
      • Satish Dhawan Space Centre SHAR
      • True Anomaly
      • Turion Space

Variation of Parameters

Author: Tianjiang Shuo

Website: https://cislunarspace.cn

Definition

The variation of parameters (also known as the variation of constants) is an orbital perturbation analysis method first proposed by Euler in 1748 and later refined and published by Lagrange in 1808. This method uses orbital elements (or canonical parameters) as variables, treating the osculating orbit as continuously changing. The osculating orbit is tangent to the actual orbit at any instant with identical velocity. By numerically integrating the variation equations, the osculating orbit at any given instant can be obtained, thereby determining the actual trajectory. The fundamental equation is:

σ˙=f(σ,t,ε)\dot{\boldsymbol{\sigma}} = \mathrm{f}(\boldsymbol{\sigma}, t, \varepsilon) σ˙=f(σ,t,ε)

where σ\boldsymbol{\sigma}σ is the vector of six orbital elements or canonical parameters, and ε\varepsilonε is a small parameter of the same order of magnitude as the perturbation acceleration.

Core Elements

Comparison with Other Special Perturbation Methods

MethodVariablesCharacteristics
Cowell's methodPosition and velocityDirect integration, simplest principle, wide applicability
Encke's methodPosition and velocity deviationsUses osculating orbit as reference, deviations are small quantities
Variation of parametersOrbital elementsOsculating orbit changes continuously, larger step sizes permissible

Core Principles

  • The osculating orbit is tangent to the actual orbit at any instant, with identical velocity
  • Orbital elements or canonical parameters change much more slowly than position and velocity, allowing larger integration step sizes
  • The choice of variation parameters should avoid singularities in the equations of motion as much as possible

Typical Variation Equations

Equation TypeApplicable ScenarioCharacteristics
Gaussian Type IArbitrary perturbation forcesPerturbation force decomposed into radial, transverse, and normal components
Gaussian Type IIAtmospheric drag analysisPerturbation force decomposed into tangential, principal normal, and binormal components
Lagrangian typeConservative perturbation forcesExpressed via partial derivatives of the disturbing function

Historical Significance

The variation of parameters was the only successful method for handling orbital perturbation problems before the Cowell and Encke methods appeared. It occupies a central position in celestial mechanics perturbation theory and remains one of the most important methods for studying spacecraft orbital perturbations to this day.

Application Value

The variation of parameters is the foundational method for orbital perturbation analysis. Through this method, Gaussian-type and Lagrangian-type perturbation equations can be established to analyze the effects of various perturbation forces on orbits, including Earth's oblateness, atmospheric drag, solar radiation pressure, and three-body gravitational attraction. In general perturbation methods, both the series expansion method and the mean element method perform analytical solutions based on the variation of parameters equations.

References

  • 郑伟, 安雪滢, 周祥, 何睿智. 空天飞行力学[M]. 国防科技大学, 2026.
  • 刘林. 航天器轨道理论[M]. 国防工业出版社.
Improve this page
Last Updated: 6/5/26, 9:08 AM
Contributors: Ou Yang Jiahong
Prev
Coordinated Universal Time
Next
Velocity Frame
地月空间入门指南
Cislunar Space Beginner's GuideYour guide to cislunar space
View on GitHub

Navigate

  • Home
  • About
  • Space News
  • Glossary

Content

  • Cislunar Orbits
  • Research
  • Resources

English

  • Home
  • About
  • Space News
  • Glossary

Follow Us

© 2026 Cislunar Space Beginner's Guide  |  湘ICP备2026006405号-1
Related:智慧学习助手 UStudy航天任务工具箱 ATK
微信公众号
欢迎关注天疆说扫码关注,手机获取航天资讯