lpiinterceptfactor.m

用于低截获概率雷达设计

% function LPIinterceptFactor
% 根据截获接收机的情况计算是否具有低截获概率性能，即计算截获因子是否小于1。
% 只考虑对方侦察接收机在导引头雷达作战目标的平台上的情况，即雷达的主瓣对准目标，侦察机的主瓣对准雷达
% 
% fuxiongjun @ 2008-03-08
close all;
clear all;

Rr0=50e3;   % 雷达到目标/侦察机的距离，km，初始值
Pt=10;     % 40  20 10
duty=0.3;       % 占空比
 PRF=1.5e3;
 PRI=1/PRF;
RCS=3900;       % 假设目标的RCS,m2
fc=17e9;        % carrier frequency
renta=3e8/fc;   % wavelength

Gt=10^(28/10);  % 雷达在目标方向上的发射天线功率增益（发射天线主瓣增益）
Gr=Gt;          % 雷达在目标方向上的接收天线功率增益（接收天线主瓣增益）
    SidBmi=-80;
    SidBi=SidBmi-30;
Si=10^(SidBi/10);% 截获接收机灵敏度


Gi=10^(12/10);   % 截获接收机在雷达发射天线方向的增益(按SLQ-32)
G =Gt;           % 雷达发射天线在截获接收机方向的增益.(雷达主瓣对准侦察机时为Gt，雷达副瓣对准侦察机时G是雷达旁瓣增益）

% Sr=;             % 雷达接收机灵敏度
% Sr也可按下式计算（如果不是直接提供的话）:
 Bn=5e6;             % band width, Hz－－实际与此无关
 T0=290;                  %  absoulte temperature
 k=1.38e-23;            % Boltzmann constant
 Fn=10^(4/10);        % Receiver noise figure,  dB->original data
 
%  SNRmin=10^(-3/10);   % 最小可检测信噪比
     % SNR可按下式计算
     SNR0=10^(19/10);% 单脉冲不脉压时进行检测所需要的信噪比(Minimum detectable signal-to-noise ratio, ->original data)
     M=128;          % coherent integration pulse number, if only use single pulse, then M equals 1
     PW=200e-6;      % 脉宽
     Pcr=PW*Bn;      % 脉压比,即“时宽带宽积”
     SNRmin=SNR0/(M*Pcr);
     
Sr= k*T0*Bn*Fn*SNRmin;


% 雷达系统损耗
Lr1=10^(10/10);            % 雷达系统损耗（不含双程传播损耗）
    att1=0.07;  %大气衰减系数：晴天：0.07；4mm/h：0.39；7.5mm/h：0.8－－已经考虑了双程
    Lr2=10^(att1*(Rr0/1e3)/10);  % 初步假设的可达到距离为Rr0，以计算大气衰减
Lr=Lr1*Lr2;

% 上面分析损耗用的是假设的一个距离，实际上分析截获因子，应该按最大的雷达作用距离和最大的截获距离相比，因此需要首先计算在pt功率下的雷达最大作用距离
% 在不同距离下，大气损耗相差很大（Ku波段），因此需要精确一些
Pa=Pt*(PW/PRI);  %  平均功率
Rmax=(Pa*PRI*M*RCS*G^2*renta^2/(4*pi)^3/SNR0/k/T0/Fn/Lr)^(1/4)  % non-coherent integration
while abs(Rmax-Rr0)>=0.0001   % 通过循环得到一个真实值
    Rr0=(Rmax+Rr0)/2;
    Lr2=10^(att1*(Rr0/1e3)/10);  % 初步假设的可达到距离为Rr0，以计算大气衰减
    Lr=Lr1*Lr2;
    Rmax=(Pa*PRI*M*RCS*G^2*renta^2/(4*pi)^3/SNR0/k/T0/Fn/Lr)^(1/4);  % non-coherent integration
end
Rmax

% 侦察机系统损耗
Li1=10^(10/10);   % 截获接收机的接收系统损耗(不含大气衰减）
Li2=Lr2/2;       % 单程大气衰减损耗
Li=Li1*Li2;


alpha=((10*Gi^2*G^2*Sr*Lr*renta^2)/(4*pi*RCS*Gt*Gr*Si^2*Li^2))^(1/4)   % 截获因子。小于1则表明该雷达为LPI雷达（相对于程序中设定的侦察接收机）