ODE Runge Kutta MATLAB 错误
ODE Runge Kutta MATLAB error
所以我正在尝试创建一个 Runge Kutta 函数,这是我的代码:
function [t,U] = RK(f, n, eta, interv)
h = (interv(2)-interv(1))/n;
t = interv(1):h:interv(2);
v(1) = eta(1);
w(1) = eta(2);
for i=1:n
k1 = f([v(i),w(i)]);
k2 = f([v(i),w(i)]+h*k1/2); %f(t(i)+h/2, u(:,i)+h*k1/2);
k3 = f([v(i),w(i)]+h*k2/2);
k4 = f([v(i),w(i)]+h*k3);
v(i+1) = v(i) + h*(k1(1)+2*k2(1)+2*k3(1)+k4(1))/6;
w(i+1) = w(i) + h*(k1(2)+2*k2(2)+2*k3(2)+k4(2))/6;
end
U = [v;w];
end
其中 U 是 2 行 n+1 列的矩阵,这是我尝试执行此函数时出现的问题,例如:
RK(sin, 10, [0,1], [5,15])
我收到错误 not enough input arguments
但是当我尝试将代码作为脚本执行并将 f 替换为 sin 时,一切正常,我得到了 U 矩阵
谁能告诉我解决办法是什么?
您正在使用 sin
作为函数句柄。只需添加 @
符号即可。
RK(@sin, 10, [0,1], [5,15])
所以我正在尝试创建一个 Runge Kutta 函数,这是我的代码:
function [t,U] = RK(f, n, eta, interv)
h = (interv(2)-interv(1))/n;
t = interv(1):h:interv(2);
v(1) = eta(1);
w(1) = eta(2);
for i=1:n
k1 = f([v(i),w(i)]);
k2 = f([v(i),w(i)]+h*k1/2); %f(t(i)+h/2, u(:,i)+h*k1/2);
k3 = f([v(i),w(i)]+h*k2/2);
k4 = f([v(i),w(i)]+h*k3);
v(i+1) = v(i) + h*(k1(1)+2*k2(1)+2*k3(1)+k4(1))/6;
w(i+1) = w(i) + h*(k1(2)+2*k2(2)+2*k3(2)+k4(2))/6;
end
U = [v;w];
end
其中 U 是 2 行 n+1 列的矩阵,这是我尝试执行此函数时出现的问题,例如:
RK(sin, 10, [0,1], [5,15])
我收到错误 not enough input arguments
但是当我尝试将代码作为脚本执行并将 f 替换为 sin 时,一切正常,我得到了 U 矩阵
谁能告诉我解决办法是什么?
您正在使用 sin
作为函数句柄。只需添加 @
符号即可。
RK(@sin, 10, [0,1], [5,15])