mti杂波抑制仿真matlab-MTI+MTD+CFAR简单仿真.pdf

mti杂波抑制仿真matlab,MTI+MTD+CFAR简单仿真 ⽂章 ⽬录 代码实现 运⾏效果 ⼏点注意 代码实现 %% 清空⼀切 c ear a ;c ose a ;c c; %% 仿真设置 USE_COMPLEX_SIG = 0; USE_MTI = 1; USE_PRINT_INFO = 1; USE_CFAR_Method = 3; USE_CFAR_Custom_Thd = 0; %% 参数 fft_num = 128; % FFT运算点数 pu se_num = fft_num + 2; % 总的脉冲数 prf = 1e3; % 脉冲重复频率(Hz) pri = 1/prf; % 脉冲重复间隔(s) fs = 1e6; % 快时间维采样频率(Hz)(⼤于等于发射脉冲带宽，发射脉冲带宽⼏乎总是⼩于等于载波频率的10%，通常为载波频率的1%) ts = 1/fs; % 快时间维采样间隔(s) fc = 100e6; % 载波频率(Hz) vr = 300; % 动 ⽬标相对雷达的径向速度(m/s) c = 3e8; % 光速(m/s) amada = c/fc; % 载波波长(m) fd = 2*vr/ amada; % 动 ⽬标多普勒频率(Hz)(或多普勒移动，即发射频率与接收频率之差) np_fast = fs/prf; % 1个脉冲重复间隔内的采样点数(快时间维的距离点数) np_s ow = pu se_num * np_fast; % 多个脉冲重复间隔内的采样点数 target_start_index = 200; % ⽬标起始距离单元 target_end_index = 203; % ⽬标结束距离单元 %% ⽣成动 ⽬标回波+杂波 m = 1:np_s ow; target = 20*exp(1i*2*pi*fd*m/fs); % ⽣成 动 ⽬标回波(此处直接⽣成包含多普勒信息的信号,不含载频)(注意式中采⽤fs,⽽⾮fr) rng defau t; c utteri = 5 * randn(1, ength(target)); c utterq = 5 * randn(1, ength(target)); c utter = comp ex (c utteri,c utterq); % ⽣成 杂波 points 1 = zeros(1,np_fast); points 1(1,target_start_index:target_end_index) = 1; % 1个脉冲重复间隔内np_fast个距离点 points_tota = repmat(points 1, 1, pu se_num); % pu se_num个脉冲重复间隔内pu se_num*np_fast个距离点 s = points_tota .* target + c utter; % ⽣成 动 ⽬标与杂波的混合信号 if USE_COMPLEX_SIG == 0 s = rea (s); end %% MTI杂波抑制(三脉冲对消处理) if USE_MTI == 1 mti_in = reshape(s, ength(s)/pu se_num, pu se_num); % 形成由多个脉冲组成的⼆维数据矩阵(慢时间维pu se_num个脉冲，快时间 维1000个距离点) b = [1, -2, 1]; % 三脉冲对消器系统函数分⼦多项式系数向量 a = [1]; % 三脉冲对消器系统函数分母多项式系数向量 mti_out = fi ter(b, a, mti_in); % 三脉冲对消运算(同⼀距离点在3个相邻脉冲重复间隔内的数据做相减运算) mti_out = mti_out(3:end,:); e se mti_in = reshape(s, ength(s)/pu se_num, pu se_num); % 形成由10个脉冲组成的⼆维数据矩阵(慢时间维10个脉冲，快时间维1000 个距离点) mti_out = mti_in(3:end,:); end %% 脉冲多普勒处理 win = repmat(hamming(fft_num),1,np_fast); % FFT采样点数即多普勒谱分析样本点数亦即MTD滤波器组滤波器个数;hamming(K)是 ⼀个列向量 mtd_out = fft(mti_out.*win); % 对同⼀距离点不同脉冲重复间隔的数据做K点FFT加窗运算(对矩阵进⾏FFT时,是分别计算各列的FFT) figure;