dbtex8l.m

几个关于阵列天线的仿真程序

lambda=1;
nElem=8;
elem21=defant('isotropElem');
elem2=cell(1,nElem);
for n=1:nElem,elem2{n}=elem21;end%for n
radius=0.5*lambda;
angle=linspace(0,2*pi,nElem+1);
angle=angle(1:(length(angle)-1));
xPos=radius*cos(angle);
yPos=radius*sin(angle);
elemPos=[xPos;yPos];
plot(elemPos(1,:),elemPos(2,:),'o')
axis('square')
xlabel('x-position')
ylabel('y-position')
title('Element positions')
for n=1:nElem
    text(elemPos(1,n),elemPos(2,n),[' ',num2str(n)])
end%for n
ant2=defant('array',elemPos,[],elem2);
smplPoints=linspace(-pi,pi,321);
figure,spantpat3(ant2,...
   {'lambda',lambda,...
    'mainPointDoa',d2r(0),...
    'smplPoints',linspace(-pi,pi,361),...
    'taperType','uniform','taperParam',[1 1 1 1 1 1 1 1].'})
figure,spantpat3(ant2,...
   {'lambda',lambda,...
    'mainPointDoa',d2r(0),...
    'smplPoints',smplPoints,...
    'taperType','uniform','taperParam',[1 1 1 1 1 1 1 1].',...
    'plotType','polar'})
theta=d2r([77 100])';
phi=zeros(size(theta));
SNR=[10 10]';
alpha=d2r([0 -150])';
dalpha=d2r([34 -66])';
dist=Inf*ones(size(theta));
sig=compsim4(ant2,lambda,100,'rndnw',[theta,...
   phi,SNR,alpha,dalpha,dist,...
   eye(size(theta,1))],'rndnw',[]);
spect1=sdoaspc('music',sig,smplPoints,3);
spect2=sdoaspc('cbf',sig,smplPoints);
      figure,splot2(spect1,'',spect2)
legend('Music','Conventional beamforming')