写出MATLAB模拟lotka-volterra竞争模型的代码

Lotka-Volterra竞争模型是用来描述两种物种在一个生态系统中相互竞争的模型，它可以用一组微分方程来表示。MATLAB代码如下：

% 参数定义
alpha = 1;
beta = 0.5;
gamma = 0.5;
delta = 2;

% 定义初始值
x0 = 2;
y0 = 1;

% 定义时间范围，步长和时间点个数
tspan = [0 100];
dt = 0.01;
t = tspan(1):dt:tspan(2);

% 定义微分方程
f = @(t, xy) [alpha*xy(1) - beta*xy(1)*xy(2); -gamma*xy(2) + delta*xy(1)*xy(2)];

% 解微分方程
[~, xy] = ode45(f, t, [x0,y0]);

% 绘图
plot(t, xy(:,1), 'r', t, xy(:,2), 'b')
xlabel('时间')
ylabel('物种数量')
legend('物种1', '物种2')

其中，参数`alpha`、`beta`、`gamma`、`delta`分别表示物种1的出生率、死亡率和与物种2的相互作用强度，物种2的死亡率和与物种1的相互作用强度。初始值`x0`和`y0`表示物种1和2的初始数量。`tspan`定义了时间范围，`dt`为时间步长，`t`为时间点序列。`f`定义了微分方程，`ode45`函数用来解微分方程，返回时间序列`t`和物种数量序列`xy`。最后，用`plot`函数将两个物种的数量随时间的变化绘制出来。