.. note:: :class: sphx-glr-download-link-note Click :ref:`here ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_intro_summary-exercises_auto_examples_plot_cumulative_wind_speed_prediction.py: Cumulative wind speed prediction ================================ Generate the image cumulative-wind-speed-prediction.png for the interpolate section of scipy.rst. .. image:: /intro/summary-exercises/auto_examples/images/sphx_glr_plot_cumulative_wind_speed_prediction_001.png :class: sphx-glr-single-img .. code-block:: python import numpy as np from scipy.interpolate import UnivariateSpline import matplotlib.pyplot as plt max_speeds = np.load('max-speeds.npy') years_nb = max_speeds.shape[0] cprob = (np.arange(years_nb, dtype=np.float32) + 1)/(years_nb + 1) sorted_max_speeds = np.sort(max_speeds) speed_spline = UnivariateSpline(cprob, sorted_max_speeds) nprob = np.linspace(0, 1, 1e2) fitted_max_speeds = speed_spline(nprob) fifty_prob = 1. - 0.02 fifty_wind = speed_spline(fifty_prob) plt.figure() plt.plot(sorted_max_speeds, cprob, 'o') plt.plot(fitted_max_speeds, nprob, 'g--') plt.plot([fifty_wind], [fifty_prob], 'o', ms=8., mfc='y', mec='y') plt.text(30, 0.05, '$V_{50} = %.2f \, m/s$' % fifty_wind) plt.plot([fifty_wind, fifty_wind], [plt.axis()[2], fifty_prob], 'k--') plt.xlabel('Annual wind speed maxima [$m/s$]') plt.ylabel('Cumulative probability') **Total running time of the script:** ( 0 minutes 0.076 seconds) .. _sphx_glr_download_intro_summary-exercises_auto_examples_plot_cumulative_wind_speed_prediction.py: .. only :: html .. container:: sphx-glr-footer :class: sphx-glr-footer-example .. container:: sphx-glr-download :download:`Download Python source code: plot_cumulative_wind_speed_prediction.py ` .. container:: sphx-glr-download :download:`Download Jupyter notebook: plot_cumulative_wind_speed_prediction.ipynb ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_