Discussion :

[didilegeek]

Bonjour,

je cherche à déterminer tous les zéros de courbes (c'est à dire les x tels que f(x)=0 ) (zéros d'une courbe, ce n'est pas mathématiques, mais c'est pour souligner le fait que je n'ai pas d'expression de la fonction de cette courbe)

Cependant , je ne suis pas en présence de fonctions analytiques, mais de tableaux de 10 points, j'ai donc quelque chose de discret!

Je voudrais donc trouver les zéros de la courbe que je trace à l'aide de ces 10 points, mais chaque zéros n'est pas forcement exactement l'un des points de mon tableau, alors que je veux le zéro exact...

Par exemple si deux points de mon tableau sont (1, -3) et (2, 3) , le zéro se situe environ en x=1,5..

Si vous avez des idées...

Merci d'avance

[FR119492]

Salut!
Le plus simple est de faire un spline et d'en chercher les zéros.
Jean-Marc Blanc

[tanguy4724]

tout dépend de ce qu'il se passe entre les points donnés par le tableau. Tu les relies par une droite ou fais une interpolation linéaire?

si on considère que tu relies les points de ton tableau par des segments de droites, tu peux simplement faire ça :

je ne sais pas sous quelle forme est ton tableau, mais supposons qu'il soit comme tel :

Code :

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

tab = [x1 x2 x3 x4 x5 x6 x7 x8 x9 x10; y1 y2 y3 y4 y5 y6 y7 y8 y9 y10];

tu peux faire :

Code :

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

1
2
3
4
5
6
7
8
9
10
x0=0;
ind = 1;
for i = 1 : 9
    if tab(2,i+1)*tab(2,i)<0 %on regarde si la courbe entre les points i et i+1 coupe l'axe des abcisses
        a = ((tab(2,i+1) - tab(2,i))/(tab(1,i+1)-tab(1,i))); %calcul du coefficient directeur
        b = tab(2,i)-a*tab(1,i); %calcul de l'ordonnée à l'origine
        x0(ind)=-b/a;
        ind=ind+1;
    end
end

x0 te renvoie la liste des x pour lesquels y = 0.

Si ton tableau est d'une forme différente, tu peux adapter facilement pour lire les abscisses et ordonnées.

[didilegeek]

merci du coup de main,

à la base j'avais fait une espèce de dichotmoie, mais le signe du produit, c'est nettement plus efficace !

bonne journée

[tanguy4724]

ravi d'avoir pu t'être utile!