spectrum

Спектр выходной мощности моделей временных рядов

Синтаксис

spectrum(sys)
spectrum(sys,{wmin, wmax})
spectrum(sys,w)
spectrum(sys1,...,sysN,w)
ps = spectrum(sys,w)
[ps,w] = spectrum(sys)
[ps,w,sdps] = spectrum(sys)

Описание

spectrum(sys) создает график спектра выходной мощности идентифицированной модели sys временных рядов. Частотный диапазон и число точек выбраны автоматически.

sys модель временных рядов, которая представляет систему:

y(t)=He(t)

Где, e(t) Гауссов белый шум и y(t) наблюдаемый выходной сигнал.

spectrum графики abs(H'H), масштабируемый отклонением e(t) и шаг расчета.

Если sys модель ввода - вывода, она представляет систему:

y(t)=Gu(t)+He(t)

Где, u(t) измеренный вход, e(t) Гауссов белый шум и y(t) наблюдаемый выходной сигнал.

В этом случае, spectrum строит спектр компонента воздействия He(t).

spectrum(sys,{wmin, wmax}) создает график спектра для частот в пределах от wmin к wmax.

spectrum(sys,w) создает график спектра с помощью частот, заданных в векторном w.

spectrum(sys1,...,sysN,w) создает график спектра нескольких идентифицированных моделей на одном графике. w аргумент является дополнительным.

Можно задать цвет, стиль линии и маркер для каждой модели. Например:

spectrum(sys1,'r',sys2,'y--',sys3,'gx');

ps = spectrum(sys,w) возвращает амплитуду спектра мощности sys для заданных частот, w. Никакой график не построен на экране.

[ps,w] = spectrum(sys) возвращает вектор частоты, w, для которого построен спектр выходной мощности.

[ps,w,sdps] = spectrum(sys) возвращает предполагаемые стандартные отклонения спектра мощности.

Для моделей дискретного времени с шагом расчета Ts, spectrum использует преобразование z = exp(j*w*Ts) сопоставлять модульный круг с действительной осью частоты. Спектр только построен для частот, меньших, чем частота Найквиста pi/Ts, и значение по умолчанию 1 (единица измерения времени) принято, когда Ts не задан.

Входные параметры

sys

Модель Identified.

Если sys модель временных рядов, она представляет систему:

y(t)=He(t)

Где, e(t) Гауссов белый шум и y(t) наблюдаемый выходной сигнал.

Если sys модель ввода - вывода, она представляет систему:

y(t)=Gu(t)+He(t)

Где, u(t) измеренный вход, e(t) Гауссов белый шум и y(t) наблюдаемый выходной сигнал.

wmin

Минимальная частота частотного диапазона, для которого построен спектр выходной мощности.

Задайте wmin в rad/TimeUnit, где TimeUnit sys.TimeUnit.

wmax

Максимальная частота частотного диапазона, для которого построен спектр выходной мощности.

Задайте wmax в rad/TimeUnit, где TimeUnit sys.TimeUnit.

w

Частоты, для которых построен спектр выходной мощности.

Задайте w в rad/TimeUnit, где TimeUnit sys.TimeUnit.

sys1,...,sysN

Идентифицированные системы, для которых построен спектр выходной мощности.

Выходные аргументы

ps

Амплитуда спектра мощности.

Если sys имеет Ny выходные параметры, затем ps массив размера [Ny Ny length(w)]. Где ps(:,:,k) соответствует спектру мощности для частоты в w(k).

Для амплитудных значений в дБ введите psdb = 10*log10(ps).

w

Вектор частоты, для которого построен спектр выходной мощности.

sdps

Предполагаемое стандартное отклонение спектра мощности.

Примеры

свернуть все

Загрузите данные об оценке.

load iddata1 z1;

Оцените одно вход модель в пространстве состояний одно выхода.

sys = n4sid(z1,2);

Постройте шумовой спектр для модели.

spectrum(sys);

Загрузите данные об оценке timeseries.

load iddata9 z9

Оцените четвертый порядок модель AR с помощью подхода наименьших квадратов.

sys = ar(z9,4,'ls');

Постройте выходной спектр модели.

spectrum(sys);

Создайте вход, состоящий из пяти синусоид, распространенных по целому интервалу частоты. Сравните спектр этого сигнала с тем из его квадрата. Разделение частоты (квадрат, имеющий спектральную поддержку на других частотах), показывает включенную нелинейность.

u = idinput([100 1 20],'sine',[],[],[5 10 1]);
u = iddata([],u,1,'per',100);
u2 = u.u.^2;
u2 = iddata([],u2,1,'per',100);
spectrum(etfe(u),'r*',etfe(u2),'+')

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

| | | | | |

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