Преобразование идентифицированной линейной модели с шумовыми каналами в модель только с измеренными каналами
mod1 = noisecnv(mod)
mod2 = noisecnv(mod,'normalize')
mod1 = noisecnv(mod) и mod2 = noisecnv(mod,'normalize') преобразование идентифицированной линейной модели с шумовыми каналами в модель только с измеряемыми каналами.
mod - любая линейная идентифицированная модель, idproc, idtf, idgrey, idpoly, или idss.
Входные шумовые каналы в mod преобразуются следующим образом: Рассмотрим модель с измеренными входными каналами u (nu каналов) и шумовыми каналами e (ny каналов) с ковариационной матрицей Λ:
= Λ = LL
где L - нижняя треугольная матрица. Обратите внимание, что mod.NoiseVariance = Λ. Модель также может быть описана с единичной дисперсией, используя нормализованный источник шума v:
(v) = I
mod1 = noisecnv(mod) преобразует модель в представление системы [G
H] с nu + ny входами и ny выходами. Все входы обрабатываются как измеренные, и mod1 не имеет шумовой модели. Прежние шумовые входные каналы имеют имена e@yname, где yname - имя соответствующего выходного документа.
mod2 = noisecnv(mod,'norm') преобразует модель в представление системы [G
HL] с nu + ny входами и ny выходами. Все входы обрабатываются как измеренные, и mod2 не имеет шумовой модели. Прежние шумовые входные каналы имеют имена v@yname, где yname - имя соответствующего выходного документа. Заметим, что коэффициент L матрицы дисперсии шума обычно является неопределенным (имеет ненулевую ковариацию). Это учитывается в описании неопределенности mod2.
Если mod - временной ряд, то есть nu = 0, mod1 - модель, описывающая передаточную функцию H с измеренными входными каналами. Аналогично, mod2 описывает передаточную функцию HL.
Обратите внимание на различие между вложенными ссылками:
mod(:,[]) дает описание характеристик шумовой модели как модели временного ряда, то есть описывает Н, а также ковариацию е. В противоположность, noisecnv(m(:,[])) или noise2meas(m) описывать только передаточную функцию H. Чтобы получить описание нормализованной передаточной функции HL, используйте noisecnv(m(:,[]),'normalize') или noise2meas('normalize').
Преобразование шумовых каналов в измеренные входы полезно для изучения свойств отдельных передаточных функций от шума к выходу. Он также полезен для преобразования идентифицированных линейных моделей в представления, которые явно не обрабатывают описания возмущений.
Идентификация модели с измеряемым компонентом (G) и нетривиальным шумовым компонентом (H). Сравните амплитуду частотной характеристики измеряемой составляющей с амплитудой спектра шумовой составляющей. Необходимо преобразовать компонент шума в измеренный с помощью noisecnv если вы хотите сравнить его поведение с действительно измеренным компонентом.
load iddata2 z2 sys1 = armax(z2,[2 2 2 1]); % model with noise component sys2 = tfest(z2,3); % model with a trivial noise component sys1 = noisecnv(sys1); sys2 = noisecnv(sys2); bodemag(sys1,sys2)
idssdata | noise2meas | tfdata | zpkdata