现代数字信号处理II课程笔记 | Lecute 1: 阵列信号处理 (1)

Modern Digital Signal Processing (II)

笔记原文链接：https://mp.weixin.qq.com/s/y6ONgpQlPuw2Seoit_JaOg

张颢老师现代数字信号处理II课程笔记
Lecute 1: 阵列信号处理 (1)

本文目录

• Modern Digital Signal Processing (II)
  • 课程概述
  • Part 1: Array Signal Processing (Non-Linear)
    • 1.1 相位补偿的基本原理
    • 1.2 阵列处理问题
    • 1.3 多信号下的角度估计：Angle FFT
    • 1.4 超分辨角度估计基础：子空间、信号子空间和噪声子空间
      • 数据模型：PCA
      • 信号模型：MIMO
      • 数据模型和信号模型的联系
      • 子空间的意义

本文主要内容

• 相控阵
• 波束成型
• 角度FFT
• 子空间的概念：信号子空间、噪声子空间

课程概述

1和2的差异

• 1：1980年以前的事情
• 2：1980年以后的事情

• 上学期的主要内容：4个概念

  • 概念1：线性 Linear
  • 概念2：正交 Orthogonal
  • 概念3：平稳 Stationary
  • 概念4：高斯 Gaussian
  • 研究的主要模型是线性模型
  • 研究的方法总是正交化
  • 研究的信号总是假设是平稳的
  • 研究的信号的统计特性主要是考虑高斯的

• 本学期的主要内容

  • 非线性：阵列信号处理 (Array Processing)

    ✅ 兴起于1980-1990s

    🚩 如：MUSIC算法

  • 非正交：稀疏性 (Sparisty)

    ✅ 兴起于2006-2012s (2012年之后：Deep learning)

    🚩 如：LASSO算法

  • 非平稳: 时频分析 (Time-Frequency Analysis)

    ✅ 兴起于1995-2000s

    🚩 如：Wavelet算法

  • 非高斯: 贝叶斯方法 (Bayesian Method)

    ✅ 兴起于2000-2005s

    🚩 如：MCMC算法

  • “非”：意味着要脱离经典的限制，要处理更加广泛的信号问题

Part 1: Array Signal Processing (Non-Linear)

1980s, 思考的问题：

• Wireless (Communication, Radar, Navigation, Sonar(机械波))
• 目的：Transmission / Recieve Signal
• $\Rightarrow$ 需要 Antenna

• 天线形式：

  • 一般都是孔径天线：形状是固定的

  • 固定天线的优点：High Gain (因为事先可以设计天线，使得能量很集中)

  • 固定天线的缺点：没法移动，需要进行反复调整 (需要伺服系统， Servo )

    ✅ 天线座需要转的快、停的稳、指的准，太难了

• 改进思路：

  • 如何改进天线座的问题？
  • 之前是基于机械的系统
  • 改成Electronic Scanning $\Rightarrow$ 如何实现电扫？
  • 相控阵！(Phased Array)

• 相控阵的特点：

  • 1 Array ：意味着Multiple Antenna Elements，有很多个天线的阵元

  • 2 何谓 相控 ？

    ✅ 假如想要接收来自某特定方向($\theta =30°$)，在不机械旋转的时候如何做到？

    ✅ 信号到达各个阵元的时间是有差别的: $s(t),s(t-{\Delta }_1),...,,s(t-{\Delta }_n)$

    ✅ 如何有效地利用起这些信号？ $\Rightarrow$ 这些信号之间存在时延(Time Difference) $\Rightarrow$ 如何进行补偿(Compensation)？ $\Rightarrow$ 补偿之后，再加到一起即可

    ✅ 所以问题是：如何对这些信号进行补偿？

1.1 相位补偿的基本原理

• 信号：
  • $s(t)=A(t)exp(2\pi f_{0}t)$
  • 其中$A(t)$: Baseband 基带
  • $exp(2\pi f_{0}t)$: 载频
• 信号时延：
  • 同时体现在基带和载频两个部分
  • $s(t-\Delta )=A(t-\Delta )exp(2\pi f_0(t-\Delta ))$
• $\Delta$的来源： 波传递的路程差导致的时延
  • $\Delta =\frac{dsin(\theta )}{c}$
  • $\Rightarrow$ $exp(2\pi f_0(t-\Delta ))=exp(2\pi f_{0}t)exp(-2\pi f_0\Delta )$ = $exp(2\pi f_{0}t)exp(-2\pi f_0\frac{dsin\theta }{c})$ = $exp(2\pi f_{0}t)exp(-\frac{d}{{\lambda }_0}sin\theta )$
  • 即：$s(t-\Delta )=A(t-\Delta )exp(2\pi f_{0}t)exp(-2\pi$ $\frac{d}{{\lambda }_0}sin\theta$ )
  • 其中： $\frac{d}{{\lambda }_0}sin\theta$ 为 空间相位 spatial phase
