Визуализируйте находящиеся на диске запасы устойчивости
diskmarginplot( строит находящиеся на диске запасы по амплитуде и фазе для SISO или цикла отрицательной обратной связи MIMO L)feedback(L,eye(N)), где N количество вводов и выводов в ответе разомкнутого цикла L.
Для ответов MIMO, diskmarginplot строит многоконтурные дисковые поля. Находящийся на диске запас по амплитуде на каждой частоте является ±GM, где GM является значением, показанным в графике в дБ. Точно так же находящийся на диске запас по фазе является ±PM степенями, где PM является значением, показанным на графике. Для получения дополнительной информации о находящихся на диске запасах по амплитуде и фазе, смотрите diskmargin.
diskmarginplot(___, строит вычисленное использование находящихся на диске запасов по амплитуде и фазе эксцентриситета E)E смещать изменение усиления к увеличению усиления (E > 0) или уменьшение усиления (E <0). Если вы использовали diskmargin получить находящиеся на диске поля с некоторым конкретным E, можно использовать этот синтаксис, чтобы видеть зависимость частоты полей в том E значение. Для E ≠ 0, нанесенным на график значением является GM = min(gmax,1/max(0,gmin)). Другими словами, график показывает, что самая большая сумма усиления изменяет [1/GM,GM] это соответствует в находящемся на диске запасе по амплитуде [gmin,gmax] из системы в заданном E.
diskmarginplot(___, строит поля на частотах, заданных w)w.
Если w массив ячеек формы {wmin,wmax}, затем график показывает поля на частотах, располагающихся между wmin и wmax.
Если w вектор частот, затем график показывает поля на каждой заданной частоте.
diskmarginplot( строит область значений одновременного усиления и изменений фазы, соответствующих находящемуся на диске запасу по амплитуде. График также показывает максимальные изменения только для фазы и только для усиления (это находящиеся на диске запасы по амплитуде и фазе). DGM)DGM вектор формы [gmin,gmax]. Чтобы построить области значений для нескольких находящихся на диске запасов по амплитуде целиком, используйте матрицу 2D столбца формы [gmin1,gmax1;...;gminN,gmaxN]. Для получения дополнительной информации о находящихся на диске запасах по амплитуде, смотрите diskmargin.
diskmarginplot( строит диск с комплексным знаком смоделированного усиления и изменений фазы, соответствующих находящемуся на диске запасу по амплитуде DGM,'disk')DGM. Для получения дополнительной информации, о как граничное диском аналитическое усиление моделей и изменения фазы, смотрите, что Анализ Устойчивости Использует Дисковые Поля.
diskmarginplot( строит область исключения в соответствии плоскости Найквиста дисковому полю DGM,'nyquist')DGM. Требование, чтобы система с обратной связью осталась устойчивой для усиления или изменений фазы в диске, соответствующем DGM суммы к требованию, чтобы ответ разомкнутого цикла остался за пределами дискообразной области исключения в плоскости Найквиста. Для получения дополнительной информации смотрите, что Анализ Устойчивости Использует Дисковые Поля.
diskmarginplot( строит область значений усиления и изменений фазы, соответствующих размеру диска alpha,E,___)alpha и эксцентриситет E. Если любой alpha или E вектор, затем график включает области значений для всех заданных значений. Если оба alpha и E векторы, затем график включает области значений для пар alpha1,E1;...;alphaN,EN.
Построить смоделированный диск неопределенности, соответствующий alphaE, используйте этот синтаксис с 'disk' флаг. Чтобы построить соответствующий диск исключения в плоскости Найквиста, используйте этот синтаксис с 'nyquist' флаг.
diskmarginplot( строит график на осях, заданных указателем осей (AX,___)axes или UIaxes объект. Используйте этот аргумент, чтобы задать оси, когда создавание приложений в Разрабатывает Приложения Используя App Designer (MATLAB). Можно использовать этот аргумент с любым из предыдущих синтаксисов, чтобы построить любой из графиков что diskmarginplot может сгенерировать.
Постройте находящиеся на диске запасы по амплитуде и фазе системы с ответом разомкнутого цикла L и ответ с обратной связью feedback(L,1).
L = tf(25,[1 10 10 10]); diskmarginplot(L)
Находящимся на диске запасом по амплитуде на каждой частоте является, по крайней мере, ±GM, где GM значение, показанное в графике в дБ. Точно так же находящимся на диске запасом по фазе является ±PM степени. Находящиеся на диске поля возвращены diskmargin команда является самыми маленькими полями по частоте. (Щелкните правой кнопкой по графику и выберите Characteristics> Minimum Disk Margin для всплывающей подсказки, содержащей информацию об этих минимальных полях.)
Сравните находящиеся на диске запасы по амплитуде и фазе двух ответов разомкнутого цикла на том же графике.
L1 = tf(25,[1 10 10 10]); L2 = tf([1 100],[1 10 20 50]); diskmarginplot(L1,L2)
Можно настроить внешний вид графиков с помощью LineSpec аргумент. Постройте поля снова, с помощью синей пунктирной линии в L1 и красные точки для L2.
diskmarginplot(L1,'b--',L2,'r.')
Находящиеся на диске запасы по амплитуде возвращены diskmargin не примите изменение фазы, и запасы по фазе не принимают изменения усиления. На практике системы испытывают и усиление и изменение фазы одновременно. Граничный диском анализ составляет такое одновременное усиление и изменение фазы.
Вычислите находящиеся на диске запасы по амплитуде и фазе системы. Затем используйте diskmarginplot визуализировать соответствующие позволенные области значений одновременного усиления и изменения фазы.
L = tf(25,[1 10 10 10]); DM = diskmargin(L); DGM = DM.GainMargin
DGM = 1×2
0.6273 1.5942
DPM = DM.PhaseMargin
DPM = 1×2
-25.8017 25.8017
diskmarginplot(DGM)
Теневая область в графике показывает область значений одновременного усиления и изменений фазы, которые сохраняют устойчивость системы с обратной связью feedback(L,1). Когда нет никакого изменения фазы, система может терпеть полную область значений изменения усиления DGM, от-4 дБ до 4 дБ. Если фазе позволяют варьироваться приблизительно ±20 градусами, допустимые спады области значений изменения усиления приблизительно до-2.5 дБ к 2 дБ. В полной области значений изменения фазы ±26 градусов система не может терпеть изменение усиления.
Чтобы визуализировать несколько областей значений на том же графике, объедините их в 2D вектор-столбец. Например, вычислите дисковые поля L с положительным и отрицательным эксцентриситетом и постройте все три области значений изменения вместе.
DMn = diskmargin(L,-2); DGMn = DMn.GainMargin; DMp = diskmargin(L,2); DGMp = DMp.GainMargin; DGMall = [DGMn;DGM;DGMp]; diskmarginplot(DGMall)
Этот график показывает, что обратная связь может терпеть большее усиление и изменения фазы, когда усиление уменьшается. Другими словами, устойчивость цикла более чувствительна, чтобы получить увеличение. Для получения дополнительной информации о том, как различный эксцентриситет влияет на находящиеся на диске оценки запаса по амплитуде и фазе, смотрите, что Анализ Устойчивости Использует Дисковые Поля.
Находящееся на диске аналитическое усиление моделей запаса по амплитуде и изменение фазы как неопределенный фактор F, умножая коэффициент усиления разомкнутого контура L. (Для получения дополнительной информации этой модели, смотрите, что Анализ Устойчивости Использует Дисковые Поля.) Устойчивая область значений усиления и изменений фазы, возвращенных diskmargin команда эквивалентна диску значений F, для которых замкнутый цикл устойчив. Когда вы получаете находящиеся на диске запасы по амплитуде и фазе системы, можно использовать diskmarginplot визуализировать соответствующий диск значений F.
L = tf(25,[1 10 10 10]); DM = diskmargin(L); DGM = DM.GainMargin
DGM = 1×2
0.6273 1.5942
diskmarginplot(DGM,'disk')
График показывает значения F в комплексной плоскости, где ось X является действительной частью, и ось Y является мнимой частью. Находящийся на диске запас по амплитуде DGM исключительно определяет этот диск и соответствующий находящийся на диске запас по фазе, DPM.
Диск F неопределенности параметрируется двумя значениями: α, который устанавливает размер диска и E, который смещает изменение усиления к увеличению усиления или уменьшению. Этой параметризацией дают:
,
где δ является нормированной неопределенностью. (Для получения дополнительной информации смотрите, что Анализ Устойчивости Использует Дисковые Поля.) Для данной области значений изменения усиления можно использовать gm2dm преобразовывать находящийся на диске запас по амплитуде DGM к α и значениям E, которые описывают соответствующий диск. diskmarginplot может построить диск F для данного α, E пара.
DGM = [0.8,2]
DGM = 1×2
0.8000 2.0000
[alpha,E] = gm2dm(DGM);
diskmarginplot(alpha,E,'disk')
В этом случае E больше нуля потому что находящаяся на диске область значений усиления DGM = [0.8,2] включает больше увеличения усиления, чем уменьшение. E = 0 представляет усиление, которое может увеличиться так, как оно может уменьшиться. E <0 представляет область значений с большим количеством уменьшения, чем увеличение. Постройте диск F для различных значений E, чтобы видеть, как область значений усиления (диаметр диска) меняется в зависимости от E.
E = [-3,0,3];
diskmarginplot(alpha,E,'disk')
Для фиксированного E α управляет размером диска, и следовательно суммой смоделированной неопределенности. Постройте диск для нескольких значений α в E = 0.
E = 0;
alpha = [0.1,0.333,0.5];
diskmarginplot(alpha,E,'disk')
Как описано в Анализе Устойчивости Используя Дисковые Поля, для данного E, самого большого диска F неопределенности, для который система с обратной связью feedback(L*F,1) остается устойчивым, может быть интерпретирован как область исключения, которую не может ввести кривая Найквиста L. Для любого значения E диск исключения содержит критическую точку (–1,0) и является касательной к кривой Найквиста. Эксцентриситет настраивает размер и положение дисков касательной. Можно использовать diskmarginplot визуализировать эти диски исключения и накладывать их на кривую Найквиста L.
Вычислите находящиеся на диске запасы по амплитуде для системы с помощью трех различных значений эксцентриситета, одно соответствие большему количеству увеличения усиления, чем уменьшение (E> 0), одно соответствие большему количеству уменьшения усиления, чем увеличение (E <0), и один сбалансированный (E = 0).
L = tf(25,[1 10 10 10]); DMdec = diskmargin(L,-2); DGMdec = DMdec.GainMargin; DM = diskmargin(L,0); DGM = DM.GainMargin; DMinc = diskmargin(L,2); DGMinc= DMinc.GainMargin;
Чтобы просмотреть соответствующие области исключения, постройте годографы Найквиста L и содержите фигуру. Затем используйте diskmarginplot с 'nyquist' отметьте, чтобы добавить области исключения в график.
nyquist(L) hold on diskmarginplot([DGMdec;DGM;DGMinc],'nyquist') hold off
Как E увеличения от –2 до 2, дисковое перемещение направо и каждый диск обеспечивает более низкие оценки классических запасов по амплитуде и фазе.
Постройте дисковые поля как функцию частоты системы со следующим ответом разомкнутого цикла.
L = tf(25,[1 10 10 10]);
Для графика задайте следующие атрибуты:
Единицы частоты: Гц
Запасы по амплитуде на логарифмической шкале, в абсолютных единицах
Сетка на
opts = diskmarginoptions; opts.FreqUnits = 'Hz'; opts.MagScale = 'log'; opts.MagUnits = 'abs'; opts.grid = 'on'; diskmarginplot(L,opts)
L — Ответ разомкнутого циклаОтвет разомкнутого цикла в виде модели динамической системы. L может быть SISO или MIMO, пока это имеет то же количество вводов и выводов. diskmarginplot строит находящиеся на диске запасы по амплитуде и фазе для отрицательной обратной связи система с обратной связью feedback(L,eye(N)).
Построить поля системы положительной обратной связи feedback(L,eye(N),+1), используйте diskmargin(-L).
Если L неопределенная модель в пространстве состояний (uss или genss с неопределенными блоками), затем diskmarginplot строит поля случайных выборок L. Чтобы визуализировать запасы устойчивости худшего случая неопределенной системы, используйте wcdiskmarginplot.
Если L модель данных частотной характеристики (такая как frd), затем diskmarginplot строит поля на каждой частоте, представленной в модели.
Если L массив моделей, затем diskmarginplot строит поля для всех моделей в массиве на той же оси и в том же стиле линии.
E — ЭксцентриситетЭксцентриситет области неопределенности использовался для расчета запасов устойчивости в виде действительного скаляра или вектора (для diskmarginplot(alpha,E) только графики).
Этот параметр смещает неопределенность, используемую, чтобы смоделировать усиление и изменения фазы к увеличению усиления или уменьшению усиления.
E по умолчанию = 0 использования сбалансированная модель изменения усиления области значений [gmin,gmax], с gmin = 1/gmax.
Положительный E использует модель с большим количеством увеличения усиления, чем уменьшение (gmin > 1/gmax).
Отрицательный E использует модель с большим количеством уменьшения усиления, чем увеличение (gmin < 1/gmax).
Для более подробной информации, о как выбор E влияет на граничный расчет, смотрите, что Анализ Устойчивости Использует Дисковые Поля.
diskmargin(L,E) СинтаксисПри графическом выводе запасов по амплитуде динамической системы по сравнению с частотой используйте E по умолчанию = 0, чтобы получить объективные оценки запасов по амплитуде и фазе. Для E = 0, находящийся на диске запас по амплитуде на каждой частоте является ±GM, где GM является значением, показанным в графике в дБ.
Если вы использовали diskmargin получить находящиеся на диске поля с некоторым конкретным E, можно использовать этот синтаксис, чтобы видеть зависимость частоты полей в том E значение. Для E ≠ 0, нанесенным на график значением является GM = min(gmax,1/max(0,gmin)). Другими словами, график показывает, что самая большая сумма усиления изменяет [1/GM,GM] это соответствует в находящемся на диске запасе по амплитуде [gmin,gmax] из системы в заданном E.
Для синтаксиса diskmarginplot(L,E), эксцентриситет E должен быть скаляр.
diskmargin(alpha,E) СинтаксисВ синтаксисе diskmargin(alpha,E), графики функций диск неопределенности, параметрированный значениями alpha и E (см., что Анализ Устойчивости Использует Дисковые Поля). Преобразовывать между находящимися на диске запасами по амплитуде и alphaE параметризация, используйте dm2gm и gm2dm. Для этого синтаксиса, diskmarginplot осуществляет отношение alpha*abs(1+E) < 2.
Для этого синтаксиса, эксцентриситет E может быть вектор, позволив вам сравнить несколько дисков на том же графике. Если alpha скаляр и E вектор, затем график показывает диски, соответствующие парам alpha,Ek для каждой записи в E. Если оба alpha и E векторы, затем график показывает диски для пар alpha1,E1;...;alphaN,EN.
w — Частоты{wmin,wmax} | векторЧастоты, на которых можно построить запасы устойчивости в виде массива ячеек {wmin,wmax} или как вектор значений частоты.
Если w массив ячеек формы {wmin,wmax}, затем график показывает поля на частотах между wmin и wmax.
Если w вектор частот, затем график показывает поля на каждой заданной частоте. Например, используйте logspace сгенерировать вектор-строку с логарифмически расположенными с интервалами значениями частоты.
Задайте частоты в модулях rad/TimeUnit, где TimeUnit TimeUnit свойство L.
opts — Постройте опцииdiskmarginplot опции установленыПостройте опции в виде diskmarginplot набор опций, который вы создаете с diskmarginoptions. Элементы, которые можно настроить, включают заголовок графика, метки осей и сетки.
DGM — Находящийся на диске запас по амплитудеНаходящийся на диске запас по амплитуде в виде двухэлементного вектора формы [gmin,gmax], или матрица 2D столбца формы [gmin1,gmax1; ...; gminN,gmaxN].
Векторный [gmin,gmax], находящийся на диске запас по амплитуде, такой как запас по амплитуде, который вы вычисляете с diskmargin. Это представляет допустимую область значений относительного изменения коэффициента усиления разомкнутого контура, таким образом, что система с обратной связью остается устойчивой, не принимая изменения фазы. В DGM, gmin <1 и gmax > 1.
Используя DGM = [gmin,gmax] строит допустимую область значений одновременного усиления и изменений фазы, соответствующих находящемуся на диске запасу по амплитуде [gmin,gmax]. Чтобы показать области значений нескольких находящихся на диске запасов по амплитуде на том же графике, используйте DGM = [gmin1,gmax1;...;gminN,gmaxN].
Можно получить [gmin,gmax] в GainMargin поле структур output diskmargin команда. Можно также получить [gmin,gmax] от классических запасов по амплитуде и фазе с помощью getDGM.
alpha — Размер усиления и изменения фазыРазмер смоделированного усиления и изменения фазы в виде скаляра или вектора. Находящееся на диске аналитическое усиление моделей запаса по амплитуде и изменение фазы как мультипликативная неопределенность F, который является диском значений, содержащих F = 1, соответствуя номинальной стоимости системы. Диск параметрируется alpha, который устанавливает размер диска и E, который смещает изменение усиления к увеличению усиления или уменьшению. (Для получения дополнительной информации этой модели, смотрите, что Анализ Устойчивости Использует Дисковые Поля.)
diskmarginplot позволяет вам построить диск F или область усиления и изменений фазы, представленных конкретным alphaE пары. Для синтаксиса diskmarginplot(alpha,E), размер диска alpha может быть вектор. Если alpha вектор и E скаляр, затем график показывает области, соответствующие парам alpha_k,E для каждой записи в alpha. Если оба alpha и E векторы, затем график показывает области для пар alpha1,E1;...;alphaN,EN.
diskmarginplot осуществляет отношение alpha*abs(1+E) <2. Преобразовывать между находящимися на диске запасами по амплитуде и alphaE параметризация, используйте dm2gm и gm2dm.
AX — Указатель осейaxes возразите | UIAxes объектОси обрабатывают в виде axes возразите или UIaxes объект. Используйте этот аргумент, чтобы задать оси, когда создавание приложений в Разрабатывает Приложения Используя App Designer (MATLAB).
diskmargin | diskmarginoptions | getDGM | gm2dm | umargin | wcdiskmarginplot