天線線列陣方向圖

1、.word格式,陣列方向圖及MATLA昉真1、線陣的方向圖2R =2+2cos -n cos®)MATLA隈序如下(2元)：clear;a=0:0.1:2*pi;y=sqrt(2+2*cos(pi-pi*cos(a);polar(a,y);圖形如下：專業(yè).專注(g1ej'，,gMeja )加以組合的話,若陣元間距為半波長的M個(gè)陣元的輸出用方向向量權(quán)重陣列的方向圖為R( ) ' gmejmWoS( ) m 1MATLA翼序如下(10個(gè)陣元)： clear;f=3e10;lamda=(3e8)/f; beta=2.*pi/lamda;n=10;t=0:0.01:2*pi;

2、d=lamda/4;W=beta.*d.*cos(t);z1=(n/2).*W)-n/2*beta* d; z2=(1/2).*W)-1/2*beta* d; F1=sin(z1)./(n.*sin(z2);i K1=abs(F1);polar(t,K1);方向圖如下：2、圓陣方向圖程序如下：clc;clear all ;close all;M = 16; %行陣元數(shù)k = 0.8090; % k = r/lambdaDOA_theta = 90; % 方位角DOA_fi = 0;% 俯仰角%形成方位角為theta ,俯仰角位fi的波束的權(quán)值m = 0 : M-1;w = exp(-j*2*p

3、i*k*cos(2*pi*m7M-DOA_theta*pi/180)*cos(DOA_fi*pi/180);%wexp(-j*2*pi*k*(cos(2*pi*m7M)*cos(DOA_theta*pi/180)*cos(DOA_fi*pi/180)+sin(2*pi*m'/M)*sin(DOA_fi*pi/180); %豎直放置% w = chebwin(M, 20) .* w; %行加切比雪夫權(quán)%繪制水平面放置的均勻圓陣的方向圖theta = linspace(0,180,360);fi = linspace(0,90,180);for i_theta = 1 : length(t

4、heta)for i_fi = 1 : length(fi) aexp(-j*2*pi*k*cos(2*pi*m7M-theta(i_theta)*pi/180)*cos(fi(i_fi)*pi/180);%a=exp(-j*2*pi*k*(cos(2*pi*m7M)*cos(theta(i_theta)*pi/180)*cos(fi(i_fi)*pi/180)+sin(2*pi*m'/M)*sin(fi(i_fi)*pi/180); %豎直放置Y(i_theta,i_fi) = w'*a;endendY= abs(Y); Y = Y/max(max(Y);Y = 20*log

5、10(Y);% Y = (Y+20) .* (Y+20)>0) - 20; %切圖Z = Y + 20;Z = Z .* (Z > 0);Y = Z - 20;figure; mesh(fi, theta, Y);view(66, 33);title('水平放置時(shí)的均勻圓陣方向圖');% title('豎面放置時(shí)的均勻圓陣方向圖');%豎直放置axis(0 90 0 180 -20 0);xlabel(' 俯仰角 /(circ)'); ylabel('方位角 /(circ)'); zlabel('P/dB

6、9;);figure; contour(fi, theta, Y);方向圖如下：水平放置時(shí)的均勻圓陣方向圖方位角用3、平面陣方向圖:clc;clear all ；close all;Row_N = 16; % 行陣元數(shù)Col_N = 16; % 列陣元數(shù) k = 0.5; % k = d/lambda DOA_theta = 90; % 方位角 DOA_fi = 0;% 俯仰角%形成方位角為theta ,俯仰角位fi的波束的權(quán)值Row_n = 0 : Row_N-1; Col_n = 0 : Col_N-1;W_Row = exp(-j*2*pi*k*Row_n'*cos(DOA_th

7、eta*pi/180)*cos(DOA_fi*pi/180);W_Col = exp(-j*2*pi*k*Col_n'*sin(DOA_theta*pi/180)*cos(DOA_fi*pi/180);行加切比雪夫權(quán)列加切比雪夫權(quán)合成的權(quán)值 N*N x 1% W_Col = exp(-j*2*pi*k*Col_n'*sin(DOA_fi*pi/180); %豎直放置W_Row = chebwin(Row_N, 20) .* W_Row; %W_Col = chebwin(Col_N, 30) .* W_Col; %W = kron(W_Row, W_Col);%繪制水平面放置的

8、平面陣的方向圖theta = linspace(0,180,180);fi = linspace(0,90,90);for i_theta = 1 : length(theta)for i_fi = 1 : length(fi)row_temp行導(dǎo)向列導(dǎo)向豎直放置合成的導(dǎo)向矢量exp(-j*2*pi*k*Row_n'*cos(theta(i_theta)*pi/180)*cos(fi(i_fi)*pi/180); % 矢量N x 1col_tempexp(-j*2*pi*k*Col_n'*sin(theta(i_theta)*pi/180)*cos(fi(i_fi)*pi/18

9、0); % 矢量N x 1% col_temp = exp(-j*2*pi*k*Col_n'*sin(fi(i_fi)*pi/180); % Y(i_theta,i_fi) = W*kron(row_temp, col_temp);%N*N x 1endendY= abs(Y); Y = Y/max(max(Y);Y = 20*log10(Y);Y = (Y+60) .* (Y+60)>0) - 60; %切圖% Z = Y + 60;% Z = Z .* (Z > 0);% Y = Z - 60;figure; mesh(fi, theta, Y);view(66, 33

10、);title('水平面放置時(shí)的面陣方向圖');axis(0 90 0 180 -60 0);xlabel('俯仰角 /(circ)'); ylabel('方位角(circ)'); zlabel('P/dB');figure; contour(fi, theta, Y);方向圖如下：水平面放置時(shí)的司庫方向圖4、CAPONS法波束形成MATLA翼序如下(陣元16,信號源3,快拍數(shù)1024): clear alli=sqrt(-1);j=i;M=16;%均勻線陣列數(shù)目P=3;%信號源數(shù)目fO=10;f1=50;f2=100;%信號頻率

11、nn=1024;%快拍數(shù)angle1=-15;angle2=15;angle3=30;%the signal angleth=angle1;angle2;angle3'SN1=10;SN2=10;SN3=10;%言噪比sn=SN1;SN2;SN3;degrad=pi/180;tt=0:.001:1024;x0=exp(-j*2*pi*f0*tt); %3 個(gè)信號 x0、 x1、 x2x1=exp(-j*2*pi*f1*tt); %x2=exp(-j*2*pi*f2*tt); %t=1:nn;S=x0(t);x1(t);x2(t);nr=randn(M,nn);ni=randn(M,nn

12、);u=nr+j*ni; % 復(fù)高斯白噪聲Ps=S*S'./nn; %信號能量ps=diag(Ps); refp=2*10.A(sn/10); tmp=sqrt(refp./ps);S2=diag(tmp)*S; % 加入噪聲 tmp=-j*pi*sin(th*degrad); tmp2=0:M-1'a2=tmp2*tmp;A=exp(a2);X=A*S2+.1*u; % 接收到的信號 Rxx=X*X'./nn; %相關(guān)矩陣invRxx=inv(Rxx);%t尋信號 th2=-90:90'tmp=-j*pi*sin(th2'*degrad);tmp2=0:M-1'a2=tmp2*tmp;A2=exp(a2); den=A2'*i