0

我试图让 matlab 显示猎物与捕食者的图表

function [ output_args ] = Untitled( input_args )

    options = odeset('RelTol', 1e-4, 'NonNegative', [1 2]);

    [t,x] = ode45('lotka_volterra', [0 30], [2 1], options);

    plot(t,x);
    legend('prey', 'predators');

end


function dxdt = lotka_volterra(t,x)

    a     = 1.2;
    b     = 0.6;
    d     = 0.3;
    gamma = 0.8;

    dxdt    = [0;0];
    dxdt(1) = a * x(1) - b * x(1) * x(2);
    dxdt(2) = d * x(1) * x(2) - gamma * x(2);

end

重点是展示生活在同一环境中的两个物种的种群之间的动态关系。

两个导数(dxdt(1) 和 dxdt(2))模拟每个总体的变化率

我尝试使用 Runge-Kutta 积分方法来获得两个物种的种群与时间的关系图

错误只是说:

Error in Untitled (line 3)
[t,x] = ode45(@(t,x) lotka_volterra, [0 30], [2 1]);

我很难过任何帮助都会很棒

4

1 回答 1

1

您可能想要使用函数句柄:

[t,x] = ode45(@lotka_volterra, [0 30], [2;1], options);

或者,更慢但可能更清晰:

[t,x] = ode45(@(t,x) lotka_volterra(t,x), [0 30], [2;1], options);
于 2016-03-31T01:54:06.597 回答