双采样率下的非均匀采样
假定一个函数
g
(
t
)
g(t)
g(t)的傅里叶变换为
G
(
f
)
=
∫
−
∞
∞
g
(
t
)
e
−
j
2
π
f
t
d
t
G(f)=\int_{-\infty}^{\infty}g(t)e^{-j2\pi ft}dt
G(f)=∫−∞∞g(t)e−j2πftdt
定义一个门函数
h
1
(
t
)
h_1(t)
h1(t),再定义一个与之互补的门函数
h
2
(
t
)
h_2(t)
h2(t),有
h
2
(
t
)
=
1
−
h
1
(
t
)
h_2(t)=1-h_1(t)
h2(t)=1−h1(t)
可见
h
1
(
t
)
h
2
(
t
)
=
0
h_1(t)h_2(t)=0
h1(t)h2(t)=0,其傅里叶变换分别为:
h
1
(
t
)
↔
H
1
(
f
)
h
2
(
t
)
↔
H
2
(
f
)
=
δ
(
f
)
−
H
1
(
f
)
h_1(t)\leftrightarrow H_1(f)\\ h_2(t)\leftrightarrow H_2(f)=\delta(f)-H_1(f)
h1(t)↔H1(f)h2(t)↔H2(f)=δ(f)−H1(f)
将它们分别与
g
(
t
)
g(t)
g(t)相乘,得到两个子信号:
w
1
(
t
)
=
g
(
t
)
h
1
(
t
)
w
2
(
t
)
=
g
(
t
)
h
2
(
t
)
w_1(t)=g(t)h_1(t)\\ w_2(t)=g(t)h_2(t)
w1(t)=g(t)h1(t)w2(t)=g(t)h2(t)
现在,以采样率
f
s
1
f_{s1}
fs1对
w
1
(
t
)
w_1(t)
w1(t)进行采样,以采样率
f
s
2
f_{s2}
fs2对
w
2
(
t
)
w_2(t)
w2(t)进行采样,得到:
w
s
1
(
t
)
=
w
1
(
t
)
∑
n
1
δ
(
t
−
n
t
s
1
)
w
s
2
(
t
)
=
w
2
(
t
)
∑
n
2
δ
(
t
−
n
t
s
2
)
w_{s1}(t)=w_1(t)\sum_{n1}\delta(t-nt_{s1})\\ w_{s2}(t)=w_2(t)\sum_{n2}\delta(t-nt_{s2})
ws1(t)=w1(t)n1∑δ(t−nts1)ws2(t)=w2(t)n2∑δ(t−nts2)
其中,
n
1
,
n
2
n_1,n_2
n1,n2为对应的采样点数,
t
s
1
=
1
f
s
1
t_{s1}=\frac{1}{f_{s1}}
ts1=fs11,
t
s
2
=
1
f
s
2
t_{s2}=\frac{1}{f_{s2}}
ts2=fs21.
我们知道:
∑
n
δ
(
t
−
n
t
s
)
↔
f
s
∑
u
δ
(
f
−
u
f
s
)
\sum_{n}\delta(t-nt_{s})\leftrightarrow f_s\sum_{u}\delta(f-uf_s)
n∑δ(t−nts)↔fsu∑δ(f−ufs)
两部分采样子信号的傅里叶变换可以分别写为:
W
s
1
(
f
)
=
W
1
(
f
)
∗
f
s
1
∑
u
1
δ
(
f
−
u
f
s
1
)
=
G
(
f
)
∗
H
1
(
f
)
∗
f
s
1
∑
u
1
δ
(
f
−
u
f
s
1
)
W
s
2
(
f
)
=
W
2
(
f
)
∗
f
s
2
∑
u
2
δ
(
f
−
u
f
s
2
)
=
G
(
f
)
∗
H
2
(
f
)
∗
f
s
2
∑
u
2
δ
(
f
−
u
f
s
2
)
W_{s1}(f)=W_1(f)*f_{s1}\sum_{u_1}\delta(f-uf_{s1})=G(f)*H_1(f)*f_{s1}\sum_{u_1}\delta(f-uf_{s1})\\ W_{s2}(f)=W_2(f)*f_{s2}\sum_{u_2}\delta(f-uf_{s2})=G(f)*H_2(f)*f_{s2}\sum_{u_2}\delta(f-uf_{s2})
Ws1(f)=W1(f)∗fs1u1∑δ(f−ufs1)=G(f)∗H1(f)∗fs1u1∑δ(f−ufs1)Ws2(f)=W2(f)∗fs2u2∑δ(f−ufs2)=G(f)∗H2(f)∗fs2u2∑δ(f−ufs2)
现在想恢复出
G
(
f
)
G(f)
G(f),进而恢复出
g
(
t
)
g(t)
g(t),我们只保留其基带成分:
W
s
1
′
(
f
)
=
G
(
f
)
∗
H
1
(
f
)
∗
f
s
1
δ
(
f
)
=
G
(
f
)
∗
f
s
1
H
1
(
f
)
W
s
2
′
(
f
)
=
G
(
f
)
∗
H
2
(
f
)
∗
f
s
2
δ
(
f
)
=
G
(
f
)
∗
f
s
2
H
2
(
f
)
W'_{s1}(f)=G(f)*H_1(f)*f_{s1}\delta(f)=G(f)*f_{s1}H_1(f)\\ W'_{s2}(f)=G(f)*H_2(f)*f_{s2}\delta(f)=G(f)*f_{s2}H_2(f)
Ws1′(f)=G(f)∗H1(f)∗fs1δ(f)=G(f)∗fs1H1(f)Ws2′(f)=G(f)∗H2(f)∗fs2δ(f)=G(f)∗fs2H2(f)
现构造一个
W
s
1
′
(
f
)
W'_{s1}(f)
Ws1′(f)和
W
s
2
′
(
f
)
W'_{s2}(f)
Ws2′(f)的线性组合,令
G
(
f
)
=
k
1
W
s
1
′
(
f
)
+
k
2
W
s
2
′
(
f
)
G(f)=k_1W'_{s1}(f)+k_2W'_{s2}(f)
G(f)=k1Ws1′(f)+k2Ws2′(f)
可得
G
(
f
)
=
G
(
f
)
∗
(
k
1
f
s
1
H
1
(
f
)
+
k
2
f
s
2
H
2
(
f
)
)
=
G
(
f
)
∗
(
k
2
f
s
2
δ
(
f
)
+
(
k
1
f
s
1
−
k
2
f
s
2
)
H
1
(
f
)
)
=
k
2
f
s
2
G
(
f
)
+
(
k
1
f
s
1
−
k
2
f
s
2
)
H
1
(
f
)
∗
G
(
f
)
G(f)=G(f)*(k_1f_{s1}H_1(f)+k_2f_{s2}H_2(f))\\=G(f)*(k_2f_{s2}\delta(f)+(k_1f_{s1}-k_2f_{s2})H_1(f))\\ =k_2f_{s2}G(f)+(k_1f_{s1}-k_2f_{s2})H_1(f)*G(f)
G(f)=G(f)∗(k1fs1H1(f)+k2fs2H2(f))=G(f)∗(k2fs2δ(f)+(k1fs1−k2fs2)H1(f))=k2fs2G(f)+(k1fs1−k2fs2)H1(f)∗G(f)
可见,只要合理选择
f
s
1
f_{s1}
fs1和
f
s
2
f_{s2}
fs2,即可保证
H
1
(
f
)
H_1(f)
H1(f)的任意性,且恢复出
G
(
f
)
G(f)
G(f).
令
k
2
f
s
2
=
1
k
1
f
s
1
−
k
2
f
s
2
=
0
k_2f_{s2}=1\\ k_1f_{s1}-k_2f_{s2}=0
k2fs2=1k1fs1−k2fs2=0
可得
k
1
=
1
f
s
1
k_1 = \frac{1}{f_{s1}}
k1=fs11,
k
2
=
1
f
s
2
k_2 = \frac{1}{f_{s2}}
k2=fs21.
G
(
f
)
=
W
s
1
′
(
f
)
f
s
1
+
W
s
2
′
(
f
)
f
s
2
G(f)=\frac{W'_{s1}(f)}{f_{s1}}+\frac{W'_{s2}(f)}{f_{s2}}
G(f)=fs1Ws1′(f)+fs2Ws2′(f)