Match Filter for PAM Signal

Match Filter for PAM Signal

1. 接收信号$e(t)$：

Tx端发送信号为：
$$s(t)=\sum _{k=-\infty }^{\infty }{a}_{m[k]}g(t-kT)\text{\hspace{1em}(1)}$$
RX端接收信号（暂不考虑衰减）：
$$e(t)=e^{j\phi }\cdot (\sum _{k=-\infty }^{\infty }{a}_{m[k]}g(t-kT-\Delta T))+n(t)\text{\hspace{1em}(2)}$$
其中$e^{j\phi }$为接收端将调制信号解调到基频时，由于发送端与接收端晶振器细微差距而产生的相差，$\Delta T$则为信号经过信道传输后的时延，$n(t)$则为高斯白噪声。

现在我们希望在接收端根据接收信号$e(t)$，恢复出原始二进制码流，具体流程如下图。本文只对接收端Match Filter设计进行讨论。

2. Optimal detection of impulse signal $g(t)$:
为简化分析流程，使用脉冲信号$g(t)$取代原有的发送信号$s(t)$

接收信号在$t=T_v$时的信噪比为：
$${\text{SNR}}_d=\frac{|\widetilde{d(T_v)}{|}^2}{{\sigma }_{n_d}^2}\text{\hspace{1em}(3)}$$
Match Filter设计的核心思路就是使得（3）达到最大值。且必须注意的是，Match filter需要接收端事先知道发送信号作为先验信息（局限性）。
设计矛盾点：
a）信号完整性：
$$\tilde{d}(t)=g(t)\ast h_E(t)\text{\hspace{1em}(4)}$$
$$\tilde{d}(T_v)=g(t)\ast h_E(t){|}_{t=T_v}\text{\hspace{1em}(5)}$$
由上式可知，$H_E(f)$设计的带宽越大，约不容易导致信号被截断，从而导致信号失真。

b）引入噪声：
$${\sigma }_{n_d}^2=\frac{N_o}{2}{\int }_{-\infty }^{\infty }|H_E(f){|}^{2}df\text{\hspace{1em}(6)}$$
由此可知，$H_E(f)$设计的带宽越大，引入的噪声越多。

因此$H_E(f)$设计需要综合考虑这两个因素：

$${\text{SNR}}_d=\frac{|\tilde{d}(T_V){|}^2}{{\sigma }_{n_d}^2}=\frac{{\left|{\int }_{-\infty }^{+\infty }G(f)H_E(f)e^{j2\pi f{T}_V}df\right|}^2}{\frac{N_0}{2}{\int }_{-\infty }^{+\infty }|H_E(f){|}^{2}df}\text{\hspace{1em}(7)}$$
$|\tilde{d}(T_V){|}^2={\left|{\int }_{-\infty }^{+\infty }G(f)H_E(f)e^{j2\pi f{T}_V}df\right|}^2$原因如下：
$$\begin{array}{rl}\tilde{D}(f)& =G(f)\cdot H_E(f)\\ \tilde{d}(T_V)& =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }G(w)\cdot H_E(w)\cdot e^{jw{T}_V}dw\\ & ={\int }_{-\infty }^{+\infty }G(f)\cdot H_E(f)\cdot e^{j2\pi f{T}_V}df\end{array}\text{\hspace{1em}(8)}$$

