Вычислите рабочую точку для нелинейной модели ARX
[
возвращает значения состояния рабочей точки, X
,U
] =
findop(sys
,'steady',InputLevel
,OutputLevel
)X
, и входные значения, U
, дляidnlarx
модель, sys
, использование установившихся технических требований ввода и вывода.
Оцените нелинейную модель ARX.
load iddata6;
M = nlarx(z6,[4 3 1]);
Найдите установившуюся рабочую точку, где уровень на входе фиксируется к 1
и выход неизвестен.
[X,U] = findop(M,'steady',1,NaN);
Оцените нелинейную модель ARX.
load iddata7;
M = nlarx(z7,[4 3*ones(1,2) 2*ones(1,2)]);
Создайте объект спецификации рабочей точки по умолчанию.
spec = operspec(M);
Установите значения для входных сигналов.
spec.Input.Value(1) = -1; spec.Input.Value(2) = 1;
Установите максимальные и минимальные значения для выходного сигнала.
spec.Output.Max = 10; spec.Output.Min = -10;
Найдите установившуюся рабочую точку с помощью данных технических требований.
[X,U] = findop(M,spec);
Оцените нелинейную модель ARX.
load iddata6;
M = nlarx(z6,[4 3 2]);
Создайте findopOptions
по умолчанию опция установлена.
opt = findopOptions(M);
Измените набор опции, чтобы задать метод поиска градиента наискорейшего спуска, имеющий до 50 итераций.
opt.SearchMethod = 'grad';
opt.SearchOptions.MaxIterations = 50;
Найдите установившуюся рабочую точку с помощью заданных опций.
[X,U] = findop(M,'steady',1,1,opt);
Оцените нелинейную модель ARX.
load iddata7;
M = nlarx(z7,[4 3*ones(1,2) 2*ones(1,2)]);
Найдите, что установившаяся рабочая точка, где введено 1 установлена в 1
и вход 2 неограничен. Исходным предположением для выходного значения является 2
.
[X,U,R] = findop(M,'steady',[1 NaN],2);
Отобразите сводный отчет.
disp(R);
SearchMethod: 'auto' WhyStop: 'Near (local) minimum, (norm(g) < tol).' Iterations: 11 FinalCost: 0 FirstOrderOptimality: 0 SignalLevels: [1x1 struct]
Загрузите данные об оценке и оцените нелинейную модель ARX.
load twotankdata;
z = iddata(y,u,1);
M = nlarx(z,[4 3 1]);
Найдите снимок состояния симуляции после 10 секунд, приняв начальные состояния нуля.
[X,U] = findop(M,'snapshot',10,z);
Загрузите данные об оценке и оцените нелинейную модель ARX.
load twotankdata;
z = iddata(y,u,1);
M = nlarx(z,[4 3 1]);
Создайте вектор начального состояния. Первые четыре состояния соответствуют задержанным выходным значениям, и итоговые три состояния соответствуют задержанным входным параметрам.
X0 = [2;2;2;2;5;5;5];
Найдите снимок состояния симуляции после 10 секунд с помощью заданных начальных состояний.
[X,U] = findop(M,'snapshot',10,z,X0);
sys
— Нелинейная модель ARXidnlarx
объектНелинейная модель ARX в виде idnlarx
объект.
InputLevel
— Установившийся уровень на входеУстановившийся уровень на входе для вычисления рабочей точки в виде вектора. Длина InputLevel
должен равняться количеству входных параметров, заданных в sys
.
Алгоритм оптимизации принимает что конечные значения в InputLevel
зафиксированные входные значения. Используйте NaN
задавать неизвестные входные сигналы с исходными предположениями 0
. Минимальные и максимальные границы для всех входных параметров имеют значения по умолчанию -Inf
и +Inf
соответственно.
OutputLevel
— Установившийся уровень на выходеУстановившийся уровень на выходе для вычисления рабочей точки в виде вектора. Длина OutputLevel
должен равняться количеству выходных параметров, заданных в sys
.
Значения в OutputLevel
укажите на исходные предположения для алгоритма оптимизации. Используйте NaN
задавать неизвестные выходные сигналы с исходными предположениями 0
. Минимальные и максимальные границы для всех выходных параметров имеют значения по умолчанию -Inf
и +Inf
соответственно.
spec
— Технические требования рабочей точкиoperspec
объектТехнические требования рабочей точки, такие как минимальные и максимальные ограничения ввода/вывода и известные входные параметры в видеoperspec
объект.
T
— Время снимка состояния рабочей точкиВремя снимка состояния рабочей точки в виде положительной скалярной величины. Значение T
должен быть в области значений [T0, N *Ts], где N является количеством входных выборок, Ts является шагом расчета, и T0 является входным временем начала (Uin.Tstart
).
Uin
— Создайте снимки вход симуляцииiddata
возразите | матрицаСоздайте снимки вход симуляции в виде одного из следующего:
Временной интервал iddata
объект с шагом расчета и входным размером, который совпадает с sys
.
Матрица со столькими же столбцов сколько там вводится каналы. Если матрица имеет строки N, входные данные принят, чтобы соответствовать временному вектору (1:N)*sys.Ts
.
X0
— Начальные состоянияНачальные состояния симуляции в виде вектор-столбца с размером равняются количеству состояний в sys
x0
обеспечивает начальные условия в то время соответствие первой входной выборке (Uin.Start
, если Uin
iddata
объект или sys.Ts
если Uin
двойная матрица).
Для получения дополнительной информации о состояниях idnlarx
модель, см. Определение idnlarx состояний.
Options
— Параметры поиска рабочей точкиfindopOptions
опция установленаПараметры поиска рабочей точки в виде a findopOptions
опция установлена.
X
— Значения состояния рабочей точкиЗначения состояния рабочей точки, возвращенные как вектор-столбец длины, равняются количеству состояний модели.
U
— Входные значения рабочей точкиВходные значения рабочей точки, возвращенные как вектор-столбец длины, равняются количеству входных параметров.
Report
— Сводные данные результата поискаСводный отчет результата поиска, возвращенный как структура со следующими полями:
Поле | Описание |
---|---|
SearchMethod | Метод поиска используется для итеративной оценки параметра. Смотрите SearchMethod \in findopOptions для получения дополнительной информации. |
WhyStop | Условие завершения алгоритма поиска. |
Iterations | Количество итераций оценки выполняется. |
FinalCost | Окончательное значение целевой функции минимизации (сумма квадратичных невязок). |
FirstOrderOptimality | - норма поискового вектора градиента, когда алгоритм поиска останавливается. |
SignalLevels | Структура, содержащая поля Input и Output , которые являются уровнями сигнала ввода и вывода рабочей точки соответственно. |
findop
вычисляет рабочую точку из установившихся технических требований рабочей точки или в снимке состояния симуляции.
Чтобы вычислить установившуюся рабочую точку, вызвать findop
использование любого из следующих синтаксисов:
[X,U] = findop(sys,'steady',InputLevel,OutputLevel) [X,U] = findop(sys,spec)
Вычислить состояния, X
, и вход, U
, из установившейся рабочей точки, findop
минимизирует норму ошибки e (t) = y (t)-f (x (t), u (t)), где:
f является средством оценки нелинейности.
u (t) является входом.
x (t) является состоянием модели.
y (t) является выход модели.
Можно задать алгоритм поиска и параметры поиска с помощью findopOptions
опция установлена.
Алгоритм использует следующие независимые переменные для минимизации:
Неизвестные (незаданные) уровни входного сигнала
Уровни выходного сигнала
Поскольку idnlarx
состояния модели являются задержанными выборками переменных ввода и вывода, значения состояния являются постоянными значениями соответствующих установившихся вводов и выводов. Для получения дополнительной информации об определении нелинейных состояний модели ARX, см. Определение idnlarx состояний.
Когда вы используете синтаксис [X,U] = findop(sys,'snapshot',T,Uin,X0)
, алгоритм симулирует выход модели до времени снимка состояния, T
. Во время снимка состояния алгоритм передает выборки ввода и вывода data2state
команда, чтобы сопоставить эти значения с вектором текущего состояния.
Примечание
Для основанных на снимке состояния расчетов, findop
не выполняет числовую оптимизацию.
Поддержка параллельных вычислений доступна для оценки с помощью lsqnonlin
метод поиска (требует Optimization Toolbox™). Чтобы включить параллельные вычисления, использовать findopOptions
, установите SearchMethod
к 'lsqnonlin'
, и набор SearchOptions.Advanced.UseParallel
к true
.
Например:
opt = findopOptions(idnlarx);
opt.SearchMethod = 'lsqnonlin';
opt.SearchOptions.Advanced.UseParallel = true;
idnlarx
| findopOptions
| idnlarx/operspec
| data2state
| sim
| idnlhw/findop
У вас есть модифицированная версия этого примера. Вы хотите открыть этот пример со своими редактированиями?
1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.
2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.
3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.
4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.
5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.