1、8.1 匹配滤波器匹配滤波器8.2 最小差错概率接收准则最小差错概率接收准则8.3 确知信号的最佳接收机确知信号的最佳接收机8.4 随相信号的最佳接收机随相信号的最佳接收机8.5 最佳接收机性能比较最佳接收机性能比较8.6 最佳基带传输系统最佳基带传输系统 第第 8 章章 数字信号的最佳接收数字信号的最佳接收返回主目录第第8章章 数字信号的最佳接收数字信号的最佳接收8.1匹匹 配配 滤滤 波波 器器 在数字通信系统中，滤波器是其中重要部件之一， 滤波器特性的选择直接影响数字信号的恢复。在数字信号接收中， 滤波器的作用有两个方面，第一是使滤波器输出有用信号成分尽可能强；第二是抑制信号带外噪声，使
2、滤波器输出噪声成分尽可能小，减小噪声对信号判决的影响。 通常对最佳线性滤波器的设计有两种准则：一种是使滤波器输出的信号波形与发送信号波形之间的均方误差最小，由此而导出的最佳线性滤波器称为维纳滤波器；另一种是使滤波器输出信噪比在某一特定时刻达到最大，由此而导出的最佳线性滤波器称为匹配滤波器。在数字通信中，匹配滤波器具有更广泛的应用。 由第7章分析的数字信号解调过程我们知道，解调器中抽样判决以前各部分电路可以用一个线性滤波器来等效，接收过程等效原理图如图 8 - 1 所示。图中，s(t)为输入数字信号， 信道特性为加性高斯白噪声信道，n(t)为加性高斯白噪声， H()为滤波器传输函数。 由数字信号
3、的判决原理我们知道，抽样判决器输出数据正确与否，与滤波器输出信号波形和发送信号波形之间的相似程度无关，也即与滤波器输出信号波形的失真程度无关， 而只取决于抽样时刻信号的瞬时功率与噪声平均功率之比， 即信噪比。信噪比越大，错误判决的概率就越小；反之，信噪比越小，错误判决概率就越大。 图 8 1 数字信号接收等效原理图H()判决s(t)n(t)r(t)y(t) t t0输出SN( )o 因此，为了使错误判决概率尽可能小，就要选择滤波器传输特性使滤波器输出信噪比尽可能大的滤波器。当选择的滤波器传输特性使输出信噪比达到最大值时，该滤波器就称为输出信噪比最大的最佳线性滤波器。下面就来分析当滤波器具有什么
4、样的特性时才能使输出信噪比达到最大。 分析模型如图 8 - 1 所示。 设输出信噪比最大的最佳线性滤波器的传输函数为H(), 滤波器输入信号与噪声的合成波为 r(t)=s(t)+n(t) 式中， s(t)为输入数字信号， 其频谱函数为S()。 n(t)为高斯白噪声， 其双边功率谱密度为 。20n 由于该滤波器是线性滤波器，满足线性叠加原理，因此滤波器输出也由输出信号和输出噪声两部分组成，即 y(t)=so(t)+no(t) (8.1 - 2)式中输出信号的频谱函数为So()，其对应的时域信号为 so(t)= dwewHwsdwewsjwtjwt)()(21)(210滤波器输出噪声的平均功率为d
5、wwHwpdwwPNinn20)()(21)(210dwwHndwwHwn2020)(4)()(221 在抽样时刻t0，线性滤波器输出信号的瞬时功率与噪声平均功率之比为 ro= 由式(8.1 - 5)可见， 滤波器输出信噪比ro与输入信号的频谱函数S()和滤波器的传输函数H()有关。在输入信号给定的情况下，输出信噪比ro只与滤波器的传输函数H()有关。使输出信噪比ro达到最大的传输函数H()就是我们所要求的最佳滤波器的传输函数。式(8.1 - 5)是一个泛函求极值的问题，采用施瓦兹(Schwartz)不等式可以容易地解决该问题。 施瓦兹不等式为dwwHndwewswHNtsjwt200200)
6、(4)()(21)(0dwwYdwwXdwwYwX222)(21)(21)()(21式中， X()和Y()都是实变量的复函数。 当且仅当 X()=KY*() 时式(8.1 - 6)中等式才能成立。式(8.1 - 7)中K为任意常数。 将施瓦兹不等式用于式(8.1 - 5)， 并令 X()=H() (8.1 - 8) Y()=S()ejt0 (8.1 - 9) 可得 ro= dwwHndwewSwHjwt20)(4)()(2102)(21)(4)()(41020220ndwwsdwwHndwewsdwwHjwrt根据帕塞瓦尔定理有 EdttSdwwS)()(2122式中E为输入信号的能量。 代入
7、式(8.1 - 10)有 ro 式(8.1 - 12)说明，线性滤波器所能给出的最大输出信噪比为 r0max= (8.1 - 13)根据施瓦兹不等式中等号成立的条件X()=KY*()， 可得不等式(8.1 - 10)中等号成立的条件为 H()=KS*()e-jt0 式中，K为常数，通常可选择为K=1。S*()是输入信号频谱函数S()的复共轭。式(8.1 - 14)就是我们所要求的最佳线性滤波器的传输函数，该滤波器在给定时刻t0能获得最大输出信噪比 。 2E02nE02nE 这种滤波器的传输函数除相乘因子Ke-jt0外，与信号频谱的复共轭相一致，所以称该滤波器为匹配滤波器。 从匹配滤波器传输函数
8、H()所满足的条件，我们也可以得到匹配滤波器的单位冲激响应h(t)：dweewKSdwewHthjwtjwtjwt0)(21)(21)(dsdwekdeSKttjwjwt)(21)(2)(0)()()(00ttKdttskS即匹配滤波器的单位冲激响应为 h(t)=Ks(t0-t) 式(8.1 - 16)表明，匹配滤波器的单位冲激响应h(t)是输入信号s(t)的镜像函数，t0为输出最大信噪比时刻。其形成原理如图 8 - 2 所示。 对于因果系统， 匹配滤波器的单位冲激响应h(t)应满足: h(t)= Ks(t0-t), t00, t0 为了满足式(8.1 - 17)的条件， 必须有：图8-2 匹
9、配滤波器单位冲激响应原理s(t)OTth(t)Ott0 s(t0-t)=0, t0 (8.1 - 18) s(t)=0, t0-t0 或 tt0 (8.1 - 19)上式条件说明，对于一个物理可实现的匹配滤波器，其输入信号s(t)必须在它输出最大信噪比的时刻t0之前结束。也就是说，若输入信号在T时刻结束，则对物理可实现的匹配滤波器， 其输出最大信噪比时刻t0必须在输入信号结束之后，即t0T。 对于接收机来说，t0是时间延迟，通常总是希望时间延迟尽可能小，因此一般情况可取t0=T。 若输入信号为s(t)， 则匹配滤波器的输出信号为 so(t)=s(t)*h(t)= dhts)()(dtKtsS)
10、()(0 式中， R(t)为输入信号s(t)的自相关函数。 上式表明， 匹配滤波器的输出波形是输入信号s(t)的自相关函数的K倍。因此， 匹配滤波器可以看成是一个计算输入信号自相关函数的相关器，其在t0时刻得到最大输出信噪比romax= 。 由于输出信噪比与常数K无关，所以通常取K=1。 例 8 - 1设输入信号如图 8 - 3(a)所示，试求该信号的匹配滤波器传输函数和输出信号波形。 02nEs(t)= 1, 0t 0, 其他2T输入信号s(t)的频谱函数为图8-3 信号时间波形(a)(b)so(t)OTt(c)2T2T2T23Ts(t)0Tt2T1h(t)0Tt1 S()= 20)(Tjw
11、tjwtdtedtets0)1 (12jwtwTjeejw匹配滤波器的传输函数为H()=S*()e-jt0= jwtwTjeejw ) 1(12匹配滤波器的单位冲激响应为 h(t)=s(t0-t)取t0=T，则有TjwtwTjeejwWH) 1(1)(2 h(t)=s(T-t)匹配滤波器的单位冲激响应如图 8 - 3(b)所示。 (2) 由式(8.1 - 21)可得匹配滤波器的输出为 so(t)=R(t-t0)= dxttxsxs)()(0,2tT,23tT, 0TtT223TtT其他= 匹配滤波器的输出波形如图 8 - 3(c)所示。可见，匹配滤波器的输出在t=T时刻得到最大的能量E= 。
12、2T8.2 最小差错概率接收准则最小差错概率接收准则 8.2.1数字信号接收的统计模型数字信号接收的统计模型 在数字信号的最佳接收分析中，我们不是采用先给出接收机模型然后分析其性能的分析方法，而是从数字信号接收统计模型出发，依据某种最佳接收准则，推导出相应的最佳接收机结构，然后再分析其性能。 数字通信系统的统计模型如图 8 - 4 所示。 图中消息空间、 信号空间、噪声空间、观察空间及判决空间分别代表消息、 发送信号、噪声、接收信号波形及判决结果的所有可能状态的集合。各个空间的状态用它们的统计特性来描述。 图 8 4 数字通信系统的统计模型XSY判决规则Rn消息空间信号空间观察空间判决空间噪声
13、空间 在数字通信系统中， 消息是离散的状态， 设消息的状态集合为X=x1, x2, , xm (8.2 - 1)若消息集合中每一状态的发送是统计独立的， 第i个状态xi的出现概率为P(xi)， 则消息X的一维概率分布为X1 x2 xmP(x1) P(x2) P(xm) 根据概率的性质有 1)(1niixp 若消息各状态x1, x2, , xm出现的概率相等，则有 P(x1)=P(x2)=P(xm)= (8.2 - 3) m1 消息是各种物理量， 本身不能直接在数字通信系统中进行传输，因此需要将消息变换为相应的电信号s(t)，用参数S来表示。将消息变换为信号可以有各种不同的变换关系，通常最直接的
14、方法是建立消息与信号之间一一对应的关系，即消息xi与信号si(i=1, 2， ， m)相对应。 这样，信号集合S也由m个状态所组成，即 S=s1, s2, , sm 并且信号集合各状态出现概率与消息集合各状态出现概率相等，即 P(s1)=P(x1) P(s2)=P(x2) P(sm)=P(xm) 同时也有1)(1niisp 若消息各状态出现的概率相等， 则有 P(s1)=P(s2)=P(sm)= (8.2 - 6)P(si)是描述信号发送概率的参数，通常称为先验概率， 它是信号统计检测的第一数据。 信道特性是加性高斯噪声信道，噪声空间n是加性高斯噪声。在前面各章分析系统抗噪声性能时，用噪声的一
15、维概率密度函数来描述噪声的统计特性, 在本章最佳接收中，为了更全面地描述噪声的统计特性，采用噪声的多维联合概率密度函数。噪声n的k维联合概率密度函数为 f(n)=f(n1, n2, , nk) (8.2 - 7)式中，n1, n2, , nk为噪声n在各时刻的可能取值。 根据随机信号分析理论我们知道，若噪声是高斯白噪声， 则它在任意两个时刻上得到的样值都是互不相关的，同时也是统计独立的；若噪声是带限高斯型的，按抽样定理对其抽样，则它在抽样时刻上的样值也是互不相关的， 同时也是统计独立的。根据随机信号分析，若随机信号各样值是统计独立的，则其k维联合概率密度函数等于其k个一维概率密度函数的乘积，即
16、 f(n1, n2, , nk)=f(n1)f(n2)f(nk) 式中， f(ni)是噪声n在ti时刻的取值ni的一维概率密度函数， 若ni的均值为零，方差为2n，则其一维概率密度函数为222exp21)(nininnf噪声n的k维联合概率密度函数为kiinikninnnf12222exp)2(1)(根据帕塞瓦尔定理， 当k很大时有dttnnnTkiin)(11020122 信号通过信道叠加噪声后到达观察空间， 观察空间的观察波形为 y=n+s 由于在一个码元期间T内， 信号集合中各状态s1, s2, , sm 中之一被发送，因此在观察期间T内观察波形为 y(t)=n(t)+si(t) (i=
17、1, 2, , m) 由于n(t)是均值为零， 方差为2n的高斯过程，则当出现信号si(t)时, y(t)的概率密度函数fsi(y)可表示为),.,2 , 1()()(1exp)2(1020)(midttstynfTiknysi fsi(y)称为似然函数，它是信号统计检测的第二数据。 根据y(t)的统计特性，按照某种准则，即可对y(t)作出判决， 判决空间中可能出现的状态r1, r2, , rm与信号空间中的各状态s1, s2, , sm相对应。 8.2.2最佳接收准则最佳接收准则 在数字通信系统中，最直观且最合理的准则是“最小差错概率”准则。由于在传输过程中，信号会受到畸变和噪声的干扰，发送
18、信号si(t)时不一定能判为ri出现，而是判决空间的所有状态都可能出现。这样将会造成错误接收，我们期望错误接收的概率愈小愈好。 在噪声干扰环境中，按照何种方法接收信号才能使得错误概率最小？我们以二进制数字通信系统为例分析其原理。 在二进制数字通信系统中，发送信号只有两种状态，假设发送信号s1(t)和s2(t)的先验概率分别为P(s1)和P(s2)，s1(t)和s2(t)在观察时刻的取值分别为a1和a2，出现s1(t)信号时y(t)的概率密度函数fs1(y)为 fs1(y)= dtatynTkn2010)(1exp)2(1同理，出现s2(t)信号时y(t)的概率密度函数fs2(y)为dtatyn
19、yfTkns2020)(1exp)2(1)(2 fs1(y)和fs2(y)的曲线如图 8 - 5 所示。 若在观察时刻得到的观察值为yi，可依概率将yi判为r1或r2。在yi附近取一小区间a，yi在区间a内属于r1的概率为 图 8- 5 fs1(y)和fs2(y)的曲线图fs1(y)fs2(y)ya2aa1yiasasdyyfqdyyf)()(212yi在相同区间a内属于r2的概率为dyyfqas)(22可以看出， q1= 即yi属于r1的概率大于yi属于r2的概率。因此，依大概率应将yi判为r1出现。 由于fs1(y)和fs2(y)的单调性质，图 8 - 5 所示的判决过程可以简化为图 8
20、- 6 所示的判决过程。 dyyfqas)(12图 8 6 判决过程示意图 ya2y 0a1PPfs1( y)fs2( y)r1r2 根据fs1(y)和fs2(y)的单调性质， 在图 8 - 6 中y坐标上可以找到一个划分点y0。在区间(-, y0， q1q2；在区间(y0, )， q1q2。 根据图 8 - 6 所分析的判决原理，当观察时刻得到的观察值yi(-, y0)时，判为r1出现；若观察时刻得到的观察值yi(y0, )时，判为r2出现。 如果发送的是s1(t)，但是观察时刻得到的观察值yi落在(y0,)区间, 被判为r2出现，这时将造成错误判决，其错误概率为 Ps1(s2)= dyyf
21、ys)(01 同理， 如果发送的是s2(t)， 但是观察时刻得到的观察值yi落在(-, y0)区间, 被判为r1出现，这时也将造成错误判决，其错误概率为 Ps2(s1)= dyyfys)(02此时系统总的误码率为Pe= p(s1)ps1(s2)+p(s2)ps2(s1) =p(s1)020)()()(21ysysdyyfspdyyf 由式(8.2 - 21)可以看出， 系统总的误码率与先验概率、 似然函数及划分点 有关，0y 在先验概率和似然函数一定的情况下，系统总的误码率Pe是划分点y0的函数。不同的y0将有不同的Pe，我们希望选择一个划分点y0使误码率Pe达到最小。使误码率Pe达到最小的划
22、分点y0称为最佳划分点。y0可以通过求Pe的最小值得到。 即00ype-P(s1)fs1(y0)+P(s2)fs2(y0)=0 (8.2 - 23) 由此可得最佳划分点将满足如下方程:)()()()(120021spspyfyfss式中y0即为最佳划分点。 如果观察时刻得到的观察值y小于最佳划分点y0，应判为r1出现，此时式(8.2 -24)左边大于右边；如果观察时刻得到的观察值y大于最佳划分点y0，应判为r2出现，此时式(8.2 - 24)右边大于左边。因此，为了达到最小差错概率，可以按以下规则进行判决:,)()()()(120021spspyfyfss,)()()()(120021spsp
23、yfyfss判为 r1( 即s1)判为 r2( 即s2) 以上判决规则称为似然比准则。在加性高斯白噪声条件下，似然比准则和最小差错概率准则是等价的。 当s1(t)和s2(t)的发送概率相等时，即P(s1)=P(s2)时，则有 fs1(y)fs2(y), 判为r1(即s1) fs1(y)fs2(y), 判为r2(即s2) 上式判决规则称为最大似然准则，其物理概念是，接收到的波形y中，哪个似然函数大就判为哪个信号出现。 以上判决规则可以推广到多进制数字通信系统中，对于m个可能发送的信号，在先验概率相等时的最大似然准则为 fsi(y)fsj(y), 判为si(i=1, 2, , m; j=1, 2,
24、 , m; ij) (8.2 - 27) 最小差错概率准则是数字通信系统最常采用的准则， 除此之外，贝叶斯(Bayes)准则、尼曼-皮尔逊(Neyman-Pearson)准则、 极大极小准则等有时也被采用。 8.3确知信号的最佳接收机确知信号的最佳接收机 在数字通信系统中，接收机输入信号根据其特性的不同可以分为两大类，一类是确知信号，另一类是随参信号。所谓确知信号是指一个信号出现后，它的所有参数(如幅度、频率、 相位、到达时刻等)都是确知的。如数字信号通过恒参信道到达接收机输入端的信号。在随参信号中，根据信号中随机参量的不同又可细分为随机相位信号、随机振幅信号和随机振幅随机相位信号(又称起伏信
25、号)。本节讨论确知信号的最佳接收问题。 信号统计检测是利用概率和数理统计的工具来设计接收机。所谓最佳接收机设计是指在一组给定的假设条件下，利用信号检测理论给出满足某种最佳准则接收机的数学描述和组成原理框图，而不涉及接收机各级的具体电路。本节分析中所采用的最佳准则是最小差错概率准则。 8.3.1二进制确知信号最佳接收机结构二进制确知信号最佳接收机结构 接收端原理图如图 8 - 7 所示。设到达接收机输入端的两个确知信号分别为s1(t)和s2(t)，它们的持续时间为(0, T)，且有相等的能量，即 E=E1= tTdttsEdttS0222021)()(噪声n(t)是高斯白噪声，均值为零，单边功率
26、谱密度为n0。 要求设计的接收机能在噪声干扰下以最小的错误概率检测信号。 根据上一节的分析我们知道，在加性高斯白噪声条件下， 最小差错概率准则与似然比准则是等价的。因此，我们可以直接利用式(8.2 - 25)似然比准则对确知信号作出判决。 在观察时间(0，T)内，接收机输入端的信号为s1(t)和s2(t)， 合成波为图 8 7 接收端原理最佳接收机s(t)n(t)输出y(t)= s1(t)+n(t),发送s1(t)时, s2(t)+n(t),发送s2(t)时, 发送 s1(t)时发送 s2(t)时 由上一节分析可知，当出现s1(t)或s2(t)时观察空间的似然函数分别为 fs1(y)= dtt
27、stynTkn0210)()(1exp)2(1dttstynyfTkns0220)()(1exp)2(1)(2其似然比判决规则为TknTknssdttstyndttstynyfyf0220021000)()(1exp)2(1)()(1exp)2(1)()(21 则判为s2(t)出现。式中，P(s1)和P(s2)分别为发送s1(t)和s2(t)的先验概率。整理式(8.3 - 5)和(8.3 - 6)可得TTdttstyUdttstyU022011)()()()(判为s1(t)出现，而TTdttstyUdttstyU022011)()()()(则判为s2(t)出现。 式中： U1= lnP(s1)
28、U2= lnP(s2)20n20n 在先验概率P(s1)和P(s2)给定的情况下，U1和U2都为常数。 根据式(8.3 - 7)和式(8.3 - 8)所描述的判决规则， 可得到最佳接收机的结构如图8 - 8 所示，其中比较器是比较抽样时刻t=T时上下两个支路样值的大小。这种最佳接收机的结构是按比较观察波形y(t)与s1(t)和s2(t)的相关性而构成的，因而称为相关接收机。其中相乘器与积分器构成相关器。接收过程是分别计算观察波形y(t)与s1(t)和s2(t)的相关函数，在抽样时刻t=T，y(t)与哪个发送信号的相关值大就判为哪个信号出现。 图 8 8 二进制确知信号最佳接收机结构积分器y(t
29、)s1(t)输出积分器s2(t)U1U2比较器 如果发送信号s1(t)和s2(t)的出现概率相等，即P(s1)=P(s2)，由式(8.3 - 9)可得U1=U2。此时，图 8 - 8 中的两个相加器可以省去，则先验等概率情况下的二进制确知信号最佳接收机简化结构如图 8 - 9 所示。 由 8.1 节匹配滤波器分析我们知道，匹配滤波器可以看成是一个计算输入信号自相关函数的相关器。设发送信号为s(t)， 则匹配滤波器的单位冲激响应为 h(t)=s(T-t) 若匹配滤波器输入合成波为 y(t)=s(t)+n(t) 图 8- 9 二进制确知信号最佳接收机简化结构 积分器y(t)s1(t)输出积分器s2
30、(t)比较器则匹配滤波器的输出在抽样时刻t=T时的样值为 u0(t)= y(t)s(t)dt T0 由式(8.3 - 12)可以看出匹配滤波器在抽样时刻t=T时的输出样值与最佳接收机中相关器在t=T时的输出样值相等,因此， 可以用匹配滤波器代替相关器构成最佳接收机，其结构如图 8 - 10 所示。 在最小差错概率准则下，相关器形式的最佳接收机与匹配滤波器形式的最佳接收机是等价的。另外，无论是相关器还是匹配滤波器形式的最佳接收机, 它们的比较器都是在t=T时刻才作出判决，也即在码元结束时刻才能给出最佳判决结果。因此，判决时刻的任何偏差都将影响接收机的性能。 图 8- 10 匹配滤波器形式的最佳接
31、收机h1(t) s1(T t)(0 t T)y(t)输出比较器h2(t) s2(T t)(0 t T) 8.3.2二进制确知信号最佳接收机误码性能二进制确知信号最佳接收机误码性能 由上一节分析可知，相关器形式的最佳接收机与匹配滤波器形式的最佳接收机是等价的，因此可以从两者中的任一个出发来分析最佳接收机的误码性能。下面从相关器形式的最佳接收机角度来分析这个问题。 最佳接收机结构如图 8 - 8 所示，输出总的误码率为 Pe=P(s1)Ps1(s2)+P(s2)Ps2(s1)其中, P(s1)和P(s2)是发送信号的先验概率。Ps1(s2)是发送s1(t)信号时错误判决为s2(t)信号出现的概率；
32、Ps2(s1)是发送s2(t)信号时错误判决为s1(t)信号出现的概率。分析Ps1(s2)与Ps2(s1)的方法相同，我们以分析Ps1(s2)为例。 设发送信号为s1(t)，接收机输入端合成波为 y(t)=s1(t)+n(t) (8.3 - 14)其中, n(t)是高斯白噪声，其均值为零，方差为2n。若 U1+ y(t)s1(t)dtU2+ y(t)s2(t)dtT0T0则判为s1(t)出现，是正确判决。若 U1+ dttstyUdttstyTT02021)()()()(则判为s2(t)出现，是错误判决。 将y(t)=s1(t)+n(t)代入式(8.3 - 16)可得U1+ dttstntSU
33、dttstntSTT)()()()()()(2012101代入U1= lnP(s1)和U2= lnP(s2)， 并利用s1(t)和s2(t)能量相等的条件可得20n20ndttstsspspndttstStnTT220110201)()(21)()(ln2)()()(式(8.3 - 18)左边是随机变量，令为，即 = (8.3 - 19) 式(8.3 - 18)右边是常数， 令为a, 即dttststnT)()()(201dttstsspspnaT2021120 )()(21)()(ln2式(8.3 - 18)可简化为 a (8.3 - 21) 判为s2(t)出现，产生错误判决。则发送s1(t
34、)将其错误判决为s2(t)的条件简化为a事件，相应的错误概率为 Ps1(s2)=P(a) 只要求出随机变量的概率密度函数， 即可计算出式(8.3 - 22)的数值。 根据假设条件， n(t)是高斯随机过程， 其均值为零，方差为2n。根据随机过程理论可知，高斯型随机过程的积分是一个高斯型随机变量。所以是一个高斯随机变量，只要求出的数学期望和方差，就可以得到的概率密度函数。 的数学期望为 E=E 0)()()()()( )(201210dttststnEdttstSTnTT的方差为2=D=E2=E 0)()()()( )(212100dttssntstSTnTT 式中En(t)n()为高斯白噪声n
35、(t)的自相关函数，由第 2 章随机信号分析可知 En(t)n()= (t-)= 20n20n(0) 0t=t将上式代入式(8.3 - 24)可得2=20ndttstsT2201)()(于是可以写出的概率密度函数为 f()= 222exp21至此, 可得发送s1(t)将其错误判决为s2(t)的概率为Ps1(s2)=P(a)=dxxb2exp212 利用相同的分析方法，可以得到发送s2(t)将其错误判决为s1(t)的概率为dxxspbs2exp21)(212系统总的误码率为Pe=P(s1)Ps1(s2)+P(s2)Ps2(s1) =P(s1)2exp(21)()2exp(21222dxxspdx
36、xbb式中b和b分别为 TTdttstsnspspdttstsnb022101202210)()(212)()(ln)()(21TTdttstsnspspdttstsnb022101202210)()(212)()(ln)()(21 由式(8.3 - 30)、式(8.3 - 31)和式(8.3 - 32)可以看出， 最佳接收机的误码性能与先验概率P(s1)和P(s2)、噪声功率谱密度n0及s1(t)和s2(t)之差的能量有关，而与s1(t)和s2(t)本身的具体结构无关。 一般情况下先验概率是不容易确定的，通常选择先验等概的假设设计最佳接收机。在发送s1(t)和s2(t)的先验概率相等时， 误
37、码率Pe还与s1(t)和s2(t)之差的能量有关，如何设计s1(t)和s2(t)使误码率Pe达到最小，是我们需要解决的另一个问题。 比较式(8.3 - 31)和式(8.3 - 32)可以看出， 当发送信号先验概率相等时，b=b，此时误码率可表示为Pe= )2(212exp212AerfcdxxA式中:TdttstsnA02210)()(21 为了分析方便，我们定义s1(t)和s2(t)之间的互相关系数为EdttstsT021)()(式中, E是信号s1(t)和s2(t)在0tT 期间的平均能量。当s1(t)和s2(t)具有相等的能量时，有 E=E1=E2=Eb (8.3 - 36)将Eb和代入
38、式(8.3 - 34)可得:A= (8.3 - 37)此时， 式(8.3 - 33)可表示为 Pe=12erfc0)1 (nEb上式即为二进制确知信号最佳接收机误码率的一般表示式。 它与信噪比 及发送信号之间的互相关系数有关。 0nEb2)1 (0nEb 上式即为发送信号先验概率相等时，二进制确知信号最佳接收机所能达到的最小误码率，此时相应的发送信号s1(t)和s2(t)之间的互相关系数=-1。也就是说，当发送二进制信号s1(t)和s2(t)之间的互相关系数=-1时的波形就称为是最佳波形。 当互相关系数=0时， 误码率为 Pe= 210nEerfcb 若互相关系数=1， 则误码率为 Pe=21
39、 若发送信号s1(t)和s2(t)是不等能量信号，如E1=0，E2=Eb，=0，发送信号s1(t)和s2(t)的平均能量为E=Eb2，在这种情况下，误码率表示式(8.3 - 40)变为 Pe= 4210nEerfcb 根据式(8.3 - 39)、 式(8.3 - 40)和式(8.3 - 41)画出的Pe 关系曲线如图 8 - 11 中所示。 在第 5 章数字基带传输系统误码率性能分析中我们知道， 双极性信号的误码率低于单极性信号，其原因之一就是双极性信号之间的互相关系数=-1，而单极性信号之间的互相关系数=0。 0nEb图 8 11 二进制最佳接收机误码率曲线 1.01011021031041
40、0514 10620261014Ebn0/ dBPe 在第 7 章数字频带传输系统误码性能分析中，2PSK信号能使互相关系数=-1，因此2PSK信号是最佳信号波形；2FSK和2ASK信号对应的互相关系数=0, 因此2PSK系统的误码率性能优于2FSK和2ASK系统；2FSK信号是等能量信号， 而2ASK信号是不等能量信号，因此2FSK系统的误码率性能优于2ASK系统。 8.4随相信号的最佳接收机随相信号的最佳接收机 确知信号最佳接收是信号检测中的一种理想情况。 实际中， 由于种种原因， 接收信号的各分量参数或多或少带有随机因素，因而在检测时除了不可避免的噪声会造成判决错误外，信号参量的未知性使
41、检测错误又增加了一个因素。因为这些参量并不携带有关假设的信息，其作用仅仅是妨碍检测的进行。造成随参信号的原因很多，主要有：发射机振荡器频率不稳定，信号在随参信道中传输引起的畸变， 雷达目标信号反射等。 随机相位信号简称随相信号，是一种典型且简单的随参信号，其特点是接收信号的相位具有随机性质，如具有随机相位的2FSK信号和具有随机相位的2ASK信号都属于随相信号。 对于随相信号最佳接收问题的分析，与确知信号最佳接收的分析思路是一致的。但是，由于随相信号具有随机相位，使得问题的分析显得更复杂一些，最佳接收机结构形式也比确知信号最佳接收机复杂。 8.4.1二进制随相信号最佳接收机结构二进制随相信号最
42、佳接收机结构 二进制随相信号具有多种形式，我们以具有随机相位的2FSK信号为例展开分析。设发送的两个随相信号为 s1(t, 1)= Acos(w1t+0)1s2(t,2)= Acos(2t+2),0tT 0 0 t T , 其他 0 t T , 其他 式中，1和2为满足正交条件的两个载波角频率；1和2是每一个信号的随机相位参数，它们的取值在区间0,2上服从均匀分布，即)(1f,210 012 其他 f(2)= ,210012 其他 s1(t, 1)和s2(t, 2)持续时间为(0, T)，且能量相等，即Eb=E1= dttSEdttsTT),(),(202221021假设信道是加性高斯白噪声信
43、道，则接收机输入端合成波为 y(t)= s1(t, 1)+n(t),发送s1(t, 1)时 s2(t, 2)+n(t),发送s2(t, 2)时 在确知信号的最佳接收中，通过似然比准则可以得到最佳接收机的结构。然而在随相信号的最佳接收中，接收机输入端合成波y(t)中除了加性高斯白噪声之外，还有随机相位， 因此不能直接给出似然函数fs1(y)和fs2(y)。此时，可以先求出在给定相位1和2的条件下关于y(t)的条件似然函数fs1(y/1)和fs2(y/2)，即 fs1(y/1)=dttstyT2101),()(01exp)2(1nkn01exp)2(1nkndttstyT2202),()(fs2(
44、y/2)=由概率论知识可得fs1(y/1)=1210120),()(1exp)2(1ddttstynTkn11111)/()(),(1111dyffdyfssdttynnETbkn)(1exp)2(120020010110)cos()(2ddttwtAynT1201100)cos()(2exp2ddttwtAynkT式中knTbdttynnEK)2()(1exp0200为常数。 令随机变量(1)为dttwtwtynAdttwtAynTT)sinsincos(cos)()cos()(2)(1111000110111000110sin(sin)(2coscos)(2tdtwtynAtwtynATT
45、)arctancos(2)sincos(21112121011110XYYXnAYXnA)cos(20110MnA式中: tdtwtyXT011cos)(tdtwtyYT011sin)(21211YXM于是， 式(8.4 - 9)可表示为)2()cos(2exp21001011200MnAKIdMnAkfs1(y) 式中， K为常数， 为零阶修正贝塞尔函数。 同理可得，出现s2(t)时y(t)的似然函数fs2(y)为)2(100MnAI fs2(y)= )2(200MnAKITtdtwtyX022cos)(TtdtwtyY022cos)(22222YXM代入M1和M2的具体表示式可得: 210
46、212101sin)(cos)(TTtdtwtytdtwtyM210222202sin)(cos)(TTtdtwtytdtwtyM 假设发送信号s1(t, 1)和s2(t, 2)的先验概率相等，采用最大似然准则对观察空间样值作出判决，即fs1(y)fs2(y), 判为s1 fs1(y)fs2(y),判为s2 判为s2将式(8.4 - 15)和式(8.4 - 16)代入上式可得:),2()2(200100MnAKIMnAKI判为S2判决式两边约去常数K后有),2()2(200100MnAIMnAI判为S1),2()2(200100MnAIMnAI判为S2 根据零阶修正贝塞尔函数的性质可知，I0(
47、x)是严格单调增加函数， 若函数I0(x2)I0(x1)，则有x2x1。因此，式(8.4 - 26)和式(8.4 - 27)中，根据比较零阶修正贝塞尔函数大小作出判决，可以简化为根据比较零阶修正贝塞尔函数自变量的大小作出判决。 此时判决规则简化为 M1 M2, 判为S1 M1 M2, 判为s2 02nA02nA02nA02nA 判决式两边约去常数并代入M1和M2的具体表示式后有 M1M2, 判为s1 M1M2, 判为s2 即 21201202sin)(cos)(tdtwtytdtwtyTt21202202sin)(cos)(tdtwtytdtwtyTt判为s1, 而21201201sin)(c
48、os)(tdtwtytdtwtyTt21202202sin)(cos)(tdtwtytdtwtyTt 式(8.4 - 32)和式(8.4 - 33)就是对二进制随相信号进行判决的数学关系式，根据以上二式可构成二进制随相信号最佳接收机结构如图 8 - 12 所示。 上述最佳接收机结构形式是相关器结构形式。 可以看出， 二进制随相信号最佳接收机结构比二进制确知信号最佳接收机结构复杂很多，实际中实现也较复杂。 与二进制确知信号最佳接收机分析相类似，可以采用匹配滤波器对二进制随相信号最佳接收机结构进行简化。 由于接收机输入信号s1(t, 1)和s2(t, 2)包含有随机相位1和2，因此无法实现与输入信
49、号s1(t, 1)和s2(t, 2)完全匹配的匹配滤波器。我们可以设计一种匹配滤波器， 它只与输入信号的频率匹配，而不匹配到相位。与输入信号s1(t, 1)频率相匹配的匹配滤波器单位冲激响应为图 8 12 二进制随相信号最佳接收机结构相关器平方器cos1t相关器平方器sin1t相关器平方器cos2t相关器平方器sin2t输入相加器开方器相加器开方器比较器输出M1M2 h1(t)= cos1(T-t), 0tT 当输入y(t)时， 该滤波器的输出为 e1(t)=y(t)*h1(t)= dTwyt)1(cos)(01)(sinsin)() 1(coscos)(101101tTwdwyTwdwytt
50、)(cossin)()(cos)(11212010201tTwdwytTydwyttt式中ttdwydwy01011cos)(sin)(arctan式(8.4 - 35)在t=T时刻的取值为e1(t)=11121201201coscossin)(cos)(MdwydwYTT 可以看出， 滤波器输出信号在t=T时刻的包络与图 8 - 12 所示的二进制随相信号最佳接收机中的参数M1相等。这表明， 采用一个与输入随相信号频率相匹配的匹配滤波器，再级联一个包络检波器，就能得到判决器所需要的参数M1。 同理，选择与输入信号s2(t, 2)的频率相匹配的匹配滤波器的单位冲激响应为 h2(t)= cos2
