为了对(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');