Создайте Диаграмму Боде следующей непрерывно-разовой динамической системы SISO.

H = tf([1 0.1 7.5],[1 0.12 9 0 0]);
bode(H)

bode автоматически выбирает область значений графика на основе системной динамики.

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

H = tf([-0.1,-2.4,-181,-1950],[1,3.3,990,2600]);
bode(H,{1,100})
grid on

Массив ячеек {1,100} задает минимальные и максимальные значения частоты в Диаграмме Боде. Когда вы обеспечиваете границы частоты таким образом, функция выбирает промежуточные точки для данных о частотной характеристике.

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

w = [1 5 10 15 20 23 31 40 44 50 85 100];
bode(H,w,'.-')
grid on

bode строит частотную характеристику на заданных частотах только.

Сравните частотную характеристику непрерывно-разовой системы к эквивалентной дискретизированной системе на той же Диаграмме Боде.

Создайте динамические системы непрерывно-разового и дискретного времени.

H = tf([1 0.1 7.5],[1 0.12 9 0 0]);
Hd = c2d(H,0.5,'zoh');

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

bode(H,Hd)

Диаграмма Боде системы дискретного времени включает вертикальную разметку частота Найквиста системы.

Задайте стиль линии, цвет или маркер для каждой системы в Диаграмме Боде с помощью входного параметра LineSpec.

H = tf([1 0.1 7.5],[1 0.12 9 0 0]);
Hd = c2d(H,0.5,'zoh');
bode(H,'r',Hd,'b--')

Первый LineSpec, 'r', задает твердую красную линию для ответа H. Второй LineSpec, 'b--', задает пунктирную синюю строку для ответа Hd.

Вычислите значение и фазу частотной характеристики системы SISO.

Если вы не задаете частоты, bode выбирает частоты на основе системной динамики и возвращает их в третьем выходном аргументе.

H = tf([1 0.1 7.5],[1 0.12 9 0 0]);
[mag,phase,wout] = bode(H);

Поскольку H является моделью SISO, первые две размерности mag и phase оба 1. Третья размерность является количеством частот в wout.

size(mag)

ans = 1×3

     1     1    42

length(wout)

ans = 42

Таким образом каждая запись по третьему измерению mag дает значение ответа на соответствующей частоте в wout.

В данном примере создайте систему с 3 входами, с 2 выводами.

rng(0,'twister'); % For reproducibility
H = rss(4,2,3);

Для этой системы bode строит частотные характеристики каждого канала ввода-вывода в отдельном графике в одной фигуре.

bode(H)

Вычислите значение и фазу этих ответов на 20 частотах между 1 и 10 радианами.

w = logspace(0,1,20);
[mag,phase] = bode(H,w);

mag и phase являются 3D массивами, в которых первые две размерности соответствуют выводу и вводят размерности H, и третья размерность является количеством частот. Например, исследуйте размерности mag.

size(mag)

ans = 1×3

     2     3    20

Таким образом, например, mag(1,3,10) является значением ответа от третьего входа до первого вывода, вычисленного на 10-й частоте в w. Точно так же phase(1,3,10) содержит фазу того же ответа.

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

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

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 является параметрической идентифицированной моделью.

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

bode(sys_np,sys_p,w);
legend('sys-np','sys-p')

Можно отобразить область уверенности на Диаграмме Боде путем щелчка правой кнопкой по графику и выбора Characteristics> Confidence Region.

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

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

load iddata2 z2;
sys_p = tfest(z2,2);
w = linspace(0,10*pi,128);
[mag,ph,w,sdmag,sdphase] = bode(sys_p,w);

Используя tfest команда требует программного обеспечения System Identification Toolbox™.

sys_p является идентифицированной моделью передаточной функции. sdmag и sdphase содержат данные о стандартном отклонении для значения и фазы частотной характеристики, соответственно.

Используйте данные о стандартном отклонении, чтобы создать 3σ график, соответствующий области уверенности.

mag = squeeze(mag);
sdmag = squeeze(sdmag);
semilogx(w,mag,'b',w,mag+3*sdmag,'k:',w,mag-3*sdmag,'k:');