Matlab - Problema de Convergencia de Cobwebbing

   
Vista:

Problema de Convergencia de Cobwebbing

Publicado por Joel (1 intervención) el 08/11/2011 02:24:29
Hola amigos mi problema consiste en lo siguiente:

Quiero utilizar el meto de cobwebbing par la siguiente serie: x(n+1)=1/(2-x(n))
pero cuando utilizo el codigo(lo pondre en la parte inferior) no me permite q el valor converja al punto de equilibrio pues los limite que se le dan llegan hasta cierto punto y luego se mantiene en una sola zona. Quisiera que me den una sugerencia para poder corregir la convergencia de mi programa (traté de cambiar los limites de la funcion identidad para que avance y converja pero no me funciona). gracias por su atencion.

%inicializo la cantidad de iteraciones
x0 = 0.2;
n = 80;

% Inicializando
x=zeros(n+1,1);
t=zeros(n+1,1);
x(1)=x0;
tt(1)=0;

%Iteracion de la ecuacion diferencial logistica
for i=1:n
t(i)=i-1;
x(i+1)=-x(i)^2*(1-x(i));
end

t(n+1)=n;
nn=100;
del=1/nn;
xstart=0;
yy=zeros(nn+1,1);
lin=zeros(nn+1,1);
xx=zeros(nn+1,1);

%Aqui esta la recta identidad
for i=1:nn+1
xx(i)=xstart+(i-1)*del;
lin(i)=xx(i);
yy(i)=-xx(i)^2*(1-xx(i));
end

%Secuencia de las verticales y horizontales para la convergencia
plot(xx,lin,xx,yy),pause
xc=zeros(24,1);
yc=zeros(24,1);
xc(1)=x0;
yc(1)=0;
xc(2)=x0;
yc(2)=-x0^2*(1-x0);
yc(3)=yc(2);
xc(3)=yc(2);

plot(xx,lin,xx,yy,xc,yc),pause
for j=3:20;
jj=2*j-4;
xc(jj)=xc(jj-1);
yc(jj)=-xc(jj)^2*(1-xc(jj));
xc(jj+1)=yc(jj);
yc(jj+1)=yc(jj);
plot(xx,lin,xx,yy,xc,yc),pause
end

plot(t,x,t,x,'o');
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder