AoA.rar_AOA_AOA matlab_AOA 估计 MATLAB_天线波束AOA_线性 天线 信道入射角估计

clear clc theta1=10*pi/180;theta2=25*pi/180;theta3=60*pi/180; N=1000; bit1=2*rem(unidrnd(2,1,N),2)-1; bit2=2*rem(unidrnd(2,1,N),2)-1; bit3=2*rem(unidrnd(2,1,N),2)-1; sig=[bit1;bit2;bit3]; SNRdB=10; n0=1/(10^(SNRdB/10)); sigma=sqrt(n0); noise00=randn(1,N); noise0=sigma*noise00; noise11=randn(1,N); noise1=sigma*noise11; noise22=randn(1,N); noise2=sigma*noise22; noise33=randn(1,N); noise3=sigma*noise33; noise44=randn(1,N); noise4=sigma*noise44; noise55=randn(1,N); noise5=sigma*noise55; noise66=randn(1,N); noise6=sigma*noise66; noise77=randn(1,N); noise7=sigma*noise77; noise=[noise0;noise1;noise2;noise3;noise4;noise5;noise6;noise7]; a1=[1; exp(-j*pi*sin(theta1)); exp(-j*2*pi*sin(theta1)); exp(-j*3*pi*sin(theta1)); exp(-j*4*pi*sin(theta1)); exp(-j*5*pi*sin(theta1)); exp(-j*6*pi*sin(theta1)); exp(-j*7*pi*sin(theta1))]; a2=[1; exp(-j*pi*sin(theta2)); exp(-j*2*pi*sin(theta2)); exp(-j*3*pi*sin(theta2)); exp(-j*4*pi*sin(theta2)); exp(-j*5*pi*sin(theta2)); exp(-j*6*pi*sin(theta2)); exp(-j*7*pi*sin(theta2))]; a3=[1; exp(-j*pi*sin(theta3)); exp(-j*2*pi*sin(theta3)); exp(-j*3*pi*sin(theta3)); exp(-j*4*pi*sin(theta3)); exp(-j*5*pi*sin(theta3)); exp(-j*6*pi*sin(theta3)); exp(-j*7*pi*sin(theta3))]; Rx=zeros(8,8); p=zeros(1,90); for m=1:1000 n=diag([noise0(m)*noise0(m),noise1(m)*noise1(m),noise2(m)*noise2(m),noise3(m)*noise3(m),... noise4(m)*noise4(m),noise5(m)*noise5(m),noise6(m)*noise6(m),noise7(m)*noise7(m)]); x=a1*sig(1,m)+a2*sig(2,m)+a3*sig(3,m)+noise(:,m); xx=x*x'; Rx=((m-1)*Rx+xx)/m; end angle=0:0.1:90; for i=1:length(angle) theta=((i-1)/10)*pi/180; Atheta=[1,exp(-j*pi*sin(theta)),exp(-j*2*pi*sin(theta)),exp(-j*3*pi*sin(theta)),... exp(-j*4*pi*sin(theta)),exp(-j*5*pi*sin(theta)),exp(-j*6*pi*sin(theta)),exp(-j*7*pi*sin(theta))].'; Pcapon=Atheta'*Rx*Atheta; p(i)=sqrt(Pcapon*Pcapon'); end p=p/max(p); hold on plot(angle,p) xlabel('Angle') ylabel('spectrum')