Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
NCTM or The Math Forum.


Math Forum
»
Discussions
»
Software
»
comp.softsys.matlab
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



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,N1,'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:N1)'*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!



