Travaux Pratiques Les Antennes DI

FACULTE POLYDISCIPLINAIRE DE LARACHE
Master : Systèmes Électroniques Et Energies Renouvelables

Module : Antennes et Rayonnement
Réaliser par : Daouairi Ilyass ilyass.daouairi@etu.uae.ac.ma
TP n°1 : Diagramme de rayonnement d’un dipôle rayonnant par MATLAB
1. Objectif de ce TP :
L’objectif est de générer le diagramme de rayonnement d'un dipôle rayonnant d'une longueur
totale de 2L, alimenté en son milieu. Cette tâche sera accomplie en utilisant le logiciel
MATLAB.
Schéma d’un dipôle rayonnant
2. Préparation :
Étant donné que la fonction caractéristique d'un dipôle rayonnant d'une longueur2 L≫ λ
alimenté en son milieu est déterminée par la relation :
|cos ( KLcos ( θ ) ) −cos ⁡(KL)|

F ( θ )=
|( 1−cosKL )∗sin ⁡(θ)|
( 2m+1 )∗λ
 Quand la longueur 2 L du dipôle est égale à (dipôle demi-onde), la fonction
2
F(θ) prend la forme :
( 2m+1 )∗λ
On sait que 2 L=
2
( 2 m+1 )∗λ
Donc L=
4
On remplace avec L par son expression dans F ( θ ):
| (
F ( θ )=
cos
K∗( 2 m+1 )∗λ
4 )
cos ( θ ) −cos ⁡(
K∗( 2 m+1 )∗λ
4
) |
|( 1−cos ⁡(
K∗( 2m+1 )∗λ
4 )
) ∗sin ⁡(θ) |
 Quand la longueur 2 L du dipôle est égale à ((2 m+ 1)∗λ ¿ (dipôle à onde entière), la
fonction F(θ) prend la forme :
On sait que 2 L=((2 m+ 1)∗λ)
( 2 m+1 )∗λ
Donc L=
2
| (
F ( θ )=
cos
K∗( 2 m+1 )∗λ
2 )
K∗( 2 m+1 )∗λ
2
) |
(| 1−cos ⁡( K∗( 2m+1
2
)∗λ
) )∗sin ⁡(θ)|
3. Simulations :
1-a le traçage de diagramme de rayonnement du dipôle rayonnant en coordonnées polaires
puis en coordonnés sphériques pour différentes longueurs :
 Dipôle à onde entière :

+Pour m=0 :
F ( θ )=
| (
cos
Kλ
2 ) Kλ
cos ( θ ) −cos ⁡( )
2 |
|( 1−cos ⁡(
Kλ
2 )
) ∗sin ⁡(θ)
|
En coordonnées polaires :
% Paramètres
K = 2*pi; % Nombre d'onde
lambda = 1; % Longueur d'onde
% Angle theta
theta = linspace(0, pi, 1000);
% Fonction caractéristique F(theta)

F_theta = abs(cos((K*lambda)/2 * cos(theta)) -
cos((K*lambda)/2)) ./ (abs((1 -
cos((K*lambda)/2)) .* sin(theta)) + eps);
% Tracé du diagramme de rayonnement en coordonnées

