Discussion :

[louftansa]

Bonjour,

Mon problème est le suivant:

J'ai une liste de points qui définissent ma grille de coordonnées ainsi que l'image de ces derniers. Je me retrouve ainsi avec une matrice à 3 colonnes (X,Y,Z). Mon problème est que mes coordonnées ne sont pas réguliérement espacées, ce qui fait que je n'arrive pas à tracer la courbe. Quelqu'un peut il m'aider??? . Ma matrice est la suivante:

A =

8,00E+01 1,40E+02 6,61145E-07
8,00E+01 1,56E+02 3,42572E-06
9,00E+01 1,40E+02 1,88792E-05
9,00E+01 1,50E+02 0,09923999
9,00E+01 1,50E+02 0,099508324
9,00E+01 1,50E+02 0,09923999
9,00E+01 1,56E+02 0,000578976
9,25E+01 1,50E+02 0,168913024
9,25E+01 1,50E+02 0,168897025
9,25E+01 1,50E+02 0,168913024
9,50E+01 1,50E+02 0,191911451
9,50E+01 1,50E+02 0,192461431
9,50E+01 1,50E+02 0,191911451
9,75E+01 1,50E+02 0,193273939
9,75E+01 1,50E+02 0,193436723
9,75E+01 1,50E+02 0,193273939
1,00E+02 1,40E+02 3,58739E-05
1,00E+02 1,50E+02 0,192378185
1,00E+02 1,50E+02 0,192930907
1,00E+02 1,50E+02 0,192378185
1,00E+02 1,56E+02 0,00113179
1,03E+02 1,50E+02 0,193273939
1,03E+02 1,50E+02 0,193436723
1,03E+02 1,50E+02 0,193273939
1,05E+02 1,50E+02 0,191911451
1,05E+02 1,50E+02 0,192461431
1,05E+02 1,50E+02 0,191911451
1,08E+02 1,50E+02 0,168913024
1,08E+02 1,50E+02 0,168897025
1,08E+02 1,50E+02 0,168913024
1,10E+02 1,40E+02 1,88792E-05
1,10E+02 1,50E+02 0,09923999
1,10E+02 1,50E+02 0,099508324
1,10E+02 1,50E+02 0,09923999
1,10E+02 1,56E+02 0,000578976

[Jerome Briot]

Mon problème est que mes coordonnées ne sont pas réguliérement espacées, ce qui fait que je n'arrive pas à tracer la courbe.

Peux-tu nous montrer ce que tu as essayé de faire sous MATLAB ?

[louftansa]

J'ai créé un petit programme pour représenter la fonction:

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
47
48
49
50
51
52
53
54
55
56

%ce programme construit une matrice Z(i,j) représentant l'amplitude du
%champ en fonction des coordonnées du capteur
%les lignes représentent les valeurs de x et les colonnes les valeurs de y

format long e


%lecture des valeurs initiales de x et y 
xini=A(1,1);
yini=A(1,2);


%lecture du nombre de lignes de la matrice
n=size (A,1);

%lecture valeurs finales x et y
xfinal=A(n,1);
yfinal=A(n,2);

%lecture du pas de temps

 pasx = input('quel est le pas de temps pour x? ');
 pasy = input('quel est le pas de temps pour y? ');

%initialisation des indices du tableau
i=1;
j=1;

% lecture des différentes valeurs de x
for m=xini:pasx:xfinal
 
  
    %balayage de toules lignes du tableau
    for l=1:n
        if A(l,1)==m
            Z(i,j)=A(l,5);
            j=j+1; 
        end
    end
     j=1;i=i+1;
end

%creation de la grille de définition de la courbe à dessiner
% pensez à l'adapter en fonction des cordonnées qui définissent Z(i,j)
% pour le meshgrid inverser l'axe des x et l'axe des Y

[X,Y]= meshgrid(yini:pasy:yfinal,xini:pasx:xfinal);


%dessin de la courbe 
surf(X,Y,Z);
xlabel ('coupe selon la largeur du défaut');
ylabel('coupe selon la longueur du défaut');
zlabel ('Amplitude');

[mr_samurai]

Salut,

ton problème est un peu délicat. J'ai joins à mon poste une fonction qui réalise une régression plus élaborée sur un ensemble de points ( pas nécessairement alignés sur une grille).

Pour l'utiliser:

Code :

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

regression(A)

avec A la matrice de ton premier poste.

PS :
tu utilises des ',' pour les nombres décimaux alors que c'est des '.' sous Matlab.

++ bonne chance.

[louftansa]

Merci beaucoup, ça marche!!!