FMCW MIMO雷达距离角度扇形谱图生成仿真代码

clc;clear;close all;
%% 目标和系统参数更新disp(strcat(['=====','设置目标和系统参数','====='])); % 单板雷达模式tarOut = ConfigureTarget;     % 生成目标信息cfgOut = ConfigureParameter(tarOut);  % 生成毫米波雷达系统参数
%% 雷达采集数据disp(strcat(['=====','雷达数据生成','=====']));
[RawData] = GenerateSigIQ(tarOut, cfgOut);
%% 物理参数更新c = physconst('LightSpeed'); % 光速fc = cfgOut.fc; % 载频 Hzlambda = c / fc; % 波长
ADCNum = cfgOut.ADCNum; % ADC采样数目ChirpNum = cfgOut.ChirpNum; % 每帧发射Chirp数目
numTx = cfgOut.numTx; % 发射天线数目numRx = cfgOut.numRx; % 接收天线数目% arr = cfgOut.array; % 阵元排列[初阶版本的阵列]virtual_array = cfgOut.virtual_array; % 虚拟阵列structarr = virtual_array.virtual_arr;
arrNum = numTx * numRx; % 阵元数目arrDx = cfgOut.arrdx; % 方位向阵元间距arrDy = cfgOut.arrdx; % 俯仰向阵元间距
validB = cfgOut.validB; % 有效带宽range_res = c / (2 * validB); % 距离分辨率
TF = cfgOut.Tc * (numTx * ChirpNum);doppler_res = lambda / (2 * TF); % 多普勒分辨率Frame = cfgOut.Frame; % 帧数
% 距离索引和速度索引velocityIndex = (-ChirpNum / 2 : 1 : ChirpNum / 2 - 1) * doppler_res;rangeIndex = (0 : ADCNum - 1) * range_res;

RawData = reshape(RawData, ADCNum, ChirpNum, [], arrNum);disp(strcat(['雷达信号的维度为：',num2str(size(RawData))]));
for frame_id = 1 : Frame    adcData = squeeze(RawData(:, :, frame_id, :));    %% 距离维FFT    disp(strcat(['=====','Range-Profile生成','=====']));    fftOut1 = rangeFFT(adcData);    rangeFFTOut = fftOut1.rangeFFT;    rangeFFTOutX = squeeze(rangeFFTOut(:,1,1));    [~,rangeFFTOutIndex] = max(rangeFFTOutX);
    figure(1);    set(gcf,'unit','centimeters','position',[20,12,10,10])    plot(rangeIndex, 10*log10(abs(rangeFFTOut(:,1,1))));    xlabel('Range(m)');ylabel('Amplitude(dB)');    title(strcat(['第',num2str(frame_id),'帧-目标距离分布']));grid minor;        scale = -90:0.5:90;   %扫面范围    %% 距离-方位图生成    X = [squeeze(rangeFFTOut(:,1,1:4)),squeeze(rangeFFTOut(:,1,5:8))];    for m=1:size(rangeFFTOut,1)        RAM(m,:) = fftshift(fft(X(m,:),361));    end    figure(2);    subplot(121);    set(gcf,'unit','centimeters','position',[30,12,10,10])    [theta,rho] = meshgrid(scale, rangeIndex); % 网格化    xaxis = rho .* cosd(theta); % 横坐标    yaxis = rho .* sind(theta); % 纵坐标    surf(yaxis,xaxis,abs(RAM)*1e3,'EdgeColor','none');    view(2);colormap('jet');    xlabel('X-Range(m)','fontsize',15,'fontname','Times New Roman');ylabel('Y-Range(m)','fontsize',15,'fontname','Times New Roman');grid on; axis xy    title(strcat(['第',num2str(frame_id),'帧-距离方位FFT谱图']));grid minor;    set(gca,'GridLineStyle','- -');    set(gca,'GridAlpha',0.2);    set(gca,'LineWidth',1.5);    set(gca,'xminortick','on');    set(gca,'ygrid','on','GridColor',[0 0 0]);    colorbar;
    %%多目标定位：CaponBF测角算法    M=8;           %阵元单元    d = 0:lambda/2:(M-1)*lambda/2;%阵列天线    N = size(rangeFFTOut,1);    Rx  =rangeFFTOut(:,:,1:8);    for r =1:size(Rx,1)        Y =squeeze( Rx(r,:,:)).';        R = (Y*Y')/size(Y,2);        R_inv = pinv(R);        idx = 1;        for i=scale            steeringVector  = (exp(1j*2*pi/lambda*d*sind(i))).'; %转向矢量转置            P_Capon(idx,r) = 1/(steeringVector'*R_inv*steeringVector); %Capon算法            idx = idx + 1;        end    end    subplot(122);    set(gcf,'unit','centimeters','position',[30,12,10,10])    [theta,rho] = meshgrid(scale, rangeIndex); % 网格化    xaxis = rho .* cosd(theta); % 横坐标    yaxis = rho .* sind(theta); % 纵坐标    surf(yaxis,xaxis,10*log10(abs(P_Capon.')),'EdgeColor','none');    view(2);colormap('jet');    xlabel('X-Range(m)','fontsize',15,'fontname','Times New Roman');ylabel('Y-Range(m)','fontsize',15,'fontname','Times New Roman');grid on; axis xy    title(strcat(['第',num2str(frame_id),'帧-距离方位CaponBF谱图']));grid minor;    set(gca,'GridLineStyle','- -');    set(gca,'GridAlpha',0.2);    set(gca,'LineWidth',1.5);    set(gca,'xminortick','on');    set(gca,'ygrid','on','GridColor',[0 0 0]);    colorbar;end