Discussion :

[NihadS]

Je voudrais dessiner une ligne verticale (À x = 0.4) et un contour, puis trouver les points d'intersection entre eux
Le contour est une matrice (2*N) sous la forme [x1 y1, ...,xN yN]
J'ai trouvé cette fonction "intersections" qui trouve l'intersection entre deux courbes

http://fr.mathworks.com/matlabcentra...-intersections

Qui peut être utilisé comme suit:

Code :

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

1
2
3
4
5
6
7
8
9
t= 0:pi/64:3*pi;
y = sin(t);
y2 = 0.6*sin(t-0.7)+0.1*randn(size(y));
[xout,yout] = intersections(t,y,t,y2,1);
plot(t,y,'linewidth',2)
set(gca,'xlim',[min(t) max(t)],'ylim',[-1.1 1.1])
hold on
plot(t,y2,'g','linewidth',2)
plot(xout,yout,'r.','markersize',18)

Comment puis-je l'utiliser pour mon cas?
Un exemple de contour est joint

[Jerome Briot]

Peux-tu nous fournir les valeurs des coordonnées du contour dans un fichier .mat (mis dans une archive zip ou rar) ?

Une solution consiste à tester le changement de signe de xi-0.4 pour deux points successifs.

[phryte]

Bonjour,
Avec la function "intersections" :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
clear
x = 0:0.1:10;
y = gaussmf(x,[2 5]);
plot(-y,x)
yy=0:length(x)-1;
xx=-0.4*ones(1,length(x));
[xout,yout] = intersections(-y,x,xx,yy,1);
hold on
plot(xout,yout,'.r','MarkerSize',23)
grid
hold on
line([-0.4 -0.4],[0 max(x)],'Color','r')

Tu peux faire aussi comme te l'as dit Jérome.

[NihadS]

Bonjour,

Merci Jerome et phryte

Jerome,je vous donne le code complet:

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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
clear, close all; clc,
%Mold dimensions and mesh elements
Lx = 0.600; Nx = 60; hx = Lx / (Nx);
Ly = 0.300; Ny = 30; hy = Ly / (Ny);
 
 
%%Loading A text file containing 8 columns: number of the node, node coordinates
%%(x,y, and z respectively), pressure, flow rate, time (8th column) and fill factor in all
%%the nodes
 
B1 = load('case_99_8.txt'); T1 = B1(:,7);
 
 
%Rearranging the vector which contains the time at each node to a matrix
%with the mold dimensions
for j = 1:Ny+1
    for i = 1:Nx+1
        k = (j-1)*(Nx+1) + i;
        TIME1(i,j) = T1(k);
    end
end
TIME1_physical = TIME1'; TIME1_physical = TIME1_physical(end:-1:1,:);
TIME1_physical(TIME1_physical==0)=max(max(TIME1_physical));
 
 
%Meshgrid
x=     0:hx:Lx       ;  y=   Ly:-hy:0;
[X,Y]= meshgrid(x, y);
 
 
%Plotting contours
figure(1)
colormap('jet');
[C, H]=      contour(X,Y,TIME1_physical);        clabel(C,H);
title('Flow Front Contours for a partial Race-tracking scenario'), xlabel('x [m]')
ylabel('y [m]'), grid
hold on
set(gca,'xcolor','k','LineWidth',2,'FontSize',14,'FontWeight','bold') 
set(gca,'ycolor','k','LineWidth',2,'FontSize',14,'FontWeight','bold')
saveas(gcf,'1.pdf')
%Storing in Ni the numbers of the levels of all the curves:
Ni=get(H,'levellist'); %Niveau
NC=Ni(end);% Choix du niveau
 
