原始CSI模型
信道状态信息 (CSI) 是来自所有路径的信号的叠加。在数学上,CSI 可以表示为:
H
(
f
,
t
)
=
∑
i
=
1
L
A
i
e
−
j
2
π
d
i
(
t
)
λ
H(f, t)=\sum_{i=1}^{L} A_{i} e^{-j 2 \pi \frac{d_{i}(t)}{\lambda}}
H(f,t)=i=1∑LAie−j2πλdi(t)
根据[1],路径可以分为静态路径和动态路径。不失一般性,假设只有一条由人体反射的动态路径,而静态分量是由环境中静态物体的视距传播和其他反射路径组成。因此,CSI 可以重写为:
H
(
f
,
t
)
=
H
s
(
f
,
t
)
+
H
d
(
f
,
t
)
=
H
s
(
f
,
t
)
+
A
(
f
,
t
)
e
−
j
2
π
d
(
t
)
λ
H(f, t)=H_{s}(f, t)+H_{d}(f, t)=H_{s}(f, t)+A(f, t) e^{-j 2 \pi \frac{d(t)}{\lambda}}
H(f,t)=Hs(f,t)+Hd(f,t)=Hs(f,t)+A(f,t)e−j2πλd(t)
其中
H
s
(
f
,
t
)
H_{s}(f, t)
Hs(f,t) 是静态分量,
A
(
f
,
t
)
,
e
−
j
2
π
d
(
t
)
λ
A(f, t), e^{-j 2 \pi \frac{d(t)}{\lambda}}
A(f,t),e−j2πλd(t)和
d
(
t
)
d(t)
d(t)分别是信号幅度,信号相位和动态路径的长度变化。
当人体目标移动短距离时,动态分量
A
(
f
,
t
)
A(f, t)
A(f,t)的信号幅度可以认为是一个常数。这是因为信号幅度是由路径长度决定的(信号幅度与动态路径的长度成反比)。理想情况下,当
d
(
t
)
d(t)
d(t)长度增加一个波长, CSI
(
H
(
f
,
t
)
(H(f, t)
(H(f,t)) 顺时针旋转
2
π
2 \pi
2π,如图1所示。
但对于商品 WiFi 设备,由于发送器和接收器不是时间同步的,因此每个 CSI 样本中都有一个随时间变化的随机相位偏移
e
−
j
θ
offset
e^{-j \theta_{\text {offset } }}
e−jθoffset ,如下所示:
H
(
f
,
t
)
=
e
−
j
θ
o
f
f
s
e
t
(
H
s
(
f
,
t
)
+
A
(
f
,
t
)
e
−
j
2
π
d
(
t
)
λ
)
H(f, t)=e^{-j \theta_{o f f s e t}}\left(H_{s}(f, t)+A(f, t) e^{-j 2 \pi \frac{d(t)}{\lambda}}\right)
H(f,t)=e−jθoffset(Hs(f,t)+A(f,t)e−j2πλd(t))
有了这个随机相位偏移,复平面中运动引起的 CSI 变化不再是一个圆,此时CSI如图2所示:
关于scaling、rotation和translation的定义
scaling ( H ( f , t ) ↦ α H ( f , t ) , α ∈ R ) (H(f, t) \mapsto \alpha H(f, t), \alpha \in \mathbb{R}) (H(f,t)↦αH(f,t),α∈R),
rotation ( H ( f , t ) ↦ e i θ H ( f , t ) , θ ∈ R ) \left(H(f, t) \mapsto e^{i \theta} H(f, t), \theta \in \mathbb{R}\right) (H(f,t)↦eiθH(f,t),θ∈R)
translation ( H ( f , t ) ↦ H ( f , t ) + β , β ∈ C ) (H(f, t) \mapsto H(f, t)+\beta, \beta \in \mathbb{C}) (H(f,t)↦H(f,t)+β,β∈C)
当对
H
(
f
,
t
)
H(f, t)
H(f,t) 进行scaling、rotation、 translation时,不会改变
H
(
f
,
t
)
H(f, t)
H(f,t)的几何形状和旋转方向 ,如图3所示:
CSI商模型
首先明确两个点:
①对于商品 WiFi 卡,例如广泛使用的 Intel 5300,Wi-Fi 卡上不同天线的时变相位偏移相同,因为它们共享相同的 RF 振荡器 [2, 3]。
②当目标移动短距离(几厘米)时,两个附近天线处的两个反射路径长度之差 d 2 ( t ) − d 1 ( t ) d_{2}(t)-d_{1}(t) d2(t)−d1(t) 可以被认为是一个常数 Δ d \Delta d Δd [4]。
对第一根天线和第二根天线的CSI做比值得到CSI商:
H
1
(
f
,
t
)
H
2
(
f
,
t
)
=
e
−
j
θ
o
f
f
s
e
t
(
H
s
,
1
+
A
1
e
−
j
2
π
d
1
(
t
)
λ
)
e
−
j
θ
o
f
f
s
e
t
(
H
s
,
2
+
A
2
e
−
j
2
π
d
2
(
t
)
λ
)
=
A
1
e
−
j
2
π
d
1
(
t
)
λ
+
H
s
,
1
A
2
e
−
j
2
π
d
1
(
t
)
+
Δ
d
λ
+
H
s
,
2
=
A
1
e
−
j
2
π
d
1
(
t
)
λ
+
H
s
,
1
A
2
e
−
j
2
π
Δ
d
λ
e
−
j
2
π
d
1
(
t
)
λ
+
H
s
,
2
\begin{aligned} \frac{H_{1}(f, t)}{H_{2}(f, t)} &=\frac{e^{-j \theta_{o f f s e t}}\left(H_{s, 1}+A_{1} e^{-j 2 \pi \frac{d_{1}(t)}{\lambda}}\right)}{e^{-j \theta_{o f f s e t}}\left(H_{s, 2}+A_{2} e^{-j 2 \pi \frac{d_{2}(t)}{\lambda}}\right)} \\ &=\frac{A_{1} e^{-j 2 \pi \frac{d_{1}(t)}{\lambda}}+H_{s, 1}}{A_{2} e^{-j 2 \pi \frac{d_{1}(t)+\Delta d}{\lambda}}+H_{s, 2}} \\ &=\frac{A_{1} e^{-j 2 \pi \frac{d_{1}(t)}{\lambda}}+H_{s, 1}}{A_{2} e^{-j 2 \pi \frac{\Delta d}{\lambda}} e^{-j 2 \pi \frac{d_{1}(t)}{\lambda}}+H_{s, 2}} \end{aligned}
H2(f,t)H1(f,t)=e−jθoffset(Hs,2+A2e−j2πλd2(t))e−jθoffset(Hs,1+A1e−j2πλd1(t))=A2e−j2πλd1(t)+Δd+Hs,2A1e−j2πλd1(t)+Hs,1=A2e−j2πλΔde−j2πλd1(t)+Hs,2A1e−j2πλd1(t)+Hs,1
其中,
H
1
(
f
,
t
)
H_{1}(f, t)
H1(f,t) 是第一根天线的 CSI,
H
2
(
f
,
t
)
H_{2}(f, t)
H2(f,t) 是第二根天线的 CSI。为了简化方程以便于说明,我们以下表述:
A
1
=
A
,
H
s
,
1
=
B
A_{1}=\mathcal{A}, H_{s, 1}=\mathcal{B}
A1=A,Hs,1=B,
A
2
e
−
j
2
π
Δ
d
λ
=
C
A_{2} e^{-j 2 \pi \frac{\Delta d}{\lambda}}=C
A2e−j2πλΔd=C and
H
s
,
2
=
D
H_{s, 2}=\mathcal{D}
Hs,2=D ;
其中,
e
−
j
2
π
d
1
(
t
)
λ
=
Z
e^{-j 2 \pi \frac{d_{1}(t)}{\lambda}}=Z
e−j2πλd1(t)=Z 表示当
d
1
(
t
)
d_{1}(t)
d1(t)增加时,顺时针旋转的单位圆 ,然后我们可以简化方程式:
H
1
(
f
,
t
)
H
2
(
f
,
t
)
=
A
Z
+
B
C
Z
+
D
\frac{H_{1}(f, t)}{H_{2}(f, t)}=\frac{\mathcal{A Z}+\mathcal{B}}{C Z+\mathcal{D}}
H2(f,t)H1(f,t)=CZ+DAZ+B
这正是 Mobius 变换[5]的形式,前提是
B
C
−
A
D
≠
0
\mathcal{B} C-\mathcal{A} \mathcal{D} \neq 0
BC−AD=0。我们进一步将其分解为以下形式:
H
1
(
f
,
t
)
H
2
(
f
,
t
)
=
B
C
−
A
D
C
2
⋅
1
Z
+
D
C
+
A
C
\frac{H_{1}(f, t)}{H_{2}(f, t)}=\frac{\mathcal{B C}-\mathcal{A} \mathcal{D}}{C^{2}} \cdot \frac{1}{Z+\frac{\mathcal{D}}{C}}+\frac{\mathcal{A}}{C}
H2(f,t)H1(f,t)=C2BC−AD⋅Z+CD1+CA
因此,方程式中的映射由以下变换组成[6]:
(i) Z ↦ Z + D C Z \mapsto \mathcal{Z}+\frac{\mathcal{D}}{\mathcal{C}} Z↦Z+CD, 是加 D C \frac{\mathcal{D}}{\mathcal{C}} CD的translation
(ii) Z ↦ 1 Z Z \mapsto \frac{1}{Z} Z↦Z1, 是复数反演(复数的倒数)
(iii) Z ↦ B C − A D C 2 Z \mathcal{Z} \mapsto \frac{\mathcal{B C}-\mathcal{A D}}{C^{2}} \mathcal{Z} Z↦C2BC−ADZ, 是与复数 B C − A D C 2 \frac{\mathcal{B C}-\mathcal{A D}}{C^{2}} C2BC−AD的相乘
(iv) Z ↦ Z + A C \mathcal{Z} \mapsto \mathcal{Z}+\frac{\mathcal{A}}{C} Z↦Z+CA, 是加 A C \frac{\mathcal{A}}{\mathcal{C}} CA的translation
在操作(iii)中的scaling和rotation ,操作(i)和(iv) 中的translation不会改变 H ( f , t ) H(f, t) H(f,t)的几何形状和旋转方向。因此,重点分析操作(ii)示,即复数反演运算 ( Z ↦ 1 Z ) \left(Z \mapsto \frac{1}{Z}\right) (Z↦Z1)是理解 CSI 商(莫比乌斯变换)效果的关键。
在数学上,复数
z
=
r
e
i
θ
z=r e^{i \theta}
z=reiθ的反演(倒数)为
1
r
e
−
i
θ
\frac{1}{r} e^{-i \theta}
r1e−iθ :新模长是原模长的倒数,新相位与原相位相反。图4展示了如何将单位圆外的一点
z
z
z映射到单位圆内的点
1
z
\frac{1}{z}
z1 。(1)将
z
=
r
e
i
θ
z=r e^{i \theta}
z=reiθ移动到与
z
z
z方向相同但新模长是原模长的倒数的点
1
r
e
i
θ
\frac{1}{r} e^{i \theta}
r1eiθ (2)应用共轭运算(即相对于实轴对称)到点
1
r
e
−
i
θ
\frac{1}{r} e^{-i \theta}
r1e−iθ 。
接下来,我们研究复数反演(记为 Z ~ \tilde{\mathcal{Z}} Z~)对顺时针圆 Z \mathcal{Z} Z的影响,如图 6和图7所示。复数反演不会改变 Z \mathcal{Z} Z 的形状[6] 。但是,它可能会改变 Z \mathcal{Z} Z 的旋转方向。 Z ~ \tilde{\mathcal{Z}} Z~是否与 Z \mathcal{Z} Z方向相同取决于原点 (0, 0) 是否在 Z \mathcal{Z} Z内部。 如图6所示,如果 (0, 0) 不在 Z \mathcal{Z} Z内,则 Z ~ \tilde{\mathcal{Z}} Z~与 Z \mathcal{Z} Z具有相同的旋转方向。 如图7所示,如果 (0, 0) 在 Z \mathcal{Z} Z内,则 Z ~ \tilde{\mathcal{Z}} Z~ 与 Z \mathcal{Z} Z具有相反的旋转方向。
将复数反演应用到 Z + D C \mathcal{Z}+\frac{\mathcal{D}}{C} Z+CD,其中 Z = e − j 2 π d 1 ( t ) λ Z=e^{-j 2 \pi \frac{d_{1}(t)}{\lambda}} Z=e−j2πλd1(t)表示单位圆, D C \frac{\mathcal{D}}{\mathcal{C}} CD是translation。因此,如果 ∣ D C ∣ = ∣ H s , 2 A 2 ∣ > 1 \left|\frac{\mathcal{D}}{C}\right|=\left|\frac{H_{s, 2}}{A_{2}}\right|>1 ∣ ∣CD∣ ∣=∣ ∣A2Hs,2∣ ∣>1(当静态成分大于动态成分时,此式一定成立), Z + D C \mathcal{Z}+\frac{\mathcal{D}}{C} Z+CD则不包含(0,0),所以 1 Z + D C \frac{1}{Z+\frac{D}{C}} Z+CD1与 Z + D C \mathcal{Z}+\frac{\mathcal{D}}{C} Z+CD 具有相同的旋转方向。
由于scaling, rotation and translation 不会改变旋转方向,因此CSI商 H 1 ( f , t ) H 2 ( f , t ) = B C − A D C 2 ⋅ 1 Z + D C + A C \frac{H_{1}(f, t)}{H_{2}(f, t)}=\frac{\mathcal{B C}-\mathcal{A D}}{C^{2}} \cdot \frac{1}{Z+\frac{\mathcal{D}}{C}}+\frac{\mathcal{A}}{C} H2(f,t)H1(f,t)=C2BC−AD⋅Z+CD1+CA具有与 1 Z + D C \frac{1}{Z+\frac{D}{C}} Z+CD1和( Z + D C \mathcal{Z}+\frac{\mathcal{D}}{C} Z+CD )相同的旋转方向。
总而言之:
- 当 LoS 路径信号未被遮挡时,CSI 商具有与 Z + D C \mathcal{Z}+\frac{\mathcal{D}}{C} Z+CD 相同的旋转方向,即动态路径长度增大,CSI商顺时针旋转
- 当 LoS 路径信号被物体衰减时,CSI 商具有与 Z + D C \mathcal{Z}+\frac{\mathcal{D}}{C} Z+CD 相反的旋转方向,即动态路径长度增大,CSI商逆时针旋转
参考文献
[1] Understanding and modeling of wifi signal based human activity recognition
[2] Spotfi: Decimeter level localization using wifi
[3] Dynamic-music: accurate device-free indoor localization
[4] FullBreathe: Full Human Respiration Detection Exploiting Complementarity of CSI Phase and Amplitude of WiFi Signals
[5] Geometry of complex numbers: circle geometry, Moebius transformation, non-euclidean geometry
[6] Visual complex analysis