共享雷达信号检测包括MTI脉冲多普勒CFAR处理-radar_signal_processing.rar

共享雷达信号检测包括MTI脉冲多普勒CFAR处理-radar_signal_processing.rar
本帖最后由 雷达信号处理 于 2013-11-2 22:13 编辑    本人近段时间对雷达信号检测与判决做了一些学习和研究，还处在入门阶段，现应用一个仿真实例对动目标信号的检测判决算法做了分析和验证。其中包括MTI，脉冲多普勒处理，CFAR自适应门限检测等。这里分享一下，大家一起研究研究，还请各位前辈指点。谢谢。

%% 雷达信号检测与判决算法研究与仿真实例

clear all;
close all;
clc;

%% 动目标

fr = 1e3;                                                                  % 脉冲重复频率（Hz）
tr = 1 / fr;                                                               % 脉冲重复间隔（s）
fs = 1e6;                                                                  % 快时间维采样频率（Hz）（大于等于发射脉冲带宽，发射脉冲带宽几乎总是小于等于载波频率的10%，通常为载波频率的1%）
ts = 1 / fs;                                                               % 快时间维采样间隔（s）
fc = 100e6;                                                                % 载波频率（Hz）
vr = 100;                                                                  % 动目标相对雷达的径向速度（m/s）                          
c = 3e8;                                                                   % 电磁波速（m/s）                                         
lamada = c / fc;                                                           % 载波波长（m）
fd = 2 * vr / lamada;                                                      % 动目标多普勒频率（Hz）（或多普勒移动，即发射频率与接收频率之差）
M = fs / fr;                                                               % 1个脉冲重复间隔内的采样点数（快时间维的距离点数）
M10 = 10 * M;                                                              % 10个脉冲重复间隔内的采样点数

m = 1 : M10;
target = 20 * exp（j * （2 * m * pi * fd / fr））;                             % 动目标回波

figure（1）;
subplot（211）;
plot（real（target））;
title（‘动目标——实部‘）;
xlabel（‘距离点‘）;
grid on;
subplot（212）;
plot（imag（target））;
title（‘动目标——虚部‘）;
xlabel（‘距离点‘）;
grid on;
f1 = linspace（0 fr length（target））;
figure（2）;
plot（f1 abs（fft（target）））;
title（‘动目标幅频特性‘）;
xlabel（‘频率（Hz）‘）;
ylabel（‘幅度‘）;
grid on;

%% 杂波

randn（‘state‘ sum（100 * clock））;
clutteri = 5 * randn（1 length（target））;
randn（‘state‘ sum（10 * clock））;
clutterq = 5 * randn（1 length（target））;
clutter = clutteri + j * clutterq;                                         % 杂波

figure（3）;
subplot（211）;
plot（real（clutter））;
title（‘杂波——实部‘）;
xlabel（‘距离点‘）;
grid on;
subplot（212）;
plot（imag（clutter））;
title（‘杂波——虚部‘）;
xlabel（‘距离点‘）;
grid on;
f2 = linspace（0 fr length（clutter））;
figure（4）;
plot（f2 abs（fft（clutter）））;
title（‘杂波幅频特性‘）;
xlabel（‘频率（Hz）‘）;
ylabel（‘幅度‘）;
grid on;

%% 动目标与杂波的混合信号

s = target + clutter;                                                      % 动目标与杂波的混合信号

figure（5）;
subplot（211）;
plot（real（s））;
title（‘动目标与杂波的混合信号——实部‘）;
xlabel（‘距离点‘）;
grid on;
subplot（212）;
plot（imag（s））;
title（‘动目标与杂波的混合信号——虚部‘）;
xlabel（‘距离点‘）;
grid on;
f3 = linspace（0 fr length（s））;
figure（6）;
plot（f3 abs（fft（s）））;
title（‘混合信号幅频特性‘）;
xlabel（‘频率（Hz）‘）;
ylabel（‘幅度‘）;
grid on;

%% MTI杂波抑制（三脉冲对消处理）

%————————————————对消运算—————————————————
x = ones（1 4）;                                                            % 动目标所在的距离点
n = length（x）;
points1 = [zeros（1 300） x zeros（1 M - 300 - n）];                       % 1个脉冲重复间隔内M个距离点（区分动目标所在的距离点301、302、303、304和杂波所在距离点）
points10 = repmat（points1 1 10）;                                         % 10个脉冲重复间隔内

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      43700  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection.html

     文件       2365  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection.png

     文件       3956  2013-11-02 21:29  radar_signal_processing\html\radar_signal_detection_01.png

     文件       6768  2013-11-02 21:29  radar_signal_processing\html\radar_signal_detection_02.png

     文件       9163  2013-11-02 21:29  radar_signal_processing\html\radar_signal_detection_03.png

     文件       7174  2013-11-02 21:29  radar_signal_processing\html\radar_signal_detection_04.png

     文件       9020  2013-11-02 21:29  radar_signal_processing\html\radar_signal_detection_05.png

     文件       6980  2013-11-02 21:29  radar_signal_processing\html\radar_signal_detection_06.png

     文件       8780  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_07.png

     文件       6613  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_08.png

     文件       7978  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_09.png

     文件       8495  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_10.png

     文件       8838  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_11.png

     文件       8837  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_12.png

     文件       8353  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_13.png

     文件       8908  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_14.png

     文件       9245  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_15.png

     文件       9132  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_16.png

     文件       8950  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_17.png

     文件      11615  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_18.png

     文件       9800  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_19.png

     文件       6785  2013-11-02 21:30  radar_signal_processing\html\radar_signal_detection_20.png

     文件      11533  2013-11-02 21:10  radar_signal_processing\radar_signal_detection.m

     目录          0  2013-11-02 21:30  radar_signal_processing\html

     目录          0  2013-11-02 21:37  radar_signal_processing

----------- ---------  ---------- -----  ----

               222988                    25