Matlab code for the Poincare reachability.
This code was used to generate the results for the 2015 AAS conference and also extended for use in the submission to Acta Astronautica.
I also forked this repo and it's now being used for the 2017 JAS submission.
The results were used in the following places
There are several large files which are not tracked in
git, instead these are backed up to Google Drive.
Make sure you do not overwrite/delete files on Google drive, but rather only pull.
Code used to regenerate the plots
L1 to Moon orbit
.matfiles are located in
poincare_intersect.mto regenerate some of the plots
Geo to L1 transfer
geo_transfer_driver.mwill generate/load the shooting trajectories from the associate mat files
Variational Integrator Comparison
The plots showing the variational integrator as compared to the ODE45 are generated by running
integrator_compare- which simulates a given trajectory for 200 non-dimensional units. It then computes the energy during this time and saves the arrays to a mat file
plot_integrator_compare('E_comparison.mat')- this function then loads the data and generates several plots
The data for the paper is stored in