Оцените частотную характеристику и спектр с помощью спектрального анализа с зависимым частотой разрешением
g = spafdr(data)
g = spafdr(data,Resol,w)
g = spafdr(data)
оценивает частотную характеристику входа к выходу G (ω) и шумовой спектр Φυ общей линейной модели
где Φυ (ω) является спектром υ (t). data
содержит выходные входные данные как iddata
объект. Данные могут быть комплексные оцененный, и или время или частотный диапазон. Это может также быть idfrd
объект, содержащий данные частотной характеристики. g
idfrd
объект с оценкой на частотах ω задан вектором-строкой w
G
также включает информацию об оценке спектра Φυ (ω) на тех же частотах. Оба результата возвращены с предполагаемыми ковариациями, включенными в g
. Нормализация спектра эквивалентна описанный в spa
.
Информация о результатах оценки и используемых опциях хранится в Report
модели свойство.
Report
имеет следующие поля:
Status
— Сводные данные состояния модели, которое указывает, была ли модель создана конструкцией или получена оценкой.
Method
— Команда оценки используется.
WindowSize
— Разрешение частоты.
DataUsed
— Атрибуты данных используются для оценки. Структура со следующими полями:
Name
— Имя набора данных.
Type
тип данных.
Length
— Количество выборок данных.
Ts
Размер шага.
InterSample
— Введите междемонстрационное поведение.
InputOffset
— Возместите удаленный из входных данных временного интервала во время оценки.
OutputOffset
— Возместите удаленный из выходных данных временного интервала во время оценки.
g = spafdr(data,Resol,w)
задает разрешение частоты и частоты.
Переменная w
частоты или задан как вектор-строка из частот в rad/
TimeUnit
, где TimeUnit
относится к TimeUnit
свойство данных, или как массив ячеек {wmin,wmax}
. В последнем случае охваченные частоты будут 50 логарифмически расположенными с интервалами точками от wmin
к wmax
. Можно изменить число точек в NP
путем ввода {wmin,wmax,NP}
.
Исключение w
или вводя его, когда пустая матрица дает значение по умолчанию, которое является 100 логарифмически расположенными с интервалами частотами между наименьшей и самой большой частотой в данных. Для данных временного интервала область значений по умолчанию идет от к , где Ts является шагом расчета данных, и N является количеством точек данных.
Аргумент Resol
задает разрешение частоты оценок. Разрешение (измеренный в rad/TimeUnit
) размер самой маленькой детали в функции частоты и спектре, который разрешен оценкой. Разрешение является компромиссом между получением оценок с прекрасными, надежными деталями и страданием от побочных, случайных эффектов: Чем более прекрасный разрешение, тем выше отклонение в оценке. Resol
может быть введен как скаляр (измеренный в rad/TimeUnit
), который задает разрешение на целом интервале частоты. Это может также быть введено как вектор-строка из той же длины как w
. Затем Resol(k)
локальное, зависимое частотой разрешение вокруг частоты w(k)
.
Значение по умолчанию Resol
, полученный путем исключения его или ввода его как пустой матрицы, Resol(k) = 2(w(k+1)-w(k))
, увеличенный, так, чтобы обоснованная оценка была гарантирована. Во всех случаях разрешение возвращено в переменной g.Report.WindowSize
.
Если данные даны во временном интервале, они сначала преобразованы в частотный диапазон. Затем средние значения Y(w)Conj(U(w))
и U(w)Conj(U(w))
формируются по частотным диапазонам w
, соответствие нужному разрешению вокруг рассматриваемой частоты. Отношение этих средних значений затем формируется для функциональной частотой оценки, и соответствующие выражения задают шумовую оценку спектра.