Оцените эмпирические передаточные функции и периодограммы
g = etfe(data)
data содержит время - или данные ввода - вывода частотного диапазона или данные timeseries:
Если data сигналы ввода - вывода временного интервала, g отношение выходного преобразования Фурье к входному преобразованию Фурье для данных.
Для апериодических данных передаточная функция оценивается на 128 равномерно распределенных частотах [1:128]/128*pi/Ts.
Для периодических данных, которые содержат целое число периодов (data.Period = integer ), ответ вычисляется на частотах k*2*pi/period для k = 0 до частоты Найквиста.
Если data сигналы ввода - вывода частотного диапазона, g отношение выхода, чтобы ввести на всех частотах, где вход является ненулевым.
Если data данные timeseries (никакие входные каналы), g периодограмма, которая является normed абсолютным квадратом преобразования Фурье данных. Соответствующая спектральная оценка нормирована, как описано в Нормализации Спектра и отличается от spectrum нормализация в продукте Signal Processing Toolbox™.
g = etfe(data,M)pi/M. Эффект M похоже на эффект M \in spaM проигнорирован для периодических данных. Используйте этот синтаксис в качестве альтернативы 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')
