Полиномиальные размеры и порядки Мультивыхода Полинома моделей

Для модели с Ny (Ny > 1) выходы и Nu входы, полиномы A, B, C, D, и F определены как массивы ячеек векторов - строк. Каждая запись в массиве ячеек содержит коэффициенты конкретного полинома, который связывает входные, выходные и шумовые значения. Порядки являются матрицами целых чисел, используемых в качестве входных параметров в команды оценки.

ПолиномРазмерностьОтношение, описанноеПорядки
ANy -by - Ny массив векторов-строкA{i,j} содержит коэффициенты зависимости между выходным yi и выходным yjna: Ny -by - Ny матрицу таким образом, чтобы каждая запись содержала степень соответствующего A полинома.
BNy -by - Nu массив векторов-строкB{i,j} содержат коэффициенты отношений между выходными yi и входными uj

nk: Ny -by - Nu матрицу таким образом, чтобы каждая запись содержала количество начальных фиксированных нулей соответствующего B полинома (входная задержка).

nb: Ny -by - Nu матрица, такая nb(i,j) = length(B{i,j})- nk(i,j).

C,DNy массив -by-1 из векторов-строкC{i} и D{i} содержат коэффициенты отношений между выходным yi и шумом ei

nc и nd Ny матрицы -by-1, такие что каждая запись содержит степень соответствующего C и D полинома, соответственно.

FNy -by - Nu массив векторов-строкF{i,j} содержит коэффициенты отношений между выходными yi и входными ujnf: Ny -by - Nu матрицу таким образом, чтобы каждая запись содержала степень соответствующего F полинома.

Для получения дополнительной информации см. idpoly.

Для примера рассмотрим набор уравнений ARMAX для модели 2 выхода, 1 входа:

y1(t) + 05 y1(t-1) + 0.9 y2(t-1) + 0.1 y2(t-2) = u (t ) + 5 u (t-1 ) + 2 u (t-2 ) + e1(t) + 0.01 e1(t-1) y2(t) + 00.5 y2(t-1)  + 03 y2(t-2)  =  10 u (t-2) + e2(t) + 0.1 e2(t-1) + 0.02 e2(t-2)

y1 и y2 представляют два выхода, а u представляет входную переменную. e1 и e2 представляют нарушениям порядка белого шума на выходах, y1 и y2, соответственно. Чтобы представить эти уравнения как полином формы ARMAX, используяidpoly, сконфигурируйте A, B и C полиномы следующим образом:

A = cell(2,2);
A{1,1} = [1 0.5];
A{1,2} = [0 0.9 0.1];
A{2,1} = [0];
A{2,2} = [1 0.05 0.3];

B = cell(2,1);
B{1,1} = [1 5 2];
B{2,1}  = [0 0 10];

C = cell(2,1);
C{1} = [1 0.01];
C{2} = [1 0.1 0.02];

model = idpoly(A,B,C)
model =
Discrete-time ARMAX model:                                                      
  Model for output number 1: A(z)y_1(t) = - A_i(z)y_i(t) + B(z)u(t) + C(z)e_1(t)
    A(z) = 1 + 0.5 z^-1                                                         
                                                                                
    A_2(z) = 0.9 z^-1 + 0.1 z^-2                                                
                                                                                
    B(z) = 1 + 5 z^-1 + 2 z^-2                                                  
                                                                                
    C(z) = 1 + 0.01 z^-1                                                        
                                                                                
  Model for output number 2: A(z)y_2(t) = B(z)u(t) + C(z)e_2(t)
    A(z) = 1 + 0.05 z^-1 + 0.3 z^-2                            
                                                               
    B(z) = 10 z^-2                                             
                                                               
    C(z) = 1 + 0.1 z^-1 + 0.02 z^-2                            
                                                               
Sample time: unspecified
  
Parameterization:
   Polynomial orders:   na=[1 2;0 2]   nb=[3;1]   nc=[1;2]
   nk=[0;2]
   Number of free coefficients: 12
   Use "polydata", "getpvec", "getcov" for parameters and their uncertainties.

Status:                                                         
Created by direct construction or transformation. Not estimated.

model является моделью ARMAX в дискретном времени с неопределенным шагом расчета. При оценке таких моделей необходимо задать порядки этих полиномов в качестве входных параметров.

В приложении Системы идентификации можно ввести матрицы непосредственно в поле Orders.

В командной строке задайте переменные, которые хранят матрицы порядка модели и задают эти переменные в команде model-estimation.

Совет

Чтобы упростить ввод порядков больших матриц в приложении Система Идентификации, задайте переменную NN=[NA NB NK] в командной строке. Можно задать эту переменную в поле Orders.

См. также

| | | | | |

Похожие примеры

Подробнее о