雷达盲区盲速计算matlab,线性调频雷达系统的建模，实现对目标的测距和测速

- 日理万妓 2023-01-15 06:59 183阅读 0赞

在网上找了一个程序，实现了发射波形，回波的产生，以及通过了匹配滤波器。。现在想实现测距和速度     希望各位帮忙看看接下去怎么做，本人是新手。。。。希望说的详细点

C=3.0e8;  %光速(m/s)

RF=3.140e9/2;  %雷达射频 1.57GHz

Lambda=C/RF;%雷达工作波长

PulseNumber=16;   %回波脉冲数

BandWidth=2.0e6;  %发射信号带宽 带宽B=1/τ，τ是脉冲宽度

TimeWidth=42.0e-6; %发射信号时宽

PRT=240e-6;   % 雷达发射脉冲重复周期(s),240us对应1/2\*240\*300=36000米最大无模糊距离

PRF=1/PRT;

Fs=2.0e6;  %采样频率

NoisePower=-12;     %(dB);%噪声功率(目标为0dB)

% ---------------------------------------------------------------%

SampleNumber=fix(Fs\*PRT);%计算一个脉冲周期的采样点数480；

TotalNumber=SampleNumber\*PulseNumber;%总的采样点数480\*16=7680；

BlindNumber=fix(Fs\*TimeWidth);%计算一个脉冲周期的盲区-遮挡样点数；

%===================================================================================%

%目标参数%

%===================================================================================%

TargetNumber=4;%目标个数

SigPower(1:TargetNumber)=\[1 1 1 0.25\];%目标功率,无量纲

TargetDistance(1:TargetNumber)=\[3000 8025 15800 8025\];%目标距离,单位m   距离参数为\[3000 8025 9000+(Y\*10+Z)\*200 8025\]

DelayNumber(1:TargetNumber)=fix(Fs\*2\*TargetDistance(1:TargetNumber)/C);% 把目标距离换算成采样点(距离门) fix函数向0靠拢取整 fix(1.3)=1(取整)

TargetVelocity(1:TargetNumber)=\[50 0 204 100\];%目标径向速度 单位m/s   速度参数为\[50 0 (Y\*10+X+Z)\*6 100\]

TargetFd(1:TargetNumber)=2\*TargetVelocity(1:TargetNumber)/Lambda; %计算目标多卜勒频移2v/λ

%====================================================================================%

%产生线性调频信号%

%====================================================================================%

number=fix(Fs\*TimeWidth);%回波的采样点数=脉压系数长度=暂态点数目+1

if rem(number,2)~=0  %rem求余

number=number+1;

end   %把number变为偶数

Chirp=zeros(1,SampleNumber\*16);

for m=0:16

for i=-fix(number/2):fix(number/2)-1

Chirp(i+m\*SampleNumber+fix(number/2)+1)=exp(j\*(pi\*(BandWidth/TimeWidth)\*((i)/Fs)^2));%exp(j\*fi)\*，产生复数矩阵Chirp

end

coeff=conj(fliplr(Chirp));%把Chirp矩阵翻转并把复数共轭，产生脉压系数

figure(1);%脉压系数的实部

%subplot(2,2,1);

plot(real(Chirp));

axis(\[0 8000 -1.5 1.5\]);

title('调频信号实部');

%-------------------------产生目标回波串-----------------------------------------------------------------------------------------%

fi=pi/3;

SignalTemp=zeros(1,SampleNumber);% 一个脉冲

SignalTemp(DelayNumber(4)+1:DelayNumber(4)+number)=sqrt(SigPower(4))\*exp(j\*fi)\*Chirp(1:number);%一个脉冲的1个目标(未加多普勒速度)

Signal=zeros(1,TotalNumber);

SignalAll=zeros(1,TotalNumber);

%Signal=zeros(1,SampleNumber);

for i=1:PulseNumber %可以换作 PulseNumber

Signal((i-1)\*SampleNumber+1:i\*SampleNumber)=SignalTemp;

end

FreqMove=exp(j\*2\*pi\*TargetFd(4)\*(0:TotalNumber-1)/Fs);%目标的多普勒速度\*时间=目标的多普勒相移

Signal=Signal.\*FreqMove;

SignalAll = SignalAll+Signal;

figure(2);

subplot(2,1,1);

plot(real(SignalAll),'r-');

title('目标信号的实部');

% grid on;

% zoom on;

subplot(2,1,2);

plot(imag(SignalAll));

title('目标信号的虚部');

% grid on;

% zoom on;

%====================================================================================%

%                                   产生系统噪声信号                                  %

%====================================================================================%

%SystemNoise=normrnd(0,10^(NoisePower/10),1,TotalNumber)+j\*normrnd(0,10^(NoisePower/10),1,TotalNumber);

%-------------------------通过高斯信道-----------------------------------------------------------------------------------------%

SNR=10;%信噪比    可作为传送参数

y=awgn(SignalAll,SNR,'measured');

figure(3);

subplot(2,1,1);

plot(real(y));

axis(\[0 8000 -1.5 1.5\]);

title('加入高斯白噪声  实部');

subplot(2,1,2);

plot(imag(y));

axis(\[0 8000 -1.5 1.5\]);

title('加入高斯白噪声  虚部');

%均值为0，标准差为10^(NoisePower/10)的噪声

%====================================================================================%

%                                   总的回波信号                                     %

%====================================================================================%

Echo=y;% +SeaClutter+TerraClutter，加噪声之后的回波

for i=1:PulseNumber   %在接收机闭锁期,接收的回波为0

Echo((i-1)\*SampleNumber+1:(i-1)\*SampleNumber+number)=0; %发射时接收为0

end

figure(4);%加噪声之后的总回波信号

subplot(2,1,1);plot(real(Echo),'r-');title('总回波信号的实部,闭锁期为0');

subplot(2,1,2);plot(imag(Echo));title('总回波信号的虚部,闭锁期为0');

%====================================================================================%

%                                 匹配滤波器                                    %

%====================================================================================%

n=fix(number/2);

Ht=exp(-j\*(pi\*(BandWidth/TimeWidth)\*((n)/Fs)^2));

Sot=conv(Echo,Ht);

figure(5);

% subplot(2,1,1);plot(real(Sot),'r-');title('通过滤波器后的信号   实部');

% subplot(2,1,2);plot(imag(Sot));title('通过滤波器后的信号   虚部');

Z=abs(Sot);

Z=Z/max(Z);%normalize正常化

Z=20\*log10(Z+1e-6);

plot(Z);

雷达盲区盲速计算matlab,线性调频雷达系统的建模，实现对目标的测距和测速

发表评论取消回复

还没有评论，来说两句吧...

相关阅读

相关 ROS2机器人激光雷达测距的C++代码

相关【基于Matlab的RDA和CSA雷达正侧视点目标仿真】——雷达图像生成与分析

相关线性调频脉冲压缩雷达仿真及Matlab源代码

相关雷达系统脉冲压缩的MATLAB实现

相关基于MATLAB的雷达信号处理仿真

相关雷达盲区盲速计算matlab,线性调频雷达系统的建模，实现对目标的测距和测速

相关激光雷达测距原理和点云畸变

相关基于echarts对雷达图的封装

相关基于ROS的Mako相机，Basler相机，禾赛激光雷达，Velodyne激光雷达，速腾激光雷达，Delphi毫米波雷达安装驱动，运行与调试教程

相关线性调频连续波测距、测速原理