跳转到内容

信号模型

当前 OpenISAC 系统包含三条链路:BS→UE 下行通信链路、UE→BS 上行通信链路,以及 BS 单站感知链路。BS 的下行波形同时承担通信信号和单站感知照射信号;UE 对同一下行波形的接收既用于通信解调,也构成双站感知观测。

x{DL,UL}x\in\{\mathrm{DL},\mathrm{UL}\} 表示下行或上行,并令 qxq_x 表示链路 xx 的接收端,即 qDL=UEq_\mathrm{DL}=\mathrm{UE}qUL=BSq_\mathrm{UL}=\mathrm{BS}。第 ll 条路径的时延分为三个层次。

真实无线传播时延记为

τl,propx(t),\tau_{l,\mathrm{prop}}^x(t),

它只描述电磁波在空间和散射环境中的传播。链路 xx 的固定 RF 群时延记为 τxRF\tau_x^\mathrm{RF},包含对应发射链和接收链的合成延迟。因此物理总链路时延为

τl,linkx(t)=τl,propx(t)+τxRF.\tau_{l,\mathrm{link}}^x(t) =\tau_{l,\mathrm{prop}}^x(t)+\tau_x^\mathrm{RF}.

接收机以当前本地解调窗口切分 OFDM 符号。令 τdqx(t)\tau_d^{q_x}(t) 表示接收端 qxq_x 当前解调窗口相对于链路 xx 发射端帧边界的时变偏移,即解调窗口起始时刻减去发射端帧边界时刻;正值表示解调窗口位于发射端帧边界之后。采样频率偏差会使解调窗口缓慢偏移,且初始同步或后续整数采样校正会直接更新当前解调窗口,因此该偏移是时变的。

定义本地延迟轴上的 TO 为真实传播时延到本地观测路径时延的公共偏移。对下行和上行分别有

τTOUE(t)τDLRFτdUE(t),τTOBS(t)τULRFτdBS(t).\tau_\mathrm{TO}^\mathrm{UE}(t) \triangleq\tau_\mathrm{DL}^\mathrm{RF}-\tau_d^\mathrm{UE}(t), \qquad \tau_\mathrm{TO}^\mathrm{BS}(t) \triangleq\tau_\mathrm{UL}^\mathrm{RF}-\tau_d^\mathrm{BS}(t).

因此两端在自己的本地延迟轴上观测到的第 ll 条路径时延分别为

τlUE(t)=τl,prop(t)+τTOUE(t),\tau_l^\mathrm{UE}(t) =\tau_{l,\mathrm{prop}}(t) +\tau_\mathrm{TO}^\mathrm{UE}(t), τlBS(t)=τl,prop(t)+τTOBS(t).\tau_l^\mathrm{BS}(t) =\tau_{l,\mathrm{prop}}(t) +\tau_\mathrm{TO}^\mathrm{BS}(t).

这里的 τlDL\tau_l^\mathrm{DL}τlUL\tau_l^\mathrm{UL} 是包含 RF 群时延的物理总链路时延,τlUE\tau_l^\mathrm{UE}τlBS\tau_l^\mathrm{BS} 则是两个接收机在本地延迟轴上的直接观测。后者与真实传播时延之差即为本文的 TO。

接收端 qxq_x 观察到的等效时变基带冲激响应写为

hx(t,τ)=l=0Lx1αlx(t)ej2π(fD,lx+Δfcx)tδ ⁣(ττl,propx(t)τxRF+τdqx(t)).h_x(t,\tau)= \sum_{l=0}^{L_x-1} \alpha_l^x(t) e^{j2\pi(f_{D,l}^x+\Delta f_c^x)t} \delta\!\left( \tau-\tau_{l,\mathrm{prop}}^x(t) -\tau_x^\mathrm{RF} +\tau_d^{q_x}(t) \right).

其中,LxL_x 是可分辨路径数;αlx(t)\alpha_l^x(t)τl,propx(t)\tau_{l,\mathrm{prop}}^x(t)fD,lxf_{D,l}^x 分别是第 ll 条路径的复散射系数、真实传播时延和多普勒频移;Δfcx\Delta f_c^x 是链路残余载波频偏。下行与上行接收信号分别满足

