Quick demo of constructing a Spectrum1D object

spectrum1d.py

import numpy as np
from astropy import units as u
from astropy.modeling import models
from specutils import Spectrum1D

mean = 5
frequencies = np.linspace(1, mean*2, 100) * u.GHz
g1 = models.Gaussian1D(amplitude=5*u.Jy, mean=mean*u.GHz, stddev=0.8*u.GHz)
spectrum = Spectrum1D(spectral_axis=frequencies, flux=g1(frequencies))