[K,CL,GAM,INFO]=loopsyn(G,Gd) [K,CL,GAM,INFO]=loopsyn(G,Gd,RANGE)
loopsyn
H ∞ оптимальный метод для синтеза управления loopshaping. Это вычисляет стабилизировавшийся H ∞controller K для объекта G, чтобы сформировать график sigma
передаточной функции цикла GK, чтобы желать, чтобы цикл сформировал Gd с точностью γ = GAM
в том смысле, что, если ω 0 является перекрестной частотой на 0 дб графика sigma
Gd (j ω), то, примерно,
(1) |
(2) |
INFO
Массива структур возвращает дополнительную информацию о проекте, включая MIMO стабильная фаза min, формирующая предварительный фильтр W, имеющий форму объект Gs = GW, контроллер для имеющего форму объекта Ks = WK, а также частотный диапазон {min ω, ω макс.}, по которому достигается формирование цикла
Входной параметр | Описание |
---|---|
G | Объект LTI |
Gd | Желаемая форма цикла (модель LTI) |
RANGE | (дополнительный, |
Выходной аргумент | Описание |
---|---|
K | Контроллер LTI |
CL= G*K/(I+GK) | LTI система с обратной связью |
GAM | Формирующая цикл точность ( |
INFO | Дополнительная выходная информация |
INFO.W | Предварительный фильтр LTI W, удовлетворяющий σ (Gd) = σ (GW) для всего ω; W всегда является минимальной фазой. |
INFO.Gs | Объект, имеющий форму LTI: Gs = GW. |
INFO.Ks | Контроллер LTI для имеющего форму объекта: Ks = WK. |
INFO.range | {min ω, ω макс.} массив ячеек, содержащий аппроксимированный частотный диапазон, по которому формирование цикла могло быть точно достигнуто к с точностью |
Объект G должен быть stabilizable и обнаруживаемым, должен иметь, по крайней мере, столько же входных параметров сколько выходные параметры и должен быть полным рангом; т.е. ,
size(G,2)
≥ size(G,1)
rank(freqresp(G,w)) = size(G,1)
для некоторой частоты w.
Порядок контроллера K может быть большим. В общем, когда Gd дан как LTI SISO, затем порядок NK контроллера, K удовлетворяет
NK = NGs + NW
= Ny NGd + NRHP + NW
= Ny NGd + NRHP + NG
где
Ny обозначает количество выходных параметров объекта G.
NRHP обозначает общее количество нестабильных полюсов и нули неминимальной фазы объекта G, включая тех на контуре устойчивости и в бесконечности.
NG, NGs, NGd и NW обозначают соответствующие порядки G, Gs , Gd и W.
Снижение сложности модели может помочь уменьшать порядок K — смотрите reduce
и ncfmr
.
Используя формулу GCD Ле и Сафонова [1], loopsyn
сначала вычисляет формирование цикла стабильной минимальной фазы, придавая предварительному фильтру квадратную форму вниз W, таким образом, что имеющий форму объект Gs = GW является квадратным, и желаемая форма, Gd достигается с хорошей точностью в частотном диапазоне {min ω, ω макс.} имеющим форму объектом; т.е. .,
σ (Gd) ≈ σ (Gs) для всего ω ∊ {min ω, ω макс.}.
Затем loopsyn
использует Перчаточника-McFarlane [2], синтез управления "нормировал взаимно-простую факторную" теорию вычислить оптимальный “формирующий цикл” контроллер для имеющего форму объекта через Ks=ncfsyn(Gs), and
, возвращает K=W*Ks
.
Если объект G является непрерывным LTI времени и
G имеет D-матрицу полного ранга, и
никакие конечные нули на ω j - ось, и
{min ω, ω макс.} = [0, ∞],
затем GW теоретически достигает совершенной подгонки точности σ (Gd) = σ (GW) для всей частоты ω. В противном случае loopsyn
использует билинейный билинейный сдвиг полюса, преобразовывают [3] из формы
Gshifted=bilin(G,-1,'S_Tust',[ωmin,ωmax]),
который приводит к совершенному пригодному для преобразованного редактора Gshift
и аппроксимированной подгонки по меньшему частотному диапазону [min ω, ω макс.] для исходного непереключенного G при условии, что ω макс.>> min ω. Для лучших результатов необходимо выбрать ω макс., чтобы быть по крайней мере в 100 раз больше, чем min ω. В некоторых случаях вычисление оптимального W для Gshifted
может быть сингулярно или плохо обусловлено для области значений [min ω, ω макс.], как тогда, когда Gshifted
имеет незатухающие нули или в непрерывно-разовом случае только, Gshifted
имеет D - матрица, которая имеет неполный ранг); в таких случаях loopsyn
автоматически уменьшает частотный диапазон далее и возвращает уменьшаемую область значений [min ω, ω макс.] как массив ячеек в выводе INFO.range=
{min ω, ω макс.}
[1] Le, V.X., и М.Г. Сафонов. Рациональный матричный GCD и проект обработки на квадрат вниз компенсаторам — теория пространства состояний. Сделка IEEE Autom. Управление, AC-36 (3):384–392, март 1992.
[2] Перчаточник, К., и Д. Макфарлэйн. Устойчивая стабилизация нормированных взаимно-простых факторных описаний объекта с H ∞-bounded неуверенность. Сделка IEEE Autom. Управление, AC-34 (8):821–830, август 1992.
[3] Чанг, R.Y., и М.Г. Сафонов. H ∞ синтез с помощью билинейного переключающего полюс преобразования. AIAA J. Руководство, Управление и Динамика, 15 (5):1111–1115, сентябрь-октябрь 1992.