帮我利用MATLAB敲出马尔萨斯模型

马尔萨斯模型是一种描述人口增长的理论模型，它假设资源有限而人口可以无限增长，导致人口过剩。在MATLAB中，你可以创建一个简单的函数来模拟这个模型。下面是一个基本的示例，我们将使用指数增长来表示人口：

% 定义常数
r = 0.05; % 自然增长率
M = 1000; % 资源上限

% 初始化人口和资源
population = 1;
resources = M;

% 设置时间步长和迭代次数
tspan = [0 100]; % 时间范围，这里设为100年
num_steps = 1000; % 迭代次数

% 创建时间向量
time = linspace(tspan(1), tspan(2), num_steps);

% 创建空矩阵存储结果
results = zeros(num_steps, 2); % 人口和资源矩阵

for i = 1:num_steps
    % 人口按自然增长率增长，直到达到资源限制
    if population <= resources
        population_next = population * (1 + r);
    else
        population_next = resources;
    end
    
    % 更新资源，假设除了人口消耗外没有其他消耗
    resources_next = resources - population;
    
    results(i, :) = [population_next, resources_next];
    population = population_next;
    resources = resources_next;
end

% 绘制人口和资源随时间变化的趋势图
figure;
plot(time, results(:, 1), 'b', time, results(:, 2), 'r');
xlabel('时间（年）');
ylabel('数量');
legend('人口', '资源');
title('马尔萨斯模型');