%Plotting the selected contour level(Putting the selected level in Ni in the form [x x]:
figure(2)
[C1, H1]=contour(X,Y,TIME1_physical,[400 400]);     %clabel(C1,H1);%%%%%%%%%%%%%%%%%%%%%%% Choix du niveau 500
title('SCENARIO 12 FLOW FRONT CURVE'), xlabel('x [m]')
ylabel('y [m]'), grid
set(gca,'xcolor','k','LineWidth',2,'FontSize',14,'FontWeight','bold') 
set(gca,'ycolor','k','LineWidth',2,'FontSize',14,'FontWeight','bold')
saveas(gcf,'2.pdf')
% C and C1 contain erroneous values.To calculate the distance they must be eliminated.
[I,J,V] = find(C1>1);% Searching for erroneous values
if length(J)>0
    deb=J(end)+1;% Beginning of calculations
else
    deb=2;
end
 
grid
set(H1,'Color','k')
axis([0 Lx 0 Ly])
 
% Plotting the discretization of the chosen contour
figure(3)
plot(C1(1,deb:end),C1(2,deb:end),'.r');% Plotting the selected level points
title('Discretization of the chosen contour'), xlabel('x [m]')
ylabel('y [m]'), grid
hold on
plot(0.6,0.15,'.k','MarkerSize',23)% The vent node
grid
axis([0 Lx 0 Ly])
hold off
set(gca,'xcolor','k','LineWidth',2,'FontSize',14,'FontWeight','bold') 
set(gca,'ycolor','k','LineWidth',2,'FontSize',14,'FontWeight','bold')
saveas(gcf,'3.pdf')
% Preparation of the indices to calculate the distances:
i=deb:length(C1);% Indice

Les coordonnées du contour sont dans C1, le fichier case_99_8.txt est attaché

[NihadS]

Bonjour phryte,

C'est le meme code de la dernier fois, alors est-ce que y est le contour C1?
Si oui, s'il vous plait expliquer moi comment je peut l'utliser, C1 est une matrice 2*N

[phryte]

Bonjour,

Avec C1 :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
clear
x = 0:0.1:10;
y = gaussmf(x,[2 5]);
plot(-y,x)
yy=0:length(x)-1;
xx=-0.4*ones(1,length(x));
C1(:,1)=-y';
C1(:,2)=x';
[xout,yout] = intersections(C1(:,1),C1(:,2),xx,yy,1);
hold on
plot(xout,yout,'.r','MarkerSize',23)
grid
hold on
line([-0.4 -0.4],[0 max(x)],'Color','r')

[NihadS]

Il m'a donne cet message:

Subscripted assignment dimension mismatch.
dans la ligne:

Code :

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

C1(:,1)=-y';

J'ai seulement change le x:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
x = 0:0.01:0.6;
% y = gaussmf(x,[2 5]);
% plot(-y,x)
yy=0:length(x)-1;
xx=-0.4*ones(1,length(x));
C1(:,1)=-y';
C1(:,2)=x';
[xout,yout] = intersections(C1(:,1),C1(:,2),xx,yy,1);
hold on
plot(xout,yout,'.r','MarkerSize',23)
grid
hold on
line([-0.4 -0.4],[0 Lx],'Color','r')

Bonjour,

Avec C1 :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
clear
x = 0:0.1:10;
y = gaussmf(x,[2 5]);
plot(-y,x)
yy=0:length(x)-1;
xx=-0.4*ones(1,length(x));
C1(:,1)=-y';
C1(:,2)=x';
[xout,yout] = intersections(C1(:,1),C1(:,2),xx,yy,1);
hold on
plot(xout,yout,'.r','MarkerSize',23)
grid
hold on
line([-0.4 -0.4],[0 max(x)],'Color','r')

[phryte]

Tu as supprimé y.
Donc erreur en ligne 6.

[NihadS]

Ok, L y sert a quoi? J'ai eu l'erreur meme avec l y
ca c'est que j'obtienne

[Jerome Briot]

Une solution consiste à tester le changement de signe de xi-0.4 pour deux points successifs.

Pour tester cette solution, ajoute le code suivant à la fin du tien :

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
x = C1(1,deb:end);
y = C1(2,deb:end);
 
xi = 0.4;
 
s = sign(x-xi);
idx(1) = strfind(s, [1 -1]);
idx(2) = strfind(s, [-1 1]);
 
yi(1) = interp1(x(idx(1):idx(1)+1), y(idx(1):idx(1)+1), xi);
yi(2) = interp1(x(idx(2):idx(2)+1), y(idx(2):idx(2)+1), xi);
 
figure(4)
plot(x, y, 'r*-', [xi xi], yi, 'gs', [xi xi], [min(y) max(y)], 'b:')
axis([0 Lx 0 Ly])

[NihadS]

Merci,

Est-ce que le code est specifiquement pour ce contour et la ligne a 0.4? parce que quand j'ai change la distance il ne marchait pas, egalement quand j'ai change le fichier (le contour C1)

Pour tester cette solution, ajoute le code suivant à la fin du tien :

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
x = C1(1,deb:end);
y = C1(2,deb:end);
 
xi = 0.4;
 
s = sign(x-xi);
idx(1) = strfind(s, [1 -1]);
idx(2) = strfind(s, [-1 1]);
 
yi(1) = interp1(x(idx(1):idx(1)+1), y(idx(1):idx(1)+1), xi);
yi(2) = interp1(x(idx(2):idx(2)+1), y(idx(2):idx(2)+1), xi);
 
figure(4)
plot(x, y, 'r*-', [xi xi], yi, 'gs', [xi xi], [min(y) max(y)], 'b:')
axis([0 Lx 0 Ly])