Этот пример показывает, как идентифицировать передаточную функцию для соответствия заданным данным частотной характеристики (FRD), содержащим дополнительный фазовый накат, вызванный входной задержкой.
В этом примере требуется лицензия на Toolbox™ системы управления.
Получение данных частотной характеристики.
Для этого примера используйте bode для получения данных о величине и фазовой характеристике для следующей системы:
2s2 + s + 1
Для получения данных частотной характеристики используются 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 модель, содержащая идентифицированную передаточную функцию.