Leslie模型matlab

Leslie模型是解决捕食问题的一类模型，本压缩包内含Leslie模型的MATLAB代码。

format long;
p=[2.12e-5 0.3965 0.6 0.8 0.6387];
for a=6:16
    p（a）=0.5688;
end
display（p）;
f=[0 0 0 0 0 80110 162700 212700 267900 326400 386000 444500 499700 549600 592200 592200]
display（‘初始各年龄组鱼数为：‘）
x0=[5.21e7 1100 443 266 213 136 77 44 25 14 8 5 3 1 1 1]‘
display（‘L矩阵为：‘）;
L=zeros（16）;
for b=1:16
    L（1b）=f（b）;
    if b>1
        L（bb-1）=p（b-1）;
    end
end
display（L）;
display（‘最大特征值c为:‘）;
c=max（eig（L））
for d=1:16
    e=d;
    x（d）=1;
    while（e-1>0）
        x（d）=p（e-1）*x（d）;
        e=e-1;
    end
    x（d）=x（d）*c^（1-d）;
end
display（‘c所对应特征向量为：‘）;
x=x‘
k=1;
x0=L*x0;
x1=x0-（x0（1）*x）;
s=max（abs（max（x1））abs（min（x1）））;
while s>=0.01
    k=k+1;
    x0=L*x0;
    x1=x0 - （x0（1）*x）;
    s=max（abs（max（x1））abs（min（x1）））;
end
display（‘生态条件不变时，鲈鱼达到稳定经过年数为：‘）
display（k）;
display（‘当经受工业污染之后成活率p变为：‘）;
p（1）=p（1）*0.75;
p（2）=p（2）*0.85;
for o=6:16
    p（o）=p（o）*0.9;
end
display（p）;
display（‘初始各年龄组鱼数为：‘）
x0=[5.21e7 1100 443 266 213 136 77 44 25 14 8 5 3 1 1 1]‘
display（‘L矩阵变为：‘）;
L=zeros（16）;
for b=1:16
    L（1b）=f（b）;
    if b>1
        L（bb-1）=p（b-1）;
    end
end
display（L）;
display（‘最大特征值c为:‘）;
c=max（eig（L））
display（‘由c值知道种群相比原来会减少得越快‘）;
for d=1:16
    e=d;
    x（d）=1;
    while（e-1>0）
        x（d）=p（e-1）*x（d）;
        e=e-1;
    end
    x（d）=x（d）*c^（1-d）;
end
display（‘c所对应特征向量（即种群最终年龄结构分布）为：‘）;
display（x）;
display（‘初始可捕捞鱼数:‘）
z0=0;
for q=4:16
z0=x0（q）+z0;
end
display（z0）;
g=1;
x0=L*x0;
z1=0;
for q=4:16
z1=x0（q）+z1;
end
while（z1>z0/2）
    z1=0;
    g=g+1;
    x0=L*x0;
for q=4:16
   z1=x0（q）+z1;
end
end
display（‘可捕捞鱼数减半的年数为：‘）
display（g）;
clear;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1806  2008-12-11 21:40  Leslie.m

     文件       3899  2017-06-24 22:30  leslie.txt

     文件      13079  2017-10-13 17:18  Leslie人口预测模型程序.docx

----------- ---------  ---------- -----  ----

                18784                    3