В этом примере показано, как получить числовые значения нескольких характеристик частотного диапазона модели динамической системы SISO, включая пиковое усиление, dc усиление, системная пропускная способность и частоты, на которых системное усиление пересекает заданную частоту.
Создайте модель передаточной функции и постройте ее частотную характеристику.
H = tf([10,21],[1,1.4,26]); bodeplot(H)
Графический вывод частотной характеристики дает общее представление о характеристиках частотного диапазона системы. H
включает явный резонирующий пик и прокручивается прочь на уровне 20 дБ/десятилетие в высокой частоте. Часто желательно получить определенные числовые значения для таких характеристик.
Вычислите пиковое усиление и частоту резонанса.
[gpeak,fpeak] = getPeakGain(H); gpeak_dB = mag2db(gpeak)
gpeak_dB = 17.7596
getPeakGain
возвращает обоих пиковое местоположение fpeak
и пиковое усиление gpeak
в абсолютных единицах. Используя mag2db
преобразовывать gpeak
к децибелам показывает, что усиление достигает максимума на уровне почти 18 дБ.
Найдите полосу, в которой системное усиление превышает 0 дБ, или 1 в абсолютных единицах.
wc = getGainCrossover(H,1)
wc = 2×1
1.2582
12.1843
getGainCrossover
возвращает вектор частот, на которых отклик системы пересекает заданное усиление. Получившийся wc
вектор показывает, что системное усиление превышает 0 дБ приблизительно между 1.3 и 12,2 рад/с.
Найдите dc усиление H
.
Предвещать график отклика показывает что усиление H
стремится к конечному значению как нуль подходов частоты. dcgain
команда находит это значение в абсолютных единицах.
k = dcgain(H);
Найдите частоту в который ответ H
списки прочь к-10 дБ относительно его dc значения.
fb = bandwidth(H,-10);
bandwidth
возвращает первую частоту, на которой отклик системы опускается ниже усиления dc заданным значением в дБ.
bandwidth
| getGainCrossover
| getPeakGain