spafdr

Оцените частотную характеристику и спектр с помощью спектрального анализа с частотно-зависимым разрешением

Синтаксис

g = spafdr(data)
g = spafdr(data,Resol,w)

Описание

g = spafdr(data) оценивает частотную характеристику входа к выходу G (ω) и шумовой спектр Φυ общей линейной модели

y(t)=G(q)u(t)+v(t)

где Φυ (ω) является спектром υ (t). data содержит выходные входные данные как iddata объект. Данные могут быть комплексные оцененный, и или время или частотный диапазон. Это может также быть idfrd объект, содержащий данные частотной характеристики. g idfrd объект с оценкой G(eiω) на частотах ω задан вектором-строкой wG также включает информацию об оценке спектра Φυ (ω) на тех же частотах. Оба результата возвращены с предполагаемыми ковариациями, включенными в 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 логарифмически расположенными с интервалами частотами между наименьшей и самой большой частотой в данных. Для данных временного интервала область значений по умолчанию идет от 2πNTs к πTs, где 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, соответствие нужному разрешению вокруг рассматриваемой частоты. Отношение этих средних значений затем формируется для функциональной частотой оценки, и соответствующие выражения задают шумовую оценку спектра.

Смотрите также

| | | | | |

Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте