种群多样性计算matlab

种群多样性可以通过多种指标进行计算，常用的有Shannon指数、Simpson指数、Pielou指数等。以下是通过matlab计算这些指数的示例代码：

1. Shannon指数

function shannon = shannon_index(population)
    % 计算Shannon指数
    [~, n] = size(population);
    shannon = 0;
    for i = 1:n
        p = sum(population(:,i)) / sum(sum(population));
        if p ~= 0
            shannon = shannon - p * log2(p);
        end
    end
end

2. Simpson指数

function simpson = simpson_index(population)
    % 计算Simpson指数
    [~, n] = size(population);
    simpson = 0;
    for i = 1:n
        p = sum(population(:,i)) / sum(sum(population));
        simpson = simpson + p^2;
    end
    simpson = 1 / simpson;
end

3. Pielou指数

function pielou = pielou_index(population)
    % 计算Pielou指数
    [~, n] = size(population);
    pielou = 0;
    for i = 1:n
        p = sum(population(:,i)) / sum(sum(population));
        if p ~= 0
            pielou = pielou + p * log2(p);
        end
    end
    pielou = -pielou / log2(n);
end

这些函数接受一个种群矩阵作为输入，其中每一列代表一个物种在不同样本中的个体数量。可以通过将这些函数结合起来，来评估种群多样性的不同方面。