MATLAB: Matlab ode function help!!!

odeode45

Hi. I am using the Matlab ode45 function for my program:

[t,X]=ode15s('odefunc',[tstart:tspan:tfinal],Xinit); %ode function say tstart=0;tspan=100;tfinal=40,000

Now, my problem is I want to divide the calculations within this ode function into several cycles, say one cycle=2000seconds. How do I do that in Matlab?

Best Answer

It seems, like you are not using ODE45, but ODE15S.

Currently you divide the simulation in 401 parts by using [tstart:tspan:tfinal]. (btw. you can omit the square brackets. a:b:c is a vector already). Now you want to divide it to 2000 seconds. If your tfinal is measured in seconds, why not simply use tspan=2000 ?

Related Solutions

MATLAB: Probability plot HELP!!!

To fit a gaussian contour to your 2d data, you could try fitgmdist.

To plot traces of the mean & max pdfs, compute the mean and std of the measurement of interest (I can't tell from this figure if its supposed to be pressure or temperature) then plot the pdf.

mew = mean(x); sig = std(x);
y = 1/sqrt(2*pi*sig^2)*exp(-(x-mew).^2/2/sig^2);
plot([y y(1)],[x x(1)],'k--')

MATLAB: 4th order Runge-Kutta code that can solve for several intial conditions

Here's a version:

% Initial conditions
t0=0;
x0=0:0.2:1;
y0=1:0.2:2;
dt=0.1;
tz=1;
t_range= t0:dt:tz;
X = zeros(numel(x0),numel(t_range));
Y = zeros(numel(y0),numel(t_range));
x_rk=zeros(1,numel(t_range));
y_rk= zeros(1,numel(t_range));
for n = 1:numel(x0)
    x_rk(1)=x0(n);
    y_rk(1)=y0(n);
    for i=1:length(t_range)-1
        [x_rk(i+1), y_rk(i+1)]=rk_method(x_rk(i),y_rk(i),dt);
    end
    X(n,:) = x_rk;
    Y(n,:) = y_rk;
    %% PLots
    subplot(numel(x0)/2,2,n)
    plot (t_range, x_rk,'b-o',t_range, y_rk,'y-o','MarkerSize',5);
    grid
    title([num2str(x0(n)),'   ',num2str(y0(n))])
    legend('x','y')
end
disp('X = '), disp(X)
disp('Y = '), disp(Y)
%% Functions
function [dxdt, dydt]= f(x,y)
dxdt =  5*x-3*y;
dydt = -6*x+2*y;
end 
function [x1, y1]=rk_method(x0,y0,dt)
[k11, k12]=f(x0,y0);     
[k21, k22]=f(x0+0.5*dt*k11,y0+0.5*dt*k12); 
[k31, k32]=f(x0+0.5*dt*k21,y0+0.5*dt*k22); 
[k41, k42]=f(x0+dt*k31, y0+dt*k32); 
x1=x0+dt*((k11/6)+((k21+k31)/3)+(k41/6));
y1=y0+dt*((k12/6)+((k22+k32)/3)+(k42/6));
end

Best Answer

Related Solutions

MATLAB: Probability plot HELP!!!

MATLAB: 4th order Runge-Kutta code that can solve for several intial conditions

Related Question