Please see the following two links in order to have an idea:

https://ch.mathworks.com/help/phased/examples/high-resolution-direction-of-arrival-estimation.html?s_tid=gn_loc_drop

https://ch.mathworks.com/help/comm/ref/rainpl.html

hi

Maybe i will be able to help you, please send your problem description . I worked on DOA estimation.

You will find a very good documentation in the help of Matlab. But if you have more specific questions, you can send it to me by email.

MUSIC and ESPRIT based on Rainfall model work as follows:

% MUSIC algorithm

[Rxx,~] = corrmtx(X',N-1,'modified'); % Correlation matrix

[V,D] = eig(Rxx); % Eigen decomposition

[d,I] = sort(diag(D),'descend');

V = V(:,I);

noise_space = V(:,M+1:N);

theta_music = (-90:0.5:90);

for i = 1:length(theta_music)

a = exp(1j*2*pi*d*sin(theta_music(i)*pi/180)*[0:N-1]');

P_music(i) = 1/(a'*noise_space*noise_space'*a);

end

P_music = 10*log10(abs(P_music)/max(abs(P_music)));

% ESPRIT algorithm

[U,S,V] = svd(X);

E = U(:,1:M)*U(:,2*M+1:2*M+N)';

R = E*E'/T;

[V,D] = eig(R);

[d,I] = sort(diag(D),'descend');

V = V(:,I);

theta_esprit = (-90:0.5:90);

for i = 1:length(theta_esprit)

a = exp(1j*2*pi*d*sin(theta_esprit(i)*pi/180)*[0:N-1]');

P_esprit(i) = 1/(a'*V(:,1:M)*V(:,1:M)'*a);

end

P_esprit = 10*log10(abs(P_esprit)/max(abs(P_esprit)));

% Plot results

figure;

plot(theta_music,P_music,'-b');

hold on;

plot(theta_esprit,P_esprit,'-r');

xlabel('Angle (degrees)');

ylabel('Normalized Power (dB)');

title('DOA Estimation with MUSIC and ESPRIT');

legend('MUSIC','ESPRIT');