Discussion :

[Hayabusa]

Bonjour,

Tout d'abord, désolé si ce titre de veut rien dire (un polynôme de degré négatif n'est pas un polynôme je pense), mais bon soit, ...

J'ai un vecteur ligne E, de 101 colonnes (ordonnées) et un vecteur ligne m=0:0.1:10. Je voudrais montrer que l'énergie est liée à la masse selon règle en 1/m^{1/3} ("1 sur racine cubique de la masse")

Le problème, c'est que la fonction polyfit ne marche que pour les degrés entiers.
Donc je ne sais pas comment faire ? J'ai usé d'une astuce en inversant la relation, puis en inversant les données, mais ... ça ne marche pas

Merci d'avance pour d'éventuelles pistes de réflexion =)

[Invité]

Bonsoir,

Petite remarque avant de commencer: pour m = 0, 1/m^(1/3) =
Ensuite, cela dépend si une constante ko (de la forme ko + k1/m^(1/3) ) intervient, si c'est le cas, cela ne marchera pas:

Code :

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

1
2
3
4
m = 0.1:0.1:10;
E = 2.3./(m.^(1/3)); % coef à évaluer: 2.3
y1 = E.*(m.^(1/3)); % on ramène la puissance à 0
p = polyfit(m,y1,0)

[Hayabusa]

merci pour ta réponse !

Oui en fait j'ai menti, E a 91 colonnes et m=1:0.1:10 ;-)
normalement il n'y a pas de constante ko ;-)

Je voudrais bien essayer ton code ... mais en fait je ne le comprends pas

sinon une tierce personne vient de me montrer la fonction "fit" dans l'aide matlab. Ca peut pas m'aider ça ?

[Invité]

Tu peux en effet aussi utiliser la fonction fit pour résoudre ton problème.
Quelle partie de mon code ne comprends-tu pas?

[Hayabusa]

Je suis en train de lire la doc de fit, ça à l'air puissant mais chaud :-/

Je ne comprends pas pourquoi tu multiplies par 2 et 3. A moins que ça soit "2 virgule 3", qui serait une valeur de k_1 à tester c'est ça ?

Et puis surtout si tu fais

Code :

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

1
2
3
4
5
6
E = 2.3./(m.^(1/3));
y1 = E.*(m.^(1/3));
ça revient quand même au meme de faire ça :
y1 = 2.3./(m.^(1/3)).*(m.^(1/3));
c'est à dire
y1 = 2.3

car les m.^(1/3) se simplifient ... non ?
et je vois pas le sens de faire un polyfit de degré 0... enfin bref je comprends pas =D

[Invité]

A moins que ça soit "2 virgule 3", qui serait une valeur de k_1 à tester c'est ça ?
car les m.^(1/3) se simplifient ... non ?

Oui c'est bien 2 "virgule" 3, et c'est bien pour simplifier

et je vois pas le sens de faire un polyfit de degré 0... enfin bref je comprends pas

Maintenant que tu le dis, c'est pas faux, un simple mean(y1) suffit.

[Hayabusa]

mais alors j'ai juste y1 = 2.3 et p = 2.3 ^^ et je vois pas trop en quoi ça m'aide :p

[Invité]

2.3 c'etait juste pour l'exemple, la ligne E = 2.3./(m.^(1/3)); est à retirer bien entendu.