My Kerbal RPC scripts

analysis.rmd 2.1KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657
  1. ---
  2. title: "Flight analysis"
  3. ---
  4. This file creates graphs of recorded flight parameters
  5. ```{r opts, include=FALSE}
  6. library(knitr)
  7. opts_chunk$set(tidy=TRUE)
  8. ```
  9. ```{r init}
  10. library(ggplot2)
  11. library(scales)
  12. flight <- read.csv("flight.csv", header=T)
  13. flight$mission_time <- as.POSIXct(flight$mission_time, origin=Sys.Date(), tz="GMT")
  14. ```
  15. ```{r graph}
  16. baseg <- ggplot(flight, aes(x=mission_time)) + xlab("Time (H:M:S)") + theme_classic() + scale_x_datetime(date_labels="%H:%M:%S") + scale_y_continuous(labels=comma)
  17. alt <- baseg + geom_line(aes(y=current_altitude)) + ylab("Altitude (m)") + ggtitle("Altitude")
  18. vspd <- baseg + geom_line(aes(y=v_speed)) + ylab("Vertical speed (m/s)") + ggtitle("Vertical speed")
  19. hspd <- baseg + geom_line(aes(y=h_speed)) + ylab("Horizontal speed (m/s)") + ggtitle("Horizontal speed")
  20. mass <- baseg + geom_line(aes(y=mass)) + ylab("Vessel mass (kg)") + ggtitle("Mass")
  21. ec <- baseg + geom_line(aes(y=electric_charge)) + ylab("Electric charge") + ggtitle("Electric charge")
  22. lf <- baseg + geom_line(aes(y=liquid_fuel)) + ylab("Liquid fuel (litres)") + ggtitle("Liquid fuel")
  23. ox <- baseg + geom_line(aes(y=oxidizer)) + ylab("Oxidizer (litres)") + ggtitle("Oxidizer")
  24. at <- baseg + geom_line(aes(y=available_thrust)) + ylab("Available thrust (Newtons)") + ggtitle("Available thrust")
  25. thrust <- baseg + geom_line(aes(y=current_thrust)) + ylab("Current thrust (Newtons)") + ggtitle("Current thrust")
  26. g_force <- baseg + geom_line(aes(y=g_force)) + ylab("G-force (m/s²)") + ggtitle("G-force")
  27. ap <- baseg + geom_line(aes(y=apoapsis)) + ylab("Apoapsis height (m)") + ggtitle("Apoapsis")
  28. pa <- baseg + geom_line(aes(y=periapsis)) + ylab("Periapsis height (m)") + ggtitle("Periapsis")
  29. orad <- baseg + geom_line(aes(y=orbital_radius)) + ylab("Orbital radius (m)") + ggtitle("Orbital radius")
  30. os <- baseg + geom_line(aes(y=orbit_speed)) + ylab("Orbital speed (m/s)") + ggtitle("Orbital speed")
  31. op <- baseg + geom_line(aes(y=orbit_period)) + ylab("Orbital period (s)") + ggtitle("Orbital period")
  32. ```
  33. ```{r drawgraphs, echo=FALSE}
  34. alt
  35. vspd
  36. hspd
  37. mass
  38. ec
  39. lf
  40. ox
  41. at
  42. thrust
  43. g_force
  44. ap
  45. pa
  46. orad
  47. os
  48. op
  49. ```