一种干扰高分辨率雷达的方式

为了对(HRR)雷达和其他高分辨成像雷达进行假目标欺骗，需要在其脉冲回波中进行RCS调制，生成具有高保真干扰信号，下面的方法它可以在频域通过数字信号处理生成假目标，数字射频存储调制器通过许多具有个别调制和可信背景的单反射器，可以生成若干复杂假目标。该调制器可以合成高分辨距离剖面图，可以将高分辨雷达假目标欺骗系统，设计为可以存储、调制和转发雷达信号的中继器。

转发一个完美复制的雷达信号，可以在雷达上产生点目标响应。不过，要在高分辨雷达图像上生成假目标信息，必须生成一个目标或具有空间范围的场景。

一般可以通过目标RCS模型同采样信号进行时域卷积，来生成一个具有距离向调制的目标信号矢量和，时域卷积其实就是一个数字有限脉冲滤波器。
如图所示，在时域中基本可以通过数字有限脉冲(FIR)滤波器实现。

通过离散傅里叶变换(DFT)实现快速傅里叶变换(FFT)，采用快速卷积处理(FCP)算法，可以计算雷达信号和点目标两者中任意一个卷积。在计算需求方面，快速傅里叶变换的卷积效率要比标准卷积效率高得多(尤其当序列长度较长时)。

这样做的好处是，采用FCP卷积可以快速直接计算出总响应，而不用单独进行每个RCS回波点硬件的延迟和调制。可以将雷达信号和目标转换到频域再进行这样的计算，然后再将雷达信号和目标转换回时域，采用这种方法有助于生成更多数量的雷达分辨单元。而直接通过时域实现时，目标或分辨单元的数量受可编程逻辑器件内存的限制，在进行电子对抗计算时，受到系统动态范围、采样位数和计算速度等因素的限制。

如果在模拟距离向回波时，也可以同步来自激励信号的视频检波信号，只采集和调制有用的脉冲信号，进一步减少计算量。

以下给出的Matlab代码仿真流程：

clear;clc;close all
T   =  1e-4;       %最大探测距离 c*T/2=15000m
tao =  10e-6;
B   =  150e6;
fs  =  2*B;
c   =  3e8;  %光速
stepm = 3e8/fs; %m/采样间隔 7.5m/p

K   =  B/tao; %调制斜率
R   =  5000;
delay = 2*R./c;
t=-tao/2:1/fs:tao/2-1/fs;
x_linear=exp(1j*pi*(K.*(t.^2)));  %发射信号

figure
subplot(211);plot(t,real(x_linear));title('时域波形')
X=fftshift(fft(x_linear));
f=linspace(0,fs,length(t))-fs/2;
subplot(212);plot(f,abs(X));title('频域')

tt = 0:1/fs:T-1/fs;
echo = zeros(1,length(tt));
signrect = rectpuls(tt-delay(1)-tao/2,tao);
MKL =signrect .*exp(1j*pi*K*(tt-delay(1)-tao/2).^2)+randn(1,length(tt))*0.2 ;

%按照一定步进生成RCS
miniLen = 150*2;%模拟总长度
miniStep = 50*2;%模拟步进
allCount = round(miniLen/stepm); %总的模拟点个数
stepCount = round(miniStep/stepm);%模拟点步进
logic_array = zeros(1, allCount);%
index = round([1:stepCount:allCount]);
logic_array(index) = 5; %设置回波点的RCS值
RCSP =   logic_array;
%将采集到的雷达信号同RCS信号进行卷积，生成回波信号
rcsfft = fft(RCSP,length(echo));
echo1=ifft(conj(rcsfft).*fft(MKL,length(MKL)));

%RCS调制后的回波信号和原始信号进行卷积进行脉冲压缩
orgonFFt = fft(x_linear,length(echo));%原始线调信号
output=ifft(conj(orgonFFt).*fft(echo1,length(echo1)));  %h和echo做fft的长度要一样

disLabel = tt*c/2;
figure;
subplot(411);plot(disLabel, abs(output)/length(output)/length(index), ...
    disLabel,real(echo1)/length(RCSP));
title('回波');xlabel('distence');
subplot(412);plot(c/2*(0:1/fs:(length(RCSP)-1)/fs),abs(RCSP));title('原始RCS');xlabel('distence');
subplot(413);plot(c/2*(0:1/fs:T-1/fs),abs(output));title('脉冲压缩输出波形');xlabel('distence');
subplot(414);plot((0:1/fs:T-1/fs)*c/2,db(abs(output)/max(abs(output))),'r');
title('脉压结果分贝图');ylabel('dB');xlabel('m');