文章目录
第九章 数字信号的最佳接收
研究目的:在噪声干扰下,如何最佳地接收有用信号,使接受性能达到最优。
最佳:相对概念,在某种准则下的最佳(差错概率准则、输出信噪比最大准则)。
9.1数字信号的统计特性
假设在数字通信系统中,发送二进制码元1和0。若噪声是高斯白噪声,则其在任意两个时刻上的抽样值都是互相独立的。若噪声是限带高斯白噪声,其最高频率分量小于
f
H
f_H
fH,则在以不小于奈奎斯特速率抽样时,各抽样值之间也是相互独立的。其k维联合概率密度函数为
KaTeX parse error: Expected 'EOF', got '&' at position 25: …n_2,\cdots,n_k)&̲=&f(n_1)f(n_2)\…
其可以改写为
f
k
(
n
1
,
n
2
,
⋯
,
n
k
)
=
1
(
2
π
σ
n
)
k
e
x
p
[
−
1
n
0
∫
0
T
s
n
2
(
t
)
d
t
]
f_k(n_1,n_2,\cdots,n_k)=\frac1{(\sqrt{2\pi}\sigma_n)^k}exp[-\frac1{n_0}\int_0^{T_s}n^2(t)dt]
fk(n1,n2,⋯,nk)=(2πσn)k1exp[−n01∫0Tsn2(t)dt]
因为接收电压为信号电压与噪声电压之和,所以噪声电压可以用接收电压减去信号电压表示,带入到上式中,得到接收电压的k维联合概率密度函数。
9.2数字信号的最佳接收
-
将错误概率最小作为最佳准则。
-
错位概率函数的表达式为
-
P e = P ( 1 ) ∫ − ∞ r 0 ′ f 1 ( r ) d r + P ( 0 ) ∫ r 0 ′ ∞ f 0 ( r ) d r P_e=P(1)\int_{-\infty}^{r_0'}f_1(\bf r)\rm d\bf r\rm+P(0)\int_{r_0'}^{\infty}f_0(\bf r)\rm d\bf r Pe=P(1)∫−∞r0′f1(r)dr+P(0)∫r0′∞f0(r)dr
-
对其求导数可以求得最佳分界点r0:
-
P ( 1 ) P ( 0 ) = f 0 ( r 0 ) f 1 ( r 0 ) \frac{P(1)}{P(0)}=\frac{f_0(\bf r_0\rm)}{f_1(\bf r_0\rm)} P(0)P(1)=f1(r0)f0(r0)
-
当先验概率相等时,最佳分界点为两条曲线的交点。
-
最大似然判决准则:
-
{ P ( 1 ) P ( 0 ) < f 0 ( r 0 ) f 1 ( r 0 ) , 判为0 P ( 1 ) P ( 0 ) > f 0 ( r 0 ) f 1 ( r 0 ) , 判为 1 \begin{cases} \frac{P(1)}{P(0)}<\frac{f_0(\bf r_0\rm)}{f_1(\bf r_0\rm)},\text{判为0}\\ \frac{P(1)}{P(0)}>\frac{f_0(\bf r_0\rm)}{f_1(\bf r_0\rm)},判为1 \end{cases} {P(0)P(1)<f1(r0)f0(r0),判为0P(0)P(1)>f1(r0)f0(r0),判为1
-
最大后验概率准则:
-
{ f r ( 0 ) > f r ( 1 ) , 判为 0 f r ( 0 ) < f r ( 1 ) , 判为 1 \begin{cases} f_r(0)>f_r(1),判为0\\ f_r(0)<f_r(1),判为1 \end{cases} {fr(0)>fr(1),判为0fr(0)<fr(1),判为1
其中 f r ( 1 ) f_r(1) fr(1)为收到r后发送1的条件概率, f r ( 0 ) f_r(0) fr(0)为收到r后发送0的条件概率。
9.3确知数字信号的最佳接收机
假设一个二进制数字通信系统中,两种接收码元的波形的持续时间相同,且能量相同,则若
W
1
+
∫
0
T
B
r
(
t
)
s
1
(
t
)
d
t
<
W
0
+
∫
0
T
B
r
(
t
)
s
0
(
t
)
d
t
s
.
t
.
{
W
0
=
n
0
2
l
n
P
(
0
)
W
1
=
n
0
2
l
n
P
(
1
)
W_1+\int_0^{T_B}r(t)s_1(t)dt<W_0+\int_0^{T_B}r(t)s_0(t)dt\\ s.t. \begin{cases} W_0=\frac{n_0}2lnP(0)\\ W_1=\frac{n_0}2lnP(1) \end{cases}
W1+∫0TBr(t)s1(t)dt<W0+∫0TBr(t)s0(t)dts.t.{W0=2n0lnP(0)W1=2n0lnP(1)
则判定位发送码元是
s
0
(
t
)
s_0(t)
s0(t),反之则为
s
1
(
t
)
s_1(t)
s1(t)。
W
0
W_0
W0和
W
1
W_1
W1可以看作是由先验概率决定的加权因子。
二进制最佳接收机的原理框图为
若二进制信号的先验概率相等,则不需要加上加权因子,积分后直接比较判决。上面的最佳接收机的核心是由相乘和积分构成的相关运算,所以称这种算法为相关接收法。
9.4确知数字信号最佳接收的误码率
- 总误码率为:
KaTeX parse error: Expected 'EOF', got '&' at position 5: P_e&̲=&P(1)P(0/1)+P(…
其中
a
=
n
0
2
l
n
P
(
0
)
P
(
1
)
−
1
2
∫
0
T
B
[
s
1
(
t
)
−
s
0
(
t
)
]
2
d
t
b
=
n
0
2
l
n
P
(
1
)
P
(
0
)
−
1
2
∫
0
T
B
[
s
0
(
t
)
−
s
1
(
t
)
]
2
d
t
a=\frac{n_0}2ln{\frac{P(0)}{P(1)}}-\frac12\int_0^{T_B}[s_1(t)-s_0(t)]^2dt\\ b=\frac{n_0}2ln{\frac{P(1)}{P(0)}}-\frac12\int_0^{T_B}[s_0(t)-s_1(t)]^2dt
a=2n0lnP(1)P(0)−21∫0TB[s1(t)−s0(t)]2dtb=2n0lnP(0)P(1)−21∫0TB[s0(t)−s1(t)]2dt
- 若先验概率相等,a=b,在误码率表达式可以化简为
P e = 1 2 π σ ξ ∫ − ∞ c e − x 2 2 σ ξ 2 d x s . t . c = − 1 2 ∫ 0 T B [ s 0 ( t ) − s 1 ( t ) ] 2 d t P_e=\frac{1}{\sqrt{2\pi}\sigma_\xi}\int_{-\infty}^{c}e^{-\frac{x^2}{2\sigma_\xi^2}}dx\\ s.t.~~~c=-\frac12\int_0^{T_B}[s_0(t)-s_1(t)]^2dt Pe=2πσξ1∫−∞ce−2σξ2x2dxs.t. c=−21∫0TB[s0(t)−s1(t)]2dt
-
当先验概率相等时,对于给定的噪声功率 σ ξ 2 \sigma_\xi^2 σξ2,误码率仅和两种码元波形之差 [ s 0 ( t ) − s 1 ( t ) ] [s_0(t)-s_1(t)] [s0(t)−s1(t)]的能量有关,而与波形本身无关,差别越大,c越小,误码率也越小。(考研复试问过!!!)先验概率不等时的误码率小于相等时的误码率,就误码率而言,先验概率相等是最坏的情况。
-
因为误码率完全取决于信号码元的区别,所以我们可以用码元的相关系数 ρ \rho ρ来定量地描述码元的区别。
ρ = ∫ 0 T B s 0 ( t ) s 1 ( t ) d t E 0 E 1 \rho=\frac{\int_0^{T_B}s_0(t)s_1(t)dt}{\sqrt{E_0E_1}} ρ=E0E1∫0TBs0(t)s1(t)dt
其中, E 0 、 E 1 E_0、E_1 E0、E1为码元的能量。当两种信号码元波形相等时,码元相关系数最大,为1。波形相反时,码元相关系数最小,为-1。
- 当两种码元的能量相等时, E 0 = E 1 = E b E_0=E_1=E_b E0=E1=Eb。
ρ = ∫ 0 T B s 0 ( t ) s 1 ( t ) d t E b P e = 1 2 { 1 − e r f [ E b ( 1 − ρ ) 2 σ ξ ] } = 1 2 e r f c [ E b ( 1 − ρ ) 2 n 0 ] s . t . e r f ( x ) = 2 π ∫ 0 x e − z 2 d z \rho=\frac{\int_0^{T_B}s_0(t)s_1(t)dt}{E_b}\\ P_e=\frac12\{1-erf[\frac{E_b(1-\rho)}{\sqrt2\sigma_\xi}]\}=\frac12erfc[\sqrt{\frac{E_b(1-\rho)}{2n_0}}]\\ s.t.~~~erf(x)=\frac2{\sqrt\pi}\int_0^xe^{-z^2}dz ρ=Eb∫0TBs0(t)s1(t)dtPe=21{1−erf[2σξEb(1−ρ)]}=21erfc[2n0Eb(1−ρ)]s.t. erf(x)=π2∫0xe−z2dz
该公式给出了理论上 二进制 等能量 数字信号误码率的最佳值(最小可能值)。所以,误码率仅与 E b / n 0 E_b/n_0 Eb/n0和相关系数 ρ \rho ρ,与信号波形和噪声功率无关。 E b / n 0 E_b/n_0 Eb/n0实际上等于信号噪声功率比 P s / P n P_s/P_n Ps/Pn。
*9.5随相数字信号的最佳接收
经过信道传输后码元相位带有随机性的信号称为随相信号。
上图是随相信号最佳接收机结构图,随相信号最佳接收机得到的误码率为
P
e
=
1
2
e
x
p
(
−
E
b
/
2
n
0
)
P_e=\frac12exp(-E_b/2n_0)
Pe=21exp(−Eb/2n0)
因为随相信号的相位带有由信道引入的随机变化,所以在接收端不可能采用相干接收方法。
*9.6起伏数字信号的最佳接收
起伏信号是包络随机起伏,相位也随机变化的信号。
起伏信号最佳接收机的结构和随相信号相同,但误码率为
P
e
=
1
2
+
(
E
‾
/
n
0
)
P_e=\frac1{2+(\overline E/n_0)}
Pe=2+(E/n0)1
其中
E
‾
\overline E
E为接收码元的统计平均能量。
9.7实际接收机和最佳接收机的性能比较
接收方式 | 实际接收机的 P e P_e Pe | 最佳接收机的 P e P_e Pe |
---|---|---|
相干接收的2ASK | 1 2 e r f c ( r / 4 ) \frac12erfc(\sqrt r/4) 21erfc(r/4) | 1 2 e r f c ( E b / 4 n 0 ) \frac12erfc(\sqrt E_b/4n_0) 21erfc(Eb/4n0) |
非相干接收的2ASK | 1 2 e x p ( − r / 4 ) \frac12exp(-r/4) 21exp(−r/4) | 1 2 e x p ( − E b / 4 n 0 ) \frac12exp(-E_b/4n_0) 21exp(−Eb/4n0) |
相干接收的2FSK | 1 2 e r f c r / 2 \frac12erfc\sqrt{r/2} 21erfcr/2 | 1 2 e r f c E b / 2 n 0 \frac12erfc\sqrt{E_b/2n_0} 21erfcEb/2n0 |
非相干接收的2FSK | 1 2 e x p ( − r / 2 ) \frac12exp(-r/2) 21exp(−r/2) | 1 2 e x p ( − E b / 2 n 0 ) \frac12exp(-E_b/2n_0) 21exp(−Eb/2n0) |
相干接收的2PSK | 1 2 e r f c r \frac12erfc\sqrt r 21erfcr | 1 2 e r f c E b / n 0 \frac12erfc\sqrt {E_b/n_0} 21erfcEb/n0 |
差分相干接收的2DPSK | 1 2 e x p ( − r ) \frac12exp(-r) 21exp(−r) | 1 2 e x p ( − E b / n 0 ) \frac12exp(-E_b/n_0) 21exp(−Eb/n0) |
相干接收的2DPSK | e r f c r ( 1 − 1 2 e r f c r ) erfc\sqrt r(1-\frac12erfc\sqrt r) erfcr(1−21erfcr) | e r f c E b / n 0 ( 1 − 1 2 e r f c E b / n 0 ) erfc\sqrt {E_b/n_0}(1-\frac12erfc\sqrt {E_b/n_0}) erfcEb/n0(1−21erfcEb/n0) |
- 在实际接收机中的信号噪声功率比r相当于最佳接收机的码元能量与噪声功率谱密度之比 E b / n 0 E_b/n_0 Eb/n0。
- 当系统恰好带宽满足奈奎斯特准则时, E b / n 0 E_b/n_0 Eb/n0就等于信号噪声功率比。
9.8数字信号的匹配滤波接收法
滤波器的作用:使有用信号尽可能强,抑制带外噪声。
最佳线性滤波器设计准则:
1.输出信号与发送信号之间的均方误差最小(维纳滤波器)。
2.令输出信噪比最大的线性滤波器称为匹配滤波器。
- 匹配滤波器的传输特性为
H ( f ) = k S ∗ ( f ) e − j 2 π f t 0 H(f)=kS^*(f)e^{-j2\pi ft_0} H(f)=kS∗(f)e−j2πft0
它等于信号码元频谱的复共轭。
- 匹配滤波器的时域特性为
h ( t ) = k s ( T B − t ) h(t)=ks(T_B-t) h(t)=ks(TB−t)
- 匹配滤波器可以看成是一个计算输入信号自相关函数的相关器。经过匹配滤波器的输出信号为
s o ( t ) = k R ( t − T B ) s_o(t)=kR(t-T_B) so(t)=kR(t−TB)
输出信号码元波形是输入信号码元波形的自相关函数的k倍。k是一个常数,与r0的最大值无关,通常取k=1。
- 最大输出信噪比与信号波形无关,只取决于信号能量E和噪声功率谱密度n0,所以匹配滤波器适用于系带数字信号和已调数字信号。
- 对于二进制确知信号,匹配滤波器构成的接收电路框图如图所示
哪个匹配滤波器的输出抽样值更大,就判决哪个为输出。
- 输入信号经过匹配滤波器,得到的输出波形实际上是一种相关积分运算。用相关运算代替匹配滤波器,可以得到相关接收法框图
PS:实际上,积分器或匹配滤波器具有清零端子,在每次抽样判决的同时进行清零(猝息),以便为下一个码元腾空积分器。
- 匹配滤波器还可以对随相信号和起伏信号进行接收。
9.9最佳基带传输系统
设计目的:无码间串扰的同时,误码率最小。
信道的传输特性 C ( f ) C(f) C(f)不易知,在系统设计时有两种分析方法:第一种是假设其具有理想特性, C ( f ) = 1 C(f)=1 C(f)=1。第二种则考虑信道的非理想特性。
- 最佳化的两个条件:
- H ( ω ) H(\omega) H(ω)无码间串扰,频域条件 ∑ i H ( w + 2 π i T B ) = C ( 常数 ) \sum_iH(w+\frac{2\pi i}{T_B})=C(常数) ∑iH(w+TB2πi)=C(常数)。
- 使系统输出差错概率最小。
9.9.1理想信道的最佳基带传输系统
假设信道传输函数
C
(
f
)
=
1
C(f)=1
C(f)=1,则基带系统的传输特性变为
H
(
f
)
=
G
T
(
f
)
⋅
G
R
(
f
)
H(f)=G_T(f)\cdot G_R(f)
H(f)=GT(f)⋅GR(f)
接收匹配滤波器的传输函数为
G
R
(
f
)
=
G
T
∗
(
f
)
e
−
j
2
π
f
t
0
G_R(f)=G_T^*(f)e^{-j2\pi ft_0}
GR(f)=GT∗(f)e−j2πft0
若没有对接收滤波器的相位进行要求,则
G
R
(
f
)
=
G
T
(
f
)
=
H
1
/
2
(
f
)
G_R(f)=G_T(f)=H^{1/2}(f)
GR(f)=GT(f)=H1/2(f)
- 关于误码率,设基带信号码元为M进制的多电平信号。d为相邻电平间隔的一半,则系统的误码率为
P e = ( 1 − 1 M ) e r f c ( d 2 σ ) P_e=(1-\frac1 M)erfc(\frac d{\sqrt2\sigma}) Pe=(1−M1)erfc(2σd)
平均码元能量为
E
=
d
2
3
(
M
2
−
1
)
E=\frac{d^2}3(M^2-1)
E=3d2(M2−1)
- 根据平均码元能量表达式求出d的表达式,带入到系统的误码率中,得到误码率的最终表达式为
P e = ( 1 − 1 M e r f c [ ( 3 M 2 − 1 ⋅ E n 0 ) 1 / 2 ] ) P_e=(1-\frac 1Merfc[(\frac3{M^2-1}\cdot\frac E{n_0})^{1/2}]) Pe=(1−M1erfc[(M2−13⋅n0E)1/2])
当M=2时,
P
e
=
1
2
e
r
f
c
(
E
/
n
0
)
P_e=\frac12erfc(\sqrt {E/n_0})
Pe=21erfc(E/n0)
与2PSK误码性能一样。
- 对于M进制多电平信号误码率,当误码率较低时,为保持误码率不变,M值增大到2倍,则信噪比约需要增大7dB。
*9.9.2非理想信道的最佳基带传输系统
接收信号码元的频谱为
G
T
(
f
)
⋅
C
(
f
)
G_T(f)\cdot C(f)
GT(f)⋅C(f),则匹配滤波器的传输函数要与接收信号码元的频谱匹配
G
R
′
(
f
)
=
G
T
∗
(
f
)
⋅
C
∗
(
f
)
G_R'(f)=G_T^*(f)\cdot C^*(f)
GR′(f)=GT∗(f)⋅C∗(f)
这时基带传输系统的总传输特性为
H
(
f
)
=
∣
G
T
(
f
)
∣
2
∣
C
(
f
)
∣
2
H(f)=|G_T(f)|^2|C(f)|^2
H(f)=∣GT(f)∣2∣C(f)∣2
该传输特性应满足消除码间串扰的等效理想低通滤波器条件。