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

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

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

Измените заголовок графика, чтобы считать "годограф Найквиста 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'}

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)

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

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

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

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

zoomcp(h)

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

Сравните частотные характеристики идентифицированных моделей в пространстве состояний порядка 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 показывает более высокую неопределенность в ее частотной характеристике, особенно близко к частоте Найквиста как показано графиком. Чтобы видеть это, сгенерируйте годограф Найквиста и покажите область уверенности в подмножестве точек, в которых отображен ответ Найквиста.

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

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