Discussion :

[titi_du_34]

Bonjour à tous et à toutes
pour commencer voici mon code :

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
clear all;
close all;
 
a=4;
b=3;
 
 
omega=0:2000;
t=10;
p0=4;
rho=2530;
h=0.002;
 
x1=a/2;
x2=b/2;
 
 
 
for m=1:4;
    for n=1:4;
 
        wmn=sin((n*pi*x1)/a).*sin((m*pi*x2)/b);
        A=((4*a*b*p0)./(m*n*pi.^2*rho*h));
        omegamn = 2*pi*(sin((n*pi*x1)/a).*sin((m*pi*x2)/b));
        w=((A*exp(1j*omega*t))./(-omega.^2+omegamn.^2))*wmn;
    end
end
 
 
 
figure(1)
plot(omega,20*log(abs(w)))
 
figure(2)
plot(omega,unwrap(angle(w)))

voila la question j'aimerai tracer le module et l'argument de w ,
j'obtiens des courbes qui ne correspondent pas à ce que j'attends ( je souhaite visualiser mes multiples résonances en fonction des modes que j'attribue.
Est ce qu'une personne a une idée sur ce qui peut clocher dans mon programme ?

merci d'avance

[Invité]

Bonjour,

Le problème est que tu écrases ton w à chaque itération.
Donc soit tu effectues le tracé à l'intérieur de ta boucle, soit tu stockes les résultats de chaque itération et traces le tout à la fin.