线性调频连续波雷达仿真matlab，汽车雷达测距测速雷达仿真

% %使用FMCW技术的汽车自适应巡航控制
% 该示例显示了如何使用调频连续波（FMCW）技术对汽车自适应巡航控制系统进行建模。
% 该示例执行了移动车辆的距离和多普勒估计。
% 与国防工业中常见的脉冲雷达系统不同，汽车雷达系统通常采用FMCW技术。 
% 与脉冲雷达相比，FMCW雷达体积更小，功耗更低且制造成本更低。 结果，FMCW雷达只能监视更小的距离。

%   版权所有2012-2014 The MathWorks，Inc.
clear;
clc;
close all;
%% FMCW Waveform 
% 考虑用于自动巡航控制（ACC）的汽车远程雷达（LRR）。 如[1]所示，这种雷达通常占用77 GHz附近的频带。 
% 雷达系统会不断估算其所安装的车辆与前方车辆之间的距离，并在两者之间的距离过近时向驾驶员发出警报。
%
% <>
%上图显示了ACC的示意图，还说明了使用FMCW技术进行距离测量的原理。 接收到的信号是发射信号的时延copy，
% 其中延迟 \Delta t与距离有关。 因为信号总是在整个频带中扫描，所以在扫描期间的任何时刻，频率差$ f_b $是发射信号和接收信号之间的常数。
% f_b通常称为拍频（频差）。 由于the sweep是线性的，因此可以从拍频中得出时间延迟，然后将延迟转换为距离。
% 在ACC设置中，雷达需要监视的最大范围约为200 m，并且系统需要能够区分相距1米的两个目标。 根据这些要求，可以计算波形参数。

fc = 77e9;
c = 3e8;
lambda = c/fc;%波长

%%
%可以根据信号传播明确的最大范围所需的时间来计算扫描时间。 
%通常，对于FMCW雷达系统，扫描时间至少应为往返时间的5至6倍。 本示例使用的系数为5.5。
range_max = 200;
tm = 5.5*range2time（range_maxc）;%range2time：将传播距离转换为传播时间

%%
%可以根据距离分辨率确定扫描带宽，并使用扫描带宽和扫描时间来计算扫描斜率。
range_res = 1;
bw = range2bw（range_resc）;
sweep_slope = bw/tm;   %斜率=带宽/扫描时间

%%
%由于FMCW信号通常会占用很大的带宽，因此将采样率盲目设置为带宽的两倍通常会增加A / D转换器硬件的压力。
%为了解决这个问题，通常可以选择较低的采样率。 这里可以考虑两件事：
% ＃对于复杂的采样信号，可以将采样率设置为与带宽相同。
% ＃FMCW雷达使用嵌入在去噪信号中的拍频（频差）估算目标距离。 
%雷达需要检测的最大拍频是对应于最大距离和最大多普勒频率的拍频之和。 因此，采样率仅需为最大拍频（差频）的两倍。
% 在此示例中，对应于最大距离的拍频由下式给出：

fr_max = range2beat（range_maxsweep_slopec）;
%range2beat将范围转换为拍频
% FB = range2beat（R，SLOPE）转换范围R（以米为单位）将线性FMCW信号消解到其相应范围拍频FB（以Hz为单位）。 FB与R具有相同的尺寸。

%%
% 另外，旅行车的最高速度约为230 km / h。 因此，最大多普勒频移和最大拍频可以计算为

v_max = 230*1000/3600;
fd_max = speed2dop（2*v_maxlambda）;%最大多普勒频移

fb_max = fr_max+fd_max; %最大拍频


% %
%此示例采用最大拍频和带宽两倍的采样率。
fs = max（2*fb_maxbw）;

%%
% The following table summarizes the radar parameters.
% 
%  System parameters            Value
%  ----------------------------------
%  Operating frequency （GHz）    77
%  Maximum target range （m）     200
%  Range resolution （m）         1
%  Maximum target speed （km/h）  230
%  Sweep time （microseconds）    7.33
%  Sweep bandwidth （MHz）        150
%  Maximum beat frequency （MHz） 27.30
%  Sample rate （MHz）          fs=150
% 下表总结了雷达参数。
% 
% ? 系统参数                 值
% ? ----------------------------------
% ? 工作频率（GHz）          fc=77
% ? 最大目标距离（米）       range_max=200
% ? 距离分辨率（m）          1
% ? 最高目标速度（km / h）   v_max=230km/h=63.8889m/s
% ? 扫描时间（微秒）         7.33
% ? 扫描带宽（MHz）          bw=150e6
% ? 最大拍频（MHz）         fb_max=27.30e6
% ? 采样率（MHz）           fs=150e6

%%
% 利用以上所有信息，您可以设置雷达系统中使用的FMCW波形。

hwav = phased.FMCWWaveform（‘SweepTime‘tm‘SweepBandwidth‘bw...
    ‘SampleRate‘fs）;

%%
%这是上扫线性FMCW信号，通常称为锯齿形。 可以检查所生成信号的时频图。
s = step（hwav）;
figure;
subplot（211）; plot（0:1/fs:tm-1/fsreal（s））;
xlabel（‘Time （s）‘）; ylabel（‘Amplitude （v）‘）;
title（‘FMCW signal‘）; axis tight;
subplot（212）; spectrogram（s321632fs‘yaxis‘）;
title（‘FMCW signal spectrogram‘）;
hold on;

%% Target Model目标模型
% ACC雷达的目标通常是前方的汽车。 本示例假定目标汽车在带雷达的情况下在x 轴以96 km / h的速度向前移动50 m。
% 根据图[1]，可以根据雷达与目标汽车之间的距离来计算汽车的雷达横截面。
car_dist = 50;
car_spee