二维数字波束形成-理想阵列和非理想阵列

function radar = Radar_para_init_csdn(antena_method, phase_method) c0 = 3e8; % Speed of light in vacuum %% radar parameters radar.para.slope = 55.256e6;%79*10^6;%Hz/us radar.para.numsample = 64;% num of samples per chirp if antena_method == 1 radar.para.numtx = 4;% num of tx else radar.para.numtx = 3;% num of tx end radar.para.numrx = 4;% num of rx radar.para.numantena = radar.para.numtx*radar.para.numrx;% num of Antena radar.para.samplerate = 2;% Msps radar.para.duration = radar.para.numsample/radar.para.samplerate; % us radar.para.bandwidth = radar.para.slope*radar.para.duration; radar.para.fc = 61.2e9+radar.para.bandwidth/2; radar.para.lambda = c0/radar.para.fc; %% antena para ant_d = radar.para.lambda / 2; if antena_method == 1 % 理想16阵元 RxPos = [3,2;3,3;2,2;2,3]; TxPos = [2,2;0,2;2,0;0,0]; else % 实际3T4R RxPos = [1,0;1,1;0,0;0,1]; TxPos = [0,2;2,0;0,0]; end figure(); plot(TxPos(:,1),TxPos(:,2),'r*'); hold on; plot(RxPos(:,1),RxPos(:,2),'bo'); daspect([1 1 1]) [virRxPos_X,virTxPos_X] = ndgrid(RxPos(:,1),TxPos(:,1));% 这种排列导致数据先排接收天线，再排发射天线，又因为matlab列优先，所以是4*3的矩阵 [virRxPos_Z,virTxPos_Z] = ndgrid(RxPos(:,2),TxPos(:,2)); xPhaseShift = ant_d * (virRxPos_X(:) + virTxPos_X(:));%发射天线和接收天线横坐标计算的阵元间距,即虚拟阵元间距 zPhaseShift = ant_d * (virRxPos_Z(:) + virTxPos_Z(:));%发射天线和接收天线纵坐标计算的阵元间距,即虚拟阵元间距 %% aoa scan Nazi = 128; Nele = 128; radar.aziGrid = linspace(-90,90,Nazi); % 单位是角度？ 21份,方位角，间隔为5度 radar.eleGrid = linspace(-90,90,Nele); % 单位是角度？ 29份，俯仰角 45 [aziMesh,eleMesh] = meshgrid(radar.aziGrid, radar.eleGrid); %生成对应的感兴趣的二维平面 radar.dirctionGrid = [aziMesh(:),eleMesh(:)]; n_directions = size(radar.dirctionGrid,1);%共有289个网格点 n_channels = radar.para.numantena; radar.directBfGrid = zeros(n_channels, n_directions); for p=1:n_directions%共289个网格点，二维平面 if phase_method == 1 tot_shift = zPhaseShift * sind(radar.dirctionGrid(p,2)) * cosd(radar.dirctionGrid(p,1)) + xPhaseShift * sind(radar.dirctionGrid(p,1)); end s = exp(-2j*pi/radar.para.lambda * tot_shift); %标准bf radar.directBfGrid(:,p) = s;%补偿相位 end end