спектр
Спектр выходной мощности моделей временных рядов
Синтаксис
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 является моделью временных рядов, которая представляет систему:
Где, e(t) является Гауссовым белым шумом, и y(t) является наблюдаемым выходным сигналом.
spectrum строит abs(H'H), масштабируемый отклонением e(t) и шага расчета.
Если sys является моделью ввода - вывода, он представляет систему:
Где, 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 не задан.
Входные параметры
|
Модель Identified. Если Где, Если Где, |
|
Минимальная частота частотного диапазона, для которого построен спектр выходной мощности. Задайте |
|
Максимальная частота частотного диапазона, для которого построен спектр выходной мощности. Задайте |
|
Частоты, для которых построен спектр выходной мощности. Задайте |
|
Идентифицированные системы, для которых построен спектр выходной мощности. |
Выходные аргументы
|
Амплитуда спектра мощности. Если Для амплитудных значений в дБ введите |
|
Вектор частоты, для которого построен спектр выходной мощности. |
|
Предполагаемое стандартное отклонение спектра мощности. |
Примеры
Постройте шумовой спектр линейной идентифицированной модели SISO
Загрузите данные об оценке.
load iddata1 z1;
Оцените одно вход модель в пространстве состояний одно вывода.
sys = n4sid(z1,2);
Постройте шумовой спектр для модели.
spectrum(sys);
Постройте Выходной спектр модели timeseries
Загрузите данные об оценке 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),'+')
