Writing a differential equation for iteration using for loop.

Question

Heya :) 2020-11-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/645423-writing-a-differential-equation-for-iteration-using-for-loop

评论： Heya :) 2020-11-13

采纳的回答： Alan Stevens

在 MATLAB Online 中打开

I have a differential equation:

p=@(t,x,y) (a-x^2-y^2)*x-omega*y;

Why for iteration we write this equation as ?

for i=2:50000
                x(i)=x(i-1)+((a-x(i-1)^2-y(i-1)^2)*x(i-1)-omega*y(i-1))*dt;
end

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

请先登录，再进行评论。

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

Answer 1

Alan Stevens 2020-11-12

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/645423-writing-a-differential-equation-for-iteration-using-for-loop#answer_542478

编辑：Alan Stevens 2020-11-12

在 MATLAB Online 中打开

This is known as simple Euler integration. It arises from

% dx/dt = (a - x^2 - y^2)*x - omega*y   
% Let dx/dt be approximated by (x(t+dt) - x(t))/dt
% and the right hand side by (a - x(t)^2 - y(t)^2)*x(t) - omega*y(t)
% so (x(t+dt) - x(t))/dt = (a - x(t)^2 - y(t)^2)*x(t) - omega*y(t)        approximately
% Multiply both sides by dt to get
% x(t+dt) - x(t) = ((a - x(t)^2 - y(t)^2)*x(t) - omega*y(t))*dt
% Add x(t) to both sides to get
% x(t+dt) = x(t) + ((a - x(t)^2 - y(t)^2)*x(t) - omega*y(t))*dt