Código fuente de Paladin Ethernet
% clc
Echo
m = 6; porcentaje del número de antenas
d = 0,5
k = 2 * pi
KD = k * d;
k = 100; La amplitud de un código Walsh binario con % de muestras finitas es 1.
d = 3; % número de señal
th=[-5/180*pi, 5/180*pi, 20/180 * pi];
s=sign(randn(D,K));% de muestreo de tiempo de señal
SIG 2 = 0,1;
n=sqrt(sig2)*randn(M,K) ; % ruido tiempo muestreo
θ=-π/6:0.01:π/6;
a =[];
a 1 =[] ;
a2 =[];
a3 =[];
Para ii=1:M
a 1 =[a 1; exp(j*(ii-1)*kd.* sin(th(1)))];
a2 =[a2; exp(j*(ii-1)*kd.* sin( th(2)))];
a3 =[a3;exp(j*(ii-1)*kd.* sin(th(3)))];
Fin
A=[a1, a2, a3];
%randn('estado ', 0);
Rss=s*s' /K
Clasificación (Rss)
rnn = n * n '/K
rns = n * s '/K;
rsn = s * n '/K;
rxx = A * Rss * A '+A * Rsn+Rns * A '+Rnn;
[V, Dia] = EIG (Rxx);
[Y, Índice]= sort(diag(Dia));
EN=V(:, Índice(1:M-D));
p>Para jj=1:18000
theta(jj)=-pi/2+jj/18000 * pi;
a =[];
Para kk=1:M
a =[a; exp(j *(kk-1)* KD * sin(theta(jj)))];
Fin
PB(jj)= ABS(real(1/ABS(a ' * EN * EN ' * a)));
Fin
PB _ dB = 10 * log 10(PB);
plot(theta/pi*180, PB_dB, ' r ');
Cuadrícula activada
xlabel('angle')
ylabel(' | P(\theta)|(dB)');
Title(['Uso del pseudoespectro musical promediado en el tiempo'] ) ;
[max, index]=maximun(PB,D);
eth =-9index/100;
Disp(['Música método (Señal relacionada):']);
Disp(['1, estimación de DOA:']);
disp([blanks(10),' d= ', num2str (d), 'lambda ']);
Disp([blanks(10)), 'El valor teórico es: theta _ 1 = ', num 2 str(th(1)/pi * 180) ,', theta_2 = '.
Disp([blanks(10)), 'La estimación de la simulación es: theta _ 1xzxx = ', num2str (eth (1)), ', ' theta _ 2 = ', num2str (eth (2) ))
Disp([blanks(10),'bias:',num 2 str(ABS(eth(1)-th(1)/PI * 180)),',',num2str ( ABS (ABS), ', num 2 str(ABS(eth(3)-th(3)/pi * 180), ' ']);
Disp(['2, estimación del número de señal: ']);
Disp([blanks(10),'Valor teórico: D = ',num2 str(D)]);
[V_s,Dia_s] = EIG( Rss);
Disp([blanks(10), 'La estimación de la simulación es: ed = ', num 2 str(length(Diag(dia _ s))]);