yUEDL(t)=hDL(t,τ)sDL(tτ)dτ+zUE(t),y_\mathrm{UE}^\mathrm{DL}(t) =\int h_\mathrm{DL}(t,\tau)s_\mathrm{DL}(t-\tau)\,d\tau +z_\mathrm{UE}(t), yBSUL(t)=hUL(t,τ)sUL(tτ)dτ+zBSUL(t).y_\mathrm{BS}^\mathrm{UL}(t) =\int h_\mathrm{UL}(t,\tau)s_\mathrm{UL}(t-\tau)\,d\tau +z_\mathrm{BS}^\mathrm{UL}(t).

tDLt_{\mathrm{DL}}tULt_{\mathrm{UL}} 分别表示参与 eRTM 的下行与上行参考 OFDM 符号时刻。TDD 中若二者间隔远小于信道相干时间,则

τl,prop(tDL)τl,prop(tUL)τl,prop,\tau_{l,\mathrm{prop}}(t_{\mathrm{DL}}) \approx \tau_{l,\mathrm{prop}}(t_{\mathrm{UL}}) \approx \tau_{l,\mathrm{prop}}, αlDL(tDL)αlUL(tUL)αl,\alpha_l^\mathrm{DL}(t_{\mathrm{DL}}) \approx \alpha_l^\mathrm{UL}(t_{\mathrm{UL}}) \approx \alpha_l,

TO 主要随两端时钟缓慢漂移,因此对于时间上足够接近的下行与上行参考符号,同一测量对还满足

τTOUE(tDL)τTOUE,τTOBS(tUL)τTOBS.\tau_\mathrm{TO}^\mathrm{UE}(t_{\mathrm{DL}}) \approx\tau_\mathrm{TO}^\mathrm{UE}, \qquad \tau_\mathrm{TO}^\mathrm{BS}(t_{\mathrm{UL}}) \approx\tau_\mathrm{TO}^\mathrm{BS}.

这里将同一测量对中两个相近时刻的 TO 分别简记为 τTOUE\tau_\mathrm{TO}^\mathrm{UE}τTOBS\tau_\mathrm{TO}^\mathrm{BS},并不假设二者相等。UE 端观测到的下行频域信道和 BS 端观测到的上行频域信道分别为

HUE[n]=l=0L1αlej2πκnΔf[τl,prop+τTOUE],H_{\mathrm{UE}}[n] =\sum_{l=0}^{L-1}\alpha_l e^{-j2\pi\kappa_n\Delta f[ \tau_{l,\mathrm{prop}} +\tau_\mathrm{TO}^\mathrm{UE}]}, HBS[n]=l=0L1αlej2πκnΔf[τl,prop+τTOBS].H_{\mathrm{BS}}[n] =\sum_{l=0}^{L-1}\alpha_l e^{-j2\pi\kappa_n\Delta f[ \tau_{l,\mathrm{prop}} +\tau_\mathrm{TO}^\mathrm{BS}]}.

定义两端 TO 之差

τTOBSUEτTOBSτTOUE=τlBSτlUE,\tau_\mathrm{TO}^{\mathrm{BS-UE}} \triangleq \tau_\mathrm{TO}^\mathrm{BS} -\tau_\mathrm{TO}^\mathrm{UE} =\tau_l^\mathrm{BS}-\tau_l^\mathrm{UE},

则有

HBS[n]HUE[n]ej2πκnΔfτTOBSUE.\boxed{ H_{\mathrm{BS}}[n] \approx H_{\mathrm{UE}}[n] e^{-j2\pi\kappa_n\Delta f \tau_\mathrm{TO}^{\mathrm{BS-UE}}} }.

假设 BS 端配备 RR 阵元的均匀线阵(ULA)。阵元间距为 dad_a,下行载波波长为 λ=c/fc\lambda=c/f_c,入射角 θ\theta 从阵列法向(broadside)计量。阵列导向矢量定义为

