Travaux Pratique N°2
Visualisation des champs électrostatiques pour une charge ponctuelle
Code Matlab :
clc;
close all;
clear all;
% Utiliser la commande « meshgrid » pour créer un maillage
dans le plan
% de dimension -2, 2 avec un pas de 0.3
[x, y] = meshgrid(-2:0.3:2);
% Placer une charge ponctuelle à l'origine du maillage dans
le plan cartésien
figure;
plot(0, 0, 'c*'); % Affichage de la charge ponctuelle à
l'origine
xlabel('Ex');
ylabel('Ey');
title('Visualisation des champs électrostatiques pour une
charge ponctuelle');
% Calcul du champ électrique pour une charge positive
E = 1 ./ sqrt(x.^2 + y.^2); % Module du champ électrique
Ex = E .* cos(atan2(y, x)); % Composante Ex du champ
Ey = E .* sin(atan2(y, x)); % Composante Ey du champ
% Utilisation de la commande « quiver » pour tracer les
vecteurs du champ E
hold on;
quiver(x, y, Ex, Ey, 'r'); % Tracé des vecteurs de champ en
rouge
% Utilisation de la commande « streamslice » pour tracer
les lignes de champ
streamslice(x, y, Ex, Ey);
1
clc;
close all;
clear all;
Nettoyage de l'environnement :
clc efface la fenêtre de commande.
close all ferme toutes les figures ouvertes.
clear all supprime toutes les variables en mémoire pour commencer avec un environnement vide.
[x, y] = meshgrid(-2:0.3:2);
Création du maillage dans le plan cartésien :
• meshgrid génère une grille 2D dans le plan xxx-yyy avec des valeurs de -2 à 2, et un pas
de 0.3 entre chaque point.
• Cela crée un maillage régulier pour les calculs futurs du champ électrique en chaque
point.
figure;
plot(0, 0, 'c*');
xlabel('Ex');
ylabel('Ey');
2
title('Visualisation des champs électrostatiques pour une charge ponctuelle');
Affichage de la charge ponctuelle :
figure ouvre une nouvelle fenêtre de figure.
plot(0, 0, 'c*') place une charge ponctuelle à l'origine du système de coordonnées, représentée par
une étoile cyan (c*) en position (0,0).
Les commandes xlabel, ylabel et title ajoutent des étiquettes aux axes et un titre à la figure pour
faciliter la compréhension du graphique.
E = 1 ./ sqrt(x.^2 + y.^2); % Module du champ électrique
Ex = E .* cos(atan2(y, x)); % Composante Ex du champ
Ey = E .* sin(atan2(y, x)); % Composante Ey du champ
Calcule l'intensité (module) du champ électrique E en chaque point du maillage en utilisant la
formule E=1/r , r= racine (x2+y2). Cela représente le champ d'une charge ponctuelle située à
l'origine.
Ex et Ey représentent les composantes du champ électrique dans les directions x et yatan2(y, x)
calcule l'angle polaire de chaque point par rapport à l'origine.
cos(atan2(y, x)) et sin(atan2(y, x)) permettent de décomposer le champ en ses composantes Ex et
Ey.
hold on;
quiver(x, y, Ex, Ey, 'r');
racé des vecteurs du champ électrique :
• hold on permet de superposer plusieurs tracés dans la même figure.
• quiver(x, y, Ex, Ey, 'r') trace les vecteurs du champ électrique E⃗ en chaque point du
maillage, en utilisant les composantes Ex et Ey. Les vecteurs sont tracés en rouge ('r').
streamslice(x, y, Ex, Ey);
Tracé des lignes de champ :
streamslice trace les lignes de champ (ou lignes de flux) du champ électrique, montrant la
direction dans laquelle une particule de charge se déplacerait sous l’influence du champ
électrostatique.
