Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Implementation of all steps of Music algorithm
Replies: 1   Last Post: Oct 1, 2013 2:59 AM

 Messages: [ Previous | Next ]
 Ste_ee Posts: 7 Registered: 9/20/13
Implementation of all steps of Music algorithm
Posted: Sep 30, 2013 2:47 PM

HI, i'm implementing all the steps of music algorithm, but i believe there are some problems.

This is m script:

N= 50; % number of antenna
K= 6; %number of signal
f0=10; %frequency of signal
fs=100; %sampling frequency

d = 0.5; %distance/lambda ratio

t = 0:1/fs:1;

x1= sin(2*pi*3*f0*t);
x2= 2*sin(2*pi*f0*t);
x3 = cos(2*pi*f0*t);

x = x1+x2+x3 + wgn(1,length(t),0); % signal + white noise

[t_matrix Rx]= corrmtx(x,N-1,'mod'); %Rx is autocorrelation matrix NxN.

[a_vector a_value] = eig(Rx); % find eigenvalues and eigenvector

[a_val_sort index]= sort(diag(a_value),1,'descend');
a_vett_sort=a_vector (:,index);
Qn = a_vett_sort(:,K+1:N); %noise subspace

angles = -90:1:90;
steering = exp(-i*2*pi*d*(0:N-1)'*sind(angles)); %steering matrix

for m= 1:length(angles)
spectrum(m)= 1/(steering(:,m)'*Qn*Qn'*steering(:,m));
end
grid on;
plot(angles,abs(spectrum));

I'm not sure of value of K and its relation with signals x1,x2,x3.
The spectrum is very strange, i don't understand why it's simmetric!
I find always +/- angles.

If is use SVD insted eig?? Rx is always square, and i believe i don't need SVD.

Is there anyone who may help me?
Thanks very much to all!

Date Subject Author
9/30/13 Ste_ee
10/1/13 Ste_ee