Этот пример показывает, как идентифицировать передаточную функцию, чтобы соответствовать данным данным о частотной характеристике (FRD), содержащим дополнительный список фазы от вызванного входной задержкой.
Этот пример требует лицензии Control System Toolbox™.
Получите данные о частотной характеристике.
В данном примере используйте bode, чтобы получить значение и данные о фазовом отклике для следующей системы:
Используйте 100 точек частоты, в пределах от от 0,1 рад/с до 10 рад/с, чтобы получить данные о частотной характеристике. Используйте frd, чтобы создать объект данных частотной характеристики.
freq = logspace(-1,1,100);
[mag, phase] = bode(tf([1 .2],[1 2 1 1],'InputDelay',.5),freq);
data = frd(mag.*exp(1j*phase*pi/180),freq);data является объектом iddata, который содержит данные о частотной характеристике для описанной системы.
Оцените передаточную функцию с помощью data. Задайте неизвестную транспортную задержку идентифицированной передаточной функции.
np = 3; nz = 1; iodelay = NaN; sys = tfest(data,np,nz,iodelay);
np и nz задают количество полюсов и нулей в идентифицированной передаточной функции, соответственно.
iodelay задает неизвестную транспортную задержку идентифицированной передаточной функции.
sys является моделью idtf, содержащей идентифицированную передаточную функцию.