|
Ситема ДУ большая – 22 ур-ния, все 1-ой степени.
Правые части – полиномы не выше 2-й степени.
При численном решении с помощью ode23s оказывается, что некоторые переменные с самого начала ползут вниз, при том что ПЧ соответствующих уравнений в начальный момент времени равны 0 (благодаря начальным условиям).
x'(0)=0, => д.б. x(0+)=x(0).
(Возможно, некоторые таким же макаром ползут вверх, но это я просто не замечаю; но когда неотрицательные по физ. смыслу переменные уходят в минус это не заметить невозможно)
С другими odeXX примерно то же самое (набор «плохих» x-ов немножко различается).
С подобными системами ДУ меньшего порядка (до 19) такой проблемы не возникало.
С чем это может быть связано? Погрешность вычислений? От больших систем начинает глючить? Я что-то делаю не так?
Думаю, скорее последнее, но куда тогда копнуть?
Запускал и из m-файла и через симулинк.
Mатлаб R2008a и R2014a.
Если понадобится, выложу систему и нач. условия; сразу не стал загромождать.
|
29.12.2017, 15:19 