[求助]关于虫口问题

[huanghuan]

典型的虫口混沌问题。状态的递推关系为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时的情况，望高手指点。

作者: Juicy

[p,k] = size(a)应该换为k=length(a)

不是一样的吗？

[Juicy]

作者: 未注册

不是一样的吗？

恩，我好糊涂:sweat:

[pslpsl]

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)