Подгонка данных частотной характеристики с помощью модели пространства состояний
B = fitfrd(A,N) B = fitfrd(A,N,RD) B = fitfrd(A,N,RD,WT)
B = fitfrd(A,N) является объектом пространства состояний с размерностью состояний N, где A является frd объект и N является неотрицательным целым числом. Частотная характеристика B тесно соответствует D - масштабировать данные частотной характеристики в A.
A должна иметь 1 строку или 1 столбец, хотя он не должен быть 1 на 1. B будет того же размера, что и A. Во всех случаях N должен быть неотрицательным скаляром.
B = fitfrd(A,N,RD) увеличивает относительную степень B чтобы быть RD. RD должно быть неотрицательным целым числом. Значение по умолчанию для RD равен 0. Если A - строка (или столбец), затем RD может быть вектором того же размера, задавая относительную степень каждой записи B. Если RD является скаляром, затем задает относительную степень для всех записей B. Можно задать значение по умолчанию для RD путем установки RD в пустую матрицу.
B = fitfrd(A,N,RD,WT) использует величину WT для взвешивания критериев подгонки оптимизации. WT может быть double, ss или frd. Если WT является скаляром, затем используется для взвешивания всех записей критериев ошибки (A-B). Если WT является вектором, он должен быть того же размера, что и A, и каждый отдельный ввод WT действует как функция взвешивания на соответствующем элементе (A-B).
Используйте fitfrd команда для подгонки данных шкалы D.
Создайте данные частотной характеристики D-шкалы из системы пятого порядка.
sys = tf([1 2 2],[1 2.5 1.5])*tf(1,[1 0.1]);
sys = sys*tf([1 3.75 3.5],[1 2.5 13]);
omeg = logspace(-1,1);
sysg = frd(sys,omeg);
bode(sysg,'r-');
Можно попытаться подогнать данные D-шкалы частотной характеристики sysg с системой первого порядка, b1. Точно так же можно подгонять данные шкалы D к системе третьего порядка, b3.
b1 = fitfrd(sysg,1); b3 = fitfrd(sysg,3);
Сравните исходные данные шкалы D sysg с частотными характеристиками модели первого и третьего порядка, рассчитанными по fitfrd.
b1g = frd(b1,omeg); b3g = frd(b3,omeg); bode(sysg,'r-',b1g,'k:',b3g,'b-.') legend('5th order system','1st order fit','3rd order fit','Location','Southwest')
Числовые проблемы обусловленности возникают, если порядок состояний подгонки N выбран таким образом, чтобы он был выше, чем требуется динамикой A.