Hi Vörös Gergely,
I understand that you are trying to switch between two similar functions, 'a' and 'b', during the integration of an ODE in MATLAB. You wish to use 'a' for the initial phase of the integration and then switch to 'b' after t = 200s. This can be accomplished using state-events in MATLAB's ODE solvers.
Please follow the steps given below to achieve this:
- Define the two functions: Start by defining the two functions 'a' and 'b'.
- Implement the Event Function: Event function will trigger when the integration time reaches 200 seconds. This function will stop the current integration. Set the 'Events' as this function's handle in odeset.
- Setup the Integration with Event Handling: Perform the first part of the integration using 'a'. When the event triggers, stop the integration.
- Continue Integration with 'b': Use the final state from the first integration as the initial condition for continuing the integration with 'b'.
Here is an example for your reference:
% Initial conditions
y0 = 1;
tspan = [0, 400]; % Total time span
% Perform the first part of the integration using 'a'
options1 = odeset('Events', @eventFcn);
[t1, y1, te, ye, ~] = ode45(@funcA, tspan, y0, options1);
% Continue with 'b' from the last point
tspan2 = [te, tspan(2)];
y0_2 = ye(end, :); % Initial condition for the next part
[t2, y2] = ode45(@(t, y) funcB(t, y), tspan2, y0_2);
% Combine the results
t = [t1; t2(2:end)]; % Avoid repeating the event time
y = [y1; y2(2:end)]; % Avoid repeating the event state
% Plot the results
figure;
plot(t, y, 'LineWidth', 2);
xlabel('Time');
ylabel('y');
grid on;
% function 'a'
function dydt = funcA(~, y)
dydt = -0.01 * y;
end
% function 'b' with an additional parameter
function dydt = funcB(~, y)
dydt = -0.01 * y + 0.5;
end
% Event function to switch from 'a' to 'b' at t = 200
function [value,isterminal,direction] = eventFcn(t, ~)
value = t - 200; % Trigger at t = 200
isterminal = 1; % Stop the integration
direction = 0; % Trigger regardless of the direction of the zero crossing
end
Pleasre refer to the following documentation link to learn more about odeset.
Hope this helps!
