Discussion :

[Brian_]

Bonsoir,

je me suis inscrit sur ce forum car j'aimerais représenter une fonction sur matlab mais je coince un peu!

J'aimerais représenter la fonction suivante qui est la composante longitudinale d'un champ électrique dans un guide d'onde :

Ez = Asin(Bx/b)cos(wt-kz) où A,B,b,w,k sont des constantes connues; il n'y a pas de limites sur t ou z mais en revanche x ne peut varier qu'entre 0 et b.

Le problème c'est que je ne sais pas par où commencer, si je dois utiliser surf,plot,meshgrid etc.
Je voulais dans un premier temps représenter cette fonction pour z fixé et après j'aurais aimé la représenter pour tout z compris entre 0 et une valeur quelconque!

Quelqu'un pourrait me donner un coup de main, me donner au moins les outils que je dois utiliser ?

Merci beaucoup,

bonne soirée

[FLB]

Je voulais dans un premier temps représenter cette fonction pour z fixé

Donc t varie, et x aussi? Tu veux donc représenter un plan?

pour tout z

Représenter les variations de 3 variables va être compliqué!

Les méthodes utiles : meshgrid et surf.
Un petit tour dans la faq est toujours utile, et si tu cherches de l'inspiration sur le rendu que tu veux avoir, regardes le portfolio de Dut.

[Brian_]

Donc t varie, et x aussi? Tu veux donc représenter un plan?

Oui c'est ça mais je pense que je vais finalement tracer pour t fixé et avec x,z qui varient.

Je vais regarder un peu comment fonctionne meshgrid parce que j'ai pas trop compris pour l'instant :p Et le lien a l'air génial merci !

Bonne soirée

[FLB]

t fixé et avec x,z qui varient.

Ca semble raisonnable
Le principe du meshgrid est 2 vecteurs->2 grilles de coordonnées, auxquels tu peux appliquer ta fonction.

[Brian_]

Je crois que j'ai eu ce que je voulais :p

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
%Champ longitudinal
A = 1; %amplitude
n = input('Mode souhaité n =');
b = 1; %1 mètre de hauteur
c = 3E8;
wc = n*pi*c/b; %pulsation de coupure
w = 1.5*wc;
k = w/c;
 
%------------------t fixé----------------------------------------
x = [0:0.01:b];
z = [0:0.01:b];
 
[X,Z] = meshgrid(x,z);
t=0;
E = A.*sin(pi*n*X/b).*cos(w*t-k*Z);
figure;
surf(X,Z,E);

Ça me donne bien le champ longitudinal en fonction de x et z à t=0 c'est génial merci!

Après j'ai voulu tracer le champ à t=0 et au milieu du guide d'onde donc pour x =0.5*b puis je me suis inspiré d'une correction que j'ai trouvé sur internet qui permettait de tracer une onde progressive et j'ai fait ça :

Code :

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

1
2
3
4
5
6
figure;
for t=0:0.01:2
    Et = sin(pi*0.5/b).*cos(2*pi*(0.05*t-z));
    plot(z,Et)
    drawnow
end

Le résultat a l'air correct en tout cas ça me convient puisqu'il me faut un champ longitudinal maximal en x = 0.5*b qui suivrait une particule pour l'accélérer. Par contre j'ai du changer la pulsation et k sinon ça allait trop vite et voyait rien :p

Mais j'ai voulu faire pareil avec la surface que j'ai tracée au début :

Code :

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

1
2
3
4
5
6
figure;
for t=0:0.01:10
    Etest = A.*sin(pi*n*X/b).*cos(0.05*t-Z);
    surf(X,Z,Etest);
    drawnow
end

Mais c'est bizarre pendant l'animation les axes bougent, genre au début ils sont entre -1 et 1 puis d'un coup ils passent de 0 à 1 >< Il y a moyen de les fixer avant que je trace ?

Le désavantage de la méthode que j'ai utilisée c'est que je ne peux pas enregistrer l'animation j'ai l'impression ?

Merci !

[FLB]

Pour sauvegarder les animations regarde

Code :

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

1
2
getframe
movie

Pour déterminer tes axes, et non les avoir automatiquement, regarde

Code :

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

1
2
3
axis
ylim
xlim

[Jerome Briot]

Code :

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

1
2
x = [0:0.01:b];
z = [0:0.01:b];

Les crochets sont inutiles :

Code :

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

1
2
x = 0:0.01:b;
z = 0:0.01:b;

[Brian_]

Bonjour,

merci à vous j'ai réussi à avoir ce que je voulais !

Je trouve juste que l'aspect est pas super beau mais j'ai pas réussi à faire mieux :p j'ai essayé interp avec winter etc mais ça donne des choses bizarres

Merci encore,

bon dimanche

[Jerome Briot]

Tu peu essayer avec

Code :

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

1
2
3
s = surf(X,Z,E);
set(s,'edgecolor','none');
shading interp

Montre nous un exemple de ce que tu souhaites obtenir si le code ci-dessus ne te convient pas.

[Brian_]

Bonjour, merci pour votre réponse

En fait par exemple en regardant ceci,j'ai l'impression que le tracé est beaucoup plus "fin" que le mien, comment ça se fait ? J'ai tracé avec un mesh comme eux mais c'est pas aussi net :p

Merci !

[Jerome Briot]

Ces graphiques ont certainement été retouchés avec un autre logiciel que MATLAB