Что такое полиномиальные модели?

Полиномиальная структура модели

Полиномиальная модель использует обобщенное понятие передаточных функций, чтобы выразить отношение между входом, u (t), выход y (t), и шумовым e (t) с помощью уравнения:

A(q)y(t)=i=1nuBi(q)Fi(q)ui(tnki)+C(q)D(q)e(t)

Переменные A, B, C, D и F являются полиномами, выраженными в операторе сдвига времени q^-1. ui является входом ith, ню является общим количеством входных параметров, и nki является входной задержкой ith, которая характеризует транспортную задержку. Отклонение белого шума e (t) принято, чтобы быть λ. Для получения дополнительной информации об операторе сдвига времени, смотрите Понимание Оператора Сдвига времени q.

На практике не все полиномы одновременно активны. Часто, более простые формы, такие как ARX, ARMAX, Ошибка на выходе и Поле-Jenkins используются. У вас также есть опция представления интегратора в источнике шума так, чтобы общая модель приняла форму:

A(q)y(t)=i=1nuBi(q)Fi(q)ui(tnki)+C(q)D(q)11q1e(t)

Для получения дополнительной информации смотрите Различные Настройки Полиномиальных Моделей.

Можно оценить полиномиальные модели, использующие данные о частотном диапазоне или время.

Для оценки необходимо задать порядок модели как набор целых чисел, которые представляют количество коэффициентов для каждого полинома, который вы включаете в свою выбранную структуру — na для A, nb для B, nc для C, без обозначения даты для D и nf для F. Необходимо также задать количество выборок nk соответствие входной задержке — потеря времени — данный количеством выборок, прежде чем выход ответит на вход.

Количество коэффициентов в полиномах знаменателя равно количеству полюсов, и количество коэффициентов в полиномах числителя равно количеству нулей плюс 1. Когда движущие силы от u (t) к y (t) содержат задержку nk выборок, затем первые nk коэффициенты B являются нулем.

Для получения дополнительной информации о семействе моделей передаточной функции, смотрите соответствующий раздел в System Identification: Теория для Пользователя, Второго Выпуска, Lennart Ljung, PTR Prentice Hall, 1999.

Понимание Оператора Сдвига времени q

Общее полиномиальное уравнение записано в терминах оператора сдвига времени q–1. Чтобы изучить этот оператор сдвига времени, рассмотрите следующее разностное уравнение дискретного времени:

y(t)+a1y(tT)+a2y(t2T)=       b1u(tT)+b2u(t2T)

где y (t) является выход, u (t) является входом, и T является шагом расчета. q-1 является оператором сдвига времени, который сжато представляет такое использование разностных уравнений q1u(t)=u(tT):

y(t)+a1q1y(t)+a2q2y(t)=         b1q1u(t)+b2q2u(t)илиA(q)y(t)=B(q)u(t)

В этом случае, A(q)=1+a1q1+a2q2 и B(q)=b1q1+b2q2.

Примечание

Это q описание абсолютно эквивалентно форме Z-преобразования: q соответствует z.

Различные настройки полиномиальных моделей

Эти структуры модели являются подмножествами следующего общего полиномиального уравнения:

A(q)y(t)=i=1nuBi(q)Fi(q)ui(tnki)+C(q)D(q)e(t)

Структуры модели отличаются тем, сколько из этих полиномов включено в структуру. Таким образом различные структуры модели обеспечивают различные уровни гибкости для моделирования динамики и шумовых характеристик.

Следующая таблица обобщает общие линейные полиномиальные структуры модели, поддержанные продуктом System Identification Toolbox™. Если вы имеете определенную структуру в виду для вашего приложения, можно решить, имеют ли динамика и шум общие или различные полюса. (q) соответствует полюсам, которые характерны для динамической модели и шумовой модели. Используя общие полюса для динамики и шума полезно, когда воздействия вводят систему во входе. F i определяет полюса, уникальные для системной динамики, и D определяет полюса, уникальные для воздействий.

Структура моделиУравнениеОписание
ARX

A(q)y(t)=i=1nuBi(q)ui(tnki)+e(t)

Шумовая модель 1A и шум связывается с моделью динамики. ARX не позволяет вам шум модели и динамика независимо. Оцените, что модель ARX получает простую модель в хороших отношениях сигнал-шум.

ARIX

Ay=Bu+11q1e

Расширяет структуру ARX включением интегратора в источнике шума, e (t). Это полезно в случаях, где воздействие не является стационарным.
ARMAX

