| | ruskiy COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 20:53 |
| | calculer les exposants de lyapunov sur octave sans les toolboxes!!!
bientot le test BDS à la main sans programme extérieur! | |
| Unpseudoalacon COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 20:54 |
| | A tes souhaits | |
| associal.return COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 20:55 |
| | elle est cool ta vie!!! | |
| Fanny30 COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 20:55 |
| | Unpseudoalacon (Godlike) (EN LIGNE) Mardi 10 Juin 2008 à 20:54
Godlike - 20 ans - Montpellier (Hérault) A tes souhaits
merci d'avoir résumé mes pensées | |
| guigui30700 COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 20:56 |
| | j'ai faillit me faire un noeud a la langue a essayer de lire tes conneries !  | |
| FP.WWW.FP COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 20:56 |
| | Oula ????
French ?
English ? | |
| Nop COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 21:33 |
| | j'en suis content pour toi , néanmoins , j'ai strictement rien compris a l'éxplication de wikipédia , si tu pouvais nous éxpliquer l'utilité de la chose , tout du moins de lyapunov pour commencer | |
| Mickael.M COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 21:42 |
| | Matlab cay le mal. | |
| | Spival Régis Mardi 10 Juin 2008 à 21:55 |
| | Si tu sais utiliser sous ASP un composant COM compilé avec VB.NET, tu pourrais m'être utile. | |
| ruskiy COMPTE SUPPRIMÉ Mardi 10 Juin 2008 à 22:30 |
| | s=1;
while s>0;
r=input('input growth rate r: ')
x0=input('input initial population x0: ')
z0=input('input second initial population z0: ')
n=input('end of time interval b: ')
x=zeros(n+1,1);
z=zeros(n+1,1);
t=zeros(n+1,1);
x(1)=x0;
z(1)=z0;
tt(1)=0;
for i=1:n
t(i)=i-1;
x(i+1)=r*x(i)*(1-x(i));
z(i+1)=r*z(i)*(1-z(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);
for i=1:nn+1
xx(i)=xstart+(i-1)*del;
lin(i)=xx(i);
yy(i)=r*xx(i)*(1-xx(i));
end
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)=r*x0*(1-x0);
yc(3)=yc(2);
xc(3)=yc(2);
zc=zeros(24,1);
zyc=zeros(24,1);
zc(1)=z0;
zyc(1)=0;
zc(2)=z0;
zyc(2)=r*z0*(1-z0);
zyc(3)=zyc(2);
zc(3)=zyc(2);
plot(xx,lin,xx,yy,xc,yc,zc,zyc),pause
for j=4:12;
jj=2*j-4;
xc(jj)=xc(jj-1);
yc(jj)=r*xc(jj)*(1-xc(jj));
xc(jj+1)=yc(jj);
yc(jj+1)=yc(jj);
zc(jj)=zc(jj-1);
zyc(jj)=r*zc(jj)*(1-zc(jj));
zc(jj+1)=zyc(jj);
zyc(jj+1)=zyc(jj);
plot(xx,lin,xx,yy,xc,yc,zc,zyc),pause
end
plot(t,x,t,x,'o',t,z,t,z,';+'),pause
s=input('Do you want to stop - if so enter 0')
end
d=(x-z)^2;
e=1:500; % si iteration de 500
k=log(d/(x(1)-z(1))²)./e';
plot(k);
lyap_est=max(k); %estimateur de lyapunov
voilà!!! bien à vous!!! | |
| |
