getrom

Obtain reduced-order models when using modal truncation method

Since R2023b

Syntax

rsys = getrom(R,Name=Value)

[rsys,info] = getrom(R,___)

rsys = getrom(R)

getrom(R,'-help')

Description

Use getrom to obtain reduced-order models from a ModalTruncation or SparseModalTruncation model order reduction task. For the full workflow, see Task-Based Model Order Reduction Workflow.

rsys = getrom(R,Name=Value) returns a reduced-order model rsys based on the options specified by one or more name-value arguments.

[rsys,info] = getrom(R,___) also returns a structure info containing additional information, such as the left and right projector matrices.

rsys = getrom(R) returns a modal decomposition of the original model sys.

• For the ordinary modal truncation, this syntax returns:

  • The modal form (same as modalreal) of sys when R.Options.ModeOnly is set to false

  • The Schur form when R.Options.ModeOnly is set to true.

• For sparse modal truncation, this syntax returns a reduced modal form based on the subset of modes computed by the algorithm. This subset is controlled by the R.Options.Focus or R.Options.MaxOrder options

getrom(R,'-help') returns help specific to the model order specification object R. The returned help shows the name-value arguments and syntaxes applicable to R.

Examples

Modal Truncation of LTI Models

This example shows how to obtain a reduced-order model for a linear time-invariant (LTI) model using the modal truncation method. In this example, you reduce a high-order model with a focus on the dynamics in a particular frequency range and damping.

Load a model and examine its frequency response.

load('highOrderModel.mat','G')
bodeplot(G)

G is a 48th-order model with several large peak regions around 5.2 rad/s, 13.5 rad/s, and 24.5 rad/s, and smaller peaks scattered across many frequencies.

Create a model order reduction task.

R = reducespec(G,"modal");
figure
view(R,"damp")

Suppose that for your application you are only interested in the dynamics between 10 rad/s and 40 rad/s with damping less than 0.04. Focus the model reduction on the region of interest to obtain a good match with a low-order approximation.

rsys = getrom(R,Frequency=[10,40],Damping=[0.022 0.04]);
bp = bodeplot(G,rsys,"r--");
bp.PhaseMatchingEnabled = "on";
legend("Full order","Reduced order");

The reduced-order model provides a good approximation for the specified targets.

Modal Truncation of Sparse State-Space Model

This example shows how to obtain a reduced-order model of a structural beam using the modal truncation method. For this example, consider a SISO sparse state-space model of a cantilever beam. This example uses the linearized model from the Linear Analysis of Cantilever Beam example.

Load the beam model.

load linBeam.mat
size(sys)

Sparse second-order model with 1 outputs, 1 inputs, and 3303 degrees of freedom.

Plot the Bode response.

bode(sys,w)

Create a model order reduction task.

R = reducespec(sys,"modal");

Set the Focus option to compute the eigenvalues only in the specified frequency range.

R.Options.Focus = [0 5e5];

Analyze the model and compute the derived information.

R = process(R);
size(R.Mode)

ans = 1×2

    26     1

The original model contains 26 modes. View the DC contribution of the modal components using the view function.

view(R,"contrib")

Discard the modes with smaller contributions. For this example, discard all modes with DC contribution smaller than 1e-6. This results in a reduced model with 12 modes.

[rsys,info] = getrom(R,MinDC=1e-6);

Plot the bode responses.

bode(sys,rsys,w)

The reduced-order model provides a good approximation for the original sparse model.

Input Arguments

R — Model order reduction specification object

ModalTruncation object | SparseModalTruncation object

Model order reduction specification object created using reducespec, specified as a ModalTruncation or SparseModalTruncation object.

Name-Value Arguments

Specify optional pairs of arguments as Name1=Value1,...,NameN=ValueN, where Name is the argument name and Value is the corresponding value. Name-value arguments must appear after other arguments, but the order of the pairs does not matter.

Example: rsys = getrom(R,Frequency=[1e-2,1e6],Method="matchDC")

Selection — Selected modes

vector

Selected modes for the reduced-order model, specified as a vector. You can select modes in one of the following ways.

• A vector containing index values into R.Modes.

• A logical vector with the same size as R.Modes.

The function retains the modes selected by these vectors. For example, if Selection is [1 2 3], it retains the modes R.Modes([1 2 3],:).

Frequency — Frequency range

vector

Frequency range of interest, specified as a vector of form [F_min,F_max]. The function discards all the modes outside this range.

Damping — Damping range

vector

Damping range of interest, specified as a vector of form [ζ_min,ζ_max]. The function discards all the modes outside this range.

MinDC — Minimum DC contribution bound

nonnegative scalar

Minimum DC contribution bound for the reduced-order model, specified as a nonnegative scalar. The function discards all the modes with normalized DC contributions smaller than this value.

Method — State elimination method

"matchDC" (default) | "truncate"

State elimination method, specified as "matchDC" or "truncate". This argument specifies how the function eliminates the states with weak contribution.

• The "matchDC" method preserves the DC gain (steady-state value of step response) and tends to match the time response better.

• The "truncate" method tends to match the frequency response better. Use this method when you want to control the relative error.

For the modal decomposition G(s) = G₁(s) + G₂(s) + D, where:

• G₁(s) contains the retained modes.

• G₂(s) contains the discarded modes.

• D is the feedthrough term.

The "truncate" method returns G_r(s) = G₁(s) + D, and "matchDC" returns G_r(s) = G₁(s) + D + G₂(0) so that G_r(0) = G(0).

Output Arguments

rsys — Reduced-order model

state-space model

Reduced-order model, returned as a state-space (ss) model.

info — Additional information

structure

Additional information about the reduced-order model, returned as a structure with these fields.

Field	Description
PL	Left projector matrix, returned as a matrix with dimensions Nx-by-Nxr, where Nx is the number of states in the original full-size model and Nxr is the number of states in the reduced-order model.
PR	Right projector matrix, returned as a matrix with dimensions Nx-by-Nxr, where Nx is the number of states in the original full-size model and Nxr is the number of states in the reduced-order model.

More About

Spectral Projector Matrices

For the model-order reduction algorithms (excluding NCF balanced truncation), the function also returns the left and right spectral projector matrices in the info argument. For both "truncate" and "matchDC" methods, the matrices A_r, B_r, C_r, D_r, and E_r of the reduced-order model are related to the matrices A, B, C, D, and E of the original model by:

• For explicit state-space models

  $$\begin{array}{cccc}{A}_r=P_L^{T}A{P}_R,& B_r=P_L^{T}B,& C_r=C{P}_r,& E_r=P_L^T{P}_R=I\end{array}$$

• For descriptor state-space models

  $$\begin{array}{cccc}{A}_r=P_L^{T}A{P}_R,& B_r=P_L^{T}B,& C_r=C{P}_r,& E_r=P_L^{T}E{P}_R\end{array}$$

Version History

Introduced in R2023b