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

1、111rI-:IIUIAIII If I /II 電/II f /I/ / 丿ICapon beamforming功率譜5060708090 00O10Capon波束形成陣列 N=16,信號(hào) 0. = -30° ,干擾為 q=60°, q=19°, 2=45。, 干擾功率分別為:40dB, 35dB, 50dB。Capon波束形成后 的方向圖和功率譜如下：-60-40-20020406080100theta/o為了比較接收數(shù)據(jù)直接估計(jì)噪聲協(xié)方差矩陣和利用丁擾+噪聲估計(jì)協(xié)方差矩陣的Capon 波束形成的差異，進(jìn)行如卞仿貞:102030 40506070807 90O

2、10%陣元數(shù)%陣元間距與波長的關(guān)系%范圍%來波方向%干擾方向1%干擾方向2%干擾方向3%采樣單元數(shù)可以看出利用擾+噪聲估計(jì)協(xié)方差矩陣的方向圖性能較優(yōu)于接收數(shù)拯直接估計(jì)噪聲協(xié)方差矩陣的方向圖。代碼：clc;clear all;close all: iiiia=sqrt(-l);element_num=8; d_lamda=l/2; tlieta=-90:0.5:90; tlieta0=-30;tlietal=60;tlieta2=19;tlieta3=45;L=1000: for i=l:L;ampO=lO*randn(l);%信號(hào)的幅度隨機(jī)產(chǎn)生，保證信號(hào)之間是不相關(guān)的amp 1=1 OO*ra

3、ndn( 1):%輸入陣列的噪聲amp2=sqrt( 10A3.5 ) *raiidn( 1);%輸入陣列的噪聲amp3=sqrt( 10A5)*randii( 1 入陣列的噪聲ampn=3;% 噪聲 x(:,i)=amp0*exp(ima*2*pi*l/2*sin(tlieta0*pi/180)*0:element_num-l,)-r.an4)l*exp(ima*2*pi*l/2*sin(tlietal*pi/180)*0:element_niuii-lr)+.aiiip2*exp(iina*2*pi*l/2*sin(theta2*pi/180)*0:element_muii-lr)+.ai

4、i3*exp(iina*2*pi*l/2*sin(theta3*pi/180)*0:element_niim-lr)+. aiiipn*(raiidn(element_num, l)+iina*randii(element_nunL 1);endRx=l/*L*x* x：R=uiv(Rx);steer=exp(ima *2 *pi* 1 /2 *sin(thetaO*pi/z 180)*0: element_num-1 *);w=R*steer/(steer,*R*steer);%Capon 最優(yōu)權(quán)矢量for j=l:length(theta);anexpqimaPl產(chǎn) dlamdarinQh

5、etaGpi/lgOPeleineiHnumiy); f(j)=w，*a;p(j)=l/(廣RS);endF=20*logl O(abs(f)/(max(abs(f);P=20*logl O(abs(p)/(max(abs(p) ;%此處是功率的對數(shù)形式figure;% subplot(121)plot(tlieta.F),grid on.hold onplot(theta0?-80:0,'.')plot(theta 1 ,-80:0/.')plot(theta2,-80:0J. *)plot(theta3 ,-80:0/. *)xlabelCtheta/o1);ylabel(T/dB*);titl已Capon beainfomiing 方向圖')% axis(-90 90 -50 0);% subplot(122)figure;plot