雷达信号通信过程仿真MATLAB 实现

这里按"雷达作为通信链路"的思路来做——也就是把雷达发射→目标散射→接收→解调看成一条"通信信道"，在数字域仿真波形设计 → 发射 → 信道（含RCS/路径损耗）→ 接收 → 解调/估计的全过程。：一套传统脉冲雷达链路（测距+测速），一套雷达通信一体化 ISAC（chirp 扩频传比特）

一、系统模型

1.1 雷达"通信过程"三层含义

下面① 给完整主干 + ② 给核心段，都能直接跑。

1.2 链路框图（① 传统雷达）

信息比特(可选) → 波形产生(FMCW chirp) → 发射天线 → 自由空间+RCS → 接收天线+噪声 → 混频+ADC → 2D-FFT → 距离/速度估计

二、MATLAB：FMCW 雷达链路仿真

2.1 主脚本fmcw_radar_link_sim.m

%% FMCW 雷达链路仿真（把雷达当"通信过程"看）clear;clc;close all;%% ===== 1. 雷达参数 =====c=3e8;fc=77e9;% 77GHz 车载雷达B=500e6;% 扫频带宽 500MHzTchirp=50e-6;% chirp 周期 50uslambda=c/fc;% 目标参数R_true=80;% 真实距离 80mv_true=15;% 径向速度 15m/s（朝向雷达为正）RCS=10;% 雷达散射截面积 dBsm，转线性 10^(RCS/10)=10% ADC & 帧fs=4e6;% 采样率 4MHzNchirp=256;% 每帧 chirp 数（速度维）Nfft_r=2048;% 距离 FFT 点数Nfft_v=256;% 速度 FFT 点数fprintf('FMCW 雷达参数: fc=%.1f GHz, B=%.1f MHz, Tchirp=%.1f us\n',fc/1e9,B/1e6,Tchirp*1e6);fprintf('目标: R=%.1f m, v=%.1f m/s, RCS=%.1f dBsm\n',R_true,v_true,10*log10(RCS));%% ===== 2. 发射波形（锯齿调频）=====t=(0:1/fs:Tchirp-1/fs).';k=B/Tchirp;% 调频斜率% 单 chirp 基带（复包络）tx_baseband=exp(1j*pi*k*t.^2);% 发射复数 chirptx_real=real(exp(1j*2*pi*fc*t).*tx_baseband);% 上变频后实信号（仿真可省略射频）% 延拓到 Nchirp 帧tx_frame=repmat(tx_baseband,Nchirp,1);% Nchirp × Nsamp%% ===== 3. 目标回波（"通信信道"）=====% 单程延时 tau = 2R/c，考虑速度引起延时变化（逐 chirp 微调）Nsamp=length(t);tau0=2*R_true/c;rx_frame=zeros(Nchirp,Nsamp);form=1:Nchirp% 当前 chirp 目标距离（匀速）Rm=R_true-v_true*(m-1)*Tchirp;tau_m=2*Rm/c;% 延时 + 多普勒相位% 数字域用分数延时近似：把 tx 平移 + 乘多普勒相delay_samples=tau_m*fs;delay_int=floor(delay_samples);delay_frac=delay_samples-delay_int;% 线性插值做分数延时（简化）rx_chirp=delay_sig(tx_baseband,delay_int,delay_frac);% 多普勒相位：fd = 2v/λ, 每 chirp 附加相位 2π*fd*m*Tchirpfd=2*v_true/lambda;doppler_phase=2*pi*fd*(m-1)*Tchirp;rx_chirp=rx_chirp*exp(1j*doppler_phase);% ===== 雷达"信道"：路径损耗 + RCS =====% 雷达方程简化：Pr ∝ RCS / (4πR^2)^2 * Pt*G^2*λ^2 / (4π)^3% 仿真里我们归一化 Pt=1, G=1，只留 RCS 和距离衰减path_loss=1/(Rm^4+eps);% 双程 R^4rx_chirp=rx_chirp*sqrt(RCS*path_loss);rx_frame(m,:)=rx_chirp.';end%% ===== 4. 噪声（AWGN，"通信过程"的接收机噪声）=====SNR_dB=15;% 接收 SNR（这里指信号峰值/噪声功率）rx_power=mean(abs(rx_frame(:)).^2);noise_power=rx_power/(10^(SNR_dB/10));noise=sqrt(noise_power/2)*(randn(size(rx_frame))+1j*randn(size(rx_frame)));rx_frame_noisy=rx_frame+noise;fprintf('接收 SNR 设定: %.1f dB\n',SNR_dB);%% ===== 5. 混频 + 距离 FFT（传统雷达接收机）=====% 去斜（dechirp）：rx × conj(tx)dechirped=rx_frame_noisy.*conj(repmat(tx_baseband.',Nchirp,1));% 距离 FFT（每 chirp 独立）R_FFT=fft(dechirped,Nfft_r,2);R_FFT=fftshift(R_FFT,2);% 速度 FFT（每距离 bin 沿 chirp 维）RD_map=fft(dechirped,Nfft_v,1);RD_map=fftshift(RD_map,1);% 距离-速度 2D-FFT（经典做法：先距离后速度，或反过来）% 这里做 2D-FFT 直接rd_2d=fft2(dechirped,Nfft_v,Nfft_r);rd_2d=fftshift(rd_2d,1);rd_2d=fftshift(rd_2d,2);%% ===== 6. 距离/速度估计 =====% 频率→距离：fr = 2kR/c → R = fr * c / (2k)freq_axis_r=(-Nfft_r/2:Nfft_r/2-1)*fs/Nfft_r;range_axis=freq_axis_r*c/(2*k);% 速度轴：fv = 2v/λ → v = fv * λ / 2freq_axis_v=(-Nfft_v/2:Nfft_v/2-1)/(Nchirp*Tchirp);vel_axis=freq_axis_v*lambda/2;% 峰值搜索[~,idx_v]=max(abs(rd_2d(:)));[est_v_idx,est_r_idx]=ind2sub(size(rd_2d),idx_v);R_est=range_axis(est_r_idx);v_est=vel_axis(est_v_idx);fprintf('\n===== 估计结果 =====\n');fprintf('距离: 真值 %.1f m | 估计 %.1f m | 误差 %.3f m\n',R_true,R_est,abs(R_est-R_true));fprintf('速度: 真值 %.1f m/s | 估计 %.1f m/s | 误差 %.3f m/s\n',v_true,v_est,abs(v_est-v_true));%% ===== 7. 可视化 =====figure('Color','w','Position',[1001001200500]);% 单 chirp 时域subplot(2,4,1);plot(t*1e6,real(tx_baseband));xlabel('t (us)');ylabel('幅度');title('发射 chirp（实部）');grid on;subplot(2,4,2);plot(t*1e6,abs(rx_frame(1,:)));hold on;plot(t*1e6,abs(rx_frame_noisy(1,:)),'r');xlabel('t (us)');ylabel('幅度');title('接收回波（蓝=无噪 红=有噪）');grid on;% 距离谱（单 chirp）subplot(2,4,3);plot(range_axis,abs(R_FFT(1,:)));xlabel('距离 (m)');ylabel('|FFT|');title('距离 FFT（单 chirp）');grid on;xlim([0150]);% 距离-多普勒图subplot(2,4,[4,7,8]);% 占3格imagesc(range_axis,vel_axis,20*log10(abs(rd_2d)/max(abs(rd_2d(:)))+eps));xlabel('距离 (m)');ylabel('速度 (m/s)');title('距离-速度 2D-FFT 图');colorbar;hold on;plot(R_true,v_true,'rp','MarkerSize',14,'LineWidth',2);legend('真值');% SNR 变化对估计影响（小扫）subplot(2,4,5);SNR_scan=0:5:30;R_err=zeros(size(SNR_scan));fors=1:length(SNR_scan)np=rx_power/(10^(SNR_scan(s)/10));noi=sqrt(np/2)*(randn(size(rx_frame))+1j*randn(size(rx_frame)));rd_tmp=fft2(rx_frame+noi,Nfft_v,Nfft_r);rd_tmp=fftshift(rd_tmp,1);rd_tmp=fftshift(rd_tmp,2);[~,idx]=max(abs(rd_tmp(:)));[ev,er]=ind2sub(size(rd_tmp),idx);R_err(s)=abs(range_axis(er)-R_true);endplot(SNR_scan,R_err,'o-','LineWidth',1.5);grid on;xlabel('SNR (dB)');ylabel('距离误差 (m)');title('SNR 对测距误差影响');sgtitle('FMCW 雷达"通信过程"仿真','FontSize',14,'FontWeight','bold');

2.2 延时子函数（手写，不搜）

functiony=delay_sig(x,d_int,d_frac)% x: 列向量, d_int 整数延时, d_frac 分数(0~1)N=length(x);y=zeros(size(x));% 整数部分：平移ifd_int>=0y(d_int+1:end)=x(1:end-d_int);elsey(1:end+d_int)=x(-d_int+1:end);end% 分数部分：邻近线性插值（简化）ifd_frac>0y=y+d_frac*(circshift(y,-1)-y);% 粗糙但够仿真用endend

参考代码 对雷达信号的通信过程进行了仿真www.youwenfan.com/contentcsw/82167.html

三、ISAC 小段：chirp 扩频"雷达通信一体化"

如果你说的"通信过程"是指雷达波形上扛比特，给一段核心，嵌到上面框架里就能跑：

%% ===== 雷达通信一体化：chirp 扩频（每个 chirp 扛 1 bit）=====Nbits=Nchirp;bits=randi([01],Nbits,1);% 调制：0→+Δf, 1→-Δf（频偏键控到 chirp 斜率上）delta_f=200e3;% 200kHz 频偏tx_isac=zeros(Nchirp,Nsamp);form=1:Nchirp k_m=k+(2*bits(m)-1)*delta_f/Tchirp;% 斜率微调 = 频移tx_isac(m,:)=exp(1j*pi*k_m*t.^2).';end% 接收端：先雷达去斜 → 2D-FFT 测距测速% 再：对同一距离 bin 沿 chirp 维做"频偏估计"解调比特% 简化：对 dechirped 沿 chirp 维做 FFT，看峰值偏哪dech_isac=tx_isac.*conj(repmat(tx_baseband.',Nchirp,1));% 粗同步假设% 实际要对着目标距离 bin 抽出来，这里示意chirp_fft=fft(dech_isac(:,find(range_axis>75&range_axis<85,1)),Nfft_v,1);% 峰值位置 → 判 bit（需训练 + 同步，这里只示意结构）fprintf('ISAC: %d 个 chirp 各扛 1 bit，BER 需同步后算\n',Nbits);

这一段就是典型论文里的FMCW-Chirp Spread Spectrum思路，雷达测距+低速通信同时做。