宽带信号的DOA估计学习笔记(三):宽带信号阵列接收模型

在宽带阵列信号处理中,阵列接收到的信号包络不再恒定,相位延迟和时延不再是简单的线性关系。

宽带时域模型

MMM个阵元组成的阵列,接收空间中PPP个宽带信号,入射角度分别为θ1,θ2, ,θP\theta_1,\theta_2,\cdots,\theta_Pθ1​,θ2​,⋯,θP​,则第kkk个阵元接收到的信号可表示为xk(t)=i=1Psi[t+τk(θi)]+nk(t) x_k(t)=\sum_{i=1}^P s_i[t+\tau_k(\theta_i)]+n_k(t)xk​(t)=i=1∑P​si​[t+τk​(θi​)]+nk​(t)

宽带频域模型

傅里叶变换的时延定理:信号时延后的傅里叶变换等于信号傅里叶变换的相位延迟
如果s(t)s(t)s(t) 的傅里叶变换为s(f)s(f)s(f),即
FT[s(t)]=s(f) FT[s(t)]=s(f)FT[s(t)]=s(f)
那么 FT[s(t+τ)]=s(f)ej2πτFT[s(t+\tau)]=s(f)e^{j2\pi\tau}FT[s(t+τ)]=s(f)ej2πτ

把第k个阵元接收到的信号,两边做傅里叶变换得到
xk(f)=i=1Psi(f)ej2πfτk(θi)+nk(f)x_k(f)=\sum_{i=1}^P {s_i(f)}e^{j2\pi f\tau_k(\theta_i)}+n_k(f)xk​(f)=i=1∑P​si​(f)ej2πfτk​(θi​)+nk​(f)
写成矩阵形式
Xf=A(f,θ)S(f)+N(f)\mathbf{X}_f=\mathbf{A}(f,\theta)\mathbf{S}(f)+\mathbf{N}(f)Xf​=A(f,θ)S(f)+N(f)

其中方向矩阵为
A(f,θ)=[ej2πfτ1(θ1)ej2πfτ1(θ2)ej2πfτ1(θP)ej2πfτ2(θ1)ej2πfτ2(θ2)ej2πfτ2(θP)ej2πfτM(θ1)ej2πfτM(θ2)ej2πfτM(θP)]\mathbf{A}(f,\mathbf{\theta})= \left[ \begin{matrix} e^{j2\pi f\tau_1(\theta_1)} & e^{j2\pi f\tau_1(\theta_2)} &\cdots &e^{j2\pi f\tau_1(\theta_P)} \\e^{j2\pi f\tau_2(\theta_1)} & e^{j2\pi f\tau_2(\theta_2)} &\cdots &e^{j2\pi f\tau_2(\theta_P)} \\ \vdots &\vdots &\ddots &\vdots \\e^{j2\pi f\tau_M(\theta_1)} & e^{j2\pi f\tau_M(\theta_2)} &\cdots &e^{j2\pi f\tau_M(\theta_P)} \end{matrix} \right]A(f,θ)=⎣⎢⎢⎢⎡​ej2πfτ1​(θ1​)ej2πfτ2​(θ1​)⋮ej2πfτM​(θ1​)​ej2πfτ1​(θ2​)ej2πfτ2​(θ2​)⋮ej2πfτM​(θ2​)​⋯⋯⋱⋯​ej2πfτ1​(θP​)ej2πfτ2​(θP​)⋮ej2πfτM​(θP​)​⎦⎥⎥⎥⎤​

与窄带信号的方向矩阵不同的是,窄带模型中的频率为固定的单值,这里的频率为信号的整个频带。

当对数字信号进行JJJ点离散傅里叶变换时,频率点也为JJJ个离散点,那么
X(fj)=A(fj,θ)S(fj)+N(fj)(j=0,1, ,J)\mathbf{X}{(f_j)}=\mathbf{A}{(f_j,\mathbf{\theta})} \mathbf{S}{(f_j)}+\mathbf{N}{(f_j)}(j=0,1,\cdots,J)X(fj​)=A(fj​,θ)S(fj​)+N(fj​)(j=0,1,⋯,J)
其中,方向矩阵为
A(fj,θ)=[ej2πfjτ1(θ1)ej2πfjτ1(θ2)ej2πfjτ1(θP)ej2πfjτ2(θ1)ej2πfjτ2(θ2)ej2πfjτ2(θP)ej2πfjτM(θ1)ej2πfjτM(θ2)ej2πfjτM(θP)]=[a(θ1,fj)a(θ2,fj)a(θP,fj)](j=0,1,2, ,J)\begin{aligned} \mathbf{A}(f_j,\mathbf{\theta})&= \left[ \begin{matrix} e^{j2\pi f_j\tau_1(\theta_1)} & e^{j2\pi f_j\tau_1(\theta_2)} &\cdots &e^{j2\pi f_j\tau_1(\theta_P)} \\e^{j2\pi f_j\tau_2(\theta_1)} & e^{j2\pi f_j\tau_2(\theta_2)} &\cdots &e^{j2\pi f_j\tau_2(\theta_P)} \\ \vdots &\vdots &\ddots &\vdots \\e^{j2\pi f_j\tau_M(\theta_1)} & e^{j2\pi f_j\tau_M(\theta_2)} &\cdots &e^{j2\pi f_j\tau_M(\theta_P)} \end{matrix} \right]\\ &=[\mathbf{a}(\theta_1,f_j)\qquad \mathbf{a}(\theta_2,f_j)\qquad \cdots \qquad \mathbf{a}(\theta_P,f_j)](j=0,1,2,\cdots,J) \end{aligned} A(fj​,θ)​=⎣⎢⎢⎢⎡​ej2πfj​τ1​(θ1​)ej2πfj​τ2​(θ1​)⋮ej2πfj​τM​(θ1​)​ej2πfj​τ1​(θ2​)ej2πfj​τ2​(θ2​)⋮ej2πfj​τM​(θ2​)​⋯⋯⋱⋯​ej2πfj​τ1​(θP​)ej2πfj​τ2​(θP​)⋮ej2πfj​τM​(θP​)​⎦⎥⎥⎥⎤​=[a(θ1​,fj​)a(θ2​,fj​)⋯a(θP​,fj​)](j=0,1,2,⋯,J)​

如果再将观察时间T0T_0T0​内的阵列接收数据分为LLL个子段,每段时间为TdT_dTd​,然后对观察数据进行JJJ点离散傅里叶变换,只要子段TdT_dTd​ 相比噪声的相关时间较长(保证DFT后的数据不相关),宽带模型就变为
Xl(fj)=A(fj,θ)Sl(fj)+Nl(fj)\mathbf{X}_l(f_j)=\mathbf{A}(f_j,\theta)\mathbf{S}_l(f_j)+\mathbf{N}_l(f_j)Xl​(fj​)=A(fj​,θ)Sl​(fj​)+Nl​(fj​)

其中,JJJ是指将带宽BBB的宽带信号划分为JJJ个子带的个数,对于子带对应的不同频率点f1,f2, ,fJf_1,f_2,\cdots,f_Jf1​,f2​,⋯,fJ​,均成立。

上一篇:@ApiParam @PathVariable @RequestParam三者区别


下一篇:如何把视频顺时转90度同时添加滚动字幕?