Преобразование идентифицированной линейной модели с шумовыми каналами для моделирования только с измеренными каналами
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 каналы) с ковариационными матричными Λ:
где L - нижняя треугольная матрица. Обратите внимание, что mod.NoiseVariance = Λ. Модель также может быть описана с единичным отклонением, с помощью нормализованного источника шума v:
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(:,[]) приводит описание характеристик модели шума как модели timeseries, то есть описывает H, а также ковариацию e. Напротив, 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