Остановить интеграцию по истечении заданного промежутка времени в Matlab

Во-первых, вы не можете сделать ничего подобного тому, что пытаетесь сделать. Все, что произойдет, это то, что ode15s будет работать столько, сколько ему нужно, и вернет результат, как обычно. Вам нужен способ прерывания решателя ОДУ в середине его выполнения, когда превышено определенное заданное количество времени. К сожалению, вы не предоставили исполняемый код, поэтому мне придется дать общие рекомендации, которые вы можете адаптировать к своей проблеме.

Одним из способов сделать это является использование функции вывода. Это функция, которая вызывается после каждого (успешного) шага интеграции решателя Matlab. Существует несколько примеров функций вывода, код которых вы можете изучить, чтобы увидеть, что вы можете сделать: odeplot, odeprint, odephas2 и odephas3. Просто введите edit odeplot в командном окне, чтобы увидеть код.

Ниже приведен пример того, как вы можете использовать эту возможность, используя встроенную функцию-пример vdp1000, которая упоминается в справке для ode15s. Вам нужно создать постоянную переменную, чтобы сохранить начальное время, инициализировать ее, а затем проверить прошедшее время относительно параметра, который вы передаете. Затем вы выдаете ошибку для прерывания ode15s, если прошедшее время превышает назначенное значение. Оператор try-catch может быть используется для обнаружения ошибки.

function ode_interupt_demo

tspan = [0 3000]; y0 = [2;0];
interupt_time = 0.05;
outputFun= @(t,y,flag)interuptFun(t,y,flag,interupt_time);
opts = odeset('AbsTol',1e-8,'RelTol',1e-5,'OutputFcn',outputFun);

try
    [t,y] = ode15s(@vdp1000,tspan,y0,opts);
catch ME
    if strcmp(ME.identifier,'interuptFun:Interupt')
        disp(ME.message);
        % Do other things
    else
        rethrow(ME); % It's possible the error was due to something else
    end
end


function status = interuptFun(t,y,flag,interupt_time)   %#ok<INUSL>
persistent INIT_TIME;
status = 0;
switch(flag)
    case 'init'
        INIT_TIME = tic;
    case 'done'
        clear INIT_TIME;
    otherwise
        elapsed_time = toc(INIT_TIME);
        if elapsed_time > interupt_time
            clear INIT_TIME;
            str = sprintf('%.6f',elapsed_time);
            error('interuptFun:Interupt',...
                 ['Interupted integration. Elapsed time is ' str ' seconds.']);
        end
end

Одним из недостатков этого подхода является то, что вы не получите возвращаемого значения из ode15s. Это может быть или не быть желательным в вашем приложении. Вероятно, есть способы обойти это, но они потребуют больше кода. Одна простая вещь, которую вы можете сделать, это заставить сообщение об ошибке возвращать последнее время и состояние системы (когда она была прервана), используя значения t и y, переданные в выходную функцию. Также возможно, что для прекращения интеграция на основе прошедшего времени — это может позволить вам по-прежнему возвращать вывод из ode15s.