Оценка эмпирических передаточных функций и периодограмм
g = etfe(data)
) + v (t)
data содержит входные-выходные данные временной или частотной области или данные временного ряда:
Если data - входные-выходные сигналы временной области, g - отношение выходного преобразования Фурье к входному преобразованию Фурье для данных.
Для непериодических данных передаточная функция оценивается в 128 равноотстоящих частот. [1:128]/128*pi/Ts.
Для периодических данных, содержащих целое количество периодов (data.Period = integer ), ответ вычисляется на частотах k*2*pi/period для k = 0 вплоть до частоты Найквиста.
Если data - входно-выходные сигналы частотной области, g - отношение выхода к входу на всех частотах, где вход ненулевой.
Если data - данные временного ряда (входные каналы отсутствуют), g - периодограмма данных, т.е. нормированный абсолютный квадрат преобразования Фурье. Соответствующая спектральная оценка нормализуется, как описано в разделе Нормализация спектра, и отличается от spectrum нормализация в продукте Toolbox™ обработки сигналов.
g = etfe(data,M)pi/M. Влияние M аналогичен эффекту M в spa. M игнорируется для периодических данных. Использовать этот синтаксис в качестве альтернативы spa для узкополосных спектров и систем, которые требуют больших значений M.
g = etfe(data,M,N)
Для непериодических данных временной области N задает частотную сетку [1:N]/N*pi/Ts rad/TimeUnit. Если не указано, N составляет 128.
Для периодических данных временной области, N игнорируется.
Для данных в частотной области: N является fmin:delta_f:fmax, где [fmin fmax] - диапазон частот в data, и delta_f является (fmax-fmin)/(N-1) rad/TimeUnit. Если не указано, ответ вычисляется на частотах, содержащихся в данных, где вход ненулевой.
Данные оценки нагрузки.
load iddata1 z1;
Оценка эмпирической передаточной функции и сглаженная спектральная оценка.
ge = etfe(z1); gs = spa(z1);
Сравните две модели на графике Боде.
bode(ge,gs)
Создайте периодический вход, смоделируйте систему и сравните частотную характеристику расчетной модели с исходной системой в возбужденных частотных точках.
Формирование периодического входного сигнала и выходного сигнала с помощью моделирования.
m = idpoly([1 -1.5 0.7],[0 1 0.5]);
u = iddata([],idinput([50,1,10],'sine'));
u.Period = 50;
y = sim(m,u);Оцените эмпирическую передаточную функцию.
me = etfe([y u]);
Сравните эмпирическую передаточную функцию с исходной моделью.
bode(me,'b*',m,'r')
Выполните операцию сглаживания для необработанных спектральных оценок с помощью окна Хэмминга и сравните отклики.
Загрузить данные.
load iddata1Оцените эмпирические передаточные функции с операцией сглаживания и без нее.
ge1 = etfe(z1); ge2 = etfe(z1,32);
Сравните модели на графике Боде.
ge2 является более гладким, чем ge1 из-за эффекта операции сглаживания.
bode(ge1,ge2)
Оцените эмпирические передаточные функции с низко- и высокочастотными расстояниями и сравните отклики.
Загрузить данные.
load iddata9Оценка эмпирических передаточных функций с низкочастотными и высокочастотными расстояниями.
ge1 = etfe(z9,[],32); ge2 = etfe(z9,[],512);
Постройте график спектра выходной мощности двух моделей.
spectrum(ge1,'b.-',ge2,'g')
