| notebook.community

notebook.community

First Last - Homework 6



In [ ]:

    
%matplotlib inline
import matplotlib.pyplot as plt

import numpy as np
from astropy.table import QTable

The file `SDSS_MainBelt.csv` contains data on objects in the asteroid belt collected by the Sloan Digital Sky Survey.

The columns are:

Name - Object name
a - semi-major axis
cA - SDSS a$^*$ color
cI - SDSS i color (near infrared)
cZ - SDSS z color (infrared)
X - X position [au] - Sun-centered coords (10/29/2014)
Y - Y position [au] - Sun-centered coords (10/29/2014)
Z - Z position [au] - Sun-centered coords (10/29/2014)

Read in the file `SDSS_MainBelt.csv`



In [ ]:

The Kirkwood Gaps - Mean Motion Resonances

Orbital period [years]: P = a ** (3/2)
Orbital Period of Jupiter: Pjup = 11.8618 years
Make a histogram of: (Pjup / P)
Bins = 150
Invert the x-axis
Draw vertical lines at (Pjup / P) = 1,2,3,4
Label the axis
Make the plot look nice



In [ ]:

The Position of the Asteroids

Make a plot of X vs Y for the asteroids
Make sure that 1 unit in X = 1 unit in Y
Make the axes go through (0,0)
Make the plot look nice



In [ ]:

The Color of the Asteroids

Make three (3) plots in one row
In each panel, put a plot of all of the asteroid positions as viewed from above. Use alpha = 0.20
In the first panel, overplot the positions of C-Type asteroids using a different color and symbol
In the second panel, overplot the positions of S-Type asteroids using a different color and symbol
In the third panel, overplot the positions of V-Type asteroids using a different color and symbol
The asteroid types can be determined from their colors (see image below)
Make the plot look nice



In [ ]:

Fitting data to a function

The datafile 'FitMe.csv' contains two columns of data (x,y)
Plot the data x,y
Fit the function below to the data
Find the values of a,b,c that best fit the data
Plot the fit with the data
Make the plots look nice

$$\Large f(x) = a \sin(bx) \ e^{-cx} $$



In [ ]:

Asteroid classes - SDSS Colors

Due Tue Nov 08 - 5pm

Make sure to change the filename to your name!
Make sure to change the Title to your name!
File -> Download as -> HTML (.html)
upload your .html and .ipynb file to the class Canvas page



In [ ]: