Radar Signal Simulation: range-angle response

1 次查看（过去 30 天）

显示更早的评论

Wei-Min Chen 2018-12-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/435298-radar-signal-simulation-range-angle-response

编辑： Wei-Min Chen 2018-12-12

在 MATLAB Online 中打开

I'm going to get range-angle response of example in ACC.

I know that there's only single sensor in RX in this example, what I want to do is replacing single sensor to ULA.

Then the size of response I get should be like:4*512 (for example)

4 is number ULA elements and 512 is fast time samples.

I studied FMCW DSP before but I still don't know how to deal with so many system objects in Matlab.

I really need an example or a tutorial.

Thanks!

%% Radar System Setup
ant_aperture = 6.06e-4;                         % in square meter
ant_gain = aperture2gain(ant_aperture,lambda);  % in dB
tx_ppower = db2pow(5)*1e-3;                     % in watts
tx_gain = 9+ant_gain;                           % in dB
rx_gain = 15+ant_gain;                          % in dB
rx_nf = 4.5;                                    % in dB
transmitter = phased.Transmitter('PeakPower',tx_ppower,'Gain',tx_gain);
receiver = phased.ReceiverPreamp('Gain',rx_gain,'NoiseFigure',rx_nf,'SampleRate',fs);
sULA = phased.ULA('NumElements',4,'ElementSpacing',lambda/2);
radar_speed = 0;
radarmotion = phased.Platform('InitialPosition',[0;0;0],'Velocity',[radar_speed;0;0]);
%% Radar Signal Simulation
% specanalyzer = dsp.SpectrumAnalyzer('SampleRate',fs,...
%     'PlotAsTwoSidedSpectrum',true,...
%     'Title','Spectrum for received and dechirped signal',...
%     'ShowLegend',true);
rng(2012);
for azimuth = 1:180
xr = complex(zeros(4,waveform.SampleRate*waveform.SweepTime));
% Update radar and target positions
[radar_pos,radar_vel] = radarmotion(waveform.SweepTime);
[tgt_pos,tgt_vel] = carmotion(waveform.SweepTime);
% Transmit FMCW waveform
sig = waveform();
txsig = transmitter(sig);
% Propagate the signal and reflect off the target
txsig = channel(txsig,radar_pos,tgt_pos,radar_vel,tgt_vel);
txsig = cartarget(txsig);
% Dechirp the received radar return
txsig = collectPlaneWave(sULA,sig,[azimuth;0],fc);
% txsig = receiver(txsig);
dechirpsig = dechirp(txsig,sig);
% Visualize the spectrum
% specanalyzer([txsig dechirpsig]);
xr(:,azimuth) = dechirpsig;
azimuth
end