Exercices en Matlab
Exercice 1 : Simulation des performances d’un radar à impulsions
Un radar à impulsions fonctionne à une fréquence de 9 GHz, avec une puissance d’émission de 300 kW,
un gain d’antenne de 30 dB, une largeur de faisceau de 2°, une durée d’impulsion de 1 µs, une
fréquence de répétition d’impulsions (PRI) de 1 ms. La sensibilité du récepteur est de -115 dBW. On
considère des cibles de surface équivalente radar σ variant de 1 m² à 100 m², situées à des distances
entre 10 km et 100 km.
1. Implémenter l’équation radar pour calculer la puissance reçue en fonction de la
distance et de la SER.
2. Tracer un graphique 2D de la puissance reçue en fonction de la distance et de la SER.
3. Calculer et afficher :
La portée maximale pour une cible de σ=10 m2
La distance non ambiguë
La distance aveugle
La résolution en distance
La résolution angulaire à 60 km
La dimension de la cellule de résolution à 60 km
4. Tracer la zone de détection (où la puissance reçue > seuil) sur un graphe.
Exercice 2 : Analyse et Simulation d'un Système Radar
L'exercice est divisé en trois parties :
Partie 1 : Calcul des Paramètres Radar
1. Implémentation de l'Équation du Radar :
o Écrivez une fonction MATLAB portee_radar qui prend les entrées suivantes :
Pt : Puissance émise (W)
G : Gain de l'antenne (dB)
lambda : Longueur d'onde (m)
sigma : Surface équivalente radar de la cible (m²)
Pe_min : Puissance minimale détectable (W)
L_sys : Pertes du système (dB)
o La fonction doit calculer et retourner la portée radar maximale ( Rmax) en utilisant
l'équation du radar. N'oubliez pas de convertir les valeurs en dB en unités
linéaires si nécessaire.
o Incluez une gestion des erreurs pour vérifier les entrées invalides (par exemple,
puissance négative, longueur d'onde nulle).
Exercices en Matlab
Tester la fonction avec
Paramètre Valeurs de Test Valeurs Test 2 Valeurs Test 3
1
Pt (W) 5000 10000 2000
G (dB) 35 40 30
λ (m) 0.03 0.05 0.1
σ (m²) 1.5 3 0.8
Pe_min (W) 1e-13 5e-13 2e-13
L_sys (dB) 5 3 8
2. Analyse de la Variation des Paramètres :
Créez un script MATLAB qui :
Définit un ensemble de paramètres radar (Pt, G, lambda, sigma, Pe_min,
L_sys).
Utilise la fonction portee_radar pour calculer Rmax.
Trace des graphiques montrant la variation de Rmax en fonction de :
Puissance émise (Pt)
Gain de l'antenne (G)
Surface équivalente radar (sigma)
Gardez les autres paramètres constants tout en faisant varier un seul.
Utilisez une échelle logarithmique pour Pt et sigma.
Étiquetez correctement les axes et ajoutez des titres aux graphiques.
Ajoutez une grille aux graphiques.
Tester le script avec :
Pt_values : de 100 W à 100 000 W (50 points logarithmiques)
G_values : de 25 dB à 45 dB (pas de 1 dB)
sigma_values : de 0.1 m² à 100 m² (50 points logarithmiques)
Paramètres fixes :
o G = 35 dB
o λ = 0.03 m
o σ = 1.5 m²
o Pe_min = 1e-13 W
o L_sys = 5 dB
Questions
Commentez l’allure des trois courbes obtenues.
Que pouvez-vous dire de la relation entre Pt, G, σ et la portée radar maximale Rmax ?
Pourquoi utilise-t-on une échelle logarithmique pour Pt et σ mais pas pour G ?
Exercices en Matlab
À votre avis, quel paramètre semble avoir l’impact le plus fort sur la portée ? Justifiez.
Partie 2 : Simulation de la Poursuite de Cible
1. Modèle de Mouvement de la Cible :
o Supposons qu'une cible se déplace dans un plan 2D à vitesse constante.
o Écrivez une fonction MATLAB position_cible qui prend :
position_initiale : [x0, y0] (m)
vitesse : [vx, vy] (m/s)
temps : Un vecteur de valeurs de temps (s)
o La fonction doit retourner la position de la cible [x, y] à chaque pas de temps.
o Testez la fonction avec les paramètres proposés.
o Tracez la trajectoire de la cible pour chaque test :
Valeurs de test
Cas de test Position Initiale [x0, y0] (m) Vitesse [vx, vy] (m/s) Temps (s)
Test 1 [0, 0] [10, 5] 0:1:10
Test 2 [100, 50] [-5, 0] 0:0.5:5
Test 3 [50, 50] [0, -2] 0:2:20
Questions
1. Quel est le type de trajectoire obtenu si la vitesse est constante dans les deux directions ?
Pourquoi ?
2. Comment la trajectoire évoluerait-elle si la vitesse dans l’axe Y était nulle ?
3. Que se passe-t-il si la vitesse dans l’axe X est négative ? Interprétez graphiquement.
4. Pourquoi utilise-t-on ici un vecteur de temps et non un scalaire ?
5. Que représente le pas de temps dans l’analyse de mouvement radar ? Et quel est son
impact si on le diminue ou l’augmente ?
6. Proposez un cas de test où la cible reste immobile. Que devrait afficher le graphique ?
7. En considérant un radar fixe à (0,0), à partir de quelle distance la cible sera-t-elle hors de
portée si le radar a un Rmax de 500 m ? Calculez et commentez.
2. Simulation de la Détection Radar :
o Écrivez un script MATLAB qui :
Définit la position et les paramètres du radar (comme dans la Partie 1).
Définit la position initiale et la vitesse de la cible.
Génère un vecteur de temps.
Calcule la position de la cible au fil du temps en utilisant
position_cible.
Exercices en Matlab
Calcule la distance entre le radar et la cible à chaque pas de temps.
Simule la détection radar :
Si la distance est inférieure ou égale à Rmax (calculée à l'aide de
portee_radar), la cible est "détectée" à ce moment-là.
Introduisez un élément aléatoire pour simuler le bruit : Ajoutez une
probabilité de détection (par exemple, 0,9) - la cible n'est détectée
qu'avec cette probabilité même si elle est à portée.
Trace la trajectoire de la cible et la position du radar.
Marque les positions détectées avec un marqueur différent.
Trace la portée du radar sous forme de cercle autour du radar.
Étiquetez correctement les axes et ajoutez des titres aux graphiques.
Ajoutez une grille aux graphiques.
Ajoutez une légende aux graphiques.
Paramètres de test
Paramètre Valeur
Position du radar [0, 0] m
Puissance émise Pt 1000 W
Gain G 30 dB
Longueur d'onde λ 0.03 m
Surface radar σ 1 m²
Puissance min Pe_min 1e-10 W
Pertes système L_sys 3 dB
Position initiale cible [500, 0] m
Vitesse cible [-20, 5] m/s
Vecteur de temps 0:1:50 s
Probabilité de détection 0.9
Pd
Questions
1. Que se passe-t-il si la probabilité de détection Pd = 1 ? Et si Pd = 0 ?
2. Pourquoi introduit-on une probabilité de détection dans la simulation ? Quel phénomène
réel cela modélise-t-il ?
3. À partir de quel instant (temps) la cible entre-t-elle dans la portée radar ?
4. Pourquoi utilise-t-on axis equal dans la représentation graphique ?
5. Si la vitesse de la cible passe à [0, 10] m/s, comment cela impacte-t-il la trajectoire et les
détections ?
6. Pourquoi la portée du radar est-elle représentée par un cercle ?
7. Modifiez la position initiale de la cible à [700, -100] m : quelles sont les conséquences
sur le moment des premières détections ?
8. Ajoutez un second radar à la position [200, 200] avec les mêmes paramètres : comment
pourrait-on modifier le script pour simuler et afficher cette double détection ?
Exercices en Matlab
Partie 3 : Analyse de la Résolution
[Link] de la Résolution
Écrivez une fonction MATLAB nommée resolution_radar qui prend en entrée :
c : vitesse de la lumière (en m/s)
tau : durée de l'impulsion radar (en secondes)
theta_bw : largeur du faisceau de l’antenne (en degrés)
R : distance de la cible (en mètres)
La fonction doit retourner :
La résolution en portée delta_R
La résolution angulaire en mètres delta_theta à la distance R
[Link] des Distances Radar Caractéristiques
Implémentez deux fonctions :
distance_non_ambigue qui calcule la distance non ambiguë R_na à partir de :
o c : vitesse de la lumière
o PRI : période de répétition des impulsions (en secondes)
distance_aveugle qui calcule la distance aveugle R_b à partir de :
o c : vitesse de la lumière
o tau : durée de l'impulsion
[Link]énario à Deux Cibles
Écrivez un script MATLAB qui :
Définit les paramètres du radar (c, tau, theta_bw, PRI)
Définit les positions de deux cibles situées à la même distance R mais à des angles
d’azimut différents
Calcule la résolution en portée et en angle à l’aide de votre fonction resolution_radar
Calcule la distance non ambiguë et la distance aveugle
Détermine si les deux cibles sont résolues en portée et en angle selon leur séparation et
les résolutions calculées
Trace :
o Le radar en (0,0)
o Les deux cibles
o La cellule de résolution autour de la première cible
o La distance aveugle sous forme de cercle
o La distance non ambiguë sous forme de cercle
Étiquète correctement les axes et ajoute une légende, un titre et une grille.
Exercices en Matlab
Essayez de représenter la zone aveugle et la zone de détection non ambiguë en couleurs
différentes sur le graphique.
Proposez un cas où deux cibles sont ambiguës en portée.
Paramètres de Test :
Paramètre Valeur
C 3×10⁸ m/s
Tau 1×10⁻⁶ s
theta_bw 2°
PRI 1×10⁻³ s
Portée des 3000 m
cibles
Position cible 1 (3000, 0) m
Position cible 2 (3000, 150) m
Consignes supplémentaires aux étudiants :
Testez en changeant :
angle2 pour voir quand les cibles deviennent non résolues angulairement
PRI pour voir l’impact sur R_na
tau pour observer l’évolution de delta_R et R_b
Proposez un cas :
Où R_cible > R_na pour forcer une ambiguïté
Où abs(angle2 - angle1) < theta_bw pour une non-résolution angulaire
Questions
1. Quelle est la différence entre la distance non ambiguë et la distance aveugle ?
2. Que se passe-t-il si deux cibles sont séparées de 100 m en angle à 3 km avec un faisceau
de 2° ? Sont-elles résolues ?
3. Comment la résolution en portée évolue-t-elle si tau augmente ? Pourquoi ?
4. Pourquoi la distance non ambiguë dépend-elle du PRI et non de tau ?
5. Que représente la cellule de résolution tracée autour de la cible 1 ?
6. Que se passerait-il si le PRI est réduit à 0.5 ms ? Impact sur R_na ?
7. Pourquoi la distance aveugle limite-t-elle la détection des cibles proches ?
8. Comment améliorer la résolution angulaire sans modifier la portée du radar ?
9. Proposez un cas où deux cibles sont non résolues en portée mais résolues en angle.
10. Modifiez la valeur de theta_bw à 1° et observez. Quelles conclusions peut-on tirer ?
Exercices en Matlab