view

Plot state contributions when using balanced truncation of normalized coprime factors method

Since R2023b

Syntax

view(R,type)

view(R)

view(___,Parent=parent)

view(___,Name=Value)

hsvp = view(___)

view(R,"-help")

Description

Use view to graphically analyze the model and select a model order reduction criteria from a model order reduction task created using reducespec. For NCFBalancedTruncation objects, you can visualize the state contributions of the normalized coprime factorization of the full-order model as either Hankel singular values or normalized state energies. For the full workflow, see Task-Based Model Order Reduction Workflow.

view(R,type) creates a plot that helps you select the order of the reduced model. R specifies the model order reduction (MOR) specification object. Use the type argument to specify the plot type as either "sigma" or "energy".

view(R) plots the default plot type for the model order reduction algorithm of R. For balanced truncation methods, this syntax plots Hankel singular values and associated error bounds.

view(___,Parent=parent) creates a plot in the specified parent graphics container, such as a Figure or TiledChartLayout. Use this syntax when you want to create a plot in a specified open figure or when creating apps in App Designer. You can specify the parent container after any of the input argument combinations in the previous syntaxes.

view(___,Name=Value) specifies additional options for customizing the appearance of Hankel singular value plots using one or more name-value arguments. For example, view(R,"sigma",YScale="Linear") plots the Hankel singular values using a linear scale for the y axis.

hsvp = view(___) returns an HSVPlot object. To customize your plot, modify properties of the plot object using dot notation.

view(R,"-help") returns help specific to the model order specification object R. The returned help shows plot types and syntaxes applicable to R.

Examples

Balanced Truncation of Normalized Coprime Factors

This example shows how to obtain a reduced-order model using the balanced truncation of normalized coprime factors method.

Load a 30-state plant model G.

load ncfModel.mat G
size(G)

State-space model with 2 outputs, 3 inputs, and 30 states.

Create a model order reduction task.

R = reducespec(G,"ncf");

To help you select a suitable target reduction order, examine the plot of Hankel singular values and approximation errors.

Create the plot.

view(R)

The function generates a Hankel singular value plot, which shows the relative energy contributions of each state in the coprime factorization of G, arranged in decreasing order by energy. The plot also shows the upper bound on the error between the original and reduced-order models that you obtain by truncating the states at that point. Examine this plot to choose the target order. For instance, for a maximum error of 0.01, you can reduce the model to 13th order.

Examine the singular values of G and of the difference between G and Gred.

Gred = getrom(R,MaxError=0.01);
sigma(G,G-Gred)
legend("G","G-Gred")

The difference is small across all frequencies, showing that the reduced-order model is a good approximation of the full-order model.

Customize State Contribution Plots

This example shows how to customize the state contribution plots obtained using the view function in the model order reduction workflow.

For this example, create a model order reduction specification for an LTI model using the balanced truncation of normalized coprime factors method.

Generate a random discrete-time state-space model with 40 states.

rng(0)
sys = drss(40);

Create a specification object.

R = reducespec(sys,"ncf");

Visualize the Hankel singular values.

view(R)

To customize the plots for balanced truncation methods, you can use hsvoptions properties as input arguments.

h = view(R,"sigma",YScale="linear");

You can further customize this plot by modifying the properties of h.

h.AxesStyle.GridVisible = "off";

Input Arguments

R — Model order reduction specification object

NCFBalancedTruncation object

Model order reduction specification object created using reducespec, specified as a NCFBalancedTruncation object.

type — Plot type

"sigma" | "energy"

Plot type, specified as one of the following.

• "sigma" — Bar chart of Hankel singular values and associated error bound

• "energy" — Bar chart of normalized energies

If you do not specify this argument, the function plots the Hankel singular values and associated error bounds of the normalized coprime factorization of the original model sys.

parent — Parent graphics container

Figure object | TiledChartLayout object | UIFigure object | UIGridLayout object | UIPanel object | UITab object | Axes object | UIAxes object

Parent graphics container, specified as one of these objects:

• Figure

• TiledChartLayout

• UIFigure

• UIGridLayout

• UIPanel

• UITab

You can also specify parent as an Axes or UIAxes object, which sets the plot parent to the parent of the specified axes object.

Output Arguments

hsvp — HSV plot object

HSVPlot object

HSV plot object, returned as an HSVPlot chart object. To customize the plot, modify the properties of hsvp using dot notation.

Version History

Introduced in R2023b

R2025a: Axes=AX input argument not recommended

Specifying the axes for plotting using Axes=AX is not recommended. Specify a parent graphics container using Parent=parent instead.

For example, the following code plots the response in an axes object obtained from a figure.

f = figure;
ax = gca;
view(R,Axes=ax)

To update the code, set the parent of the plot to the figure instead.

f = figure;
view(R,Parent=f)

R2025a: Specify parent graphics container

You can now plot state contributions in a parent graphics container, such as a Figure or TiledChartLayout object, using the parent argument.

Specifying parent as an Axes or UIAxes object sets the plot parent to the parent of the specified axes object.