Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Restructure sys-mor tutorial #1141

Merged
merged 10 commits into from
Nov 30, 2020
Merged

Restructure sys-mor tutorial #1141

merged 10 commits into from
Nov 30, 2020

Conversation

pmli
Copy link
Member

@pmli pmli commented Oct 29, 2020

This is still a work in progress. The idea is to add an LTI tutorial and for the existing BT tutorial to reuse the model from there. An IRKA tutorial could be added in a new PR

@pmli pmli added the pr:change Change in existing functionality label Oct 29, 2020
@pmli pmli added this to the 2020.2 milestone Oct 29, 2020
@pmli pmli force-pushed the sys-mor-tutorials branch from db752b5 to b043f91 Compare November 26, 2020 21:42
@pmli pmli marked this pull request as ready for review November 26, 2020 21:46
@pmli pmli requested review from sdrave and lbalicki November 26, 2020 21:46
Copy link
Member

@sdrave sdrave left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Very nice!


.. jupyter-execute::

fom.gramian('c_lrcf')
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This seems a bit cryptic here. Maybe you can add a sentence about the argument and what other values are possible.

+ C^{\operatorname{T}} C
& = 0.
\end{align*}

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Maybe one sentence (and perhaps a reference), for what these Gramians are useful.

Co-authored-by: Stephan Rave <stephanrave@uni-muenster.de>
Copy link
Member

@lbalicki lbalicki left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Agreed, very nice new structure!

invertible matrix.
In particular, there exist invertible transformation matrices
:math:`T, S \in \mathbb{R}^{n \times n}` such that the realization with
:math:`\tilde{E} = S^{\operatorname{T}} E T = I`,
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This only holds if E is invertible right? Maybe it should be mentioned that we consider an LTI system with invertible E matrix at the beginning of the paragraph.

:math:`\hat{B} = W^{\operatorname{T}} B`,
:math:`\hat{C} = C V`.

It is known that the reduced-order model is asymptotically stable if
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think it should also be mentioned that we assume the full order model to be asymptotically stable. Maybe we should do this when the LTI system is introduced in this tutorial.

\sup_{u \neq 0}
\frac{\lVert y \rVert_{\mathcal{L}_\infty}}{\lVert u \rVert_{\mathcal{L}_2}}.

The computation of the :math:`\mathcal{H}_2` is based on the system Gramians
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

:math:\mathcal{H}_2 norm

@pmli
Copy link
Member Author

pmli commented Nov 29, 2020

Thanks for the comments, I made updates.

@pmli pmli requested review from sdrave and lbalicki November 29, 2020 21:39
@pmli pmli merged commit 39b872b into master Nov 30, 2020
@pmli pmli deleted the sys-mor-tutorials branch November 30, 2020 09:04
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
pr:change Change in existing functionality
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants