信号检测与估计试题——答案(不完整版)
《信号检测与估计》第九章习题解答

T x 2 (t )dt −2 A
0
Tx(t
0
)sin
(ω
0
t
+θ
)dt
+
A
2
T 0
sin
2
(ω
0
t
+θ
)dt
⎟⎞ ⎠
由于 ∫0Tsin 2 (ω0t
+ θ )dt
=
1 2
∫0T(1 − cos 2(ω0t
+ θ ))dt
=
T 2
,得到
( ) ∫ ∫ f
x A,θ
−1
= Fe N0
T x 2 (t )dt 2 A
ω0
ω0
s(t
)
=
⎪ ⎨
A
⎪
⎪ ⎪⎩
A(1
+
cos
ω0t)
− 2mπ < t ≤ 2mπ
ω0 2mπ
<t
≤
(2ωm0 +1)π
ω0
ω0
试证明时延τ
的无偏估计量的方差为
σ
2 τˆ
≥
3 + 4m
(2E / N0 )ω02
。其中 E 为信号能量。
解:略
9.4 接收信号 x(t) = s(t) + n(t),s(t)的到达有时延τ ,求时延τ 的无偏估计量τˆ 的最小方差。其中 n(t)
⎤ ⎥⎦
∫ ∫ ∫ ∫ [ ] = 4
N
2 0
T / 2 ∂s(t −τ )
−T / 2 ∂τ
T /2
E
−T / 2
n(t )n(u )
∂s(u −τ ) dudt
信号检测与估计试题——答案(不完整版)

一、概念:1. 匹配滤波器。
概念:所谓匹配滤波器是指输出判决时刻信噪比最大的最佳线性滤波器。
应用:在数字信号检测和雷达信号的检测中具有特别重要的意义。
在输出信噪比最大准则下设计一个线性滤波器是具有实际意义的。
2. 卡尔曼滤波工作原理及其基本公式(百度百科)首先,我们先要引入一个离散控制过程的系统。
该系统可用一个线性随机微分方程(Linear Stochastic Difference equation)来描述:X(k)=A X(k-1)+B U(k)+W(k)再加上系统的测量值:Z(k)=H X(k)+V(k)上两式子中,X(k)是k时刻的系统状态,U(k)是k时刻对系统的控制量。
A和B是系统参数,对于多模型系统,他们为矩阵。
Z(k)是k时刻的测量值,H是测量系统的参数,对于多测量系统,H为矩阵。
W(k)和V(k)分别表示过程和测量的噪声。
他们被假设成高斯白噪声(White Gaussian Noise),他们的covariance 分别是Q,R(这里我们假设他们不随系统状态变化而变化)。
对于满足上面的条件(线性随机微分系统,过程和测量都是高斯白噪声),卡尔曼滤波器是最优的信息处理器。
下面我们来用他们结合他们的covariances 来估算系统的最优化输出(类似上一节那个温度的例子)。
首先我们要利用系统的过程模型,来预测下一状态的系统。
假设现在的系统状态是k,根据系统的模型,可以基于系统的上一状态而预测出现在状态:X(k|k-1)=A X(k-1|k-1)+B U(k) (1)式(1)中,X(k|k-1)是利用上一状态预测的结果,X(k-1|k-1)是上一状态最优的结果,U(k)为现在状态的控制量,如果没有控制量,它可以为0。
到现在为止,我们的系统结果已经更新了,可是,对应于X(k|k-1)的covariance还没更新。
我们用P表示covariance:P(k|k-1)=A P(k-1|k-1) A’+Q (2)式(2)中,P(k|k-1)是X(k|k-1)对应的covariance,P(k-1|k-1)是X(k-1|k-1)对应的covariance,A’表示A的转置矩阵,Q是系统过程的covariance。
《信号检测与估计》第十章习题解答

《信号检测与估计》第十章习题解答10.1 设线性滤波器的输入信号为()()()t n t s t x +=,其中()[]0E =t s ,()[]0E =t n ,并且已知()ττ-e =S R ,()ττ-2e=N R ,()0=τsn R ,求因果连续维纳滤波器的传递函数。
解:连续维纳滤波器与离散维纳滤波器的形式是相同的,即()()()()+⎦⎤⎢⎣⎡−⋅⋅=s B s P s B s H xs w112opt σ 因此需要求解()t s 的复功率谱和()t x 的时间信号模型。
考虑到信号与噪声不相关,因此观测数据的功率谱就等于信号的复功率谱加上噪声的复功率谱。
对观测数据的复功率谱进行谱分解,就可以得到()t x 的时间信号模型。
()t s 的复功率谱为()()()20s -10s 1-s --121111e e e e s s s d d d s P S −=−++=+==∫∫∫∞−+∞++∞∞−τττττττ ()t n 的复功率谱为()2s -2-44e es d s P N −==∫+∞∞−τττ因此,观测数据的复功率谱为()()()()()()()()()s s s s ss ss s P s P s P N S X −+−++=−+−=+=2211-226441122 取12=w σ()()()()s s ss B +++=2126()()()()()()()()()s s s s s s s s B s P s B s P N xs +=−==1-2-262-2-1-2612--2令()()()s B s P s F xs -=,()τf 是()s F 的拉普拉斯反变换。
要求()τf 是因果的,可将s 平面右半平面的极点扔掉,()()()[]12e 61,e Re e21-s s +=−==∫τττπτs F s ds s F jf C给()τf 取因果,并做拉普拉斯变换,得到()s d s F +⋅+=⋅⋅+=∫∞++11126e e 1260s --τττ()()()()()()())()()122261112626211112opt +++=+×+×+++×=⎦⎤⎢⎣⎡−⋅⋅=+ss ss s s s B s P s B s H xs wσ10.2 设已知()()()n n n s n x +=,以及()()()z z z G S 4.014.0192.01−−=−,()1=z G N ,()0=z G sn ()n s 和()n n 不相关。
信号检测与估计答案1

信号检测与估计答案15-2 若观测方程为i i x s n =+()1,2,,i N =,其中信号()2~0,s s N σ,噪声()()2~0,1,2,,i n n N i N σ=独立同分布,且信号与噪声满足{}0i E sn =。
求s 的最大后验概率估计ˆMAP s。
解:依题意,以信号s 为条件的观测样本的概率密度函数为()()()2112221,,|exp 22N i i N Nnnx s f x x s σπσ=⎡⎤-⎢⎥⎢⎥=-⎢⎥⎢⎥⎣⎦∑信号s 的概率密度函数为()222ss f s σ⎛⎫=- ⎪⎝⎭则由上面两式可得()()()()()211222212221ln ,,|ln exp 221ln 22N i i N N nn Ni i N n n x s f x x s ss x s s σπσσπσ==⎧⎫⎡⎤⎧⎫-⎪⎪⎢⎥⎪⎪∂∂⎪⎪⎪⎪⎢⎥=-⎨⎨⎬⎬∂∂⎢⎥⎪⎪⎪⎪⎢⎥⎪⎪⎪⎪⎩⎭⎣⎦⎩⎭⎡⎤-⎢⎥∂⎢⎥=-∂⎢⎥⎢⎥⎣⎦∑∑()22222ln ln 22s s s s f s s s s s s⎧⎫⎡⎤⎛⎫∂∂⎪⎪=-⎥⎨⎬ ⎪∂∂⎥⎝⎭⎪⎪⎦⎩⎭⎡⎤∂=-⎢⎥∂⎢⎥⎣⎦=-σσσ最大后验概率准则为()ˆmax |MAP f θθθ=x ,即()ˆ|0MAPf θθθθ=∂⎡⎤=⎢⎥∂⎣⎦x ,又可表示为()()ˆln |ln 0MAPf f θθθθθθ=∂∂⎡⎤+=⎢⎥∂∂⎣⎦x ,将之前结果带入其中可得2221ˆNs MAP ii ns sx N σσσ==+∑ 。
5-4已知观测信号0()cos()()x t A t n t ωθ=++(0)t T ≤≤,式子中()n t 是零均值,功率谱为2N 的高斯白噪声,θ是在[0,2)π上均匀分布的随机变量,求A 的最大似然估计和估计量的均方误差。
解:0()cos()()x t A t n t ωθ=++()x t 的似然函数为:020002220000022000000()cos()()1(|,)exp [()cos()]1exp [()2()cos()cos ()]12exp [()()cos()2TTTTTT x t A t n t f x A F x t A t dt N F x t dt x t A t dt At dt N A A T F x t dt x t t dt N N N ωθθωθωθωθωθ=++⎧⎫=⋅--+⎨⎬⎩⎭⎧⎫=⋅--+++⎨⎬⎩⎭⎧⎫=⋅-++-⎨⎬⎩⎭⎰⎰⎰⎰⎰⎰因为1(),022f θθππ=≤≤ 所以202200000(|)(|,)()12exp{}exp [()()2Tf x A f x A f d A TAq F x t dt I N N N πθθθ=⎧⎫=⋅--⎨⎬⎩⎭⎰⎰ 其中22200002200000()sin ()cos 12ln (|)ln ()ln ()2T TT q x t tdt x t tdt A T Aqf x A F x t dt I N N N ωω⎡⎤⎡⎤=+⎢⎥⎢⎥⎣⎦⎣⎦=--+⎰⎰⎰令000ln (|)20()0f x A AT AqI A N A N ∂∂=⇒-++∂∂ (1)假设SNR,即02Aq N 足够大,则00022()Aq AqI N N ≈0022ˆ(1)0MLAT q q A N N T⇒-+=⇒=由2220000()sin ()cos T Tq x t tdt x t tdt ωω⎡⎤⎡⎤=+⎢⎥⎢⎥⎣⎦⎣⎦⎰⎰知22202221()exp(())()242T T TqA T qATf q q I σσσ=-+所以222323240001()2T T qq x x q T q E qf q dq AT e dq e AT x e dx AT σσσ=-+∞+∞+∞-⎛⎫−−−→ ⎪==←−−− ⎪⎝⎭⎰⎰⎰ 所以221ˆ()()2MLE A E q AT A T T ==⋅= (无偏估计) 200024ˆvar(),var()44T ML N T N T N q A T T σ====5-11. 假定已知信号112()cos cos 2...cos p s t a t a t a p t ωωω=+++212()sin sin 2...sin p s t b t b t b p tωωω=+++观测信号12()()()()x t s t s t n t =++,()n t 是均值为0、均方差为1的高斯白噪声。
信号检测与估计理论(复习题解)

H1)
s2 1k
s1k s0k
k 1
k 1
第3章 信号状态的统计检测理论 例题解答
Var(l |1H1) Var(l | H
)
E
N
k 1
nk
s1k
N k 1
nk
s0k
2
N
2
2
s n
1k
k1
N
s2 0k k 1
信号检测与估计理论
内容提要 例题解答
第1章 信号检测与估计概论 信号的随机性及其统计处理方法。
内容提要
第1章 信号检测与估计概论
略
例题解答
第2章 信号检测与估计理论的基础知识 内容提要
一. 离散随机信号
1. 概率密度函数p(x)及特性: 非负,全域积分等于1,落入[a,b]间的概率。
2. 统计平均量:均值,方差。
解:似然函数为
p(x
|
H0
)
1
2
2 n
N
2
exp
N
k 1
( xk
s0k
2
2 n
)2
p(x
|
H1)
1
2
2 n
N
2
exp
N
k 1
(xk s1k
2
2 n
)2
第3章 信号状态的统计检测理论 例题解答
其中,观测噪声n服从对称三角分布,如图3.1(a)所示。
若似然比检测门限 1,求最佳判决式,图示判决域,计算P(H1 | H0 )。
信号检测与估计理论(复习题解)-精选文档

a ba 0 图 2. 1 (b)
ab y
2 b y x
2 2 y 4 x
第2章 信号检测与估计理论的基础知识 例题解答
例 2 . 3 设连续随机信号 x ( t ) a cos( t ), 其振幅 a 和频率 已知 相位 在 [ , ) 范围内均匀分布。分析 该信号的广义平稳 并求其自 差函数 。 解 : 分析该信号是否满足广 义平稳的条件。 信号的均值 ( t ) E a cos( t ) a cos( t ) p ( ) d x
2 1 ( y b ) / 2 1 x p ( y ) exp 2 2 2 2 2 x x 1 2
2 1 y ( 2 b ) x exp 2 2 8 8 x x 1 2
二. 离散随机信号矢量
1. 概率密度函数描述 。 2. 统计平均量:均值矢量 , 协方差, 协方差矩阵。 3. 各分量之间的互不相关 性和相互统计独立性及 关系。 4. 高斯离散随机信号矢量 的概率密度函数及特 点: x ~ N ( μ , C ), 互不相关等价于相互统 计独立 , 独立同分布 x x
E ( x b ) b
y
2 y
2 2 22 E ( y b ) E ( x b b ) E ( x 0 ) a / 6
第2章 信号检测与估计理论的基础知识 例题解答
当 a b 2 a 时, p ( y ) 的函数曲线如图 2 . 1 (b)所示 。 p ( x) p( y ) 1/ a 1/ a
第 1章
信号检测与估计概论
信号检测与估计试题及答案

(1). 若 为常数,求 的最大似然估计。
ˆ 1 N ln xiБайду номын сангаасN i 1
(2). 判断 的最大似然估计是否是有效估计? 因为
ˆ HX B ,其中 H C M N , B C M 1
(1). 用最小均方误差准则确定矩阵 H , B 。 (用 , x 的一阶和二阶统计量表 示。 )
H cov( , x ) cov1 ( x , x ) B E ( ) cov( , x ) cov 1 ( x , x ) E ( x)
2 ) ,做 H1 判决,反之做 H 0 判决。 ln 2 3
2
4. 求解下列问题 (1). 什么是序贯检测?
A1 , D1 ( x) A0 , D0 other , more obervation
(2). 对二元检测 P D1 H 0 , P D0 H1 若,推导瓦尔特序贯检测的门
1 (2). 若是线性调频信号, 即 s1 (t ) A1 cos(1t t 2 ) 0 t T , 2 / 1 T , 2
是常数,再求 Pe 结果相同。
3. 设有两种假设分别为:
H 0 : P0 ( x)
x2 1 exp 2 2 2 1 x A, A 0 H1 : P 1 ( x) 2 A 0 x >A
(2). ˆ 是否无偏
是无偏估计。
7. 求解下列问题。 (1). 什么是卡尔曼滤波,写出卡尔曼滤波的状态方程,观测方程和滤波方程
《信号检测与估计》第十二章习题解答

《信号检测与估计》第十二章习题解答12.1 采用下式给出的有偏自相关函数的定义,并加窗,得到BT 谱估计器:()()()()()()⎪⎩⎪⎨⎧−−−−−=−+=+=∑∗1,,2,11ˆ1,,1,01ˆL L N N m m R N m m n x n x N m R X X ()⎪⎩⎪⎨⎧−≤=其它011N m m W N()()()()∑−−−=−⋅⋅=11e ˆˆN N m m j X N X m R m W G ωω证明该BT 估计器与周期图相同。
解:()()()()()()()()()()()()()()()()()211111111e 1e e 1e e 1e 1e ˆˆωωωωωωωωj N N m n m j nj N N m nj n m j N N m m j N N N m m j X N XX N m n x n x N m n x n x N m n x n x N m W m R m W G =⋅+⋅⋅=⋅⋅+=⋅⎥⎦⎤⎢⎣⎡+⋅=⋅⋅=∑∑∑∑∑∑∑−−−=+−−∗−−−=−+−∗−−−=−∗−−−=− 12.2 设自相关函数()3,2,1,0,==m m R m X ρ。
试用Levinson-Durbin 递推法求解AR (3)模型参量。
解: ()()ρ−=−=0111X X R R a 110=a()()221121101ρσ−=⋅−=X R a ()()012211122=+−=σX X R a R a ρ−=⋅+=11221121a a a a ()2212222211ρσσ−=⋅−=a因此模型为一阶 ()()[]()012322222133=⋅+−=σX X X R a R a R a021332232=⋅+=a a a aρ−=⋅+=22332131a a a a()2222332311ρσσ−=⋅−=a 所以模型为()()()n w n x n x +−=1ρ12.3 设5=N 的数据记录为:10=x ,21=x ,32=x ,43=x ,54=x ,AR 模型的阶数3=p 试用Levinson-Durbin 递推法求模型参量。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
一、概念:1. 匹配滤波器。
概念:所谓匹配滤波器是指输出判决时刻信噪比最大的最佳线性滤波器。
应用:在数字信号检测和雷达信号的检测中具有特别重要的意义。
在输出信噪比最大准则下设计一个线性滤波器是具有实际意义的。
2. 卡尔曼滤波工作原理及其基本公式(百度百科)首先,我们先要引入一个离散控制过程的系统。
该系统可用一个线性随机微分方程(Linear Stochastic Difference equation)来描述:X(k)=A X(k-1)+B U(k)+W(k)再加上系统的测量值:Z(k)=H X(k)+V(k)上两式子中,X(k)是k时刻的系统状态,U(k)是k时刻对系统的控制量。
A和B是系统参数,对于多模型系统,他们为矩阵。
Z(k)是k时刻的测量值,H是测量系统的参数,对于多测量系统,H为矩阵。
W(k)和V(k)分别表示过程和测量的噪声。
他们被假设成高斯白噪声(White Gaussian Noise),他们的covariance 分别是Q,R(这里我们假设他们不随系统状态变化而变化)。
对于满足上面的条件(线性随机微分系统,过程和测量都是高斯白噪声),卡尔曼滤波器是最优的信息处理器。
下面我们来用他们结合他们的covariances 来估算系统的最优化输出(类似上一节那个温度的例子)。
首先我们要利用系统的过程模型,来预测下一状态的系统。
假设现在的系统状态是k,根据系统的模型,可以基于系统的上一状态而预测出现在状态:X(k|k-1)=A X(k-1|k-1)+B U(k) (1)式(1)中,X(k|k-1)是利用上一状态预测的结果,X(k-1|k-1)是上一状态最优的结果,U(k)为现在状态的控制量,如果没有控制量,它可以为0。
到现在为止,我们的系统结果已经更新了,可是,对应于X(k|k-1)的covariance还没更新。
我们用P表示covariance:P(k|k-1)=A P(k-1|k-1) A’+Q (2)式(2)中,P(k|k-1)是X(k|k-1)对应的covariance,P(k-1|k-1)是X(k-1|k-1)对应的covariance,A’表示A的转置矩阵,Q是系统过程的covariance。
式子1,2就是卡尔曼滤波器5个公式当中的前两个,也就是对系统的预测。
现在我们有了现在状态的预测结果,然后我们再收集现在状态的测量值。
结合预测值和测量值,我们可以得到现在状态(k)的最优化估算值X(k|k):X(k|k)= X(k|k-1)+Kg(k) (Z(k)-H X(k|k-1)) (3)其中Kg为卡尔曼增益(Kalman Gain):Kg(k)= P(k|k-1) H’ / (H P(k|k-1) H’ + R) (4)到现在为止,我们已经得到了k 状态下最优的估算值X(k|k)。
但是为了要令卡尔曼滤波器不断的运行下去直到系统过程结束,我们还要更新k 状态下X(k|k)的covariance :P(k|k)=(I-Kg(k) H )P(k|k-1) (5)其中I 为1的矩阵,对于单模型单测量,I=1。
当系统进入k+1状态时,P(k|k)就是式子(2)的P(k-1|k-1)。
这样,算法就可以自回归的运算下去。
卡尔曼滤波器的原理基本描述了,式子1,2,3,4和5就是他的5 个基本公式。
根据这5个公式,可以很容易用计算机编程实现。
3. 白噪声的概念白噪声定义:将噪声用n(t)表示,功率谱密度 为常数,具有这种特性的噪声称为白噪声。
说明:这种称呼来源于光学。
因为光学中将包括全部可见波长的光称为白光,所以我们也将包括了全部频率成分的噪声称为白噪声。
说明:实际上,完全理想的白噪声是不存在的,但只要噪声功率谱均匀分布的范围超过电子系统工作的频率范围很多时,就可以近似认为是白噪声。
例如,热噪声功率谱密度均匀分布的部分高达1013Hz ,因此可将它看成白噪声。
4. 二元信号检测理论模型 第三章ppt 第4页开始5. 谱、功率谱、频谱的联系与区别。
区别:1、 一个信号的频谱,只是这个信号从时域表示转变为频域表示,只是同一种信号的不同的表示方式而已, 而功率谱是从能量的观点对信号进行的研究,其实频谱和功率谱的关系归根揭底还是信号和功率,能量等之间的关系。
2、 频谱是个很不严格的东西,常常指信号的Fourier 变换,是一个时间平均(time average )概念;功率谱的概念是针对功率有限信号的(能量有限信号可用能量谱分析),所表现的是单位频带内信号功率随频率的变换情况。
保留频谱的幅度信息,但是丢掉了相位信息,所以频谱不同的信号其功率谱是可能相同的。
3、功率谱是随机过程的统计平均概念,平稳随机过程的功率谱是一个确定函数;而频谱是随机过程样本的Fourier 变换,对于一个随机过程而言,频谱也是一个2)(0nP n =ω“随机过程”。
(随机的频域序列)4、功率概念和幅度概念的差别。
此外,只能对宽平稳的各态历经的二阶矩过程谈功率谱,其存在性取决于二阶局是否存在并且二阶矩的Fourier变换收敛;而频谱的存在性仅仅取决于该随机过程的该样本的Fourier变换是否收敛。
联系:1、功率谱可以从两方面来定义,一个是自相关函数的傅立叶变换,另一个是时域信号傅氏变换模平方然后除以时间长度。
第一种定义就是常说的维纳辛钦定理,而第二种其实从能量谱密度来的。
根据parseval定理,信号傅氏变换模平方被定义为能量谱,能量谱密度在时间上平均就得到了功率谱。
2、在频域分析信号分两种:(1).对确定性信号进行傅里叶变换,分析频谱信息。
(2).随机信号的傅里叶信号不存在,转向研究它的功率谱。
随机信号的功率谱和自相关函数是傅里叶变换对(即维纳辛钦定理)。
功率谱估计有很多种方法6. 二元通信系统的检测性能与那些因素有关。
(可能不准确)概率转移机构,观测空间R,最佳判决准则第三章ppt第4页开始二、计算:1. 在二元数字通信系统中,假设为H1时,信源输出为正电压A,假设为H0时,信源输出为零电平。
信号在通信信道传输过程中叠加了高斯噪声n(t);每种信源的持续时间为T,在接收端对接收到的信号x(t) 在T时间内进行N次独立采样,样本为x k(k=1,2,…,N)。
已知噪声样本n k是均值为零、方差为σn2的高斯噪声。
(1)试建立信号检测系统的信号模型;(2)若似然检测门限已知,确定似然比检验的判决表达式;(3)计算判决概率P(H1|H0)和P(H1|H1)。
解:(1) 接收信号模型为:01()(),0()(),0H x t n t t T H x t A n t t T=≤≤=+≤≤::在(0,T )内进行N 次独立采样后,接收信号模型为:01,1,2,...,,1,2,...,k k k k H x n k N H x A n k N===+=:: 其中 x k 之间相互独立。
(2)已知22()2k k n n p n σ⎛⎫=- ⎪⎝⎭在两种假设情况下,似然函数为:()202212(|)2(|)2k k n kk n x p x H x A p x H σσ⎛⎫=- ⎪⎝⎭⎛⎫- ⎪=- ⎪⎝⎭由于N 次采样的样本 x k 之间是独立同分布(iid)的,所以()210021211121(|)(|)exp()2(|)(|)exp()2NNkNk k k n NNkNk k k nxp x H p x H xA p x H p x H σσ====⎛⎫==--⎛⎫==-∑∏∑∏这样,似然比函数为()221112202221(|)()exp()(|)22exp()2λσσσσ===--==+=-∑∑∑NNk kk k nnNk k n nx A xP x H x P x H ANA x似然比函数检验(LRT )为102221exp(2H Nk k nn H ANA x ησσ=≥-<∑) 取对数:12221ln 2H Nk k nn H ANA x ησσ=≥-<∑进一步整理得1021ln 12H Nn kk H A x NNA ση=≥+<∑ (3) 检验统计量11Nkk xN=∑是N 个信号的平均值,它是x k (k=1,2,…,N)的函数,是个随机变量。
121ln 1(),()2H Nnk k H A l x x l x NNAσηγγ=≥==+∴<∑令,说明:由于N 次采样的样本x k 之间是独立同分布(iid)的,因此 l(x) 在两种假设情况下均服从高斯分布,均值和方差计算过程如下。
假设H 0情况下,均值和方差分别为:[][][]00112220001111()||011()||()|N N k k k k N Nn k k k k E l x H E x H E n N N D l x H E x H E l x H E n N NN σ====⎡⎤⎡⎤===⎢⎥⎢⎥⎣⎦⎣⎦⎡⎤⎡⎤⎛⎫⎛⎫=-==⎢⎥⎢⎥ ⎪⎪⎝⎭⎝⎭⎢⎥⎢⎥⎣⎦⎣⎦∑∑∑∑假设H 1情况下,同样的方法计算均值和方差为:[]()[][]()11112221111111()||11()||()|N N k k k k N Nn k k k k E l x H E x H E A n AN N D l x H E x H E l x H E A n A N N Nσ====⎡⎤⎡⎤==+=⎢⎥⎢⎥⎣⎦⎣⎦⎡⎤⎡⎤⎛⎫⎛⎫=-=+-=⎢⎥⎢⎥ ⎪⎪⎝⎭⎝⎭⎢⎥⎢⎥⎣⎦⎣⎦∑∑∑∑用l 表示l(x),有()()20211|~0,1|~,n n l H N N l H N A N σσ⎧⎛⎫⎪⎪⎪⎝⎭⎨⎛⎫⎪ ⎪⎪⎝⎭⎩根据判决准则,()10022|(|)2n P H H p l H dlNl dl γγσ∞∞=⎛⎫=- ⎪⎝⎭⎰⎰2ln :2n Awhere NAσηγ=+()()11122|(|)2n P H H p l H dlN l A dl γγσ∞∞=⎛⎫-=- ⎪ ⎪⎝⎭⎰⎰2. 例题3.4-2:在OOK 通信系统中,两个假设下的观测信号模型为nA x H n x H +==:10:其中,观测噪声n~N (0,σn2);信号A 是常数,且A>0。
若两个假设的先验概率P(Hj)未知,代价因子C00=C11=0,C10=C01=1,采用极小化极大准则,试确定检测门限和平均错误概率。
解:在两个假设下,观测量x 的概率密度函数分别为()202212(|)2(|)2n n x p x H x A p x H σσ⎛⎫=- ⎪⎝⎭⎛⎫-=- ⎪ ⎪⎝⎭似然比函数为: ()()21220|2(exp |22n n P x H Ax A x P x H λσσ⎡⎤∴==-⎢⎥⎣⎦)假设判决门限为η,则: 12222exp 22H n n HAx A ησσ⎡⎤≥-⎢⎥<⎣⎦ 化简得: 12ln 2H n H AxAσηγ≥+=< 显然,检验统计量l(x)=x .X 在两种假设情况均服从高斯分布,根据判决准则,有()100222|(|)22F n n P P H H p l H dll dl u du Q γγσγσ∞∞∞==⎛⎫=- ⎪⎝⎭⎛⎫=- ⎪⎝⎭⎛⎫= ⎪⎝⎭⎰⎰()()01122|(|)21M n n P P H H p l H dll A dl A Q γγσγσ-∞==⎛⎫- ⎪=- ⎪⎝⎭⎛⎫-=- ⎪⎝⎭⎰⎰因为C 00=C 11=0,C 10=C 01=1, 此时极小化极大方程 为:P F =P M ,即极小化极大方程为**1n n A Q Q γγσσ⎛⎫⎛⎫-=- ⎪ ⎪⎝⎭⎝⎭解得*2Aγ=此时平均错误概率**2Fn n A Pe P r Q Q γσσ⎡⎤⎡⎤===⎢⎥⎢⎥⎣⎦⎣⎦() 式中22nASNR σ=3. 例题3.4-3:在二元数字通信系统中,假设为H1时信源输出为1,假设为H0时信源输出为0,信号在通信信道上传输时叠加了均值为零、方差为1的高斯噪声。