ncfmargin

Вычислите нормированный взаимно-простой запас устойчивости обратной связи контроллера объекта

Синтаксис

[marg,freq] = ncfmargin(P,C)
[marg,freq] = ncfmargin(P,C,sign)
[marg,freq] = ncfmargin(___,tol)

Описание

пример

[marg,freq] = ncfmargin(P,C) возвращает нормированный взаимно-простой запас устойчивости многомерной обратной связи, состоящей из контроллера C в отрицательной обратной связи с объектом P:

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

[marg,freq] = ncfmargin(P,C,sign) задает знак связи обратной связи, принятой для граничного вычисления. По умолчанию, sign = -1. Установите sign = +1 для соединения позитивных откликов.

[marg,freq] = ncfmargin(___,tol) вычисляет нормированную взаимно-простую факторную метрику с заданной относительной точностью.

Примеры

свернуть все

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

p = tf(4,[1 -0.001]); 	
cL = 1;				
cH = 10;

Вычислите запас устойчивости системы с обратной связью с контроллером низкого усиления.

[margL,~] = ncfmargin(p,cL)
margL = 0.7069

Точно так же вычислите запас устойчивости системы с обратной связью с контроллером высокого усиления.

[margH,~] = ncfmargin(p,cH)
margH = 0.0995

Системы с обратной связью с низким усилением и контроллерами высокого усиления нормировали взаимно-простые запасы устойчивости приблизительно 0,71 и 0.1, соответственно. Этот результат показывает, что система с обратной связью с контроллером низкого усиления более устойчива к неструктурированным возмущениям, чем система с контроллером высокого усиления.

Чтобы наблюдать это различие в робастности, создайте неопределенный объект, punc, который имеет дополнительную несмоделированную динамику в высокой частоте по сравнению с номинальным объектом.

punc = p + ultidyn('uncstruc',[1 1],'Bound',1);
sigma(p,punc,'r--')

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

[stabmargL,~] = robstab(feedback(punc,cL))
stabmargL = struct with fields:
           LowerBound: 0.9980
           UpperBound: 1
    CriticalFrequency: Inf

[stabmargH,~] = robstab(feedback(punc,cH))
stabmargH = struct with fields:
           LowerBound: 0.0998
           UpperBound: 0.1000
    CriticalFrequency: Inf

Как ожидалось устойчивый анализ устойчивости показывает, что система с обратной связью с контроллером низкого усиления более надежно стабильна в присутствии несмоделированной динамики LTI. На самом деле эта система с обратной связью может терпеть почти 100% заданной неуверенности. Напротив, система с обратной связью с контроллером высокого усиления может терпеть только приблизительно 10% заданной неуверенности.

Рассмотрите объект и стабилизировавшийся контроллер.

P1 = tf([1 2],[1 5 10]);
C = tf(4.4,[1 0]);

Вычислите запас устойчивости для этого объекта и контроллера.

b1 = ncfmargin(P1,C)
b1 = 0.1961

Затем, вычислите разрыв между P1 и встревоженным объектом, P2.

P2 = tf([1 1],[1 3 10]);
[gap,nugap] = gapmetric(P1,P2)
gap = 0.1391
nugap = 0.1390

Поскольку запас устойчивости, b1 = b(P1,C) больше, чем разрыв между этими двумя объектами, C также, стабилизирует P2. Как обсуждено в Метриках Разрыва и Запасах устойчивости, запасе устойчивости b2 = b(P2,C) удовлетворяет неравенство asin(b(P2,C)) ≥ asin(b1)-asin(gap). Подтвердите этот результат.

b2 = ncfmargin(P2,C);
[asin(b2) asin(b1)-asin(gap)]
ans = 1×2

    0.0997    0.0579

Входные параметры

свернуть все

Объект, заданный как модель динамической системы. P может быть SISO или MIMO, пока P*C имеет то же количество вводов и выводов. P может быть непрерывным временем или дискретным временем. Если P является обобщенной моделью в пространстве состояний (genss или uss) затем, ncfmargin использует текущее значение или номинальную стоимость всех блоков системы управления в P.

Объект, заданный как модель динамической системы. C может быть SISO или MIMO, пока P*C имеет то же количество вводов и выводов. C может быть непрерывным временем или дискретным временем. Если C является обобщенной моделью в пространстве состояний (genss или uss) затем, ncfmargin использует текущее значение или номинальную стоимость всех блоков системы управления в P.

По умолчанию ncfmargin принимает соединение отрицательной обратной связи между P и C. Чтобы вычислить поля для системы с обратной связью с позитивными откликами, используйте [marg,freq] = ncfmargin(P,C,+1).

Знак связи обратной связи, заданной или как -1 или как +1.

Значение по умолчанию, sign = -1, задает отрицательную обратную связь. Установка sign = +1 принимает связь позитивных откликов для граничного вычисления, как в следующей схеме.

Относительная точность для вычисленного поля, заданного как значение положительной скалярной величины меньше чем 1. Значение по умолчанию 0.001, или точность на 0,1%.

Выходные аргументы

свернуть все

Нормированный взаимно-простой устойчивый запас устойчивости, возвращенный как скаляр в области значений [0,1]. Это количество, также известное как метрический запас устойчивости разрыва, является индикатором робастности с обратной связью к неструктурированным возмущениям. Для архитектуры управления отрицательной обратной связи Это задано как:

b(P,C)=[IC](I+PC)1[IP]1=[IP](I+CP)1[IC]1.

Значения, больше, чем 0,3 обычно, указывают на хорошие поля робастности. Если система с обратной связью нестабильна, то marg = 0. Количество b (P, C) –1 является усилением сигнала от воздействий на вводе и выводе объекта к вводу и выводу контроллера.

Частота, на которой происходит граничный marg, возвратилась как скаляр. Если система с обратной связью нестабильна, то freq = NaN.

Больше о

свернуть все

Запас устойчивости и метрики разрыва

b запаса устойчивости (P, C) связан с метрикой разрыва, которая дает численное значение δ (P 1, P 2) для расстояния между двумя системами LTI (см. gapmetric).

И для разрыва и для ν - метрики разрыва, следующий устойчивый результат производительности содержит:

arcsin b (P 2, C 2) ≥ arcsin b (P 1, C 1) – arcsin δ (P 1, P 2) – arcsin δ (C 1, C 2),

Чтобы интерпретировать этот результат, предположите, что номинальный объект P 1 стабилизируется контроллером C 1 с запасом устойчивости b (P 1, C 1). Затем если P 1 встревожен к P 2, и C 1 встревожен к C 2, запас устойчивости ухудшается не больше, чем вышеупомянутой формулой.

ν - разрыв всегда меньше чем или равен разрыву, таким образом, его прогнозы с помощью вышеупомянутого результата робастности более трудны.

Советы

  • В то время как ncfmargin принимает цикл отрицательной обратной связи, команда ncfsyn разрабатывает контроллер для цикла позитивных откликов. Поэтому, чтобы вычислить поле с помощью контроллера, разработанного с ncfsyn, используйте [marg,freq] = ncfmargin(P,C,+1).

Алгоритмы

Вычисление нормированного взаимно-простого запаса устойчивости как описано в Главе 16 [1].

Ссылки

[1] Чжоу, K., Дойл, J.C., основы устойчивого управления. Лондон, Великобритания: Пирсон, 1997.

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

| | | |

Представлено до R2006a