Функции изменения и уточнения структуры модели

Для изменения и уточнения структуры модели служат следующие функции.

Функция fixpar

Функция fixpar фиксирует коэффициенты в структурах тета-моделей для переменных состояния и ARX, созданных с помощью функ­ций ms2th и arx2th:

thn = fixpar(tho,matrix)

thn = fixpar(tho,matrix,elements,parval)

где:

• tho — исходная тета-модель;

• аргумент matrix (строковая переменная):

o если исходная модель задана через переменные состояния, при­нимает одно из значений: 'А', 'В', 'С', 'D', 'К' или 'х0';

o если исходная модель — ARX-типа, то параметр matrix прини­мает одно из значений: 'А1', 'А2', ..., 'В0', 'В1', ... ;

• аргумент elements представляет собой матрицу из двух столбцов, содержащую индексы фиксируемых (то есть не подлежащих ка­кому-либо последующему изменению) элементов матрицы, опре­деляемой аргументом matrix, например, запись thn = fixpar(tho,'A',[1,2;3,5]); означает, что будут зафиксированы элементы А12 и A35. По умолчанию фиксируются начальные значения параметров или их оцен­ки (последние по времени);

• parval — вектор, элементы которого явно задают некоторые фиксируемые параметры (k-й элемент этого вектора задает величину k-ro элемента в строке elements).

Возвращаемая величина thn — тета-модель с частично зафиксирован­ными параметрами.

Функция sett

Функция sett явным образом задает интервал дискретизации в мо­делях:

modn = sett(mod,T)

Аргументы:

• mod — имя (идентификатор) модели в тета-формате или в частот­ном формате (см. описание функции spa);

• Т — задаваемый интервал дискретизации.

Возвращаемая величина modn — модель объекта с новым значением Т.

Функция thinit

Функция thinit устанавливает начальные (случайные) значения ко­эффициентов, подлежащих оцениванию в модели для переменных состояния (тета-формат) с последующим применением для такого оценивания функции pem. Записывается в виде:

th = thinit(th0)

th = thinit(th0,R,pars,sp)

где

• th0 — исходная тета-модель;

• аргументы R и pars задают параметры закона случайного задания коэффициентов модели по следующему правилу: коэффициент с номером к устанавливается равным величине pars(k) + e*sqrt(R(k)), где е — гауссова случайная величина с нулевым математическим ожиданием и единичной дисперсией. По умолчанию все элемен­ты вектора R — единицы, а элементы вектора pars равны номи­нальным значениям коэффициентов модели;

• sp — строковая переменная с возможными значениями 'b' (ста­бильна как модель, так и прогнозируемый по ней выход), 's' (ста­бильна только модель) и 'р' (стабилен только прогнозируемый выход); по умолчанию sp='р'. Возвращаемая величина th — тета-модель сзаданными начальными значениями коэффициентов.

Функция unfixpar

Функция unfixpar является обратной по отношению к функции fixpar, то есть делает зафиксированные коэффициенты модели свободными:

thn = unfixpar(tho,matrix)

thn = unfixpar(tho,matrix,elements)

Аргументы и возвращаемая величина идентичны рассмотренным для функции fixpar.