2D MUSIC algorithm for range-azimuth FMCW data processing
显示 更早的评论
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
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));
回答（1 个）
Shirleyuue Jiang
2019-8-22
0 个投票
5 个评论
mk14
2019-8-24
You can add my QQ if you want to discuss 81214332
Shirleyuue Jiang
2019-8-26
Patrick Nkam
2020-6-18
编辑：Patrick Nkam
2020-6-18
Hi Shirleyuue Jiang ,
I work already with MUSIC, and i have some problems to get the Ranges values. Can you help me?
I habe problem with the steering-vector aa
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.
类别
在 帮助中心 和 File Exchange 中查找有关 Particle & Nuclear Physics 的更多信息
2019-8-22
2020-8-11
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
