2D MUSIC algorithm for range-azimuth FMCW data processing

Question

Shirleyuue Jiang 2019-8-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/477124-2d-music-algorithm-for-range-azimuth-fmcw-data-processing

评论： Emre Kurtoglu 2020-8-11

Hi, I want to draw a map of range-azimuth. I can achieve 1D music, but there is some problem with 2D music for range-azimuth.

%X_ is the signal.

J = fliplr(eye(m1*m2,m1*m2));

Cx_ = 1/(2*p1*p2)*(X_*X_' + J*(X_*X_')*J);

[EV,D] = eig(Cx_); %D for eigenvalues and EV for eigenvector

EVA_temp = diag(D)';

[EVA,I] = sort(EVA_temp);

Sort_EV = fliplr(EV(:,I)); %eig

% W_ = Sort_EV(:,K+1:end); %(m1*m2)*(m1*m2-K)

W_ = Sort_EV(:,4:end);

for r = 1:length(VecR)

for theta = 1:length(VecA)

a_1 = [];

for mm1 = 0:m1-1

tao_i(mm1+1) = 2/c*( VecR(r) + mm1*d*sin(VecA(theta)) );

aa = exp(1i*2*pi*tao_i(mm1+1)*(linspace(fc,fc+alpha*(m2-1)/fs,m2) ) ); %0-

a_1 = cat(2,a_1,aa);

end

a = a_1.';

Smusic(r,theta) = 1/(a'*W_*W_'*a);

end

S_ = Smusic/(max(max(Smusic)));

figure,imagesc(VecA*180/pi,VecR(1:50),log(abs(S_)));xlabel('Azimuth(deg)');ylabel('Range(m)');title('range-azimuth');

% figure,imagesc(abs(Smusic));

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Shirleyuue Jiang 2019-8-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/477124-2d-music-algorithm-for-range-azimuth-fmcw-data-processing#answer_388587

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"). Is there any problem with "aa"?

5 个评论
显示 3更早的评论隐藏 3更早的评论

Kun-Lin Hsieh 2020-7-23

Hi Shirleyuue Jiang

I'm working on implementing this algorithm as well, and trying to extend it into a 3D problem by considering elevation dimention. However, there are some problem I met, and I'm even not sure if I have completely realize the whole procedure of this algorithm. Is it possible to share your code kindly if you have sucessfully completed it? We could also share our results to each other.

Emre Kurtoglu 2020-8-11

Hi all,

Is there any update about this implementation? I am able to implement 1D MUSIC but cannot get results using this method and it takes infinitely long time to process since there are three consecutive 'for' loops.

请先登录，再进行评论。