双采样率下的非均匀采样

双采样率下的非均匀采样

假定一个函数 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)ej2π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)=1h1(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δ(tnts1)ws2(t)=w2(t)n2δ(tnts2)
其中, 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δ(tnts)fsuδ(fufs)
两部分采样子信号的傅里叶变换可以分别写为:
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δ(fufs1)=G(f)H1(f)fs1u1δ(fufs1)Ws2(f)=W2(f)fs2u2δ(fufs2)=G(f)H2(f)fs2u2δ(fufs2)
现在想恢复出 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)+(k1fs1k2fs2)H1(f))=k2fs2G(f)+(k1fs1k2fs2)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=1k1fs1k2fs2=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)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值