Anonymous
×
Create a new article
Write your page title here:
We currently have 106 articles on MOR Wiki. Type your article name above or click on one of the titles below and start writing!



Earth Atmosphere: Difference between revisions

remove preliminary and stub warnings
Mlinaric (talk | contribs)
Add categories, add Python code, edit math
Line 1: Line 1:
[[Category:benchmark]]
[[Category:benchmark]]
[[Category:SLICOT]]
[[Category:SLICOT]]
[[Category:linear]]
[[Category:time invariant]]
[[Category:first differential order]]
[[Category:SISO]]
[[Category:SISO]]
[[Category:Dense]]
[[Category:Dense]]
Line 8: Line 11:
This benchmark models the track of an [[wikipedia:Storm_track|atmospheric storm track]].
This benchmark models the track of an [[wikipedia:Storm_track|atmospheric storm track]].
More details can be found in <ref name="farrell95"/> and <ref name="chahlaoui02"/>, <ref name="chahlaoui05"/>.
More details can be found in <ref name="farrell95"/> and <ref name="chahlaoui02"/>, <ref name="chahlaoui05"/>.


==Origin==
==Origin==


This benchmark is part of the '''SLICOT Benchmark Examples for Model Reduction'''<ref name="chahlaoui05"/>.
This benchmark is part of the '''SLICOT Benchmark Examples for Model Reduction'''<ref name="chahlaoui05"/>.


==Data==
==Data==
Line 19: Line 20:
The system matrices <math>A</math>, <math>B</math>, <math>C</math> are available from the [http://slicot.org/20-site/126-benchmark-examples-for-model-reduction SLICOT benchmarks] page: [http://slicot.org/objects/software/shared/bench-data/eady.zip eady.zip] and are stored as MATLAB [https://www.mathworks.com/help/matlab/import_export/mat-file-versions.html .mat] file.
The system matrices <math>A</math>, <math>B</math>, <math>C</math> are available from the [http://slicot.org/20-site/126-benchmark-examples-for-model-reduction SLICOT benchmarks] page: [http://slicot.org/objects/software/shared/bench-data/eady.zip eady.zip] and are stored as MATLAB [https://www.mathworks.com/help/matlab/import_export/mat-file-versions.html .mat] file.


Here is [https://www.python.org Python] code for loading the matrices (<math>B</math> and <math>C</math> are stored as arrays of 8-bit unsigned integers):
:<syntaxhighlight lang="python">
import numpy as np
from scipy.io import loadmat
mat = loadmat('eady.mat')
A = mat['A']
B = mat['B'].astype(np.float64)
C = mat['C'].astype(np.float64)
</syntaxhighlight>


==Dimensions==
==Dimensions==
Line 25: Line 37:


:<math>
:<math>
\begin{array}{rcl}
\begin{align}
\dot{x}(t) &=& Ax(t) + Bu(t) \\
  \dot{x}(t) &= A x(t) + B u(t) \\
y(t) &=& Cx(t)
  y(t) &= C x(t)
\end{array}
\end{align}
</math>
</math>


Line 36: Line 48:
<math>B \in \mathbb{R}^{598 \times 1}</math>,
<math>B \in \mathbb{R}^{598 \times 1}</math>,
<math>C \in \mathbb{R}^{1 \times 598}</math>.
<math>C \in \mathbb{R}^{1 \times 598}</math>.


==Citation==
==Citation==
Line 64: Line 75:
   doi =          {10.1175/1520-0469(1995)052<1642:SDOTMA>2.0.CO;2}
   doi =          {10.1175/1520-0469(1995)052<1642:SDOTMA>2.0.CO;2}
  }
  }


==References==
==References==

Revision as of 23:53, 1 September 2023


Description: Model of an Atmospheric Storm Track

This benchmark models the track of an atmospheric storm track. More details can be found in [1] and [2], [3].

Origin

This benchmark is part of the SLICOT Benchmark Examples for Model Reduction[3].

Data

The system matrices A, B, C are available from the SLICOT benchmarks page: eady.zip and are stored as MATLAB .mat file.

Here is Python code for loading the matrices (B and C are stored as arrays of 8-bit unsigned integers):

import numpy as np
from scipy.io import loadmat

mat = loadmat('eady.mat')
A = mat['A']
B = mat['B'].astype(np.float64)
C = mat['C'].astype(np.float64)

Dimensions

System structure:

x˙(t)=Ax(t)+Bu(t)y(t)=Cx(t)

System dimensions:

A598×598, B598×1, C1×598.

Citation

To cite this benchmark, use the following references:

  • For the benchmark itself and its data:
Niconet e.V., SLICOT - Subroutine Library in Systems and Control Theory, http://www.slicot.org
@MANUAL{slicot_cdisc,
 title =        {{SLICOT} - Subroutine Library in Systems and Control Theory},
 organization = {Niconet e.V.}
 address =      {\url{http://www.slicot.org}},
 key =          {SLICOT}
}
  • For the background on the benchmark:
@ARTICLE{FarI95,
 author =       {B.F. Farrell and P.J. Ioannou},
 title =        {Stochastic Dynamics of the Midlatitude Atmospheric Jet},
 journal =      {Journal of the Atmospheric Sciences},
 volume =       {52},
 number =       {10},
 pages =        {1642--1656},
 year =         {1995},
 doi =          {10.1175/1520-0469(1995)052<1642:SDOTMA>2.0.CO;2}
}

References

  1. B.F. Farrell, P.J. Ioannou. Stochastic Dynamics of the Midlatitude Atmospheric Jet. Journal of the Atmospheric Sciences, 52(10): 1642--1656, 1995.
  2. Y. Chahlaoui, P. Van Dooren, A collection of Benchmark examples for model reduction of linear time invariant dynamical systems, Working Note 2002-2: 2002.
  3. 3.0 3.1 Y. Chahlaoui, P. Van Dooren, Benchmark Examples for Model Reduction of Linear Time-Invariant Dynamical Systems, Dimension Reduction of Large-Scale Systems, Lecture Notes in Computational Science and Engineering, vol 45: 379--392, 2005.