polaires
polarplot(theta, F_theta, '-b');
title('Diagramme de Rayonnement - Dipôle a onde
entière');
% Ajustements d'affichage facultatifs

grid on;
J'ai ajouté eps dans le dénominateur pour éviter une division par zéro potentiellement
lorsqu`sin(θ) est très proche de zéro.
En coordonnées sphérique :
% Paramètres
% Angles azimutaux et d'élévation
phi = linspace(0, 2*pi, 360);
% Meshgrid pour les coordonnées sphériques
[Phi, Theta] = meshgrid(phi, theta);
F_theta = abs(cos((K*lambda)/2 * cos(Theta)) -
cos((K*lambda)/2)) ./ (abs((1 - cos((K*lambda)/2)) .*
sin(Theta)) + eps);
% Conversion des coordonnées sphériques en coordonnées
cartésiennes
X = F_theta .* sin(Theta) .* cos(Phi);
Y = F_theta .* sin(Theta) .* sin(Phi);
Z = F_theta .* cos(Theta);

sphériques
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Diagramme de Rayonnement - Dipôle a onde entière en
Coordonnées Sphériques');
xlabel('X');
ylabel('Y');
zlabel('Z');
On utilise la fonction surf pour créer une représentation en 3D des coordonnées sphériques.
+Pour m=1 :
| (
F ( θ )=
cos
3 Kλ
2 )
3 Kλ
2
)
|
|( 1−cos ⁡(
3 Kλ
2 ) |
) ∗sin ⁡(θ)
% Paramètres
% Angle theta

F_theta = abs(cos((3*K*lambda)/2 * cos(theta)) -
cos(3*K*lambda/2)) ./ (abs((1 - cos(3*K*lambda/2)) .*
sin(theta)) + eps);

polaires
% Paramètres
F_theta = abs(cos((3*K*lambda)/2 * cos(Theta)) -
cos((3*K*lambda)/2)) ./ (abs((1 - cos((3*K*lambda)/2)) .*
sin(Theta)) + eps);
cartésiennes
% Tracé du diagramme de rayonnement en coordonnées sphériques

figure;
xlabel('X');
ylabel('Y');
zlabel('Z');
+Pour m=2 :
| (
F ( θ )=
cos
5 Kλ
2 )
5 Kλ
2
)
|
|( 1−cos ⁡(
5 Kλ
2 ) |
) ∗sin ⁡(θ)
% Paramètres
% Angle theta

sin(theta)) + eps);

polaires
title('Diagramme de Rayonnement - Dipôle a onde entière
');
% Ajustements d'affichage facultatifs

grid on;
% Paramètres



sin(Theta)) + eps);

cartésiennes

sphériques
figure;
xlabel('X');
ylabel('Y');
zlabel('Z');
+Pour m=3 :
| (
F ( θ )=
cos
7 Kλ
2 )
7 Kλ
2
)
|
|( 1−cos ⁡(
7 Kλ
2 ) |
) ∗sin ⁡(θ)
% Paramètres
% Angle theta

cos(7*K*lambda/2)) ./ (abs((1 - cos(7*K*lambda/2))
.* sin(theta)) + eps);
% Paramètres


sin(Theta)) + eps);

cartésiennes

sphériques
figure;
title('Diagramme de Rayonnement - Dipôle a onde
entière en Coordonnées Sphériques');
xlabel('X');
ylabel('Y');
zlabel('Z');
 Dipôle demi-onde :
+Pour m=0 :
| (
F ( θ )=
cos
Kλ
4 ) Kλ
cos ( θ ) −cos ⁡( )
4 |
|( 1−cos ⁡(
Kλ
4 ) |
) ∗sin ⁡(θ)
% Paramètres
% Angle theta

F_theta = abs(cos((K*lambda)/4 *
cos(theta)) - cos(K*lambda/4)) ./
(abs((1 - cos(K*lambda/4)) .*
% Paramètres



F_theta = abs(cos((K*lambda)/4 * cos(Theta)) -
cos(K*lambda/4)) ./ (abs((1 -
cos(K*lambda/4)) .* sin(Theta)) + eps);
% Conversion des coordonnées sphériques en

coordonnées cartésiennes
% Tracé du diagramme de rayonnement en

coordonnées sphériques
figure;
title('Diagramme de Rayonnement - Dipôle Demi-
Onde en Coordonnées Sphériques');
xlabel('X');
ylabel('Y');
zlabel('Z');
+Pour m=1 :
| (
F ( θ )=
cos
3 Kλ
4 )
3 Kλ
4
)
|
|( 1−cos ⁡(
3 Kλ
4 ) |
) ∗sin ⁡(θ)
% Paramètres
% Angle theta

cos(3*K*lambda/4)) ./
(abs((1 -
cos(3*K*lambda/4)) .* sin(theta)) + eps);
% Paramètres


F_theta = abs(cos((3*K*lambda)/4 * cos(Theta))
- cos(3*K*lambda/4)) ./ (abs((1 -
cos(3*K*lambda/4)) .* sin(Theta)) + eps);


figure;
Onde en Coordonnées Sphériques');
xlabel('X');
ylabel('Y');
zlabel('Z');
+Pour m=2 :
| (
F ( θ )=
cos
5 Kλ
4 )
5 Kλ
4
)
|
|( 1−cos ⁡(
5 Kλ
4 ) |
) ∗sin ⁡(θ)
% Paramètres
lambda = 1; % Longueur d'onde (à titre d'exemple,
vous pouvez ajuster cela)
% Angle theta

cos(5*K*lambda/4)) ./ (abs((1 -
% Paramètres





sphériques
figure;
title('Diagramme de Rayonnement - Dipôle Demi-Onde
en Coordonnées Sphériques');
xlabel('X');
ylabel('Y');
zlabel('Z');
+Pour m=3 :
| (
F ( θ )=
cos
7 Kλ
4 )
7 Kλ
4
)
|
|( 1−cos ⁡(
7 Kλ
4 ) |
) ∗sin ⁡(θ)
% Paramètres
% Angle theta

% Paramètres





sphériques
figure;
title('Diagramme de Rayonnement - Dipôle Demi-Onde
en Coordonnées Sphériques');
xlabel('X');
ylabel('Y');
zlabel('Z');
Remarque : Un dipôle demi-onde aura des lobes principaux plus larges et des lobes
secondaires plus prononcés par rapport à un dipôle onde entière de même longueur.
B- Conclusion :
 Commentaires :
-Les diagrammes en coordonnées polaires fournissent une représentation bidimensionnelle du
rayonnement dans le plan horizontal.
-Les diagrammes en coordonnées sphériques offrent une vue tridimensionnelle du
rayonnement en prenant en compte à la fois l'angle d'élévation et l'angle azimutal.
-La directivité plus élevée est souvent associée à une meilleure performance en termes de
portée et de capacité à transmettre ou recevoir des signaux dans des directions spécifiques.
 Conclusion :
L’analyse des diagrammes de rayonnement en coordonnées polaires et sphériques pour
différentes longueurs de dipôle permet de comprendre les propriétés directionnelles et la
directivité du dipôle rayonnant. Ces informations sont cruciales dans la conception d'antennes
pour des applications spécifiques, où la directionnalité du rayonnement joue un rôle essentiel.
2-a Un dipôle électrique infinitésimal horizontal de courant constant I0 est placé

symétriquement par rapport à l'origine et dirigé le long de l'axe des x.
( ) √ 1−sinθ cosφ
− jkr
k I0 Le 2 π ×120 π × I 0 × L
|E (θ , φ)|= jη
4 πr
√ 1−sinθ2 cosφ2 |E (θ , φ)|=
4 πrλ
2 2
|E (θ , φ)|= ( ηk I 0 L
4 πr ) √1−sinθ2 cosφ2 |E (θ , φ)|=
60 I 0 πL
r λ
√ 2
1−sinθ cosφ
2
60 I 0
|E (θ , φ)|= |F ( θ , φ )|
r
πL
|F (θ , φ )|=
λ
√ 2
1−sinθ cosφ
2
( 2m+1 )∗λ
 Quand la longueur 2 L du dipôle est égale à (dipôle demi-onde), la fonction
2
F(θ) prend la forme :
( 2m+1 )∗λ
On sait que 2 L=
2
( 2 m+1 )∗λ
Donc L=
4
π ( 2 m+ 1 ) λ
|F (θ , φ )|=
λ
×
4
√ 1−sinθ2 cosφ 2
( 2 m+ 1 ) π
|F (θ , φ )|=
4
√ 1−sinθ 2 cosφ2
Pour m=1 :
3π
|F (θ , φ )|=
4
√ 2
1−sinθ cosφ
2
+Pour φ=0 ° :
{cosφ=1
φ=0 ° 3π
|F (θ)|=
4
√ 1−sinθ2 |F (θ)|=
3π
4
cosθ
% Paramètres
theta = linspace(0, 2*pi, 180); % Angle
d'élévation (0 à pi pour une demi-onde)
% Calcul de la fonction caractéristique

F_theta = (3*pi/4) * cos(theta);

coordonnées polaires
figure;
En coordonnées
polarplot(theta, sphérique : '-b');
abs(F_theta),
% Onde');
Paramètres
grid =on;
theta linspace(0, pi, 180); % Angle d'élévation
(0 à pi pour une demi-onde)
phi = linspace(0, 2*pi, 360); % Angle azimutal (0
à 2*pi pour une onde complète)
% Création de la grille sphérique

[Theta, Phi] = meshgrid(theta, phi);

F_spher = (3*pi/4) * cos(Theta);

X = F_spher .* sin(Theta) .* cos(Phi);
Y = F_spher .* sin(Theta) .* sin(Phi);
Z = F_spher .* cos(Theta);

figure;
Onde');
xlabel('X');
ylabel('Y');
zlabel('Z');
axis equal;
+Pour φ=90 ° :
{cosφ=0
φ=0 °
|F (θ)|=
3π
4
√1 |F (θ)|=
3π
4
% Paramètres
theta = linspace(0, 2*pi, 180); % Angle
d'élévation (0 à pi pour une demi-onde)
F_theta = (3*pi/4) * ones(size(theta)); %
Fonction caractéristique constante

figure;
Onde');
grid on;
% Paramètres
theta = linspace(0, pi, 180); % Angle d'élévation


F_spher = (3*pi/4) * ones(size(Theta));

+Pour θ=90
coordonnées ° :
cartésiennes
{sinθ=1
Y = F_spher
|
Z = F_spher .* cos(Theta); |
F (φ) =
3π
θ=90°.* sin(Theta) .* sin(Phi);
4
√
1−cosφ
2
3π
| |
F (φ) =
coordonnées sinφ
sphériques
figure; 4
Onde');
xlabel('X');
ylabel('Y');
zlabel('Z');
axis equal;
+Pour θ=90 ° :
{sinθ=1
θ=90°
|F (φ)|=
3π
4
√ 1−sinφ
2
|F (φ)|=
3π
4
cosφ
% Paramètres
phi = linspace(0, 2*pi, 360); % Angle azimutal
(0 à 2*pi pour une onde complète)
F_phi = (3*pi/4) * sin(phi);

figure;
polarplot(phi, abs(F_phi), '-b');
Onde');
grid on;
% Paramètres
theta = linspace(0, pi, 180); % Angle d'élévation


F_spher = (3*pi/4) * cos(Phi);

Y = F_spher .* sin(Theta) .* sin(Phi);
Z = F_spher .* cos(Theta);

figure;
Onde');
xlabel('X');
ylabel('Y');
zlabel('Z');
axis equal;
-b Un dipôle électrique infinitésimal horizontal de courant constant I0 est placé

symétriquement par rapport à l'origine et dirigé le long de l'axe des y.
( 2 π ×1204 πrλπ × I × L ) √ 1−sinθ sin φ

− jkr
k I0 Le
|E (θ , φ)|= jη
4 πr
√ 1−sinθ2 sin φ2 |E (θ , φ)|= 0 2 2
|E (θ , φ)|= ( ηk4 IπrL ) √1−sinθ sin φ

0 2 2
|E (θ , φ)|=
60 I 0 πL
r λ
√ 1−sinθ2 sin φ 2
60 I 0
|E (θ , φ)|= |F ( θ , φ )|
r
πL
|F (θ , φ )|=
λ
√ 1−sinθ 2 s∈φ2
 Quand la longueur 2 L du dipôle est égale à ((2 m+ 1)∗λ ¿ (dipôle à onde entière), la
fonction F(θ) prend la forme :
On sait que 2 L=((2 m+ 1)∗λ)
( 2 m+1 )∗λ
Donc L=
2
π ( 2 m+ 1 )∗λ
|F (θ , φ )|=
λ
×
2
√ 1−sinθ 2 s∈φ2
( 2 m+ 1 )∗π
|F (θ , φ )|=
2
√ 1−sinθ 2 s∈φ 2
Pour m=1 :
3π
|F (θ , φ )|=
2
√ 2
1−sinθ cosφ
2
+Pour φ=0 ° :
{cosφ=1
φ=0 °
|F (φ)|=
3π
2
√ 1−sin φ
2
|F (φ)|=
3π
2
cos φ
% Paramètres
phi = linspace(0, 2*pi, 360); % Angle azimutal
(0 à 2*pi pour une onde complète)
F_phi = (3*pi/4) * sin(phi);

figure;
polarplot(phi, abs(F_phi), '-b');
Onde');
grid on;
+Pour φ=0 ° :
{
− jkr
φ=0 ° k I0 Le
sinφ=0
Eψ = jη √1
4 πr
− jkr
k I0 Le
Eψ = jη
4 πr
% Paramètres
f = 1e9; % Fréquence en Hertz
c = 3e8; % Vitesse de la lumière en m/s
lambda = c / f; % Longueur d'onde
% Dipôle
L = lambda / 2; % Longueur du dipôle
I0 = 1; % Courant du dipôle
% Constantes
eta = 120 * pi; % Impédance du milieu (pour l'air,
environ 120*pi ohms)
+Pour φ=90 ° :
{
− jkr − jkr
φ=90 ° k I0 Le k I0 Le
sinφ=1
Eψ = jη
4 πr
√1−sinθ2 Eψ = jη
4 πr
cosθ
% Paramètres
% Dipôle
% Constantes
eta = 120 * pi; % Impédance du milieu (pour l'air, environ
120*pi ohms)
% Angle theta (en radians)

theta = linspace(0, 2*pi, 361);
% Distance radiale
r = lambda;
% Calcul du champ électrique

k = 2 * pi / lambda;
E_psi = 1j * eta * (k * I0 * exp(-1j * k * r) / (4 * pi *
r)) * cos(theta);
% Tracer le diagramme de rayonnement

figure;
polarplot(theta, abs(E_psi));
title('Diagramme de rayonnement dun
dipôle électrique');
+Pour θ=90 ° :
{
− jkr − jkr
θ=90° k I0 Le k I0 Le
sinθ=1
Eψ = jη
4 πr
√1−sinφ 2
Eψ = jη
4 πr
cosφ
% Paramètres
% Dipôle
% Constantes
eta = 120 * pi; % Impédance du milieu (pour l'air,
environ 120*pi ohms)
% Angle phi (en radians)

% Distance radiale
r = lambda;
% Coordonnées cartésiennes (x, y, z)

x = r * cos(phi);
y = zeros(size(phi)); % Dipôle placé le long de l'axe y
z = r * sin(phi);
% Calcul du champ électrique

k = 2 * pi / lambda;
E = 1j * eta * (k * I0 * L * exp(-1j * k * r) / (4 * pi *
r)) * cos(phi);
% Tracer le diagramme de rayonnement

figure;
polarplot(phi, abs(E));
title('Diagramme de rayonnement d''un dipôle électrique');

Travaux Pratiques Les Antennes DI

Transféré par

Informations du documentcliquez pour développer les informations du document

Droits d'auteur :

Formats disponibles

Travaux Pratiques Les Antennes DI

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Travaux Pratiques Les Antennes DI

Transféré par

Droits d'auteur :

Formats disponibles

FACULTE POLYDISCIPLINAIRE DE LARACHE

Master : Systèmes Électroniques Et Energies Renouvelables

TP n°1 : Diagramme de rayonnement d’un dipôle rayonnant par MATLAB

Schéma d’un dipôle rayonnant

|cos ( KLcos ( θ ) ) −cos ⁡(KL)|

 Dipôle à onde entière :

% Fonction caractéristique F(theta)

% Tracé du diagramme de rayonnement en coordonnées

% Ajustements d'affichage facultatifs

% Tracé du diagramme de rayonnement en coordonnées

% Fonction caractéristique F(theta)

% Tracé du diagramme de rayonnement en coordonnées

% Tracé du diagramme de rayonnement en coordonnées sphériques

% Fonction caractéristique F(theta)

% Tracé du diagramme de rayonnement en coordonnées

% Ajustements d'affichage facultatifs

% Angles azimutaux et d'élévation

% Meshgrid pour les coordonnées sphériques

% Fonction caractéristique F(theta)

% Conversion des coordonnées sphériques en coordonnées

% Tracé du diagramme de rayonnement en coordonnées

% Fonction caractéristique F(theta)

% Meshgrid pour les coordonnées sphériques

% Fonction caractéristique F(theta)

% Conversion des coordonnées sphériques en coordonnées

% Tracé du diagramme de rayonnement en coordonnées

% Fonction caractéristique F(theta)

% Angles azimutaux et d'élévation

% Meshgrid pour les coordonnées sphériques

% Fonction caractéristique F(theta)

% Conversion des coordonnées sphériques en

% Tracé du diagramme de rayonnement en

% Fonction caractéristique F(theta)

% Meshgrid pour les coordonnées sphériques

% Fonction caractéristique F(theta)

% Conversion des coordonnées sphériques en

% Tracé du diagramme de rayonnement en

% Fonction caractéristique F(theta)

% Angles azimutaux et d'élévation

% Meshgrid pour les coordonnées sphériques

% Fonction caractéristique F(theta)

% Conversion des coordonnées sphériques en

% Tracé du diagramme de rayonnement en coordonnées

% Fonction caractéristique F(theta)

% Angles azimutaux et d'élévation

% Meshgrid pour les coordonnées sphériques

% Fonction caractéristique F(theta)

% Conversion des coordonnées sphériques en

% Tracé du diagramme de rayonnement en coordonnées

2-a Un dipôle électrique infinitésimal horizontal de courant constant I0 est placé

% Calcul de la fonction caractéristique

% Tracé du diagramme de rayonnement en

% Création de la grille sphérique

% Calcul de la fonction caractéristique

% Conversion des coordonnées sphériques en

% Tracé du diagramme de rayonnement en

% Tracé du diagramme de rayonnement en