hsvd

(Not recommended) Hankel singular values of dynamic system

collapse all in page

hsvdis not recommended. Usebalredinstead. For more information, seeCompatibility Considerations.

Syntax

hsv= hsvd(sys) hsv= hsvd(sys,opts) [hsv,baldata] = hsvd(___) hsvd(___)

Description

hsv= hsvd(sys)computes the Hankel singular valueshsvof thedynamic systemsys. In state coordinates that equalize the input-to-state and state-to-output energy transfers, the Hankel singular values measure the contribution of each state to the input/output behavior. Hankel singular values are to model order what singular values are to matrix rank. In particular, small Hankel singular values signal states that can be discarded to simplify the model (seebalred).

For models with unstable poles,hsvdonly computes the Hankel singular values of the stable part and entries ofhsvcorresponding to unstable modes are set toInf.

hsv= hsvd(sys,opts)computes the Hankel singular values using options that you specify usinghsvdOptions. Options include offset and tolerance options for computing the stable-unstable decompositions. The options also allow you to limit the HSV computation to energy contributions within particular time and frequency intervals. SeebalredOptionsfor details.

[hsv,baldata] = hsvd(___)returns additional data to speed up model order reduction. You can use this syntax with any of the previous combinations of input arguments.

hsvd(___)displays a Hankel singular values plot.

Examples

collapse all

Compute Hankel Singular Values of System With Near-Unstable Pole

Open Live Script

Create a system with a stable pole very near to 0, and display the Hankel singular values.

sys = zpk([1 2],[-1 -2 -3 -10 -1e-7],1); hsv = hsvd(sys)

hsv =5×110⁵× 1.6667 0.0000 0.0000 0.0000 0.0000

Notice the dominant Hankel singular value with magnitude $1 0^{5}$ , which is so much larger that the significant digits of the other modes are not displayed. This value is due to the near-unstable mode at $s = 1 0^{- 7}$ . Use the'Offset'option to treat this mode as unstable.

opts = hsvdOptions('Offset',1e-7); hsvu = hsvd(sys,opts)

hsvu =5×1Inf 0.0688 0.0138 0.0024 0.0001

The Hankel singular value of modes that are unstable, or treated as unstable, is returned asInf. Create a Hankel singular-value plot while treating this mode as unstable.

hsvd(sys,opts)

ans =5×1Inf 0.0688 0.0138 0.0024 0.0001

The unstable mode is shown in red on the plot.

By default,hsvduses a linear scale. To switch the plot to a log scale, right-click on the plot and selectY Scale > Log. For information about programmatically changing properties of HSV plots, seehsvplot.

Frequency-Limited Hankel Singular Values

Open Live Script

Compute the Hankel singular values of a model with low-frequency and high-frequency dynamics. Focus the calculation on the high-frequency modes.

Load the model and examine its frequency response.

loadmodeselectGmsbodeplot(Gms)

图包含2轴对象。坐标轴对象1ylabel Magnitude (dB) contains an object of type line. This object represents Gms. Axes object 2 with ylabel Phase (deg) contains an object of type line. This object represents Gms.

Gmshas two sets of resonances, one at relatively low frequency and the other at relatively high frequency. Compute the Hankel singular values of the high-frequency modes, excluding the energy contributions to the low-frequency dynamics. To do so, usehsvdOptionsto specify a frequency interval above 30 rad/s.

opts = hsvdOptions('FreqInterval',[30 Inf]); hsvd(Gms,opts)

ans =18×110^-4× 0.6237 0.4558 0.3183 0.2468 0.0895 0.0847 0.0243 0.0028 0.0000 0.0000 ⋮

Tips

To create a Hankel singular-value plot with more flexibility to programmatically customize the plot, usehsvplot.

Version History

Introduced before R2006a

expand all

R2021a:`hsvd`is not recommended

Starting in R2021a, use thebalredcommand to compute Hankel singular values. This approach facilitates better workflow since you can also compute reduced-order model approximations using the same command. Furthermore, the enhanced workflow also includes new options that preserve roll-off characteristics.

The following table shows some typical uses ofhsvdand how to update your code to usebalredinstead.

Not Recommended	Recommended
`hsv = hsvd(sys)`	`[~,info] = balred(sys,order)`computes the Hankel singular values and the error approximations of`sys`. For more information, see`balred`.
`hsv = hsvd(sys,opts)`	`[~,info] = balred(sys,order,opts)`computes the Hankel singular values with options specified in the option set`opts`. For more information, see`balred`.

There are no plans to removehsvdat this time.