Variable-step solvers варьируется размер шага во время симуляции, уменьшая размер шага, чтобы увеличить точность, когда состояния модели изменяются быстро и увеличивают размер шага, чтобы не делать ненужные шаги, когда состояния модели изменяются медленно. Вычисление размера шага добавляет к вычислительным издержкам на каждом шаге, но может сократить общее количество шагов, и следовательно время симуляции, требуемое обеспечить заданный уровень точности для моделей с быстрым изменением или кусочно-непрерывными состояниями.
Когда вы устанавливаете управление Типом панели настройки Решателя к Variable-step, управление Решателем позволяет вам выбирать один из решателей переменного шага. Как с решателями фиксированного шага, набор решателей переменного шага включает дискретный решатель и набор непрерывных решателей. Однако различающийся решатели фиксированного шага, размер шага варьируется динамически на основе локальной ошибки.
Выбор между двумя типами решателей переменного шага зависит от того, ли блоки в ваших состояниях модели define и, если так, тип состояний, которые они задают. Если ваша модель не задает состояний или задает только дискретные состояния, выберите дискретный решатель. Если модель имеет непрерывные состояния, непрерывные решатели используют численное интегрирование, чтобы вычислить значения непрерывных состояний на следующем временном шаге.
Если модель не имеет никаких состояний или только дискретных состояний, Simulink® использует дискретный решатель, чтобы симулировать модель, даже если вы задаете непрерывный решатель.
Решатели переменного шага динамически варьируются размер шага во время симуляции. Каждый из этих решателей увеличивает или уменьшает размер шага с помощью своего локального контроля ошибок, чтобы достигнуть допусков, которые вы задаете. Вычисление размера шага на каждом временном шаге добавляет к вычислительным издержкам. Однако это может сократить общее количество шагов и время симуляции, требуемое обеспечить заданный уровень точности.
Можно далее категоризировать непрерывные решатели с переменным шагом как один шаг или многоступенчатый, одноразрядный или переменный порядок, и явный или неявный. Смотрите Один Шаг По сравнению с Многоступенчатыми Непрерывными Решателями для получения дополнительной информации.
Переменный шаг явные решатели спроектирован для нежестких проблем. Simulink обеспечивает четыре таких решателя:
ode45
ode23
ode113
odeN
| Решатель ОДУ | Метод с одним шагом | Многоступенчатый метод | Порядок точности | Метод |
|---|---|---|---|---|
ode45 | X | Носитель | Рунге-Кутта, Dormand-принц (4,5) пара | |
ode23 | X | Низко | Рунге-Кутта (2,3) пара Bogacki & Shampine | |
ode113 | X | Переменная, низко к высоко | Реализация PECE Адамса-Бэшфорта-Маултона | |
odeN | X | Смотрите порядок точности на Фиксированном Шаге непрерывные явные решатели | Смотрите метод интегрирования на Фиксированном Шаге непрерывные явные решатели |
Выберите odeN решатель, когда скорость симуляции важна, например, когда
Модель содержит много сброса решателя и/или нулевых пересечений
Solver Profiler не обнаруживает не пройдено шагов при профилировании модели
Если ваша проблема жестка, попытайтесь использовать один из переменного шага неявные решатели:
ode15s
ode23s
ode23t
ode23tb
| Решатель ОДУ | Метод с одним шагом | Многоступенчатый метод | Порядок точности | Метод сброса решателя | 'MaxOrder' | Метод |
|---|---|---|---|---|---|---|
ode15s | X | Переменный, низкий и средний | X | X | Числовые формулы дифференцирования (NDFs) | |
ode23s | X | Низко | Модифицированная формула Розенброка второго порядка | |||
ode23t | X | Низко | X | Метод трапеций с помощью interpolant | ||
ode23tb | X | Низко | X | TR-BDF2 |
Для ode15sode23t, и ode23tb выпадающее меню для Solver reset method появляется в разделе Solver details панели Настройки. Этот параметр управляет, как решатель обрабатывает вызванный сброс, например, обнаружением пересечения нулем. Позволенными опциями является Fast и Robust. Fast указывает, что решатель не повторно вычисляет якобиан для сброса решателя, тогда как Robust указывает, что решатель делает. Следовательно, Fast установка в вычислительном отношении быстрее, но она может использовать размер небольшого шага в определенных случаях. Чтобы протестировать на такие случаи, запустите симуляцию с каждой установкой и сравните результаты. Если нет никакого различия в результатах, можно безопасно использовать Fast установка и экономит время. Если результаты значительно отличаются, попытайтесь уменьшать размер шага для быстрой симуляции.
Для ode15s решатель, можно выбрать максимальный порядок числовых формул дифференцирования (NDFs), который применяет решатель. Начиная с ode15s использование сначала - через формулы пятого порядка, Maximum order параметр позволяет вам выбирать порядки 1 - 5. Для жесткой проблемы можно хотеть запуститься с порядка 2.
Следующая таблица обеспечивает советы для приложения переменного шага неявные решатели. Для примера, сравнивающего поведение этих решателей, смотрите, что Решатели Переменного Шага Исследования Используют Жесткую Модель.
Для жесткой проблемы решения могут измениться на масштабе времени, который очень мал по сравнению с интервалом интегрирования, в то время как решение интереса изменяется на намного более длительном масштабе времени. Методы, которые не спроектированы для жестких проблем, неэффективны на интервалах, где решение медленно изменяется, потому что эти методы используют временные шаги, достаточно маленькие, чтобы разрешить самое быстрое изменение. Для получения дополнительной информации смотрите Шемпина, L. F. Числовое Решение Ordinary Differential Equations, Chapman & Hall, 1994.
Решатели переменного шага используют стандартные методы управления, чтобы контролировать локальную ошибку на каждом временном шаге. Во время каждого временного шага решатели вычисляют значения состояния в конце шага и определяют локальную ошибку — предполагаемая ошибка этих значений состояния. Они затем сравнивают локальную ошибку с приемлемой ошибкой, которая является функцией и относительного допуска (rtol) и абсолютного допуска (атолл). Если локальная ошибка больше приемлемой ошибки для какого-либо состояния, решатель уменьшает размер шага и попробовал еще раз.
Relative tolerance измеряет ошибку относительно размера каждого состояния. Относительный допуск представляет процент значения состояния. Значение по умолчанию, 1e-3, означает, что вычисленное состояние точно в пределах 0,1%.
Absolute tolerance является пороговым ошибочным значением. Этот допуск представляет приемлемую ошибку как значение измеренного нуля подходов состояния.
Решатели требуют ошибки для iсостояние th, ei, чтобы удовлетворить:
Следующий рисунок показывает график состояния и областей, в которых относительный допуск и абсолютный допуск определяют приемлемую ошибку.
Ваша модель имеет глобальный абсолютный допуск, который можно установить на панели Решателя диалогового окна Configuration Parameters. Этот допуск применяется ко всем состояниям в модели. Можно задать auto или действительный скаляр. Если вы задаете auto (значение по умолчанию), Simulink первоначально устанавливает абсолютный допуск к каждому состоянию на основе относительного допуска. Если относительный допуск является большим 1e-3, abstol инициализируется в 1e-6. Однако для reltol меньший, чем 1e-3, abstol поскольку состояние инициализируется в reltol * 1e-3. В то время как симуляция прогрессирует, абсолютный допуск к каждому состоянию сброс к максимальному значению, которое состояние приняло до сих пор, времена относительный допуск к тому состоянию. Таким образом, если изменения состояния от 0 до 1 и reltol 1e-3, abstol инициализирует в 1e-6, и к концу симуляции достигает 1e-3 также. Если состояние идет от 0 до 1 000, то abstol изменения в 1.
Теперь, если изменения состояния от 0 до 1 и reltol установлен в 1e-4, затем abstol инициализирует в 1e-7, и к концу симуляции достигает значения 1e-4.
Если вычисленное начальное значение для абсолютного допуска не подходит, можно определить соответствующее значение сами. Вам придется запустить симуляцию несколько раз, чтобы определить соответствующее значение для абсолютного допуска. Можно также задать, должен ли абсолютный допуск адаптироваться так же к его auto установка путем включения или отключения AutoScaleAbsTol параметр. Для получения дополнительной информации смотрите Автоматическую шкалу абсолютный допуск.
Несколько блоков позволяют вам указывать, что абсолютные значения допуска для решения модели утверждают, что они вычисляют или которые определяют их выход:
Абсолютные значения допуска, которые вы задаете для этих блоков, заменяют глобальные настройки в диалоговом окне Configuration Parameters. Вы можете хотеть заменить глобальную установку, если, например, глобальная установка не обеспечивает достаточный контроль ошибок для всех ваших состояний модели, потому что они значительно различаются в величине. Можно установить блок абсолютный допуск на:
auto
–1 (тот же as auto)
positive scalar
real vector (наличие размерности равняется количеству соответствующих непрерывных состояний в блоке),
Если вы действительно принимаете решение установить абсолютный допуск, иметь в виду, что слишком низко значения заставляет решатель делать слишком много шагов около значений почти нулевого состояния. В результате симуляция медленнее.
С другой стороны, если вы устанавливаете абсолютный допуск слишком высоко, ваши результаты могут быть неточными как одно или несколько непрерывных состояний в вашем нуле подхода модели.
Если симуляция завершена, можно проверить точность результатов путем сокращения абсолютного допуска и выполнения симуляции снова. Если результаты этих двух симуляций удовлетворительно близки, то можно чувствовать себя уверенно в их точности.
