In this exercise we use historical air temperature measurements from Potsdam, to be found at https://opendata.dwd.de/climate_environment/CDC/observations_germany/climate/monthly/kl/historical/monatswerte_KL_03987_18930101_20221231_hist.zip. The file 'product_klima_monat...' contains the monthly mean temperatures in column 'MO_TT'.
np.loadtxt( file, delimiter=';', skiprows=1, usecols=5 )
.plt.plot()
function.scipy.fftpack.fft()
.The yearly cycle is the dominant feature of monthly temperature data. In order to see other features, we can use yearly data by averaging data of each year. The resulting timeseries has 1/12th of the length of the original series.
cumsum(x)
). Assuming power-law decay of autocorrelations: what is the exponent?