Discussion :

[Nabuchodonosor15]

bonjour,

voila j'ai un probleme sous matlab je veux trouver la valeur en 0. Or cette valeur c'est celle que je passe en parametre de ma fonction.

U(x)=E-Ri et le but c'est de trouver U(x)-E+Ri=0 je n'ai pas solve ni fsolve donc je compare les deux cotés de l'équation en faisant varier "x" mais la précision de matlab n'etant pas bonne je n'y arrive pas.

j'ai trouvé un petit truc pour cela c'est "eps" je fais donc

eps((U(x))=E-Ri
mais ca ne marche pas non plus, qqun aurai t il une solution a ce probleme?

merci bcp

[loique]

Salut

t'as essayé

Code :

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

fzero

?

a+

[christophe_halgand]

Si une égalité ne peut pas fonctionner, il te reste l'intervalle...

Christophe

[FR119492]

Salut!

je veux trouver la valeur en 0

Ton problème est formulé de manière incompréhensible: est-ce en x=0 ou en i=0?

U(x)-E+Ri=0

Qu'est-ce qui est donné: U(x), E, R, i et sous quelle forme? qu'est-ce que tu veux calculer?
Jean-Marc Blanc

[Nabuchodonosor15]

merci pour vos réponses
j'ai réussis a me debrouiller tout seul en utilisant la fonction "roots" car je n'ai pas fsolve ou solve
c'est pas tres pratique mais je me debrouillerai avec.

par contre matlab n'est pas tres pratique coté égalité car si j'utilise la meme formule
U(x) - E+Ri = 0 je connais tout sauf x je le fait donc varier en utilisant la méthode de la dichotomie
comme cela

Code :

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

1
2
3
4
5
for x= 0:0.1:2
        if ( 0 <U(x) - E+Ri < 0.1)    % je mets une fourchette
              break;
       end
end

mais ca ne marche pas tres bien .. certainement que j'ai un probleme quelque part
merci

[Caro-Line]

Tu ne peux pas mettre un test avec 2 comparaisons comme tu le fais.
MATLAB ne vérifie alors que la dernière inégalité.
Pour preuve :

Code :

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

1
2
0<1<0.5 %donne 0 donc faux => OK
1.5<1<2%donne 1 donc vrai => KO

Il faut mettre 2 inégalités avec un && :

Code :

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

if (0 <U(x) - E+Ri) && (U(x) - E+Ri< 0.1)

Sinon loique avait proposé la fonction FZERO : l'as-tu essayé ? Si oui pourquoi ne convient-elle pas ?

[Jerome Briot]

Il faut mettre 2 inégalités avec un && :

Code :

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

if (0 <U(x) - E+Ri) && (U(x) - E+Ri< 0.1)

Je dirais plutôt :

Code :

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

if (U(x) - E+Ri)>0 && (U(x) - E+Ri)< 0.1

[Nabuchodonosor15]

oui merci ca marche tres bien !

pour fzero je ne sais pas, j'ai essayé avec j'avais écris en suivant l'help mais avec ma fonction ca ne fonctionnait pas...



j'ai une autre question si je peux abuser :
voila mon problème

mon code est fait ainsi je cherche une valeur pour laquelle mon équation est nulle pas dichotonomie est je trouve xdeb=0.973

cette valeur a ma seconde valeur donc a xdeb(2)=0.973
xdeb(1) lui est nulle =0
et ensuite je veux calculer les termes suivant :

Code :

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

1
2
3
4
5
6
7
8
9
for k=2:kmax;
 
 
          xdeb(k+1)=(xdeb(k)*2)-xdeb(k-1);
          debut=xdeb(k+1)
          variable1=-(e(k)+r*(ucalcul*(1/u(k))))
          variable2=Tension_permanente(vyt5_linear(debut),vyt20_linear(debut),vyt40_linear(debut),I,Temp,tfloat,age)
          xdeb=variable2+variable1;
end;

et quand je fais ca il me met une erreur :

??? Attempted to access xdeb(3); index out of bounds because numel(xdeb)=1.

Error in ==> batterie_cadmium_nickel at 311
xdeb(k+1)=(xdeb(k)*2)-xdeb(k-1);

apparement il aime pas le troisieme élément et effectivement quand j'affiche xdeb(2) il met qu'il est égale à :
0 0.9730

donc je pense que l'erreur vient de la au lieu d'avoir une valeur pour xdeb(2) il fait une matrice ( logique pour matlab mais moi ca marrange pas )

comment puis je faire ?? s'il vous plait

[Caro-Line]

Non le problème vient de la dernière ligne de ta boucle :

Code :

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

xdeb=variable2+variable1;

En faisant cela tu écrases xdeb (qui contient pour k=1 3 éléments xdeb(1), xdeb(2) et xdeb(3)) par un scalaire (donc tu n'as plus que un xdeb(1) qui ne correspond plus à ce que tu voulais)

[Nabuchodonosor15]

oui je venais de trouver le probleme merci
c'est vraiment une erreur bete ... dsl de t'avoir fait perdre ton temps pour cela