First of all python is free! Then, I think of python as a super cool programming language that is also a super calculator. It allows you to do higher level math calculations but it can also help find cheaper airplane tickets in the Internet (ask me if you’re interested). Here is a short tutorial to get you started with python after you have installed canopy.
Quick start:
- Install Enthought Canopy Express from https://www.enthought.com/downloads/.
- Open the Canopy editor and create a new file.
- Read the following mini tutorial and try some commands in the prompt.
- Type in the given equations from 1 to 6 and check the results with print statements. Don’t forget to hit run (Ctr + R or Cmd + R).
mini tutorial continued …
variables in python
you can basically use any typical variable name in python such as ‘perigee’ or ‘radius_earth’. Try to give descriptive variable names so your code is easily readable by others.
assign numerical values to variables
it’s as easy as you’d expect
perigee = 418.0
important note: in python 2 an integer number (such as ’18’) is treated differently than a floating number (such as ‘18.0’). To do math with floating numbers in python 2 you can’t have integers in the mix. So that’s why I usually initialize my integer variables with a ‘.0’ at the end so I won’t get strange results. You should try to see what I mean.
math calculations in python
To do math in python is also very simple and as expected
Examples:
Rp = perigee + radius_earth R = a*(1-e**2)/(1+e*cos(nu))
important note: in python the exponent operator is ‘**’ instead of ‘^’ as in other programming languages. Here is a short list of the operators you can use in python:
http://www.tutorialspoint.com/python/python_basic_operators.htm
print your results
The way you print in python 2 is different from python 3. This example is for python 2:
print 'Period = ',P,' min'
last but not the least: numpy
numpy is a python library that helps you do math. You can read more about it here: http://wiki.scipy.org/Tentative_NumPy_Tutorial
For now you just need to import this python library and use it blindly. On the first line of your script type:
from numpy import *
this will import all the functionality of numpy into your python script. For this tutorial it’s specially useful to be able to use the constant ‘pi’ and also for the function ‘cos’. Otherwise you would have to make your own ‘pi’ and ‘cos’.
Here is the source code to make the magic happen:
#iss orbital speed and period from numpy import * perigee = 418.0 #km apogee = 426.0 #km radius_earth = 6371.0 #km mu_earth = 3.986e5 # km^3/s^2 Rp = perigee + radius_earth Ra = apogee + radius_earth print 'Radius Perigee: ', Rp, 'km' print 'Radius Apogee : ', Ra, 'km' # 1) eccentricity e = (Ra - Rp)/(Ra + Rp) print 'eccentricity: ', e # 2) semimajor axis a = (Ra+Rp)/2.0 print 'a =', a, ' km' # 3) mechanical energy epsilon = -mu_earth/(2*a) # 4) altitude at nu=90deg nu = 90*pi/180 R = a*(1-e**2)/(1+e*cos(nu)) print 'R =', R, 'km' # 5) orbital speed V = sqrt(2*(mu_earth/R + epsilon)) print("V = {:2f} km/s".format(V)) # 6) orbital period P = 2*pi*sqrt(a**3/mu_earth) print("P = {:.2f} sec = {:.2f} min = {:.2f} hr".format(P, P/60, P/60/60))