?About the author: A Matlab simulation developer who loves scientific research. He cultivates his mind and improves his technology simultaneously. For cooperation on MATLAB projects, please send a private message.
Personal homepage: Matlab Research Studio
Personal credo: Investigate things to gain knowledge.
For more complete Matlab code and simulation customization content, click
Intelligent optimization algorithm Neural network prediction Radar communication Wireless sensor Power system
Signal processing Image processing Path planning Cellular automaton Drone
Content introduction
Abstract: Satellite perturbation motion simulation is an important research direction in the aerospace field. This paper presents an algorithmic procedure based on the Runge-Kutta sum force model for simulating satellite motion in a perturbed environment. Through this algorithm, the orbit and attitude of the satellite can be predicted more accurately and provide an important reference for satellite control and navigation.
Introduction: With the continuous development of aerospace technology, satellites are increasingly used in earth orbit. However, satellites are affected by various disturbances in orbit, such as gravity, atmospheric drag, solar radiation pressure, etc. These disturbances can cause changes in the satellite’s orbit and attitude. Therefore, accurately simulating the perturbed motion of a satellite is crucial for satellite control and navigation.
Runge-Kutta algorithm: The Runge-Kutta algorithm is a commonly used numerical integration method for solving ordinary differential equations. This algorithm discretizes time-continuous problems into time-discrete problems to approximate the real solution. In simulating the perturbation motion of a satellite, the Runge-Kutta algorithm can be used to solve for the position and velocity of the satellite.
Force model: Satellites are affected by various forces in orbit, such as gravity, atmospheric drag, solar radiation pressure, etc. In order to accurately simulate the perturbation motion of the satellite, a corresponding force model needs to be established. The force model can calculate the magnitude and direction of various forces based on the satellite’s motion state and environmental parameters.
Algorithm steps: The satellite perturbation motion algorithm based on the Runge-Kutta sum force model can be divided into the following steps:
-
Initialization: Set the initial position, speed and attitude of the satellite, as well as other related parameters.
-
Calculate force: Calculate the magnitude and direction of various forces based on the satellite’s motion status and environmental parameters. These forces include gravity, atmospheric drag, solar radiation pressure, etc.
-
Update status: Use the Runge-Kutta algorithm to calculate the position and velocity for the next time step based on the current position, velocity, and force.
-
Update attitude: Calculate the attitude of the next time step based on the motion state of the satellite and the effect of force.
-
Repeat steps 3 and 4 until the simulation ends.
Result analysis: Through the satellite perturbation motion algorithm based on the Runge-Kutta sum force model, the motion of the satellite in the perturbation environment can be simulated more accurately. Through this algorithm, changes in the orbit and attitude of the satellite can be predicted, and it provides an important reference for satellite control and navigation. However, this algorithm also has some limitations, such as higher accuracy requirements for the force model and higher computational complexity.
Conclusion: Satellite perturbation motion simulation is an important research direction in the aerospace field. This paper presents an algorithmic procedure based on the Runge-Kutta sum force model for simulating satellite motion in a perturbed environment. Through this algorithm, the orbit and attitude of the satellite can be predicted more accurately and provide an important reference for satellite control and navigation. However, this algorithm also has some limitations that require further research and improvement.
Part of the code
%------------------------------------------------- ----------------------------</code><code>%</code><code>% Chebyshev approximation of 3-dimensional vectors</code><code>%</code><code>% Inputs:</code><code>% N Number of coefficients</code><code>% Ta Begin interval</code><code>% Tb End interval</code><code>% Cx Coefficients of Chebyshev polyomial (x-coordinate)</code><code>% Cy Coefficients of Chebyshev polyomial (y-coordinate)</code><code>% Cz Coefficients of Chebyshev polyomial (z-coordinate)</code><code>%</code><code>% Last modified: 2018/01/27 Meysam Mahooti</code><code>% </code><code>%- -------------------------------------------------- -----------------------</code><code>function ChebApp = Cheb3D(t, N, Ta, Tb, Cx, Cy, Cz)</code><code>?</code><code>% Check validity</code><code>if ( (t<Ta) || (Tb<t) )</code><code> error(' ERROR: Time out of range in Cheb3D::Value\\ ');</code><code>end</code><code>?</code><code>% Clenshaw algorithm</code><code>tau = (2*t-Ta-Tb)/(Tb-Ta); </code><code>?</code><code>f1 = zeros(1,3);</code><code> f2 = zeros(1,3);</code><code>?</code><code>for i=N:-1:2</code><code> old_f1 = f1;</code><code> f1 = 2*tau*f1-f2 + [Cx(i),Cy(i),Cz(i)];</code><code> f2 = old_f1;</code><code>end</code><code>?</code><code>ChebApp = tau*f1-f2 + [Cx(1),Cy(1),Cz(1)];</code><code>?</code><code>?
% ---------------------------------------------- --------------------------------</code><code>%</code><code>% function finddays </code><code>%</code><code>% This function finds the fractional days through a year given the year,</code><code>% month, day, hour, minute and second.</code> ><code>%</code><code>?</code><code>% revisions</code><code>% -</code><code>%</code><code>% inputs description range / units</code><code>% year - year 1900 .. 2100</code><code>% mon - month 1 .. 12</code><code>% day - day 1 .. 28,29, 30,31</code><code>% hr - hour 0 .. 23</code><code>% min - minute 0 .. 59</code><code>% sec - second 0.0 .. 59.999</code> code><code>%</code><code>% outputs:</code><code>% days - day of year plus fraction of a</code><code>% day days</code><code> %</code><code>% locals:</code><code>% lmonth - length of months of year</code><code>% i - index</code><code>%</code><code>% coupling :</code><code>% none.</code><code>%</code><code>% references :</code><code>% vallado 2007, 207, ex 3-12 </code><code>%</code><code>% [days] = finddays (year,month,day,hr,min,sec);</code><code>% ------- -------------------------------------------------- --------------------</code><code>function [days] = finddays (year,month,day,hr,min,sec)</code> <code>?</code><code>for i= 1:12</code><code> lmonth(i) = 31;</code><code> if i == 2</code><code> lmonth(i)= 28;</code><code> end</code><code> if i == 4 | i == 6 | i == 9 | i == 11</code><code> lmonth (i)= 30;</code><code> end</code><code>end</code><code>?</code><code>if (rem(year,4) == 0)</code><code> lmonth(2)= 29;</code><code> if (rem(year,100) == 0) & amp; (rem(year,400) ~= 0)</code> <code> lmonth(2)= 28;</code><code> end</code><code>end</code><code>?</code><code>i = 1;</code><code>days= 0.0;</code><code>while (i < month) & amp; ( i < 12 )</code><code> days= days + lmonth(i);</code><code> i= i + 1;</code><code>end</code><code>?</code><code>days= days + day + hr/24.0 + min/1440.0 + sec/86400.0;</code><code>?</code><code>?
Operation results
References
Montenbruck O., Gill E.; Satellite Orbits: Models, Methods and Applications; Springer Verlag, Heidelberg; Corrected 3rd Printing (2005).
Montenbruck O., Pfleger T.; Astronomy on the Personal Computer; Springer Verlag, Heidelberg; 4th edition (2000).
Seeber G.; Satellite Geodesy; Walter de Gruyter, Berlin, New York; 2nd completely revised and extended edition (2003).
Vallado D. A; Fundamentals of Astrodynamics and Applications; McGraw-Hill, New York; 4th edition (2013).
Damn it. 2000. Department of Defense World Geodetic System 1984. NIMA-TR 8350.2, Edition 3, Amendment 1. Washington, DC: Headquarters, National Imagery and Mapping Agency. National Imagery and Mapping Agency.