FMCW雷达仿真：基于L形阵列4D点云获取

Poulen

已于 2024-04-07 21:34:37 修改

阅读量1.4k

点赞数 17

分类专栏： MMWave案例解析文章标签： 4-D点云 FMCW Radar L型阵列

于 2023-12-21 19:08:51 首次发布

本文链接：https://blog.csdn.net/Poulen/article/details/135135494

版权

MMWave案例解析专栏收录该内容

21 篇文章 75 订阅

订阅专栏

摘要：本期内容为3D点云目标获取的延续工作，在距离、速度、方位角估计的基础上，通过设计L型阵列结构，进一步实现目标俯仰角的估计，最终实现目标4-D点云的获取。首先，通过中频信号建立仿真信号模型，并设计1发N收的L型阵列结构接收目标回波信号；通过对回波信号进行相应的信号处理，最终获得目标的4-D点云。

1.L型阵列

有关于信号建模的方法，各位读者可以回顾博主往期3D点云的内容。如图1所示为L形阵列的阵列结构，有M+N-1个阵元。此L型阵列由x轴上阵元数为N的均匀线阵和y轴上阵元数为M的均匀线阵构成，阵列间距为d。假设空间有K个信源照射到此阵上。其二维波达方向为（θk,ϕk），k = 1,2,3,…,K，其中,θk和ϕk分别表示第k个信源的俯仰角和方位角。

假设入射到此阵列上的信源数为K，则x轴上N个阵元对应的导向矩阵为

y轴上M个阵元对应的导向矩阵为

根据上述的导向矩阵，我们在建立信号模型时，首先需要预先确定Nx和Ny的天线数量，其次就是根据不同的接收天线之间依次加入对应的导向信息，因此信号建模代码如下：

for frameid = 1 : numframe
        for rxid = 1 : numrx - 1
            for chirpid = 1 : numchirp
                for targetid = 1 : numtarget 
                    for  sampleid = 1 : numsamples
                        R = target(targetid,1) + ((chirpid-1)*T)*target(targetid,2) + (frameid-1)*numchirp*T*target(targetid,2); %计算目标随速度变化的距离,以忽略快时间对距离的影响
                        delay = 2 * R / parameter.c;
                        fixphase = exp(1j*2*pi*(f0*delay-0.5*slope*delay^2));         %中频信号的固定相位，末尾的相位为剩余视频相位，为了更真实还原毫米波雷达数据，本程中没有将其省略
                        fastsampdata(sampleid) = abs(randn)^2 * exp(1j*2*pi*slope*delay*t(sampleid)) * fixphase; %这里采用随机起伏   
                    end
                    %用来判断不同的天线间
                    if(rxid <= Nx)
                        azimuthphase = exp(1j*(2*pi/lambda)*cosd(target(targetid,3))*sind(target(targetid,4)) *(rxid -1) * Detarx);  %目标方位角信息
                    else
                        azimuthphase = exp(1j*(2*pi/lambda)*sind(target(targetid,3))*sind(target(targetid,4)) *(rxid - Nx) * Detarx);  %目标方位角信息
                    end
                    rawdata(:,chirpid,rxid,frameid) = rawdata(:,chirpid,rxid,frameid) + fastsampdata.* azimuthphase;
                end
            end
            %加入高斯噪声
            MixIQ = rawdata(:,:,rxid,frameid);
            xigema = std(MixIQ)/db2mag(snr);
            MixIQ = MixIQ + ((randn(size(MixIQ)) + 1i*randn(size(MixIQ)))).*xigema;
            rawdata(:,:,rxid,frameid) = MixIQ;
        end
end

信号时域图如下：