hsvdOptions

Создайте набор опции для вычислительных сингулярных значений Ганкеля и балансировки ввода/вывода

Синтаксис

opts = hsvdOptions opts = hsvdOptions(Name,Value)

Описание

opts = hsvdOptions возвращает опции по умолчанию для hsvd и balreal команды.

opts = hsvdOptions(Name,Value) возвращает набор опций с опциями, заданными одним или несколькими Name,Value парные аргументы.

Входные параметры

Аргументы в виде пар имя-значение

Задайте разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в одинарных кавычках (' 'Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

`'FreqIntervals'`	Интервалы частоты для вычисления ограниченных частотой сингулярных значений Ганкеля (`hsvd`) или сбалансированная реализация (`balreal`) в виде матрицы с двумя столбцами. Каждая строка задает интервал частоты `[fmin fmax]`, где `fmin` и `fmax` неотрицательные частоты, выраженные в единице частоты модели. Например: Ограничить расчет областью значений между 3 рад/с и 15 рад/с, принимая единицу частоты модели - rad/s, установите `FreqIntervals` к `[3 15]`. Ограничить расчет двумя интервалами частоты, 3-15 рад/с и 40-60 рад/с, `[3 15; 40 60]` использования. Задавать все частоты ниже частоты среза `fcut`, используйте `[0 fcut]`. Чтобы задать все частоты выше сокращения, используйте `[fcut Inf]` в непрерывное время или `[fcut pi/Ts]` в дискретное время, где `Ts` шаг расчета модели. Значение по умолчанию, `[]`, не налагает ограничения частоты и эквивалентен `[0 Inf]` в непрерывное время или `[0 pi/Ts]` в дискретное время. Однако, если вы задаете `TimeIntervals` значение кроме `[]`, затем этот предел заменяет `FreqIntervals = []`. Если вы задаете обоих `TimeIntervals` значение и `FreqIntervals` значение, затем расчет использует объединение этих интервалов. Значение по умолчанию: `[]`
`'TimeIntervals'`	Временные интервалы для вычисления ограниченных временем сингулярных значений Ганкеля (`hsvd`) или сбалансированная реализация (`balreal`) в виде матрицы с двумя столбцами. Каждая строка задает временной интервал `[tmin tmax]`, где `tmin` и `tmax` неотрицательные времена, выраженные в единице измерения времени модели. Программное обеспечение вычисляет вклады состояния в импульсную характеристику системы в этих временных интервалах только. Например: Ограничить расчет областью значений между 3 с и 15 с, принимая единицу измерения времени модели - секунды, установите `TimeIntervals` к `[3 15]`. Ограничить расчет двумя временными интервалами, 3-15 с и 40-60 с, `[3 15; 40 60]` использования. Задавать все случаи от нуля до времени среза `tcut`, используйте `[0 tcut]`. Чтобы задать все случаи после сокращения, используйте `[tcut Inf]`. Значение по умолчанию, `[]`, не налагает ограничения времени и эквивалентен `[0 Inf]`. Однако, если вы задаете `FreqIntervals` значение кроме `[]`, затем этот предел заменяет `Timeintervals = []`. Если вы задаете обоих `TimeIntervals` значение и `FreqIntervals` значение, затем расчет использует объединение этих интервалов.
`'AbsTol, RelTol'`	Допуск абсолютной и относительной погрешности к устойчивому/нестабильному разложению. Значения положительной скалярной величины. Для входной модели G с нестабильными полюсами, `hsvd` и `balreal` сначала извлеките устойчивую динамику путем вычисления устойчивого/нестабильного разложения G → GS + GU. `AbsTol` и `RelTol` допуски управляют точностью этого разложения путем гарантирования, что частотные характеристики G и GS + GU отличаются не больше, чем `AbsTol` + `RelTol`abs (G). Увеличение этих допусков помогает разделить поблизости устойчивые и нестабильные режимы за счет точности. Смотрите `stabsep` для получения дополнительной информации. Значение по умолчанию:* `AbsTol = 0; RelTol = 1e-8`
`'Offset'`	Возместите для устойчивого/нестабильного контура. Значение положительной скалярной величины. В устойчивом/нестабильном разложении устойчивый термин включает только удовлетворение полюсов: `Re(s) < -Offset * max(1,\|Im(s)\|)` (Непрерывное время) `\|z\| < 1 - Offset` (Дискретное время) Увеличьте значение `Offset` обрабатывать полюса близко к контуру устойчивости как нестабильные. Значение по умолчанию: `1e-8`

Для получения дополнительной информации об опциях и как они влияют на вычисление, смотрите hsvd. Ограниченные временем и ограниченные частотой вклады состояния вычисляются с помощью ограниченной временем и ограниченной частотой управляемости и наблюдаемости грамиана, как описано в gram и в [1].

Примеры

свернуть все

График сингулярного значения Ганкеля с почти нестабильным полюсом

Скрипт Open Live Script

Вычислите сингулярные значения Ганкеля системы, данной:

$s y s = \frac{(s + 0.5)}{(s + 10^{- 6}) (s + 2)}$

Используйте Offset опция, чтобы обеспечить hsvd исключить полюс в $s = 1 0^{- 6}$ от устойчивого срока устойчивого/нестабильного разложения.

sys = zpk(-.5,[-1e-6 -2],1); 
opts = hsvdOptions('Offset',.001); 
hsvd(sys,opts)

График показывает, что одно состояние обработано как нестабильное. hsvd вычисляет энергетические вклады устойчивых состояний только.

Ограниченные частотой сингулярные значения Ганкеля

Скрипт Open Live Script

Вычислите сингулярные значения Ганкеля модели с низкочастотной и высокочастотной динамикой. Фокусируйте вычисление на высокочастотных режимах.

Загрузите модель и исследуйте ее частотную характеристику.

load modeselect Gms
bodeplot(Gms)

Gms имеет два набора резонансов, один в относительно низкой частоте и другом в относительно высокой частоте. Вычислите сингулярные значения Ганкеля высокочастотных режимов, исключая энергетические вклады в низкочастотную динамику. Для этого используйте hsvdOptions задавать интервал частоты выше 30 рад/с.

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

Ссылки

[1] Гавронский, W. и Дж.Н. Джуэнг. “Снижение сложности модели в Интервалах Ограниченного времени и Частоты”. Международный журнал Системной Науки. Издание 21, Номер 2, 1990, стр 349–376.

Смотрите также

balreal | balred | gram | hsvd

Документация