登录论坛

查看完整版本 : [求助]关于虫口问题


huanghuan
2008-11-05, 18:26
典型的虫口混沌问题。状态的递推关系为Xn+1=a*(Xn-Xn^2)

Xn+1和Xna分别为次年和当年的虫口数。 x1=0.5 , n=50。用图形的方法研究系数a的变化为a=[2.0 3.2 3.5 3.8]对xn+1变化的影响。图形横、纵坐标分别为n和Xn+1。

我的程序如下:

n=50;
a=[2.0 3.2 3.5 3.8];
x(1)=0.5;
[p,k]=size(a);
for j=1:k
for i=1:n
x(i+1)=a(j)*(x(i)-x(i)^2)
plot(i,x)
end
end

此程序只能求出n=50时的情况,望高手指点。

未注册
2008-11-06, 16:11
[p,k] = size(a)应该换为k=length(a)

不是一样的吗?

Juicy
2008-11-11, 07:46
不是一样的吗?

恩,我好糊涂:sweat:

pslpsl
2008-11-13, 12:03
n=50;
a=[2.0 3.2 3.5 3.8];
x(1)=0.5;
[p,k]=size(a);
for j=1:k
for i=1:n
x(i+1)=a(j)*(x(i)-x(i)^2);
end
end
plot(1:n+1,x)