Capon-波束形成matlab仿真(附源代碼)

1、精選優(yōu)質文檔-傾情為你奉上精選優(yōu)質文檔-傾情為你奉上專心-專注-專業(yè)專心-專注-專業(yè)精選優(yōu)質文檔-傾情為你奉上專心-專注-專業(yè)Capon波束形成陣列N=16, 信號，干擾為，干擾功率分別為：40dB,35dB,50dB。Capon波束形成后的方向圖和功率譜如下：為了比較接收數據直接估計噪聲協(xié)方差矩陣和利用干擾+噪聲估計協(xié)方差矩陣的Capon波束形成的差異，進行如下仿真：可以看出利用干擾+噪聲估計協(xié)方差矩陣的方向圖性能較優(yōu)于接收數據直接估計噪聲協(xié)方差矩陣的方向圖。代碼：clc;clear all ;close all;ima=sqrt(-1);element_num=8; %陣元數d_lamda

2、=1/2; %陣元間距與波長的關系theta=-90:0.5:90; %范圍theta0=-30; %來波方向theta1=60; %干擾方向1theta2=19; %干擾方向2theta3=45; %干擾方向3L=1000; %采樣單元數for i=1:L; amp0=10*randn(1);%信號的幅度隨機產生，保證信號之間是不相關的 amp1=100*randn(1);%輸入陣列的噪聲 amp2=sqrt(103.5)*randn(1);%輸入陣列的噪聲 amp3=sqrt(105)*randn(1);%輸入陣列的噪聲 ampn=3;%噪聲 x(:,i)=amp0*exp(ima*2*p

3、i*1/2*sin(theta0*pi/180)*0:element_num-1)+. amp1*exp(ima*2*pi*1/2*sin(theta1*pi/180)*0:element_num-1)+. amp2*exp(ima*2*pi*1/2*sin(theta2*pi/180)*0:element_num-1)+. amp3*exp(ima*2*pi*1/2*sin(theta3*pi/180)*0:element_num-1)+. ampn*(randn(element_num,1)+ima*randn(element_num,1);endRx=1/L*x* x;R=inv(Rx)

4、;steer=exp(ima*2*pi*1/2*sin(theta0*pi/180)*0:element_num-1);w=R*steer/(steer*R*steer);%Capon最優(yōu)權矢量for j=1:length(theta); a=exp(ima*2*pi*d_lamda*sin(theta(j)*pi/180)*0:element_num-1); f(j)=w*a; p(j)=1/(a*R*a); endF=20*log10(abs(f)/(max(abs(f);P=20*log10(abs(p)/(max(abs(p);%此處是功率的對數形式figure;% subplot(121)plot(theta,F),grid on,hold onplot(theta0,-80:0,.)plot(theta1,-80:0,.)plot(theta2,-80:0,.)plot(theta3,-80:0,.)xlabel(theta/o);ylabel(F/dB);title(Capon beamforming方向圖)% axis(-90 90 -50 0);% subplot(122)figure;plot(