Вычислите ограничение системы по данным LTI ss массив
Примечание
ltiarray2uss будет удалено в следующем релизе. Использовать ucover вместо этого.
usys = ltiarray2uss(P,Parray,ord)
[usys,wt] = ltiarray2uss(P,Parray,ord)
[usys,wt,diffdata] = ltiarray2uss(P,Parray,ord)
[usys,wt,diffdata] = ltiarray2uss(P,Parray,ord,'InputMult')
[usys,wt,diffdata] = ltiarray2uss(P,Parray,ord,'OutputMult')
[usys,wt,diffdata] = ltiarray2uss(P,Parray,ord,'Additive')
Область команды ltiarray2ussвычисляет неопределенную систему usys с номинальным значением P, и область значений поведения включает в себя заданный массив систем, Parray.
usys = ltiarray2uss(P,Parray,ord), usys сформулирован как вход мультипликативная модель неопределенности,
usys = P*(I + wt*ultidyn('IMult',[size(P,2) size(P,2)])), где wt является стабильной скалярной системой, величина переграждает относительное различие, (P - Parray)/P. Порядок состояний функции взвешивания, используемой для связывания мультипликативного различия между P и Parray является ord. Оба P и Parray должен быть в классах ss/tf/zpk/frd. Если P является frd затем usys будет ufrd объект, в противном случае usys будет uss объект. The ultidyn атом называется на основе имени переменной Parray в вызов рабочей области.
[usys,wt] = ltiarray2uss(P,Parray,ord), возвращает значение wt используется для связывания нормы по бесконечности ((P - Parray)/P).
[usys,wt] = ltiarray2uss(P,Parray,ord,'OutputMult'), использует мультипликативную неопределенность на выходе объекта (в отличие от входной мультипликативной неопределенности). Формула для usys является
usys = (I + wt*ultidyn('Name',[size(P,1) size(P,1)])*P).
[usys,wt] = ltiarray2uss(P,Parray,ord,'Additive'), использует аддитивную неопределенность.
usys = P + wt*ultidyn('Name',[size(P,1) size(P,2)]). wt является переграждением частотного диапазона нормы по бесконечности (Parray - P).
[usys,wt] = ltiarray2uss(P,Parray,ord,'InputMult'), использует мультипликативную неопределенность на входе объекта управления (это значение по умолчанию). Формула для usys является usys = P*(I + wt*ultidyn('Name',[size(P,2) size(P,2)])) .
[usys,wt,diffdata] = ltiarray2uss(P,Parray,ord,type) возвращает норму различия (абсолютное различие для аддитива и относительное различие для мультипликативной неопределенности) между номинальной моделью P и Parray. wt удовлетворяет diffdata(w_i) < |wt(w_i)| во всех частотных точках.
Рассмотрим передаточную функцию третьего порядка с неопределенным усилением, постоянной времени фильтра и слегка демпфированным гибким режимом. Эта модель используется, чтобы представлять физическую систему, из которой получены данные частотной характеристики.
gain = ureal('gain',10,'Perc',20); tau = ureal('tau',.6,'Range',[.42 .9]); wn = 40; zeta = 0.1; usys = tf(gain,[tau 1])*tf(wn^2,[1 2*zeta*wn wn^2]); sysnom = usys.NominalValue; parray = usample(usys,30); om = logspace(-1,2,80); parrayg = frd(parray,om); bode(parrayg)
Данные частотной характеристики в parray представляет собой 30 экспериментов, выполненных в системе. Область команды ltiarray2uss используется для генерации неопределенной модели, umod, на основе данных частотной характеристики. Первоначально вход мультипликативная неопределенная модель используется, чтобы охарактеризовать набор 30 частотных характеристик. Первый и вторые порядки вычисляют входной мультипликативный вес неопределенности из данных.
[umodIn1,wtIn1,diffdataIn] = ltiarray2uss(sysnom,parrayg,1); [umodIn2,wtIn2,diffdataIn] = ltiarray2uss(sysnom,parrayg,2); bodemag(wtIn1,'b-',wtIn2,'g+',diffdataIn,'r.',om) title('Input Multiplicative Uncertainty Model Using ltiarray2uss') legend('1st order','2nd order','difference','Location','SouthEast')
Альтернативно, аддитивная неопределенная модель используется, чтобы охарактеризовать сбор 30 частотных характеристик.
[umodAdd1,wtAdd1,diffdataAdd] = ltiarray2uss(sysnom,parrayg,1,'Additive'); [umodAdd2,wtAdd2,diffdataAdd] = ltiarray2uss(sysnom,parrayg,2,'Additive'); bodemag(wtAdd1,'b-',wtAdd2,'g+',diffdataAdd,'r.',om) title('Additive Uncertainty Model Using ltiarray2uss') legend('1st order','2nd order','difference')
