Finding average step size and slope with ode45

Question

EB 2019-10-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/484460-finding-average-step-size-and-slope-with-ode45

评论： Star Strider 2019-10-12

采纳的回答： Star Strider

I want to be able to calculate the average step size and slope for each value in the loop.

I first want to plot the log values of dT and TOL.

I believe my average step size (dT) should be the average of 1/T.

t = [0, 5]; % time span

y0 = 1; % initial condition

dy0 = 2; % initial condition first derivative

dydt = @(t,y) [y(2); ((1 - y(1)^2))*y(2) - y(1)];

for j = 1:10

TOL(j) = 1*exp(-1-j);

options = odeset('AbsTol',TOL(j),'RelTol',TOL(j));

[T{j},Y{j}] = ode45(dydt, t, [y0, dy0], options);

dT = mean(1/length(T)); % average step size

plot(log(dT),log(TOL));

hold on;

grid

j = j+1;

end

hold off;

Errors I'm facing:

my plot doesn't show up
my slope receives the error 'Undefined function 'log' for input arguments of type 'cell'.'
not sure if I am thinking of average step size correctly

Side note: Would I be able to generate multiple slopes using polyfit in one graph?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2019-10-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/484460-finding-average-step-size-and-slope-with-ode45#answer_395582

在 MATLAB Online 中打开

Subscript ‘dT’ as well. The code will be more efficient if you move the plot outside the loop:

t = [0, 5];    % time span
y0 = 1;         % initial condition
dy0 = 2;     % initial condition first derivative
dydt = @(t,y) [y(2); ((1 - y(1)^2))*y(2) - y(1)];
for j = 1:10
    TOL(j) = 1*exp(-1-j);
    options = odeset('AbsTol',TOL(j),'RelTol',TOL(j));
    [T{j},Y{j}] = ode45(dydt, t, [y0, dy0], options);
    
    dT(j) = mean(1/length(T)); % average step size    
end
figure
plot(log(dT),log(TOL));
grid