Time varying force with ode23

Question

Christo van Rensburg 2020-4-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/520621-time-varying-force-with-ode23

评论： Christo van Rensburg 2020-4-25

采纳的回答： darova

在 MATLAB Online 中打开

Hi there

I have the following code:

tf = 1;
t = linspace(0,3);
Ft = zeros(length(t),1);
% Here I generate the force values at different time intervals.
for i = 1:length(t)
    if t(i) <= tf
        Ft(i) = Fo*sin(2*pi*t(i)/tf);
    elseif t(i) > tf
        Ft(i) = 0;
    end
end
ic = [0, 0]; % Initial conditions
dXdt = @(t, x)     [x(2);
                            Ft * 1/(m1+m2)];
        
[t, x] = ode23(dXdt, [0, 3], ic);

I'm struggling to incorportate this time varying force F(t) into the ode23 function.

Any help as to how I can maybe get this done?

Thanks in advance!

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

请先登录，再进行评论。

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

Answer 1

darova 2020-4-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/520621-time-varying-force-with-ode23#answer_428301

在 MATLAB Online 中打开

Use function handle

Ft = @(t) (t<=tf)*Fo*sin(2*pi*t/tf);
dXdt = @(t, x)     [x(2);
                        Ft(t) * 1/(m1+m2)];

2 个评论
显示无隐藏无

Walter Roberson 2020-4-25

That would work if the driving force was not discontinuous.

Christo van Rensburg 2020-4-25

It seems to work though.

Walter I'll try your method now, because what you're saying makes sense.

Thanks for the help!

请先登录，再进行评论。

Answer 2

Walter Roberson 2020-4-25

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/520621-time-varying-force-with-ode23#answer_428302

在 MATLAB Online 中打开

You cannot do that the way you want. The ode* routines cannot handle discontinuities in the driving functions. You need to break it up into two time periods, one with the force and one without the force.

tf = 1;
dXdt_F = @(t, x) [x(2);  Fo*sin(2*pi*t/tf) / (m1+m2)];
dXdt_NF = @(t, x) [x(2); 0];
ic = [0, 0]; % Initial conditions
[t1, x1] = ode23(dXdt_F, [0, tf], ic);
ic2 = x1(end,:);
[t2, x2] = ode23(dXdt_NF, [tf, 3], ic2);
t = [t1; t2(2:end)];
x = [x1; x2(2:end,:)];