【雷达】基本的雷达和干扰机仿真（Matlab实现）

原创已于 2024-06-17 16:39:08 修改 · 873 阅读

12 ·

CC 4.0 BY-SA版权

文章标签：

#信号处理 #matlab

于 2024-06-17 16:26:06 首次发布

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

💥1 概述

基本的雷达和干扰机仿真是对雷达系统和干扰机性能及相互作用的数字化模拟。通过建立数学模型和算法，模拟雷达发射信号、接收回波以及干扰机产生的干扰信号。在仿真中，详细考虑了雷达的工作原理，包括信号的生成、发射、传播、目标反射和接收处理等环节。对于干扰机，研究其干扰方式，如噪声干扰、欺骗干扰等，以及干扰信号的特性和对雷达性能的影响。这种仿真有助于深入理解雷达与干扰机之间的复杂关系，评估不同干扰策略下雷达的检测、跟踪和识别能力。为雷达系统的优化设计和抗干扰措施的制定提供了重要的理论依据和实验数据，能够在实际应用前，以较低成本和风险预测系统性能，从而提高雷达在复杂电磁环境中的作战效能。

📚2 运行结果

部分代码：

clear
rad=pi/180;
% RADAR PARAMETERS **************************************************************
Pt=1e3;          % radar trasmitter power
hr=0;           % radar height
Waz=3;           % radar antenna az aperture in m (sinc pattern)
Wel=1;           % radar antenna el aperture in m (sinc pattern)
TA=1000;         % antenna temperature
Te=1000;         % effective temperature of radar receiver
Ts=TA+Te;        % radar system noise temperature
er=0.8;          % radar antenna efficiency
freq=1e9;        % frequency
c=3e8;           % velocity of light
wave=c/freq;     % wavelength
k=2*pi/wave;     % wave number (propagation constant)
Gpdb=20;         % radar processing gain in db
Gp=10^(Gpdb/10);
tau=1e-6;        % radar pulse width
Bn=1/tau;        % radar noise bandwidth
boltz=1.38e-23;  % Boltzman's constant
No=boltz*Ts*Bn;  % noise power
SNRmindb=10;     % threshold SNR
SNRmin=10^(SNRmindb/10);
Smin=SNRmin*No;  % minimum detectable signal (MDS) in db
%------------- NEW ANTENNA MODEL (TAYLOR) ---------------------
% plot an antenna pattern for reference
 it=0; 
 SLLdb=-40; BLdb=-60;
 BL=10^(BLdb/20);
   for ang=-180:.5:180
       Thrad=ang*rad; Phrad=0;  it=it+1; A(it)=ang;
% elevation pattern factor (phi = 90 plane)
Fel=1;
arg1=k*sin(Thrad)*sin(Phrad)/2;
if abs(arg1)>1e-5
   Fel=abs(sin(Wel*arg1)/arg1/Wel)^2;
end
% azimuth pattern factor (phi = 90 plane)
Naz=floor(Waz/(0.5*wave));   % must be even
if Naz<6, Naz=6; end  % use at least 6 elements
dspc=Waz/(Naz-1);
arg2=k*sin(Thrad)*cos(Phrad);
if SLLdb <= -15
    nbar=max(5,floor(Naz/5));
    amp=tayl(Naz,SLLdb,nbar);