.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/plot_manova_ensemble.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_plot_manova_ensemble.py>`
        to download the full example code.

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

.. _sphx_glr_auto_examples_plot_manova_ensemble.py:


Manova Ensemble
=================

Defines Manova Ensemble random matrices

.. GENERATED FROM PYTHON SOURCE LINES 7-13

.. code-block:: Python


    # Author: Alejandro Santorum Varela
    # License: BSD 3-Clause

    from skrmt.ensemble import ManovaEnsemble








.. GENERATED FROM PYTHON SOURCE LINES 14-16

Manova Ensemble random matrices are considered to be double Wishart random
matrices.

.. GENERATED FROM PYTHON SOURCE LINES 18-28

Manova Real Ensemble (MRE)
---------------------------

Random matrices of MRE are formed by sampling two random real standard
guassian matrices (:math:`\mathbf{X}` and :math:`\mathbf{Y}`) of size
:math:`m \times n_1` and :math:`m \times n_2` respectively. Then, matrix
:math:`\mathbf{A} = \dfrac{\mathbf{X}\mathbf{X}^T}{\mathbf{X}\mathbf{X}^T  + \mathbf{Y}\mathbf{Y}^T}`
generates a matrix of the MRE.

They are also known as 1-Jacobi random matrices (beta = 1).

.. GENERATED FROM PYTHON SOURCE LINES 30-32

A random matrix of Manova Real Ensemble can be sampled using
scikit-rmt with the following code.

.. GENERATED FROM PYTHON SOURCE LINES 32-36

.. code-block:: Python


    mre = ManovaEnsemble(beta=1, m=3, n1=5, n2=5)
    print(mre.matrix)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    [[ 0.72525922 -0.11275618 -0.16722456]
     [-0.2002338   0.22270727  0.05215273]
     [-0.07550653  0.00267971  0.82476058]]




.. GENERATED FROM PYTHON SOURCE LINES 37-47

Manova Complex Ensemble (MCE)
------------------------------

Random matrices of MCE are formed by sampling two random complex standard
guassian matrices (:math:`\mathbf{X}` and :math:`\mathbf{Y}`) of size
:math:`m \times n_1` and :math:`m \times n_2` respectively. Then, matrix
:math:`\mathbf{A} = \dfrac{\mathbf{X}\mathbf{X}^T}{\mathbf{X}\mathbf{X}^T  + \mathbf{Y}\mathbf{Y}^T}`
generates a matrix of the MCE.

They are also known as 2-Jacobi random matrices (beta = 2).

.. GENERATED FROM PYTHON SOURCE LINES 49-51

A random matrix of Manova Complex Ensemble can be sampled using
scikit-rmt with the following code.

.. GENERATED FROM PYTHON SOURCE LINES 51-55

.. code-block:: Python


    mce = ManovaEnsemble(beta=2, m=3, n1=5, n2=5)
    print(mce.matrix)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    [[ 0.39670461+0.10974469j -0.18481543+0.08933746j -0.11432087-0.01577289j]
     [-0.04875815+0.07830072j  0.40153892+0.0692018j  -0.14120906-0.14290005j]
     [-0.24747441+0.19002681j -0.47625583+0.14751422j  0.45769172-0.17894649j]]




.. GENERATED FROM PYTHON SOURCE LINES 56-74

Manova Quaternion Ensemble (MQE)
---------------------------------

Random matrices of MQE are formed by sampling two random complex standard 
guassian matrices (:math:`\mathbf{X_1}` and :math:`\mathbf{X_2}`), both of
size :math:`m \times n_1`. Another two random complex standard guassian matrices
(:math:`\mathbf{Y_1}` and :math:`\mathbf{Y_2}`), both of size :math:`m \times n_2`,
are sampled. They are stacked forming matrices :math:`\mathbf{X}` and
:math:`\mathbf{Y}`:

:math:`\mathbf{X} = (\mathbf{X_1}\ \mathbf{X_2}; -\mathbf{X_2}^*\ \mathbf{X_1}^*)`

:math:`\mathbf{Y} = (\mathbf{Y_1}\ \mathbf{Y_2}; -\mathbf{Y_2}^*\ \mathbf{Y_1}^*)`

Finally, matrix :math:`\mathbf{A} = \dfrac{\mathbf{X}\mathbf{X}^T}{\mathbf{X}\mathbf{X}^T  + \mathbf{Y}\mathbf{Y}^T}`
generates a matrix of the MQE.

They are also known as 4-Jacobi random matrices (beta = 4).

.. GENERATED FROM PYTHON SOURCE LINES 76-78

A random matrix of Manova Quaternion Ensemble can be sampled using
scikit-rmt with the following code.

.. GENERATED FROM PYTHON SOURCE LINES 78-80

.. code-block:: Python


    mqe = ManovaEnsemble(beta=4, m=2, n1=5, n2=5)
    print(mqe.matrix)



.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    [[ 0.44500187-0.0030627j   0.07797313-0.00523991j  0.00550757-0.00992674j
       0.01551386-0.00407187j]
     [ 0.04305868+0.00015381j  0.59433566-0.00020158j  0.00309942-0.01652418j
      -0.00667558+0.00967728j]
     [-0.00550757-0.00992674j -0.01551386-0.00407187j  0.44500187+0.0030627j
       0.07797313+0.00523991j]
     [-0.00309942-0.01652418j  0.00667558+0.00967728j  0.04305868-0.00015381j
       0.59433566+0.00020158j]]





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

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


.. _sphx_glr_download_auto_examples_plot_manova_ensemble.py:

.. only:: html

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

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

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

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

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

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: plot_manova_ensemble.zip <plot_manova_ensemble.zip>`


.. only:: html

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

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