Оцените частотную характеристику и спектр с помощью спектрального анализа с зависимым частотой разрешением
g = spafdr(data)
g = spafdr(data,Resol,w)
g = spafdr(data) оценивает частотную характеристику входа к выводу G (ω) и шумовой спектр Φυ общей линейной модели
где Φυ (ω) является спектром υ (t). data содержит выходные входные данные как объект iddata. Данные могут быть комплексные оцененный, и или время или частотный диапазон. Это может также быть объект idfrd, содержащий данные частотной характеристики. g является объектом idfrd с оценкой на частотах ω задан вектором - строкой w. g также включает информацию об оценке спектра Φυ (ω) на тех же частотах. Оба результата возвращены с предполагаемыми ковариациями, включенными в g. Нормализация спектра эквивалентна описанный в spa.
Информация о результатах оценки и используемых опциях хранится в свойстве Report модели. Report имеет следующие поля:
Состояние Сводные данные состояния модели, которое указывает, была ли модель создана конструкцией или получена оценкой.
Method — Команда оценки используется.
WindowSize — Разрешение частоты.
DataUsed — Атрибуты данных используются для оценки. Структура со следующими полями:
Имя Имя набора данных.
Ввод Тип данных.
Длина Количество выборок данных.
Ts — 'SampleTime' .
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, соответствуя нужному разрешению вокруг рассматриваемой частоты. Отношение этих средних значений затем формируется для функциональной частотой оценки, и соответствующие выражения задают шумовую оценку спектра.
