Оценка частотной характеристики и спектра с помощью спектрального анализа с частотно-зависимым разрешением
g = spafdr(data)
g = spafdr(data,Resol,w)
g = spafdr(data) оценивает частотный отклик «вход-выход» G (λ) и шумовой спектр Фt общей линейной модели
) + v (t)
где Γ (λ) - спектр (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, что соответствует желаемому разрешению вокруг рассматриваемой частоты. Отношение этих средних затем формируется для оценки частотной функции, и соответствующие выражения определяют оценку шумового спектра.