Hello, I want to make function Runge-Kutta orde 4 with equation depend on time, but in MATLAB i get some error

Question

cindyawati cindyawati 2023-5-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1968079-hello-i-want-to-make-function-runge-kutta-orde-4-with-equation-depend-on-time-but-in-matlab-i-get

评论： Cris LaPierre 2023-5-20

采纳的回答： Cris LaPierre

在 MATLAB Online 中打开

I want to use this equation with runge-kutta orde 4

function M1 =fRK4M1(M1)
 %M1 =M1(j)+1./(1+exp(-T(j)));
 m1 = (50*M1(j+1)*(1-(M1(j+1)/75))-2*(10^-4)*M1(j+1)*M1(j+1)-M1(j+1)*((5*10^-4).*M2(j+1))-((10-3)*(10^-4)*M1(j+1)*M3(j+1))-((100-10)*0.1*M1(j+1)*O(j+1))-((10^-3)*M1(j+1)));
 M1 = dt *m1;
end

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

请先登录，再进行评论。

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

Answer 1

Cris LaPierre 2023-5-19

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1968079-hello-i-want-to-make-function-runge-kutta-orde-4-with-equation-depend-on-time-but-in-matlab-i-get#answer_1239729

在 MATLAB Online 中打开

You have not defined j, so in MATLAB, by default that letter is used to create complex numbers. As the error message states, your index must be positive integers.

j+1
ans = 1.0000 + 1.0000i
a=1:5;
a(j+1)
Array indices must be positive integers or logical values.

2 个评论
显示无隐藏无

cindyawati cindyawati 2023-5-20

编辑：Cris LaPierre 2023-5-20

在 MATLAB Online 中打开

Thank you for your response @Cris LaPierre. So, I must using looping like this? Because j is depend on time

t(1)=0;
dt=0.01; %time interval
t=0:dt:100; %time span
for j = 1:length((t)-1) 
    T(j+1)=T(j)+dt;
    M1 =M1(j)+1./(1+exp(-T(j)));
end
Unrecognized function or variable 'T'.