hsvdOptions

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

свернуть все на странице

hsvdOptions не рекомендуется. Использование balredOptions вместо этого. Для получения дополнительной информации см. Вопросы совместимости.

Синтаксис

opts = hsvdOptions
opts = hsvdOptions(Name,Value)

Описание

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

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

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

Аргументы name-value

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

FreqIntervals

Интервалы частоты для вычисления ограниченных частотой сингулярных значений Ганкеля (balred) или сбалансированная реализация (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

Временные интервалы для вычисления ограниченных временем сингулярных значений Ганкеля (balred) или сбалансированная реализация (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 значение, затем расчет использует объединение этих интервалов.

SepTol

Максимальная потеря значения точности в устойчивом и нестабильном разложении. Для моделей с нестабильными полюсами, balred первые извлечения устойчивое использование динамики stabsep. Используйте 'SepTol'чтобы управлять точностью разложения.

Для получения дополнительной информации смотрите balredOptions.

Значение по умолчанию: SepTol = 'auto'

Offset

Возместите для устойчивого/нестабильного контура. Значение положительной скалярной величины. В устойчивом/нестабильном разложении устойчивый термин включает только удовлетворение полюсов:

  • Re(s) < -Offset * max(1,|Im(s)|) (Непрерывное время)

  • |z| < 1 - Offset (Дискретное время)

Увеличьте значение Offset обрабатывать полюса близко к контуру устойчивости как нестабильные.

Значение по умолчанию: 1e-8

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

Примеры

свернуть все

Скрипт Open Live Script

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

sys=(s+0.5)(s+10-6)(s+2)

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

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

       Inf
    0.1875

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

Скрипт Open Live Script

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

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

load modeselect Gms
bodeplot(Gms)

Figure contains 2 axes objects. Axes object 1 contains an object of type line. This object represents Gms. Axes object 2 contains an object of type line. This object represents Gms.

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

opts = hsvdOptions('FreqInterval',[30 Inf]);
hsvd(Gms,opts)
ans = 18×1
10-4 ×

    0.6237
    0.4558
    0.3183
    0.2468
    0.0895
    0.0847
    0.0243
    0.0028
    0.0000
    0.0000
      ⋮

Вопросы совместимости

развернуть все

Не рекомендуемый запуск в R2021a

Ссылки

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

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

| | |

Введен в R2010a

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

Поддержка

Для просмотра документации необходимо авторизоваться на сайте