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

Неуверенность объекта

Номинальная модель объекта управления состоит из нестабильной системы первого порядка.

Pnom = tf(2,[1 -2]);

Семейство встревоженных объектов является изменениями Pnom. Все объекты имеют один нестабильный полюс, но местоположение этого полюса отличается через семейство.

p1 = Pnom*tf(1,[.06 1]);              % extra lag
p2 = Pnom*tf([-.02 1],[.02 1]);       % time delay
p3 = Pnom*tf(50^2,[1 2*.1*50 50^2]);  % high frequency resonance
p4 = Pnom*tf(70^2,[1 2*.2*70 70^2]);  % high frequency resonance
p5 = tf(2.4,[1 -2.2]);                % pole/gain migration
p6 = tf(1.6,[1 -1.8]);                % pole/gain migration

Покрытие неопределенной модели

Поскольку обратная связь разрабатывает цели, мы должны заменить этот набор моделей с одной неопределенной моделью объекта управления, чья область значений поведений включает p1 через p6. Это - одно использование команды ucover. Эта команда берет массив моделей LTI Parray и номинальная модель Pnom и моделирует различие Parray-Pnom как мультипликативная неуверенность в системной динамике.

Поскольку ucover ожидает массив моделей, используйте команду stack, чтобы собрать модели объекта управления p1 через p6 в один массив.

Parray = stack(1,p1,p2,p3,p4,p5,p6);

Затем, используйте ucover, чтобы "покрыть" область значений поведений Parray с неопределенной моделью формы

P = Pnom * (1 + Wt * Delta)

где вся неуверенность сконцентрирована в "несмоделированной динамике" Delta (объект ultidyn). Поскольку усиление Delta однородно ограничено 1 на всех частотах, фильтр "формирования", который Wt используется, чтобы получить, как относительная сумма неуверенности меняется в зависимости от частоты. Этот фильтр также упоминается как функция взвешивания неуверенности. Попробуйте фильтр 4-го порядка Wt за этот пример:

orderWt = 4;
Parrayg = frd(Parray,logspace(-1,3,60));
[P,Info] = ucover(Parrayg,Pnom,orderWt,'InputMult');

Получившаяся модель P является одно входом, объектом неопределенного пространства состояний (USS) одно вывода с номинальной стоимостью Pnom.

P

P =

  Uncertain continuous-time state-space model with 1 outputs, 1 inputs, 5 states.
  The model uncertainty consists of the following blocks:
    Parrayg_InputMultDelta: Uncertain 1x1 LTI, peak gain = 1, 1 occurrences

Type "P.NominalValue" to see the nominal value, "get(P)" to see all properties, and "P.Uncertainty" to interact with the uncertain elements.

tf(P.NominalValue)

ans =
 
    2
  -----
  s - 2
 
Continuous-time transfer function.

Предвещать график значения подтверждает, что формирующий фильтр Wt "покрывает" относительное изменение в поведении объекта. Как функция частоты, уровень неуверенности составляет 30% в 5 радах/секунда (-10dB = 0.3), 50% в 10 радах/секунда и 100% вне 29 радов/секунда.

Wt = Info.W1;
bodemag((Pnom-Parray)/Pnom,'b--',Wt,'r'); grid
title('Relative Gaps vs. Magnitude of Wt')

Создание модели проекта разомкнутого цикла

Чтобы разработать устойчивый контроллер для неопределенной модели объекта управления P, мы выбираем желаемую пропускную способность с обратной связью и минимизируем чувствительность к воздействиям на объекте вывод. Управляющую структуру показывают ниже. Сигналы d и n являются воздействием загрузки и шумом измерения. Диспетчер использует шумное измерение объекта вывод y, чтобы сгенерировать управляющий сигнал u.

Рисунок 1: управляющая структура.

Фильтры Wperf и Wnoise выбраны, чтобы осуществить желаемую пропускную способность и некоторый соответствующий спад. Передаточная функция с обратной связью от [d;n] до y

    y = [Wperf * S , Wnoise * T] [d;n]

где S=1/(1+PC) и T=PC/(1+PC) являются чувствительностью и дополнительными функциями чувствительности. Если мы разрабатываем контроллер, который сохраняет усиление с обратной связью от [d;n] до y ниже 1, то

    |S| < 1/|Wperf| ,    |T| < 1/|Wnoise|

Путем выбора соответствующего значения профилирует для Wperf и Wnoise, мы можем осуществить маленькую чувствительность (S) в пропускной способности и соответствующем спаде (T) вне пропускной способности.

Например, выберите Wperf в качестве фильтра нижних частот первого порядка с усилением DC 500 и перекрестное соединение усиления в желаемой пропускной способности desBW:

desBW = 4.5;
Wperf = makeweight(500,desBW,0.33);
tf(Wperf)