a(θ)=[1ejμ(θ)ej(R1)μ(θ)]T,μ(θ)=2πdaλsinθ.\boldsymbol a(\theta)= \begin{bmatrix} 1 & e^{j\mu(\theta)} & \cdots & e^{j(R-1)\mu(\theta)} \end{bmatrix}^{T}, \qquad \mu(\theta)=\frac{2\pi d_a}{\lambda}\sin\theta.

设场景中共有 Q=P+CQ=P+C 个可分辨反射分量,其中 PP 个为运动目标、CC 个为静态或近静态杂波。令 τs,pprop\tau_{s,p}^\mathrm{prop} 为第 pp 个分量的真实往返传播时延,τsensRF\tau_\mathrm{sens}^\mathrm{RF} 为单站收发链的固定群时延,则

τs,plink=τs,pprop+τsensRF,\tau_{s,p}^\mathrm{link} =\tau_{s,p}^\mathrm{prop}+\tau_\mathrm{sens}^\mathrm{RF}, hBSsens(t,τ)=p=1Qβpa(θp)δ(ττs,plink)ej2πfD,s,pt.\boldsymbol h_\mathrm{BS}^\mathrm{sens}(t,\tau) =\sum_{p=1}^{Q} \beta_p\boldsymbol a(\theta_p) \delta(\tau-\tau_{s,p}^\mathrm{link}) e^{j2\pi f_{D,s,p}t}.

因此多通道接收向量为

yBSsens(t)=hBSsens(t,τ)sDL(tτ)dτ+zBSsens(t),\boldsymbol y_\mathrm{BS}^\mathrm{sens}(t) =\int \boldsymbol h_\mathrm{BS}^\mathrm{sens}(t,\tau) s_\mathrm{DL}(t-\tau)\,d\tau +\boldsymbol z_\mathrm{BS}^\mathrm{sens}(t),

其中 zBSsens(t)CN(0,σ2IR)\boldsymbol z_\mathrm{BS}^\mathrm{sens}(t)\sim\mathcal{CN}(\boldsymbol 0,\sigma^2\boldsymbol I_R)。对距离 rpr_p、径向速度 vpv_p 和雷达散射截面 σRCS,p\sigma_{\mathrm{RCS},p} 的点目标,窄带远场模型给出

βp=c2σRCS,p(4π)3rp4fc2ejϕp,τs,pprop=2rpc,fD,s,p=2vpfcc.\beta_p= \sqrt{\frac{c^2\sigma_{\mathrm{RCS},p}} {(4\pi)^3r_p^4f_c^2}}e^{j\phi_p}, \qquad \tau_{s,p}^\mathrm{prop}=\frac{2r_p}{c}, \qquad f_{D,s,p}=\frac{2v_pf_c}{c}.

校准并去除 τsensRF\tau_\mathrm{sens}^\mathrm{RF} 后,单站距离和速度分别由 r=cτprop/2r=c\tau^\mathrm{prop}/2v=cfD/(2fc)v=cf_D/(2f_c) 得到;这里约定 vp>0v_p>0 表示目标沿径向接近。

名义采样间隔为 Ts=1/BT_s=1/B。若接收端 qxq_x 的实际采样间隔为 Ts,qx=TsΔTs,qxT_{s,q_x}=T_s-\Delta T_{s,q_x},则第 ll 条通信路径的离散样本可近似写为

yqx,l[k]αlxsx ⁣(kTsτl,propxτxRF+τdqx[0]kΔTs,qx)ej2π(fD,lx+Δfcx)kTs.y_{q_x,l}[k]\approx \alpha_l^x s_x\!\left( kT_s -\tau_{l,\mathrm{prop}}^x -\tau_x^\mathrm{RF} +\tau_d^{q_x}[0] -k\Delta T_{s,q_x} \right) e^{j2\pi(f_{D,l}^x+\Delta f_c^x)kT_s}.

该式忽略了 (fD,lx+Δfcx)ΔTs,qx(f_{D,l}^x+\Delta f_c^x)\Delta T_{s,q_x} 的二阶小量。固定 RF 群时延形成 TO 的静态部分,当前解调窗口相对于发射端帧边界的偏移形成其时变部分,并以负号计入 TO。CFO 产生随时间累积的公共相位,SFO 则使 TO 随时间缓慢漂移并形成随子载波变化的相位斜率。