.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/stochastic_processes/spectral/spectral_1d_1v.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_auto_examples_stochastic_processes_spectral_spectral_1d_1v.py>`
        to download the full example code or to run this example in your browser via Binder

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_auto_examples_stochastic_processes_spectral_spectral_1d_1v.py:


One-dimensional & one variable
=================================================================

In this example, the Spectral Representation Method is used to generate stochastic processes from a prescribed Power
Spectrum. This example illustrates how to use the SRM class for a one dimensional and one variable case and compare the
statistics of the generated stochastic processes with the expected values.

.. GENERATED FROM PYTHON SOURCE LINES 13-15

Import the necessary libraries. Here we import standard libraries such as numpy and matplotlib, but also need to
import the SRM class from the StochasticProcesses module of UQpy.

.. GENERATED FROM PYTHON SOURCE LINES 18-24

.. code-block:: default


    from UQpy.stochastic_process import SpectralRepresentation
    import numpy as np
    import matplotlib.pyplot as plt
    plt.style.use('seaborn')


.. GENERATED FROM PYTHON SOURCE LINES 25-26

The input parameters necessary for the generation of the stochastic processes are given below:

.. GENERATED FROM PYTHON SOURCE LINES 29-47

.. code-block:: default


    n_sim = 10000  # Num of samples

    n = 1 # Num of dimensions

    m = 1 # Num of variables

    T = 100  # Time(1 / T = dw)
    nt = 256  # Num.of Discretized Time
    F = 1 / T * nt / 2  # Frequency.(Hz)
    nw = 128  # Num of Discretized Freq.

    # # Generation of Input Data(Stationary)
    dt = T / nt
    t = np.linspace(0, T - dt, nt)
    dw = F / nw
    w = np.linspace(0, F - dw, nw)


.. GENERATED FROM PYTHON SOURCE LINES 48-49

Make sure that the input parameters are in order to prevent aliasing.

.. GENERATED FROM PYTHON SOURCE LINES 52-58

.. code-block:: default


    t_u = 2*np.pi/2/F

    if dt>t_u:
        print('Error')


.. GENERATED FROM PYTHON SOURCE LINES 59-60

Defining the Power Spectral Density Function.

.. GENERATED FROM PYTHON SOURCE LINES 63-77

.. code-block:: default


    S = 125 / 4 * w ** 2 * np.exp(-5 * w)

    SRM_object = SpectralRepresentation(n_sim, S, dt, dw, nt, nw)
    samples = SRM_object.samples

    fig, ax = plt.subplots()
    plt.title('Realisation of the Spectral Representation Method')
    plt.plot(t, samples[0, 0])
    ax.yaxis.grid(True)
    ax.xaxis.grid(True)
    plt.show()

    print('The mean of the samples is ', np.mean(samples), 'whereas the expected mean is 0.000')
    print('The variance of the samples is ', np.var(samples), 'whereas the expected variance is ', np.sum(S)*dw*2)

.. rst-class:: sphx-glr-timing

   **Total running time of the script:** ( 0 minutes  0.000 seconds)


.. _sphx_glr_download_auto_examples_stochastic_processes_spectral_spectral_1d_1v.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example


    .. container:: binder-badge

      .. image:: images/binder_badge_logo.svg
        :target: https://mybinder.org/v2/gh/SURGroup/UQpy/master?urlpath=lab/tree/notebooks/auto_examples/stochastic_processes/spectral/spectral_1d_1v.ipynb
        :alt: Launch binder
        :width: 150 px



    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: spectral_1d_1v.py <spectral_1d_1v.py>`

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: spectral_1d_1v.ipynb <spectral_1d_1v.ipynb>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_