Подходящие данные о частотной характеристике с моделью в пространстве состояний
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 or 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.