Сгенерируйте диаграмму Боде.

sys = rss(5);
h = bodeplot(sys);

Измените модули на Гц и подавьте график фазы. Для этого отредактируйте свойства указателя на график, h.

setoptions(h,'FreqUnits','Hz','PhaseVisible','off');

График автоматически обновляется при вызове setoptions.

Создайте Диаграмму Боде динамической системы.

sys = tf(1,[1 1]); 
h = bodeplot(sys);

Зафиксировать фазу на уровне 1 рад/с до 750 степеней. Для этого получите свойства графика. Затем измените свойства PhaseMatchingFreq и PhaseMatchingValue для согласования фазы с заданной частотой.

p = getoptions(h); 
p.PhaseMatching = 'on'; 
p.PhaseMatchingFreq = 1; 
p.PhaseMatchingValue = 750;

Обновите график.

setoptions(h,p);

Первая диаграмма Боде имеет фазу -45 степеней на частоте 1 рад/с. Установка опций соответствия фазы так, чтобы при 1 рад/с фаза была близка к 750 степеням приводит к второй диаграмме Боде. Обратите внимание, что, однако, фаза может быть только -45 + N * 360, где N является целым числом, и поэтому график установлен на ближайшую допустимую фазу, а именно 675 степени (или 2 * 360 - 45 = 675).

Сравните частотные характеристики идентифицированных моделей пространства состояний порядка 2 и 6 вместе с их 2 $$\sigma$$ доверительные области.

load iddata1
sys1 = n4sid(z1, 2); 
sys2 = n4sid(z1, 6);

Обе модели дают около 70% подгонки данным. Однако sys2 показывает более высокую неопределенность в его частотной характеристике, особенно близкую к частоте Найквиста, как показано графиком:

w = linspace(8,10*pi,256);
h = bodeplot(sys1,sys2,w);
setoptions(h,'PhaseMatching','on','ConfidenceRegionNumberSD',2);

Щелкните правой кнопкой мыши график и выберите «Характеристики» > «Доверием область», чтобы включить характеристику доверия области. Кроме того, введите showConfidence(h) для построения графика области доверия.

Сравните частотную характеристику параметрической модели, идентифицированную из входных/выходных данных, с непараметрической моделью, идентифицированной с использованием тех же данных. Идентифицируйте параметрические и непараметрические модели на основе данных.

load iddata2 z2;
w = linspace(0,10*pi,128);
sys_np = spa(z2,[],w);
sys_p = tfest(z2,2);

spa и tfest требуется программное обеспечение System Identification Toolbox™. sys_np является непараметрической идентифицированной моделью. sys_p является параметрической идентифицированной моделью.

Создайте Диаграмму Боде, которая включает обе системы.

opt = bodeoptions;  
opt.PhaseMatching = 'on';
bodeplot(sys_np,sys_p,w,opt);