nyquistplot

Годограф Найквиста с дополнительными опциями индивидуальной настройки графика

Синтаксис

h = nyquistplot(sys)
nyquistplot(sys,{wmin,wmax})
nyquistplot(sys,w)
nyquistplot(sys1,sys2,...,w)
nyquistplot(AX,...)
nyquistplot(..., plotoptions)

Описание

h = nyquistplot(sys) рисует Годограф Найквиста модели динамической системы sys. Также возвращается указатель на график h. Можно использовать этот указатель для настройки графика с помощью getoptions и setoptions команды. Напечатать

help nyquistoptions 

для получения списка доступных опций графика.

Частотная область значений и число точек выбираются автоматически. Посмотрите bode для получения дополнительной информации о понятии частоты в дискретном времени.

nyquistplot(sys,{wmin,wmax}) рисует годограф Найквиста для частот между wmin и wmaxrad/TimeUnit, где TimeUnit - временные модули входа динамической системы, заданная в TimeUnit свойство sys).

nyquistplot(sys,w) использует пользовательский вектор w частот (в rad/TimeUnit, где TimeUnit - временные модули входа динамической системы, заданная в TimeUnit свойство sys), при котором должен быть оценен ответ Nyquist. Посмотрите logspace для генерации логарифмически разнесенных векторов частоты.

nyquistplot(sys1,sys2,...,w) рисует Годографы Найквиста нескольких моделей sys1, sys2,... на одном графике. Вектор частоты w опционально. Можно также задать цвет, стиль линии и маркер для каждой системы, как в

nyquistplot(sys1,'r',sys2,'y--',sys3,'gx')

nyquistplot(AX,...) графики в осях с указателем AX.

nyquistplot(..., plotoptions) строит график отклика Nyquist с опциями, заданными в plotoptions. Напечатать

help nyquistoptions 

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

Примеры

Настройка заголовка годографа Найквиста

Создайте Годограф Найквиста модели динамической системы и сохраните указатель на график.

sys = tf(100,[1,2,1]);
h = nyquistplot(sys);

Figure contains an axes. The axes contains 2 objects of type line. This object represents sys.

Измените заголовок графика следующим образом: «Годограф Найквиста of sys». Для этого используйте getoptions чтобы извлечь существующие опции графика из указателя на график h.

opt = getoptions(h)
opt =

                         FreqUnits: 'rad/s'
                          MagUnits: 'dB'
                        PhaseUnits: 'deg'
                   ShowFullContour: 'on'
          ConfidenceRegionNumberSD: 1
    ConfidenceRegionDisplaySpacing: 5
                        IOGrouping: 'none'
                       InputLabels: [1x1 struct]
                      OutputLabels: [1x1 struct]
                      InputVisible: {'on'}
                     OutputVisible: {'on'}
                             Title: [1x1 struct]
                            XLabel: [1x1 struct]
                            YLabel: [1x1 struct]
                         TickLabel: [1x1 struct]
                              Grid: 'off'
                         GridColor: [0.1500 0.1500 0.1500]
                              XLim: {[-20 100]}
                              YLim: {[-80 80]}
                          XLimMode: {'auto'}
                          YLimMode: {'auto'}

The Title опция является структурой с несколькими полями.

opt.Title
ans = struct with fields:
         String: 'Nyquist Diagram'
       FontSize: 11
     FontWeight: 'bold'
      FontAngle: 'normal'
          Color: [0 0 0]
    Interpreter: 'tex'

Измените String поле Title структура и использование setoptions для применения изменения к графику.

opt.Title.String = 'Nyquist Plot of sys';
setoptions(h,opt)

Figure contains an axes. The axes contains 2 objects of type line. This object represents sys.

Масштабирование критической точки

Постройте график частотной характеристики Найквиста динамической системы. Присвойте имя переменной указателю на график, чтобы вы могли получить к нему доступ для дальнейшей манипуляции.

sys = tf(100,[1,2,1]);
h = nyquistplot(sys);

Figure contains an axes. The axes contains 2 objects of type line. This object represents sys.

Изменение масштаба критической точки (-1,0). Вы можете сделать это в интерактивном режиме, щелкнув правой кнопкой мыши на графике и выбрав «Масштабирование» (-1,0). Кроме того, используйте zoomcp команда на указателе на графике h.

zoomcp(h)

Figure contains an axes. The axes contains 2 objects of type line. This object represents sys.

Постройте графики идентифицированных моделей с доверительными областями в выбранных точках

Сравните частотные характеристики идентифицированных моделей пространства состояний порядка 2 и 6 вместе с их 1-std Доверие области, отображенные на каждой 50-й частотной выборке.

load iddata1
sys1 = n4sid(z1,2); % discrete-time IDSS model of order 2
sys2 = n4sid(z1,6); % discrete-time IDSS model of order 6

Обе модели дают около 76% подгонки данным. Однако sys2 показывает более высокую неопределенность в его частотной характеристике, особенно близкую к частоте Найквиста, как показано графиком. Чтобы увидеть это, сгенерируйте годограф Найквиста и отобразите доверительную область в подмножестве точек, в которых отображается реакция Nyquist.

w = linspace(10,10*pi,256);
h = nyquistplot(sys1,sys2,w);
setoptions(h,'ConfidenceRegionDisplaySpacing',50,...
             'ShowFullContour','off');

Чтобы включить отображение доверия области, щелкните правой кнопкой мыши график и выберите Characteristics > Confidence Region.

Совет

  • Можно изменить свойства графика, например модули. Для получения информации о способах изменения свойств графиков смотрите Способы настроить графики.

  • В контекстном меню доступны две опции масштаба, которые применяются специально к годографам Найквиста:

    • Full View - Клипы неограниченные ветви годографа Найквиста, но все еще включают критическую точку (-1, 0).

    • Zoom on (-1,0) - Масштабирование вокруг критической точки (-1,0). Для программного доступа к масштабу критической точки используйте zoomcp команда. См. «Изменение масштаба критической точки».

  • Чтобы активировать маркеры данных, которые отображают вещественные и мнимые значения на заданной частоте, щелкните где угодно на кривой. Следующий рисунок показывает nyquist постройте график с помощью маркера данных.

См. также

| |

Представлено до R2006a