ec2a

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

clear
close all
clc
r=[0.5; 1.5; 2.5; 3.2; 3.5; 3.6; 3.8; 4];
T=30; % período de simulação
x=0.95*ones(size(r)); % vetor com a população inicial
X=x;
for it=1:T-1
x=r.*x.*(1-x);
X=[X x];
end
tempo=(1:T)-1;
for i=1:length(r)
subplot(length(r)/2,2,i)
stem(tempo,X(i,:))
title(['r=' num2str(r(i))])
end
print -dpng simumapalogistico2 % salva a figura
fprintf('Periodo de oscilacao para r=3.5 => 4 \n')
fprintf('Periodo de oscilacao para r=3.6 => 14 \n')
fprintf('Periodo de oscilacao para r=4 => Sinal não periódico \n') %não é possível observar periodicidade no sinal