Оцените частотную характеристику и спектр с помощью спектрального анализа с частотно-зависимым разрешением
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 определяет частотное разрешение оценок. Разрешение (измеренное в рад/ TimeUnit) - размер наименьшей детали в частотной функции и спектре, который разрешается оценкой. Разрешение является компромиссом между получением оценок с тонкими, надежными деталями и страданием от ложных, случайных эффектов: Чем больше разрешение, тем выше отклонение в оценке. Resol может быть введен как скаляр (измеренный в рад/ 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, соответствующее желаемому разрешению вокруг рассматриваемой частоты. Отношение этих средних значений затем формируется для оценки частотной функции, и соответствующие выражения определяют оценку шумового спектра.
