SGP4 Orbit Propagator with GUI
- 1 file
- 50 downloads
Main function: OrbProp.m
NOTE: 1. Output generated by this program has not been validated, other than ensuring close agreement with AGI STK results. Prior to MAR 2018, expired Earth Orbit Parameters (EOP) were used for xp, yp, LOD, etc. This incurred slight variation from STK results (~0-30s offset). The program now downloads EOP data automatically from the internet, or allows selection of an EOP file the user may have.
2. Mapping toolbox and external program needed to execute distance, elevation, reckon, and azimuth functions; and plot LLA ground tracks.
This script runs the SGP4 propagator using the manual method, afspc, with
WGS72 ellipsoid. SGP4 output in ManualOut.out is transformed from TEME
to ECEF. The transformation teme2ecef is executed for each line of output in ManualOutput.out using time parameters for that iteration.
Input: TLE, UTC start, UTC stop, observer location, sector parameters. TLE input is a single two line file.
Output:
1. Access_Summary.txt: Summary of visibility times from observer's location, with az/el/slantrange
2. Access_Transit.txt: Used for plotting ground tracks for passes in Access_Summary.txt. Only visibility periods are included. Format: line number, lat, lon, alt (m geodetic), el(deg), az(degT), timevector
3. LLA_Transit.txt: Used for plotting ENTIRE ground track for period specified in the GUI. Format: line number, lat, lon, alt (m geodetic), timevector
4. PolarPlot.png: Az-El polar plot of passes in Access_Summary.txt
5. KML_Transit_Altitude.kml: Plot orbit at altitude for entire period in Google Earth
6. KML_Transit_GroundTrack.kml: Plot orbit ground track for entire period in Google Earth
--External program using mapping toolbox required for plotting ground tracks.
--A sector may be specified for the observer that constrains visibility windows.
--Geodetic horizon at observer location used for visibility when no sector specified.
Will appreciate any feedback/rating if difficulty is encountered or works as needed.
Cite As
Damon DeLuca (2020). SGP4 Orbit Propagator with GUI (https://www.mathworks.com/matlabcentral/fileexchange/55179-sgp4-orbit-propagator-with-gui), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (10)
Updates
| 1.2.0.0 | Updated program to incorporate Earth Orientation Parameters. Results now agree with AGI STK. |
|
| 1.1.0.0 | Added access period labels to polar plot. Changed from .piz file to .zip. |
|
| 1.0.0.0 | None |
|
| 1.0.0.0 | 1. Updated to include KML output files.
|
|
| 1.0.0.0 | None |
|
| 1.0.0.0 | Included additional required Matlab functions not previously included.
|
|
| 1.0.0.0 | Mapping toolbox needed. Will include those functions next update. |
|
| 1.0.0.0 | None |
|
| 1.0.0.0 | None |
|
| 1.0.0.0 | None |
|
| 1.0.0.0 | None |
|
| 1.0.0.0 | None |
|
| 1.0.0.0 | None |
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
Tags
Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
bavya lakshmi (view profile)
Michael Ernst (view profile)
Not working with R2018a:
Index in position 2 exceeds array bounds (must not exceed 1).
Error in OrbProp>run_Callback (line 657)
polar(accesstransit(:,6).*3.14159/180,90-accesstransit(:,5),'r');
Error in gui_mainfcn (line 95)
feval(varargin{:});
Error in OrbProp (line 83)
gui_mainfcn(gui_State, varargin{:});
Error in
matlab.graphics.internal.figfile.FigFile/read>@(hObject,eventdata)OrbProp('run_Callback',hObject,eventdata,guidata(hObject))
Error while evaluating UIControl Callback.
Zohair Abu Shaban (view profile)
Is there a way to obtain the velocity or range rate at a given orbit point too?
Ahsan Waqas (view profile)
Damon DeLuca (view profile)
#Bing:
KML output for access transit ground tracks is included. Thanks for the comment!
Blewis (view profile)
Bing (view profile)
Could you also please add in the result a kml file for all the access transit ground tracks?
Stephan Rediske (view profile)
Sean Kim (view profile)
christie harper (view profile)
need to state that mapping toolbox is required