Discussion :

[doinier]

salut..

je voudrais svp que vous m'aidiez un résoudre ce problème: c'est une équation d’ordre 16 que je voudrais tracer

chaque solution après une resolution numérique.

A(xp)X^16+B(xp)X^14+B(xp)X^12+C(xp)X^11+D(xp)X^1O+E(xp)X^9+F(xp)X^8+G(xp)X^7+H(xp)X^6+...+P(xp)

pour "xp" allant de 0 à 1.5 c-a-d xp=[0:0.01:1.5]

je voudrais tracer chaque solution individuellement dans cette intervalle de "xp"

c-a-dire X1,X2,...........X16.

A(xp),B(xp)...... veux dire que les coefficients sont en fonction de "xp".

Merci de m'aider SVP.

[le fab]

salut

permier problème : y a t'il une seule solution par xp ? probablement non
sinon tu peux utiliser des fonctions comme fsolve ou fzero qui te trouveront une solutions
la méthode est d'encapsuler dans une boucle for sur xp cette résolution et d'associer la(les) solution(s) à xp via une matrice / un tableau de cellules puis de tracer le résultat via plot

[doinier]

salut à tous

je voudrais S.V.P que vous m'aidiez un résoudre se problème: c'est une équation d’ordre 16 que je voudrais tracer

chaque solution après une résolution numérique.

A(xp)X^16+B(xp)X^14+B(xp)X^12+C(xp)X^11+D(xp)X^1O+E(xp)X^9+F(xp)X^8+G(xp)X^7+H(xp)X^6+...+P(xp)

pour "xp" allant de 0 à 1.5 c-a-d xp=[0:0.01:1.5]

je voudrais tracer chaque solution individuellement dans cette intervalle de "xp"

c-a-dire X1,X2,...........X16.

A(xp),B(xp)...... veux dire que les coefficients sont en fonction de "xp".

Merci de m'aider SVP.

voici l'explication sous forme de code
problème: ne trace pas les solutions individuellement en prenant différentes valeur de n "c-a-d

n=1,2,3,4.......16

"

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 i=1
for xp=0:0.001:1.5
    A=......
    B=......
    C=......
    D=......
    E=......
    ...etc
 
    polynom=[A 0 B C D ......]                           %écriture du polynôme
 
    racine=roots(polynom)                                 %différentes racines 
 
    for n=1:numel(racine)
       OMER2(i,1)=kp ;                                     % range les éléments de xp  
       OMER3(i,1)=(racine(n));                           % range la partie réelle de ONE
         i=i+1;
    end
end
plot(OMER2(:,1),OMER3(:,1),'k.')

[le fab]

salut

solution rapide, tu encadre ton plot d'une boucle for sur i pour tracer toutes tes courbes
pour ne pas écraser la courbe precedente à chaque iteration, tu utilise la fonction hold

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
figure
hold on
for i=1:16
ton code sans la première ligne
end