为了推导出如何设计$H_E(f)$使得（7）达到最大值，需要用到以下不等式：
$${\left|{\int }_{-\infty }^{+\infty }G(f)H_E(f)e^{j2\pi f{T}_V}df\right|}^2\le {\left[{\int }_{-\infty }^{+\infty }|G(f)|\cdot |H_E(f)|\cdot 1df\right]}^2\text{\hspace{1em}(9)}$$
$${\left({\int }_{-\infty }^{+\infty }|G(f)|\cdot |H_E(f)|df\right)}^2\le \left({\int }_{-\infty }^{+\infty }|G(f){|}^{2}df\right)\cdot \left({\int }_{-\infty }^{+\infty }|H_E(f){|}^{2}df\right)\text{\hspace{1em}(10)}$$
因此（7）最大值为：
$${\text{SNR}}_d\le \frac{\left({\int }_{-\infty }^{+\infty }|G(f){|}^{2}df\right)\cdot \left({\int }_{-\infty }^{+\infty }|H_E(f){|}^{2}df\right)}{\frac{N_0}{2}{\int }_{-\infty }^{+\infty }|H_E(f){|}^{2}df}=\frac{{\int }_{-\infty }^{+\infty }|G(f){|}^{2}df}{N_0/2}=\frac{E_g}{N_0/2}\text{\hspace{1em}(11)}$$
所以$H_E(f)$设计为：
$$H_E(f)=\gamma G^{\ast }(f)e^{-j2\pi f{T}_V}\text{\hspace{1em}(12)}$$
$$h_E(t)=\gamma g^{\ast }(T_V-t)\text{\hspace{1em}(13)}$$
其中$\gamma$可以作为一个放大或缩小系数，用于调整$H_E(f)$的幅度。
验证：
$${\text{SNR}}_d=\frac{{\left|{\int }_{-\infty }^{+\infty }G(f)\gamma G^{\ast }(f)e^{-j2\pi f{T}_V}{e}^{j2\pi f{T}_V}df\right|}^2}{\frac{N_0}{2}{\gamma }^2{\int }_{-\infty }^{+\infty }|G(f){|}^{2}df}=\frac{{\gamma }^2{E}_g^2}{\frac{N_0}{2}{\gamma }^2{E}_g}=\frac{E_g}{N_0/2}\text{\hspace{1em}(14)}$$
应用$H_E(f)$后的检测信号为：
$$\begin{array}{rl}\tilde{d}(t)& =\tilde{e}(t)\ast h_E(t)=\tilde{e}(t)\ast \gamma g^{\ast }(T_V-t)\\ & =\left[\sum _{k=-\infty }^{+\infty }{a}_m[k]g(t-kT)\right]\ast \gamma g^{\ast }(T_V-t)\\ & =\gamma \sum _{k=-\infty }^{+\infty }{a}_m[k]\left[g(t-kT)\ast g^{\ast }(T_V-t)\right]\end{array}\text{\hspace{1em}(15)}$$

$$\begin{array}{rl}g(t-kT)\ast g^{\ast }(T_V-t)& ={\int }_{-\infty }^{+\infty }{g}^{\ast }(T_V-t+\tau )g(\tau -kT)d\tau \\ & ={\int }_{-\infty }^{+\infty }{g}^{\ast }(T_V-t+{\tau }^{\prime })g({\tau }^{\prime }-kT)d{\tau }^{\prime }\\ \text{Substitution: }{\tau }^{\prime }=\tau -t+T_V,d\tau =d{\tau }^{\prime }:& \\ & ={\int }_{-\infty }^{+\infty }g({\tau }^{\prime }+t-T_V-kT)g^{\ast }({\tau }^{\prime })d{\tau }^{\prime }\\ & ={\phi }_{gg}(t-kT-T_V)\end{array}\text{\hspace{1em}(16)}$$

ISI Free Detection：
条件：脉冲信号$g(t)$需要满足Nyquist第一定律：
φ g g [ λ T ] = ∫ − ∞ + ∞ g ( t ′ + λ T ) g ∗ ( t ′ )   d t ′ = def { E g for  λ = 0 0    for  λ ∈ Z ∖ { 0 } = δ [ λ ]   E g (17) \varphi_{gg}[\lambda T] = \int_{-\infty}^{+\infty} g(t' + \lambda T) g^*(t') \, dt' \overset{\text{def}}{=} \begin{cases} E_g \quad \text{for } \lambda = 0 \\ 0 \quad \ \ \ \text{for } \lambda \in \mathbb{Z} \setminus \{0\} \end{cases} = \delta[\lambda] \, E_g \tag{17} φgg​[λT]=∫−∞+∞​g(t′+λT)g∗(t′)dt′=def{Eg​for λ=00   for λ∈Z∖{0}​=δ[λ]Eg​(17)

检测信号在$t=\ell T$时为:
$$\begin{array}{rl}\tilde{d}[\ell ]& \stackrel{\text{def}}{=}\tilde{d}(\ell T+T_V)\stackrel{\text{def}}{=}{a}_m[\ell ]\\ \tilde{d}(\ell T+T_V)& =\gamma \sum _{k=-\infty }^{+\infty }{a}_m[k]{\phi }_{gg}(\ell T+T_V-kT-T_V)\\ & =\gamma \sum _{k=-\infty }^{+\infty }{a}_m[k]{\phi }_{gg}((\ell -k)T)\stackrel{\text{def}}{=}{a}_m[\ell ]\end{array}\text{\hspace{1em}(18)}$$
所以$\gamma =\frac{1}{E_g}$:
$$H_E(f)=\frac{1}{E_g}{G}^{\ast }(f)e^{-j2\pi f{T}_V}\text{\hspace{1em}(19)}$$