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（‘初始各年齡組魚數(shù)為：‘）
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所對應(yīng)特征向量為：‘）;
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（‘生態(tài)條件不變時(shí)，鱸魚達(dá)到穩(wěn)定經(jīng)過年數(shù)為：‘）
display（k）;
display（‘當(dāng)經(jīng)受工業(yè)污染之后成活率p變?yōu)椋骸?
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（‘初始各年齡組魚數(shù)為：‘）
x0=[5.21e7?1100?443?266?213?136?77?44?25?14?8?5?3?1?1?1]‘
display（‘L矩陣變?yōu)椋骸?
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所對應(yīng)特征向量（即種群最終年齡結(jié)構(gòu)分布）為：‘）;
display（x）;
display（‘初始可捕撈魚數(shù):‘）
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（‘可捕撈魚數(shù)減半的年數(shù)為：‘）
display（g）;
clear;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

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

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

?????文件??????13079??2017-10-13?17:18??Leslie人口預(yù)測模型程序.docx

-----------?---------??----------?-----??----

????????????????18784????????????????????3