CobWeb Plot em Octave

Publicado por Zaraki Kenpachi (última atualização em 01/03/2012)

[ Hits: 3.672 ]

Homepage: none

Download chaos.m




Descrição gráfica da dinâmica de uma aplicação quadrática f(x) = x² + L , iniciando no ponto x0 e tomando n iteradas.

  



Esconder código-fonte

#########################################################################################
#####                                                                                                                                                              #####
#####              Gráfico CobWeb da função quadrática f(x) = x²+l                                                                 #####  
#####                                                                                                                                                              #####
#########################################################################################


function [u,v] = chaos(x0, l, n)

# x0 = ponto inicial, isto é, f^0(x0);
# l = parâmetro da função quadrática;
# n = número de iterações;

k=2*n;

x= zeros(1,k+1);
y=zeros(1,k+1);

x(1)=x0;
y(1)=0;

   for i=2:k;
      if (rem (i, 2) == 0) #i é par
      x(i) = x(i-1);
      y(i) = x(i-1).*x(i-1)+l;
      else #i é ímpar
      x(i) = y(i-1);
      y(i) = y(i-1);
      endif
   end

 p=((1+sqrt(1-4*l))./2);

 u=-p:0.01:p;      # p=((1+sqrt(1-4*l))./2) é o maior ponto fixo.
 v=u.*u +l;        # note que precisamos ter \delta = 1 - 4*l > 0  

plot(x,y,u,v,u,u)


endfunction

Scripts recomendados

Octave - Método de Euler-Modificado para solução de EDO

Algorítimo de maior e menor

Cálculo da raiz de uma função pelo Método da Bisseção - Octave

Octave - Calcular raiz pelo método da Secante

Função Fatorial Recursiva em Fortran 95


  

Comentários

Nenhum comentário foi encontrado.


Contribuir com comentário




Patrocínio

Site hospedado pelo provedor RedeHost.
Linux banner
Linux banner
Linux banner

Destaques

Artigos

Dicas

Tópicos

Top 10 do mês

Scripts