ans =
 
  0.33 s + 4.248
  --------------
   s + 0.008496
 
Continuous-time transfer function.

Точно так же выберите Wnoise как фильтр высоких частот второго порядка со значением 1 в 10*desBW. Это обеспечит коэффициент усиления разомкнутого контура PC к спаду с наклоном-2 для частот вне 10*desBW.

NF = (10*desBW)/20;  % numerator corner frequency
DF = (10*desBW)*50;  % denominator corner frequency
Wnoise = tf([1/NF^2  2*0.707/NF  1],[1/DF^2  2*0.707/DF  1]);
Wnoise = Wnoise/abs(freqresp(Wnoise,10*desBW))

Wnoise =
 
    0.1975 s^2 + 0.6284 s + 1
  ------------------------------
  7.901e-05 s^2 + 0.2514 s + 400
 
Continuous-time transfer function.

Проверьте, что границы 1/Wperf и 1/Wnoise на S и T действительно осуществляют желаемую пропускную способность и спад.

bodemag(1/Wperf,'b',1/Wnoise,'r',{1e-2,1e3}), grid
title('Performance and roll-off specifications')
legend('Bound on |S|','Bound on |T|','Location','NorthEast')

Затем используйте connect, чтобы создать соединение разомкнутого цикла (блок-схема в рисунке 1 без блока контроллера). Задайте каждый блок, появляющийся в рисунке 1, назовите сигналы входящими и из каждого блока и позвольте connect сделать проводное соединение:

P.u = 'u';   P.y = 'yp';
Wperf.u = 'd';   Wperf.y = 'Wperf';
Wnoise.u = 'n';  Wnoise.y = 'Wnoise';
S1 = sumblk('e = -ym');
S2 = sumblk('y = yp + Wperf');
S3 = sumblk('ym = y + Wnoise');
G = connect(P,Wperf,Wnoise,S1,S2,S3,{'d','n','u'},{'y','e'});

G является неопределенной системой с 2 выводами, с 3 входами, подходящей для устойчивого синтеза контроллера с dksyn.

Устойчивый синтез контроллера

Проект выполняется с автоматизированной устойчивой командой проекта dksyn. Целевая пропускная способность составляет 4,5 рад/с.

ny = 1; nu = 1;
[C,CL,muPerf] = dksyn(G,ny,nu);
muPerf

muPerf = 1.1487

Когда устойчивый показатель эффективности, muPerf близок 1, контроллер, достигает целевой пропускной способности с обратной связью и спада. Как показывает опыт, если muPerf - меньше чем 0,85, то производительность может быть улучшена, и если muPerf больше, чем 1,2, то желаемая пропускная способность с обратной связью не достижима для заданной неуверенности объекта.

Здесь muPerf - приблизительно 1, таким образом, цели достигнуты. У получившегося контроллера C есть 18 состояний:

size(C)

State-space model with 1 outputs, 1 inputs, and 18 states.

Можно использовать reduce и команды dksynperf, чтобы упростить этот контроллер. Вычислите приближения порядков 1 - 17.

Cappx = reduce(C,1:17);

Для каждого контроллера уменьшаемого порядка используйте dksynperf, чтобы вычислить устойчивый показатель эффективности и сравнить его с muPerf. Сохраните контроллер самый низкоуровневый с производительностью не хуже, чем 1,05 * muPerf, ухудшение производительности 5% или меньше.

for k=1:17
   Cr = Cappx(:,:,k);  % controller of order k
   bnd = dksynperf(lft(G,Cr));
   if bnd.UpperBound < 1.05 * muPerf
      break % abort with the first controller meeting the performance goal
   end
end

order(Cr)

ans = 7

Это приводит к контроллеру 6-го порядка Cr с сопоставимой производительностью. Сравните Cr с контроллером полного порядка C.

opt = bodeoptions;
opt.Grid = 'on';
opt.PhaseMatching = 'on';
bodeplot(C,'b',Cr,'r--',opt)
legend('18-state controller C','6-state controller Cr','Location','NorthEast')

Устойчивая валидация контроллера

Постройте ответы разомкнутого цикла моделей объекта управления p1 через p6 с упрощенным контроллером Cr.

bodeplot(Parray*Cr,'g',{1e-2,1e3},opt);

Постройте ответы на воздействие шага на объекте вывод. Они сопоставимы с желаемой пропускной способностью с обратной связью и устойчивы к изменениям объекта, как ожидалось от Устойчивого mu-значения Производительности приблизительно 1.

step(feedback(1,Parray*Cr),'g',10/desBW);

Отличаясь целевая пропускная способность с обратной связью

Тот же процесс проектирования может быть повторен для различных значений пропускной способности с обратной связью desBW. Выполнение так приводит к следующим результатам: