《高等数学》课程教学资源（matlab及案例）7Logistic方程求解与混沌

实验7Logistic方程求解与混沌

实验7 Logistic方程求解与混沌

Logistic方程与混沌在生物学中，有一个刻画生物种群个体总量增长情况的著名的方程一一逻辑斯谛（Logistic）方程：r为比例系数Xn+1 = rx,(1- xn)其中x.为某生物群体的第n代的个体总数与该群体n所能达到的最大保有量时的个体数之比选定初值和比例系数r的值后，由方程就能生成一n个数列:Xi,X2,L ,X

Logistic方程与混沌 n 在生物学中，有一个刻画生物种群个体总量增长情 况的著名的方程——逻辑斯谛（Logistic）方程： n 其中xn为某生物群体的第n代的个体总数与该群体 所能达到的最大保有量时的个体数之比。 n 选定初值和比例系数r的值后，由方程就能生成一 个数列：

考察迭代格式(Logistic方程）初值 Xx。=0.1Xn+1 = rx,(1- x,)1.当参数r取值分别为1.2，2.5，3.2，3.5，3.8考察其迭代序列的收敛情况

考察迭代格式(Logistic方程 ) 初值 1. 当参数r取值分别为1.2，2.5，3.2，3.5，3.8 考察其迭代序列的收敛情况

程序clc;clf;x=0.1; y=[];r=1.2;%改变取值得到相应的图形hold onaxis([0 100 0 1])for i=1:100x=r*x*(1-x);y=[y,x];plot(i,x,k.,markersize',10)fprintf(x(%d)=%.10f/n',i,x);endt=1:100;plot(t,y,'k-');grid

clc;clf; x=0.1; y=[ ]; r=1.2; %改变取值得到相应的图形 hold on axis([0 100 0 1]) for i=1:100 x=r*x*(1-x);y=[y,x]; plot(i,x,'k.','markersize',10) fprintf('x(%d)=%.10f\n',i,x); end t=1:100; plot(t,y,'k-'); grid 程序

.....--------------------------:---1-----11-I-:09r11-112---------------nV--0.8--------....1--1-11-1110.72-J--11111:--------.----...----.--1----:.I.10.6-11111.--.--------------:-------:-+-0.5----.----------------:-1111:--:---:--0.4L41.1>J1J....--------.-1:----------------0.3-.......-.....---..-....--------1:1----:--0.224oo00o0086中o809000o4008e80e042tooooooecolooocoecoolocooooecloooeooooothtjooo-----111--1-1-.0.11VLVJ1J11J----.-----...-.---------02050701030406080900100

10.90.8米米米米冰米米米米米临米米0.7米米米*米米米米米米米米米米米米米0.60.50.4r=30.3*0.20.1011/110203040506070809010100

10.9米y****************************+****************0.8*0.7米0.6米0.50.420.3*0.2r=3.50.10050100150200250300

0.9料米米米米米米米0.80.70.6米米米米米米米米米米米米米米米米米米咪米米米米米米米米米米米0.50.40.30.2r =3.550.10500100150200250300

0.9*楼米0.8米米0.70.6米*0.5****米?0.4米米*******米米米米楼米米森料米米米E0.30.2r =3.580.111上1111110050100150200250300350400450500

1r米米热费楼楼米冰米米梦米拌米*发楼李楼楼L1i米F0.9*米华来米菜米X萧米4米**米米米*楼楼楼*******米米米米米料*米米華*米*禁米*精***0.8淼*********米米*米米**￥米楼**类****華**米*******0.7*米米-米米米米淼米米椒米*林米米米米?*米**楼米下米****米*米米米米楼米米*类米******米0.6***米米米类米*****米米样米米米楼米***米米***米米0.5米*****米米*米?楼十**米********米*米***0.4米*棒博米**米***+****-**0.3*********菜***米******米米*米米米米擀米米米米米***米0.2s楼米米米米***米米米米X**米*米3.80.1r=111L111110050050100400450150200300350250