• 显然， 空间相位中包含了重要的信息：来波方向$\theta$ !
  • spatial phase的来源是因为有多根天线

---

窄带假设：Narrowband Assumption

• $s(t-\Delta )=A(t-\Delta )exp(2\pi f_{0}t)exp(-2\pi$ $\frac{d}{{\lambda }_0}sin\theta$ ) 可近似为：
  • $s(t-\Delta )=A(t)exp(2\pi f_{0}t)exp(-2\pi$ $\frac{d}{{\lambda }_0}sin\theta$ )
• 即窄带假设认为 $A(t)\approxeq A(t-\Delta )$
  • 原因：$A(t)$的带宽相比于载频而言很小
  • 载频：动辄上G；
  • 基带带宽：往往就几M，几十M, 上百M很特殊

---

• 从而有：

  • $s(t-\Delta) = s(t)exp(-2\pi \frac{d}{\lambda_0} sin\theta) $
  • 从而 把时间上的时延 转换为了 相移
  • 因此，只需要对每一个阵元进行相位补偿(乘以相移) $exp(-j\varphi )$
  • 即可实现信号的对齐
  • 这样，就可以对特定角度进行成像
  • 即 相控阵
  • 通过相位的处理，即可抵消阵元、阵元之间的差异

• 电扫的含义：

  • 不需要动天线，只要 动$\theta$，就可以实现扫描 的目的！！
  • 还可以 多波束 $\Rightarrow$ 因为电磁能量可以线性叠加
  • 需要移相器，会有一些工程问题
  • 但是完全没有机械动作！而且孔径天线不可能同时处理多个方向

1.2 阵列处理问题

• 有两个问题

• 问题1：波束形成(beamforming, Scanning) ，对不同的阵元馈以不同的补偿相位。但是应该馈多少？馈$\theta$的相位的话只能顾及一个方向，如何顾及所有方向？另外还希望抑制旁瓣(增加主瓣)

  • 需要求解一组系数 ${\omega }_1,{\omega }_2,...,{\omega }_n$ $\Rightarrow$ 即对这组数进行设计

  • 将这组系数加到各个阵元上：$X(t)={\Sigma }_{k=1}^n{\omega }_k{S}_k(t)$

  • 希望最终的$X(t)$能够满足各种各样的要求

    ✅ 例如：希望能够最小化X的能量: $mi{n}_{\omega }E|X{|}^2$， s.t., 在给定方向上能够正常接收: $w^{T}g({\theta }_0)=1$ $\Rightarrow$ 基于相位补偿，但比相位补偿更高级

• 问题2：DOA(Direction of Arrival)，来波方向估计

  • 为什么不把每个方向都扫一遍？哪个方向大就知道是哪个方向的 $\Rightarrow$ 已经用过了，最早的DOA(机械二战，电扫1970s)

  • 但 对性能有要求 ：

    🚩 分辨率: 2°, 1°, 甚至 <1° $\Rightarrow$ 靠转的方法，依赖于幅度，但幅度受噪声影响很大

    🚩 精度 $\Rightarrow$ 来波方向的信息全部来自相位 $\Rightarrow$ 所以所有基于幅度的方法都不会SOTA， 要使用基于相位的方法

  • 另一个问题：多信号问题(Multiple Signals) (Same Frequency)下的DOA

    ✅ 会产生干涉

    ✅ 处理多信号时，关键在于窄带假设所忽略的基带：利用随机性，在随机意义上的不相关 $\Rightarrow$ 它们之间的不相关天然形成了对信号的分离

    ✅ 分离后，再利用相位信息完成DOA

1.3 多信号下的角度估计：Angle FFT

• 有一组均匀分布的水平天线阵元

  • 阵元间间距为d
  • 均匀线阵(Uniform Linear Array, ULA)
  • 第k个阵元的接收信号：$s_k(t)=x_1^{(k)}(t)+x_2^{(k)}(t)+...+x_m^{(k)}(t)={\Sigma }_{i=1}^m{x}_m^{(k)}(t)$
  • (m: signal sources; n: elements)
  • x i ( k ) ( t ) = A i ( t − Δ i ( k ) ) e x p ( j 2 π f 0 t ) e x p ( − j 2 π f 0 Δ i ( k