Aide MatLab

[invite9e95248d]

Salut,

Alors voila, j'ai un petit algorithme à modéliser, pour illustré un théorème (théoreme de Sharkovskii pour ceux que ça interesse) mais je ne connais pas du tout matlab et je n'aurais surement pas l'occasion de m'en resservir donc j'aimerais savoir si qqun pouvait me donner le code, surtout que je pense que ça ne doit pas etre long.

Voila l'algorithme:

f(x)=x(1-x)

On ne s'interesse qu'a l'intervalle [0,1]

On choisit un point Xo et on calcule

f(Xo)

puis on prend l'intersection entre la droite y=f(Xo) et y=x

Ca donne un nouveau point X1, et on recommence le meme procédé.

Si c'est possible ça serait super si on pouvait voir la construction se faire sur un graphique ^^
Voila si qqun pouvait me filer un coup de main ça serait cool
Merci
-----

[azt]

Salut,
voilà ce que j'écrirais, j'espère que cela correspond à tes attentes

x=0:0.05:1;
y=(x.*(1-x));

newplot;
hold on;

    axis([0 1 0 0.25])
    plot(x,x);
    plot(x,y);

    x0=0.7;
    max=20;

    for i=0:max ;
        plot([x0;x0],[x0;x0.*(1-x0)],'-+');
        plot([x0; x0.*(1-x0)],[x0.*(1-x0);x0.*(1-x0)],'-+');
        x0=x0.*(1-x0);
    end;
    
hold off;

[invite9e95248d]

merci
je vais essayé ça, en tout cas ça me donne déjà une idée générale
Si je comprends ton code, la il va me calculer les 20 premieres itérations c'est ça ?

[azt]

Tu as tout compris !

[A voir en vidéo sur Futura]

[invite9e95248d]

un grand merci j'ai testé et c'est exactement ce qu'il me fallait

[invitec66e174c]

Bonjour

Quelqu'un saurait t'il comment trouver le minimun d'une fonction à valeur complexe (f->f(x); x réel et f(x) complexe.) avec matlab. un exemple svp

[invite3e5ede0a]

Bonjour

Quelqu'un saurait t'il comment trouver le minimun d'une fonction à valeur complexe (f->f(x); x réel et f(x) complexe.) avec matlab. un exemple svp

et qu'est ce que le minimum d'un complexe ? le min de la partie réelle ? de la partie imaginaire ? du module ? Par defaut, la commande min de matlab renvoie le minimum de la partie réelle.

[invitec66e174c]

en effet je cherche le x réelle qui minimise la fonction f complexe(f(x)=Re(x) + i*Im(x) )
Apparement on peut utiliser fminsearch de matlab mais je ne sais pas comment faire pour les fonction cmplxe

[invite551c2897]

Bonjour.

le x réelle qui minimise la fonction f

Si tu as les valeurs de f :
[y,i]=min(f)
Tu as un forum dédié à Matlab ici :
http://www.developpez.net/forums/

[invite3e5ede0a]

en effet je cherche le x réelle qui minimise la fonction f complexe(f(x)=Re(x) + i*Im(x) )

Comme je l'ai dis, ça n'a pas de sens de chercher le minimum d'une fonction complexe ! Par contre, chercher le min de la partie réelle, de la partie imaginaire, du module, etc. cela a du sens !

[invite551c2897]

Tu peux faire aussi :
Si tu as les valeurs de f :
[valeur,rang]=min(real(f))
ou
Si tu as les valeurs de f :
[valeur,rang]=min(imag(f))

[invitec66e174c]

Merci beaucoup

[invitec66e174c]

et si jamais le minimum( x qui minimise) la partie reelle est differente de celle qui minimise la partie imaginaire le quel prendre. Faut il prendre le x qui minimise le module?

[invite3e5ede0a]

et si jamais le minimum( x qui minimise) la partie reelle est differente de celle qui minimise la partie imaginaire le quel prendre. Faut il prendre le x qui minimise le module?

tout dépend de l'application. Il n'y a pas de "règle". Que cherches-tu à faire ?

[invitec66e174c]

En effet je dois resoudre le problème suivant avec la methode du recuit simulé j'ai dc besoins d'une solution de depart.
Probleme : Approximation par une fonction de type rationnel
il s'agit d'approcher une fonction S:w->S(w)
R->C^n*n (matrice complexe)
voici l'ennoncé


Trouver :
1. deux entiers p et r
2. p réels a0, a1,...,a(p-1)
3. r+1 matrices B0, B1,....,Br € R^n*n (matrice reelles)
qui minimise
∫||S(w)-(R(iw).R^* (iw))/|P(iw)|^2 ||dw
avec
S:w→S(w);w reel et S(w)est une matrice complexe .
S est donnée ,c^' estla fonction à approcher
R(z)=B_0+B_1 z+⋯+B_(r ) z^(r );z∊C et Bj∊R^(n×n)
P(z)=a_0 +a_1 z+⋯+ a_(p-1) z^(p-1) + z^p ;z∊C et aj∊R

Sous les contraintes :
. r<p
. toutes les racines de p(z) sont dans le demi plan gauche (P(z)=0 alors Re(z)<0)

ps:R^* c'est la transposé conjugué de R