Discussion :

[bendesarts]

Bonjour,

J'aimerais tracer une fonction qui est définie par plusieurs morceaux.

Voilà le code que j'ai réalisé mais j'ai un problème de dimensions.

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
x=[-10:0.01:20];
Eq=230769;
N=414;
L=20;
 
% pression entre -10 et -8 mm
 
x1=[-10:0.01:-7.99];
R1_1=36*sin(20*pi/180)+x1;
R2_1=83*sin(20*pi/180)-x1;
p = 0.3979*sqrt(Eq*N/(2*L)*(1./R1_1+1./R2_1));
 
% pression entre -8 et 18 mm
 
x2=[-8:0.01:18];
R1_2=36*sin(20*pi/180)+x2;
R2_2=83*sin(20*pi/180)-x2;
p = 0.3979*sqrt(Eq*N/(L)*(1./R1_2+1./R2_2));
 
% pression entre 18 et 20 mm
 
x3=[18.01:0.01:20];
R1_3=36*sin(20*pi/180)+x3;
R2_3=83*sin(20*pi/180)-x3;
p = 0.3979*sqrt(Eq*N/(2*L)*(1./R1_3+1./R2_3));
 
figure(1);hold on
title('Expression de p');
plot(x, p, 'r');

Ma méthode est surement loin d'être optimal. Est-ce vous auriez des idées pour mieux faire?

Merci d'avance pour votre aide

[phryte]

Bonjour.
Une idée :

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
clear
x=[-10:0.01:20];
Eq=230769;
N=414;
L=20;
 
% pression entre -10 et -8 mm
 
x1=[-10:0.01:-7.99];
R1_1=36*sin(20*pi/180)+x1;
R2_1=83*sin(20*pi/180)-x1;
p1 = 0.3979*sqrt(Eq*N/(2*L)*(1./R1_1+1./R2_1));
figure(1)
plot(x1,p1)
hold on
% pression entre -8 et 18 mm
 
x2=[-8:0.01:18];
R1_2=36*sin(20*pi/180)+x2;
R2_2=83*sin(20*pi/180)-x2;
p2 = 0.3979*sqrt(Eq*N/(L)*(1./R1_2+1./R2_2));
plot(x2,p2)
% pression entre 18 et 20 mm
 
x3=[18.01:0.01:20];
R1_3=36*sin(20*pi/180)+x3;
R2_3=83*sin(20*pi/180)-x3;
p3 = 0.3979*sqrt(Eq*N/(2*L)*(1./R1_3+1./R2_3));
plot(x3,p3)
title('Expression de p');
grid

[Jerome Briot]

Une version plus courte... avec une légère tricherie

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
x = -10:0.01:20;
Eq = 230769;
N = 414;
L = 20;
 
R1 = 36*sin(20*pi/180)+x;
R2 = 83*sin(20*pi/180)-x;
 
M = L*(~(x>=-8 & x<=18)+1);
M(x==-8 | x==18) = NaN; % On triche un peu
 
p = 0.3979*sqrt(Eq*N./M.*(1./R1+1./R2));
 
figure
plot(x,p)
grid on
title('Expression de p');