线性调频波形测距测速信号处理——全代码+注释

clear all

close all

clc

%% 参数设置

fs=600e6;%采样率

fc=10.45e9;% 波形发射载频

t=10e-6;%脉宽

f0=50e6;%波形中频频率

B=10e6;%带宽

u=B/(2*t);%调频斜率

Tv=100e-6;% 脉冲重复周期

Num=64;% 测速脉冲数

lamd=fs/B;% 抽取带宽

Ns=round(fs*t);

NTv=round(fs*Tv);

tt=0:1/fs:t-1/fs;

ff=0:1/t:fs-1/t;

%% 线性调频波形调制

yf=exp(1i*2*pi*((f0-B).*tt+u*tt.^2));

y_lfm=real(yf);

y_lfm_fft=fft(y_lfm);

figure

plot(tt.*1e6,y_lfm./max(y_lfm),'k')

xlabel('时间(us)')

ylabel('归一化幅度')

title('线性调频波形时域')

set(gca,'FontName','宋体','FontSize',12);

figure

plot((ff-fs/2)./1e6,abs(fftshift(y_lfm_fft))./max(abs(y_lfm_fft)),'k')

xlabel('时间(us)')

ylabel('归一化幅度')

title('线性调频波形频域')

set(gca,'FontName','宋体','FontSize',12);

%% 匹配副本处理

bellm=exp(1i*2*pi*f0.*tt);%匹配副本的下变频本振

y_match=yf.*bellm;%发射波形下变频

window_h=hamming(length(y_match));%汉明窗函数

y_match_w=y_match.*window_h';%匹配副本加窗

%% 探测目标设定

c=3e8; % 光传播速度

R=50;% 目标距离

dt=R*1e3*2/c; % 回波时延

Ndt=round(fs*dt);

v=20; % 目标径向速度

df=2*v/(c/fc);% 速度对应的多普勒频移

tt_v=0:1/fs:(Num*NTv-1)*1/fs;

fd=exp(2*j*pi*df.*tt_v);% 多普勒频率调制函数

%% 接收处理仿真

for k=1:Num;

yr(k,:)=[zeros(1,Ndt),y_lfm.*fd((k-1)*NTv+1:(k-1)*NTv+Ns)];% 回波延时

yr_n(k,:)=awgn(yr(k,:),20,'measured');% 加性噪声

tt_r=0:1/fs:(length(yr_n(k,:))-1)*1/fs;

bell_i=cos(2*pi*f0.*tt_r);% 回波下变频本振I路信号

bell_q=sin(2*pi*f0.*tt_r);% 回波下变频本振Q路信号

yr_n_i(k,:)=yr_n(k,:).*bell_i;% 正交调制

yr_n_q(k,:)=yr_n(k,:).*bell_q;% 正交调制

% yr_n_l_i(k,:)=filter(b_d,a_d,yr_n_i(k,:));% 低通滤波

% yr_n_l_q(k,:)=filter(b_d,a_d,yr_n_q(k,:));% 低通滤波

yr_n_l(k,:)=yr_n_i(k,:)+1i*yr_n_q(k,:);

yr_n_l_fft(k,:)=fft(yr_n_l(k,:));

y_match_w_fft=fft(y_match_w,length(yr_n_l_fft(k,:)));

%脉冲压缩

res_fft(k,:)=yr_n_l_fft(k,:).*conj(y_match_w_fft);

res(k,:)=ifft(res_fft(k,:));

res_c(k,:)=res(k,1:lamd:end);%等带宽抽取

end

%% 测距

res_c_1=res_c(1,:);

tt_s=0:1/fs*lamd:(length(res_c_1)-1)*1/fs*lamd;

dd_s=tt_s*c/2;

figure

plot(dd_s./1e3,-20*abs(log10(abs(res_c_1)./max(abs(res_c_1)))),'k');

ylim([-100 0])

xlabel('距离(km)')

ylabel('幅度(dB)')

title('脉冲压缩测距')

set(gca,'FontName','宋体','FontSize',12);

%% 测速

for k=1:length(res_c(1,:));

res_c_fft(:,k)=fft(res_c(:,k));

end

\~,ind_s_temp\]=max(abs(res_c(1,:))); res_v=abs(res_c_fft(:,ind_s_temp)); vv_s=0:1/(Num\*Tv)\*(c/fc)/2:(Num-1)\*1/(Num\*Tv)\*(c/fc)/2; figure plot(vv_s,res_v./max(res_v),'k'); xlabel('速度(m/s)') ylabel('归一化幅度') title('相参测速') set(gca,'FontName','宋体','FontSize',12); %% 相参处理后距离-速度三维图 res_c_fft=fft(res_c); %相参积累 figure surf(dd_s,vv_s,abs(res_c_fft)./max(max(abs(res_c_fft)))) xlabel('距离(m)') ylabel('速度(m/s)') zlabel('归一化幅度') set(gca,'FontName','宋体','FontSize',12); shading interp

相关推荐
刘小小_算法工程师10 天前
「ECG信号处理——(18)基于时空特征的心率变异性分析」2025年6月23日
信号处理
zhangfeng113311 天前
MocapApi 中文文档 和github下载地址 NeuronDataReader(以下简称 NDR)的下一代编程接口
github·信号处理
风靡晚11 天前
汽车毫米波雷达增强感知:基于相干扩展和高级 IAA 的超分辨率距离和角度估计.
算法·汽车·信息与通信·信号处理·fmcw
芳菲菲其弥章13 天前
【信号与系统二】连续时间傅里叶变换
信号处理·傅里叶分析
FPGA_ADDA13 天前
宽带中频10.4G采集卡
fpga开发·信号处理·高速数据采集·10g采集卡
棱镜研途14 天前
学习笔记丨数字信号处理(DSP)的应用——图像处理篇
图像处理·人工智能·信号处理
军训猫猫头14 天前
100.Complex[]同时储存实数和虚数两组double的数组 C#例子
算法·c#·信号处理
科研小刘带你玩学术15 天前
【跨界新视野】信号处理遇上VR/AR:下一代沉浸体验的核心技术与您的发表蓝海
人机交互·信号处理·虚拟现实·ei会议·多模态融合·论文发表·学术发表
gihigo199816 天前
RA信号处理
信号处理
FF-Studio16 天前
【DSP笔记 · 第7章】信号处理的“整形”大师:FIR滤波器与线性相位的奥秘
笔记·自动化·音视频·音频·信号处理