Discussion :

[VON PUTT KAMER]

Bonjour,
je souhaite écrire un programme sur MATLAB, qui décrit une expérience de laboratoire.
en fait, supposons que j'ai une fonction f(t)=v*cos(w*t)+exp(-t) sur un intervalle t=[0,10]s avec v=non nul par exemple, je voudrais la représenter dans cette intervalle mais qu’au bout d'un temps t=5s, v=0 (j'annule l'amplitude).
est ce possible d'obtenir une syntaxe me permettant d'écrire ce type de programme de telle sorte que j'ai une courbe traduisant cette expérience.?
en fait v est nul pour t>5s
Merci..

[phryte]

Bonjour,

Un essai :

Code :

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

1
2
3
4
5
6
t=0:0.1:10;
w=1;
v=3*ones(1,length(t));
v=v.*(t<6);
f=v.*cos(w*t)+exp(-t);
plot(t,f)

[VON PUTT KAMER]

Bonjour,

Un essai :

Code :

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

1
2
3
4
5
6
t=0:0.1:10;
w=1;
v=3*ones(1,length(t));
v=v.*(t<6);
f=v.*cos(w*t)+exp(-t);
plot(t,f)

Merci ,Bonjour
Pourquoi lorsque v est nul la fonction n'a pas l'allure de exp(-t)?

[Antoane]

Bonjour,

le code fonctionne, mais lorsque v= (i.e. pour t>6), l'exponentielle est d'amplitude négligeable par rapport au cos(w*t) (pour t<6). En effet, pour t<6, le cosinus est d'amplitude proche de v, soit environ 3 alors que pour t>6, l'exponentielle est d'amplitude max: 1/exp(6) ~ 1/(3^6) ~0.001
La fonction a donc l'air identiquement nulle pour t>6.

En changeant le paramètres, on rend les différentes contributions plus visibles :

Code :

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

1
2
3
4
5
6
t=0:0.01:10;
w=10;
v=0.1*ones(1,length(t));
v=v.*(t<6);
f=v.*cos(w*t)+exp(-0.2*t);
plot(t,f)

[VON PUTT KAMER]

en fait je revient toujours sur cette courbe
j'ai ce programme ci-dessus

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
rho=2600;
f = 10000;
z = 7e-3;
y = 20;
M = 10e9;
V = 1e-5;
b = 0;
d = 0;
g0 =0;
% % % %A la matrice contenant les éléments du coefficients de fourier an % % % % % % % 
% % % % B la matrice contenant les éléments du coefficients de fourier bn % % % % % % %
t=linspace(0,9e-3,20000);
def=V*sin(2*pi*f*t);
figure(1)
plot(t,def)
% % % % % % % % % % % % % % % % % % % % % % % % % 
A=[0.5-3*d*V^2/16 0 -0.25+d*V^2/8 0 -d*V^2/32];
B=[0 -0.25*b*V 0 b*V/12 0];
% % % % % % % % % % % % % % % % % 
w= 2*pi*f;
zy=z/y;
Ev=M*V*V/y;
N=w*zy;  
    r2=(A(2)-N*2*B(2))/(1 + (2*w*zy)^2);
    r3=(A(3)-N*3*B(3))/(1 + (3*w*zy)^2);
    r4=(A(4)-N*4*B(4))/(1 + (4*w*zy)^2);
    r5=(A(5)-N*5*B(5))/(1 + (5*w*zy)^2);
% % % % % % % % % % % % % % % % % % % % % % % %     
    p2=(B(2)+N*2*A(2))/(1 + (2*w*zy)^2);
    p3=(B(3)+N*3*A(3))/(1 + (3*w*zy)^2);
    p4=(B(4)+N*4*A(4))/(1 + (4*w*zy)^2);
    p5=(B(5)+N*5*A(5))/(1 + (5*w*zy)^2);
% % % % % % %mise sur pied de la fonction de ramollissement % % % % % % % % % % % % % % % % %  
g=(g0 - Ev*(A(1)/2 + r2+r3+r4+r5 ))*exp(-t/zy) + Ev*(A(1)/2 +r2*cos(2*w*t) + p2*sin(2*w*t)+ r3*cos(3*w*t) + p3*sin(3*w*t)+ r4*cos(4*w*t)+ p4*sin(4*w*t)+r5*cos(5*w*t)+ p5*sin(5*w*t) );
figure (2)
plot(t,(g0 - Ev*(A(1)/2 + r2+r3+r4+r5 ))*exp(-t/zy))
% % % % % % % expression de la variation relative de module élastique % % % % % % % % % % % % 
deltaM=((1-g).*(1-2*b*def-3*d*def.^2)-1).*(t<=5)+(((g0 - Ev*(A(1)/2 + r2+r3+r4+r5 ))*exp(-t/zy)).*(t>=5));
figure (3)
plot(t,deltaM)

mais la seconde parti n'apparait pas (la partie croissante) mais si je les dessine à part ca marche (les deux courbures apparaissent bien)

Code :

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

1
2
3
4
figure (3)
hold on
plot(t,deltaM)
plot(t,(g0 - Ev*(A(1)/2 + r2+r3+r4+r5 ))*exp(-t/zy))

I need help

[VON PUTT KAMER]

bonjour la communauté./
je voudrais un programme de propagation d'onde acoustique dans un milieu fissuré svp ou alors un guide ou une méthode pour tenir compte des fissures dans la résolution de l'équation d'onde.