我必须求解以下形式的常微分方程组:
dx/ds = 1/x * [y* (g + s/y) - a*x*f(x^2,y^2)]
dy/ds = 1/x * [-y * (b + y) * f()] - y/s - c
其中x和y是我需要找出的变量,s是自变量;其余的都是常数。到目前为止,我试图用 ode45 解决这个问题,但没有成功:
y = ode45(@yprime, s, [1 1]);
function dyds = yprime(s,y)
global g a v0 d
dyds_1 = 1./y(1) .*(y(2) .* (g + s ./ y(2)) - a .* y(1) .* sqrt(y(1).^2 + (v0 + y(2)).^2));
dyds_2 = - (y(2) .* (v0 + y(2)) .* sqrt(y(1).^2 + (v0 + y(2)).^2))./y(1) - y(2)./s - d;
dyds = [dyds_1; dyds_2];
return
其中@yprime 具有方程组。我收到以下错误消息:
YPRIME 返回一个长度为 0 的向量,但初始条件向量的长度为 2。YPRIME 返回的向量和初始条件向量必须具有相同数量的元素。
有任何想法吗?谢谢