DOA估计基于matlab常规波束和CAPON算法DOA估计【含Matlab译码 2661期】.zip

%MUSIC ALOGRITHM %DOA ESTIMATION BY CLASSICAL_MUSIC clear all; %close all; clc; source_number=3;%信元数 sensor_number=16;%阵元数 N_x=1024; %信号长度 snapshot_number=N_x;%快拍数 w=[pi/4 pi/6 pi/3].';%信号频率 l=sum(2*pi*3e8./w)/3;%信号波长 d=0.5*l;%阵元间距 snr=15;%信噪比 source_doa=[-45 0 60];%两个信号的入射角度 A=[exp(-1j*(0:sensor_number-1)*d*2*pi*sin(source_doa(1)*pi/180)/l);exp(-1j*(0:sensor_number-1)*d*2*pi*sin(source_doa(2)*pi/180)/l);exp(-1j*(0:sensor_number-1)*d*2*pi*sin(source_doa(3)*pi/180)/l)].';%阵列流型 s=sqrt(10.^(snr/10))*exp(1j*w*[0:N_x-1]);%仿真信号 %x=awgn(s,snr); x=A*s+(1/sqrt(2))*(randn(sensor_number,N_x)+1j*randn(sensor_number,N_x));%加了高斯白噪声后的阵列接收信号 R=x*x'/snapshot_number; iR=inv(R); %[V,D]=eig(R); %Un=V(:,1:sensor_number-source_number); %Gn=Un*Un'; [U,S,V]=svd(R); Un=U(:,source_number+1:sensor_number); Gn=Un*Un'; searching_doa=-90:0.1:90;%线阵的搜索范围为-90~90度 for i=1:length(searching_doa) a_theta=exp(-1j*(0:sensor_number-1)'*2*pi*d*sin(pi*searching_doa(i)/180)/l); Pmusic(i)=a_theta'*a_theta./abs((a_theta)'*Gn*a_theta); Pcapon(i)=1./abs((a_theta)'*iR*a_theta); end plot(searching_doa,10*log10(Pmusic/max(Pmusic)),'k-',searching_doa,10*log10(Pcapon/max(Pcapon)),'b--'); %axis([-90 90 -90 15]); xlabel('DOAs/degree'); ylabel('Normalized Spectrum/dB'); legend('Music Spectrum','Capon Spectrum'); title('Comparation of MUSIC and Capon for DOA Estimation'); grid on;