В этом примере показано, как оценить модель передаточной функции, когда структура ожидаемой модели известна, и примените ограничения к коэффициентам знаменателя и числителю.
Загрузите данные временного интервала.
load iddata1 z1; z1.y = cumsum(z1.y);
cumsum
интегрирует выходные данные z1
. Предполагаемая передаточная функция должна поэтому содержать интегратор.
Создайте модель передаточной функции с ожидаемой структурой.
init_sys = idtf([100 1500],[1 10 10 0]);
int_sys
idtf
модель с тремя полюсами и одним нулем. Коэффициент знаменателя для s^0
термин является нулем, который указывает на тот int_sys
содержит интегратор.
Задайте ограничения на числитель и коэффициенты знаменателя модели передаточной функции. Для этого сконфигурируйте поля в Structure
свойство:
init_sys.Structure.Numerator.Minimum = eps; init_sys.Structure.Denominator.Minimum = eps; init_sys.Structure.Denominator.Free(end) = false;
Ограничения указывают, что числитель и коэффициенты знаменателя являются неотрицательными. Кроме того, последний элемент коэффициентов знаменателя (сопоставленный с s^0
назовите), не допускающий оценку параметр. Это ограничение обеспечивает один из предполагаемых полюсов, чтобы быть в s = 0
.
Создайте набор опции оценки, который задает использование метода поиска Levenberg–Marquardt.
opt = tfestOptions('SearchMethod','lm');
Оцените передаточную функцию для z1
использование init_sys
и опция оценки установлена.
sys = tfest(z1,init_sys,opt);
tfest
использует коэффициенты init_sys
инициализировать оценку sys
. Кроме того, оценка ограничивается ограничениями, которые вы задаете в Structure
свойство init_sys
. Получившийся idtf
модель sys
содержит значения параметров, которые следуют из оценки.