fitmagfrd

Подходящие данные о значении частотной характеристики с моделью в пространстве состояний минимальной фазы с помощью чебышевского журналом проекта значения

свернуть все на странице

Синтаксис

B = fitmagfrd(A,N)
B = fitmagfrd(A,N,RD)
B = fitmagfrd(A,N,RD,WT)
B = fitmagfrd(A,N,RD,WT,C)

Описание

B = fitmagfrd(A,N) стабильный, объект ss минимальной фазы, с размерностью состояния N, значение частотной характеристики которого тесно совпадает с данными о значении в A. A является объектом frd 1 на 1, и N является неотрицательным целым числом.

B = fitmagfrd(A,N,RD) обеспечивает относительную степень B, чтобы быть RD. RD должен быть неотрицательным целым числом, значением по умолчанию которого является 0. Можно задать значение по умолчанию для RD установкой RD к пустой матрице.

B = fitmagfrd(A,N,RD,WT) использует значение WT, чтобы взвесить критерии подгонки оптимизации. WT может быть double, ss или frd. Если WT является скаляром, то он используется, чтобы взвесить все записи ошибочных критериев (A-B). Если WT является вектором, он должен быть одного размера как A и каждая отдельная запись действий WT как функция взвешивания на соответствующей записи (A-B). Значение по умолчанию для WT равняется 1, и можно задать его установкой WT к пустой матрице.

B = fitmagfrd(A,N,RD,WT,C) осуществляет дополнительные ограничения значения на B, заданный значениями C.LowerBound и C.UpperBound. Они могут быть пустыми, double или frd (с C.Frequency, равным A.Frequency). Если C.LowerBound непуст, то значение B ограничивается лечь выше C.LowerBound. Никакая нижняя граница не осуществляется на частотах, где C.LowerBound равен-inf. Точно так же поле UpperBound может использоваться, чтобы задать верхнюю границу на значении B. Если C является double или frd (с C.Frequency, равным A.Frequency), то ограничения верхней и нижней границы на B взяты непосредственно из A как:

если C (w) == –1, то осуществите abs (B (w)) <= abs (A (w))
если C (w) == 1, то осуществите abs (B (w))> = abs (A (w))
если C (w) == 0, то никакое дополнительное ограничение

где w обозначает частоту.

Примеры

свернуть все

Подходящие данные о частотной характеристике со стабильной моделью в пространстве состояний Минимальной Фазы

Скрипт Open Live Script

Создайте данные о значении частотной характеристики из системы пятого порядка.

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]); 
omega = logspace(-1,1); 
sysg = abs(frd(sys,omega)); 
bodemag(sysg,'r');

Соответствуйте данным о значении минимальной фазой, стабильной системой третьего порядка.

ord = 3; 
b1 = fitmagfrd(sysg,ord); 
b1g = frd(b1,omega); 
bodemag(sysg,'r',b1g,'k:');
legend('Data','3rd order fit');

Соответствуйте данным о значении системой третьего порядка, ограниченной лечь ниже и выше определенных данных.

C2.UpperBound = sysg;
C2.LowerBound = [];
b2 = fitmagfrd(sysg,ord,[],[],C2); 
b2g = frd(b2,omega); 
C3.UpperBound = [];
C3.LowerBound = sysg;
b3 = fitmagfrd(sysg,ord,[],[],C3); 
b3g = frd(b3,omega); 
bodemag(sysg,'r',b1g,'k:',b2g,'b-.',b3g,'m--') 
legend('Data','3rd order fit','3rd order fit, below data',...
       '3rd order fit, above data')

Соответствуйте данным о значении системой второго порядка, ограниченной лечь ниже и выше определенных данных.

ord = 2;
C2.UpperBound = sysg;
C2.LowerBound = [];
b2 = fitmagfrd(sysg,ord,[],sysg,C2);
b2g = frd(b2,omega);
C3.UpperBound = [];
C3.LowerBound = sysg;
b3 = fitmagfrd(sysg,ord,[],sysg,C3);
b3g = frd(b3,omega);
bgp = fitfrd(genphase(sysg),ord);
bgpg = frd(bgp,omega);
bodemag(sysg,'r',b1g,'k:',b2g,'b-.',b3g,'m--',bgpg,'r--')
legend('Data','3rd order fit','2d order fit, below data',...
       '2nd order fit, above data','bgpg')

Ограничения

Этот объект входа frd должен быть или скаляром, 1 на 1 возражают или, строка или вектор-столбец.

Алгоритмы

fitmagfrd использует версию чебышевского журналом проекта значения, решая

   min f     subject to (at every frequency point in A):  
           |d|^2 /(1+ f/WT) < |n|^2/A^2 < |d|^2*(1 + f/WT)

плюс дополнительные ограничения, наложенные с C. n, d обозначает числитель и знаменатель, соответственно, и B = n/d. n и d имеют порядки (N-RD) и N, соответственно. Проблема решена с помощью линейного программирования для фиксированного деления пополам f and, чтобы минимизировать f. Альтернативным приближенным методом, который не может осуществить ограничения, заданные C, является B = fitfrd(genphase(A),N,RD,WT).

Ссылки

Оппенхейм, A.V., и Р.В. Шаффер, Цифровая обработка сигналов, Prentice Hall, Нью-Джерси, 1975, p. 513.

Бойд, S. и Vandenberghe, L., выпуклая оптимизация, издательство Кембриджского университета, 2004.

Смотрите также

fitfrd

Документация