MATLAB Answers

0

Spatial smoothing - 2D MUSIC [FMCW MIMO radar processing]

Asked by Luca Romano on 18 Feb 2019
Latest activity Commented on by mk14
on 24 Aug 2019 at 12:59
Hi everyone,
I'm using a 2D Music algorithm for the estimation of Range-Azimuth info in a ULA FMCW MIMO radar.
The algorithm procedure is pretty simple (see "Belfiori F., Application of 2D MUSIC Algorithm to Range-Azimuth FMCW Radar data" or " Manokhin G., MUSIC-based algorithm for range-azimuth FMCW radar data processing without estimating number of targets"). I have some doubt with regard the application of the 2D FB spatial smoothing in order to decorrelate the covariance matrix of the data; I'm wondering if I do some error in the following code section (in the way I order the x_win_new vector):
x= DataV'; %Data collected by the ULA (31 Spatial samples-31Antennas), 501 time samples
N = size(x,1); % Number of antennas
M = size(x,2); % Number of samples in fast-time dimension
% window dimension for the smoothing dim[m1xm2]
m1 = 25;
m2 = 300;
p1 = N-m1+1; %positions in the spatial dimension
p2 = M-m2+1; %positions in the time dimension
N_submat = p1*p2; %number of all the possible scan
M_tot = m1*m2;
ind_tot=1;
for p1_ind = 1:p1
for p2_ind = 1:p2
x_win_new(:,ind_tot) = reshape (x(p1_ind:m1+(p1_ind-1),p2_ind:m2+(p2_ind-1)),[1,M_tot])';
ind_tot=ind_tot+1;
end
end
X_win_cor = x_win_new*x_win_new';
X_win_cor_conj = conj(X_win_cor);
J= fliplr(eye(M_tot)); %transition matrix
C=(1/2*N_submat)*(X_win_cor + J*X_win_cor_conj*J); %data smoothed covariance matrix (forward+backward)
so in my case I'm scanning the full 2D [NxM] matrix in this way:
x x x . . x x x . . . . . . . .
x x x . . x x x x x x . . x x x
x x x . . x x x x x x . . x x x
. . . . . . . . x x x . . x x x
If I select m1=M and m2=N (so no spatial smoothing, only forward/backward averaging) I get good results as shown in the attached image.
Music2D_1.png
Thanks a lot in advance,
Best regards,
Luca,

  2 Comments

Hi there, did you manage to solve the problem? I am also working on 2D music.
Hi, I have some problem with 2D music( Range-Azimuth), can you show me your whole code?

Sign in to comment.

2 Answers

Answer by mk14
on 5 Aug 2019

Hi there, can you please share your whole code. I cannot plot the data in the right way. Thanks!
Mark

  3 Comments

Hi, I can imagesc the data, but the result is not so good. Can we share the code with each other?
mk14
on 23 Aug 2019 at 12:12
I moved on to the capon beamformer instead of MUSIC.
So in your experiment, do capon be better than MUSIC?

Sign in to comment.


Answer by Shirleyuue Jiang on 24 Aug 2019 at 10:48

Hi, I have some problem with 2D music( Range-Azimuth), can you show me your whole code?
Best regards,
Shirley

  1 Comment

mk14
on 24 Aug 2019 at 12:59
For practical implementation, I think MUSIC is prone to interference. Moreover, when using capon you don't need to know number of source.

Sign in to comment.