Discussion :

[anthonystaltaro]

Bonjour tout le monde, j'ai un petit probleme.

En fait je dispose d'un ensemble de points correspondant à deux courbes, les valeurs de la courbe1 sont stockées dans le vecteur 1 et les valeurs de la courbe 2 dans le vecteur 2.

Ces deux courbes se croisent si on les traces avec excel....

Je voudrait donc savoir s'il est possible de déterminer les coordonnées des points d'intersections par interpolations entre les points que je donne avec matlab.

merci par avance

[Jerome Briot]

Si les deux courbes sont y1 et y2, tu dois chercher les zéros de y1-y2

Voir cette discussion récente : Trouver les zéros

[chris05]

Cette question semble redondante ...

J'ai donc créé une petite fonction qui interpole 2 vecteurs à partir de fonctions polynomiales d'ordre n et calcule leurs intersections.

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
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
 
function [inter] = intersections (X1,X2,F1,F2,n)
 
 
syms x  %définition de la variable x
 
%initialisation
f1b=0;
f2b=0;
 
 
%interpolation polynomiale
f1=polyfit(X1,F1,n);
f2=polyfit(X2,F2,n);
 
 
%ecriture des fonctions d'interpolation
for i=1:n+1
    f1b=f1b+ f1(i)*x^(n-i+1);
    f2b=f2b+ f2(i)*x^(n-i+1);
end
f_inter=f1b-f2b;%calcul de la différence des fonctions d'interpolation
 
 
%application numerique des fonctions d'interpolation pour le tracé
xb=min(min(X1),min(X2)):(max(max(X1),max(X2))-min(min(X1),min(X2)))/100:max(max(X1),max(X2));
f1b=subs(f1b,x,xb);
f2b=subs(f2b,x,xb);
 
 
%tracé des fonctions
figure()
hold on
plot(X1,F1,'--r')
plot(X2,F2,'--y')
plot(xb,f1b,'r')
plot(xb,f2b,'y')
 
title('tracé des fonctions et de leurs interpolations')
xlabel('x')
ylabel('f')
legend('F1','F2','interpolation de F1','interpolation de F2')
hold off
 
%calcul des intersections
inter=eval(solve(f_inter,x));

F1: vecteur contenant les f1(x)
F2: vecteur contenant les f2(x)
X1 est le vecteur x pour la fonction F1
X2 celui pour F2
n est le degré du polynome d'interpolation.

Voila, en esperant que cela rendra service a quelques personnes.

[chris05]

personne pour me dire si cela fonctionne????