Этот пример показывает, как получить числовые значения нескольких характеристик частотного диапазона модели динамической системы 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.7579
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