В этом примере показано, как оценить модель передаточной функции, когда структура ожидаемой модели известна, и примените ограничения к коэффициентам знаменателя и числителю.
Загрузите данные временного интервала.
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 содержит значения параметров, которые следуют из оценки.