A(q)y(t)=i=1nuBi(q)ui(tnki)+C(q)e(t)

Расширяет структуру ARX путем обеспечения большей гибкости для моделирования шумового использования параметров C (скользящее среднее значение белого шума). Используйте ARMAX, когда воздействия доминирования войдут во входе. Такие воздействия называются воздействиями загрузки.
ARIMAX

Ay=Bu+C11q1e

Расширяет структуру ARMAX включением интегратора в источнике шума, e (t). Это полезно в случаях, где воздействие не является стационарным.
Поле-Jenkins (BJ)

y(t)=i=1nuBi(q)Fi(q)ui(tnki)+C(q)D(q)e(t)

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

Используйте модели BJ, когда шум не входит во входе, но первичен воздействие измерения, Эта структура обеспечивает дополнительную гибкость для моделирования шума.

Ошибка на выходе (OE)

y(t)=i=1nuBi(q)Fi(q)ui(tnki)+e(t)

Используйте, когда это необходимо, чтобы параметризовать динамику, но не хотеть оценивать шумовую модель.

Примечание

В этом случае, шумовые модели H=1 в общем уравнении и белом источнике шума e (t) влияет только на выход.

Полиномиальные модели могут содержать одни или несколько выходных параметров и нуля или больше входных параметров.

Поддержки приложений System Identification прямая оценка ARX, ARMAX, моделей OE и BJ. Можно добавить шумовой интегратор в ARX, ARMAX и формы BJ. Однако можно использовать polyest оценить все пять полиномов или любое подмножество полиномов в общем уравнении. Для получения дополнительной информации о работе с pem, смотрите Используя полиоценку, чтобы Оценить Полиномиальные Модели.

Представление непрерывного времени полиномиальных моделей

В непрерывное время общее уравнение частотного диапазона записано в терминах переменной Преобразования Лапласа s, который соответствует операции дифференцирования:

A(s)Y(s)=B(s)F(s)U(s)+C(s)D(s)E(s)

В случае непрерывного времени базовая модель временного интервала является дифференциальным уравнением, и целые числа порядка модели представляют количество предполагаемого числителя и коэффициентов знаменателя. Например, na=3 и nb=2 соответствуют следующей модели:

A(s)=s4+a1s3+a2s2+a3B(s)=b1s+b2

Можно только оценить модели полинома непрерывного времени непосредственно с помощью данных частотного диапазона непрерывного времени. В этом случае необходимо установить Ts свойство данных к 0, чтобы указать, что вы имеете данные частотного диапазона непрерывного времени и используете oe команда, чтобы оценить модель полинома Ошибки на выходе. Модели непрерывного времени других структур, такие как ARMAX или BJ не могут быть оценены. Можно получить те формы только прямой конструкцией (использующий idpoly), преобразование из других типов модели, или путем преобразования модели дискретного времени в непрерывное время (d2c). Обратите внимание на то, что форма OE представляет передаточную функцию, выраженную как отношение числителя (B) и знаменатель (F) полиномы. Поскольку такие формы рассматривают использование моделей передаточной функции, представленных idtf модели. Можно оценить модели передаточной функции с помощью и времени и данных о частотном диапазоне. В дополнение к полиному числителя и полиному знаменателя можно также оценить транспортные задержки. Смотрите idtf и tfest для получения дополнительной информации.

Мультивыведите полиномиальные модели

Для модели полинома MIMO с ny выходные параметры и входные параметры nu, отношение между вводами и выводами для l th выход может быть записано как:

j=1nyAlj(q)yj(t)=i=1nuBli(q)Fli(q)ui(tnki)+Cl(q)Dl(q)el(t)

Массив полинома A (Aij; i =1:ny, j =1:ny), хранятся в A свойство idpoly объект. Диагональные полиномы (Aii; i =1:ny), monic, то есть, ведущие коэффициенты являются тем. Недиагональные полиномы (Aij; ij ), содержат задержку по крайней мере одной выборки, то есть, они запускают с нуля. Для получения дополнительной информации о порядках мультивыходных моделей смотрите Полиномиальные Размеры и Порядки Мультивыходных Моделей Полинома.

Можно создать мультивыходные модели полинома при помощи idpoly команда или оценивает их использующий ar, arx, bj, oe, armax, и polyest. В приложении можно оценить такие модели путем выбора набора мультивыходных данных и устанавливания порядков соответственно в диалоговом окне Polynomial Models.

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

| | | | | |

Связанные примеры

Больше о