自-信号检测与估计试题——答案(不完整版)
- 格式:doc
- 大小:271.50 KB
- 文档页数:11
《信号检测与估计》第十章习题解答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 不相关。
《信号检测与估计》第十二章习题解答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 递推法求模型参量。
时间:6月16日(星期一)晚上6:30-8:30 地点:六教104室(上课教室)试卷共8题,其中4题属于教材第一章内容,其余4题分别的其他章节。
请同学们对匹配滤波器,离散卡尔曼滤波,离散维纳滤波,高斯白噪声下确知信号的检测,K -L 展开,高斯白噪声信道中的单参量信号估计等内容重点关注。
1.1 (付柏成 20060150)在例1.2中,设噪声均方差电压值为σ=2v ,代价为f c =2,m c =1。
信号存在的先验概率P =0.2。
试确定贝叶斯意义下最佳门限β,并计算出相应的平均风险。
解:根据式(1-15),可以算出00.8280.21f mQc Pc ⨯Λ===⨯ 而判决门限2201ln 0.52ln88.822βσ=+Λ=+= 根据式(1-21)可知平均风险1010Pr 0.2r 0.8R Qr r =+=+01100.2(|)0.8(|)m f c P D H c P D H =+ 而011(|)(|)D P D H p x H dx =⎰1100(|)(|)D P D H p x Hdx =⎰而212(1)(|)]2x p x H σ-=-202(|)]2x p x H σ=-所以20112(1)(|)(|)]2D D x P D H p x H dx dx σ-==-⎰⎰22(1)]2x dx βσ-=-⎰=17.82()()(3.91)22β-Φ=Φ=Φ 同理1121002(|)(|)]2D D x P D H p x H dx dx σ==-⎰⎰22)2x dx βσ∞=- 8.821()1()1(4.41)22β=-Φ=-Φ=-Φ 所以0.21(3.91)0.82[1(4.41)]R =⨯⨯Φ+⨯⨯-Φ 1.2 (关瑞东 20060155)假定加性噪声()n t 服从均值为零,方差为的正态分布。
此时,两个假设为01:()():()1()H x t n t H x t n t ==+我们根据()x t 的两次独立测量值12,x x 作判断,则12,x x 是统计独立的,在假设1H 下其均值为1a =1,在假设0H 下均值为0a =0,因而在两种假设下它们的联合概率密度函数可写为22/221()(|)(2)exp()2nn i k k i x a p x H πσσ-=-=-∑ (0,1;2)k n == 于是,似然比等于22011012210()(|)()exp[](|)2n ii a a n a a p x H x x p x H σσ=--Λ==-∑如果0()x Λ≥Λ,则选择假设1H ,否则选择假设0H 。
时间:6月16日(星期一)晚上6:30-8:30 地点:六教104室(上课教室)试卷共8题,其中4题属于教材第一章内容,其余4题分别的其他章节。
请同学们对匹配滤波器,离散卡尔曼滤波,离散维纳滤波,高斯白噪声下确知信号的检测,K -L 展开,高斯白噪声信道中的单参量信号估计等内容重点关注。
1.1 (付柏成 20060150)在例1.2中,设噪声均方差电压值为σ=2v ,代价为f c =2,m c =1。
信号存在的先验概率P =0.2。
试确定贝叶斯意义下最佳门限β,并计算出相应的平均风险。
解:根据式(1-15),可以算出00.8280.21f mQc Pc ⨯Λ===⨯ 而判决门限2201ln 0.52ln88.822βσ=+Λ=+= 根据式(1-21)可知平均风险1010Pr 0.2r 0.8R Qr r =+=+01100.2(|)0.8(|)m f c P D H c P D H =+ 而011(|)(|)D P D H p x H dx =⎰1100(|)(|)D P D H p x Hdx =⎰而212(1)(|)]2x p x H σ-=-202(|)]2x p x H σ=-所以20112(1)(|)(|)]2D D x P D H p x H dx dx σ-==-⎰⎰22(1)]2x dx βσ-=-⎰=17.82()()(3.91)22β-Φ=Φ=Φ 同理1121002(|)(|)]2D D x P D H p x H dx dx σ==-⎰⎰22)2x dx βσ∞=- 8.821()1()1(4.41)22β=-Φ=-Φ=-Φ 所以0.21(3.91)0.82[1(4.41)]R =⨯⨯Φ+⨯⨯-Φ 1.2 (关瑞东 20060155)假定加性噪声()n t 服从均值为零,方差为的正态分布。
此时,两个假设为01:()():()1()H x t n t H x t n t ==+我们根据()x t 的两次独立测量值12,x x 作判断,则12,x x 是统计独立的,在假设1H 下其均值为1a =1,在假设0H 下均值为0a =0,因而在两种假设下它们的联合概率密度函数可写为22/221()(|)(2)exp()2nn i k k i x a p x H πσσ-=-=-∑ (0,1;2)k n == 于是,似然比等于22011012210()(|)()exp[](|)2n ii a a n a a p x H x x p x H σσ=--Λ==-∑如果0()x Λ≥Λ,则选择假设1H ,否则选择假设0H 。
2011《信号检测与估计》复习纲要“信号检测与估计”理论是现代信息科学的一个重要组成部分,它是把所要处理的问题,归纳为一定的“数学模型”→运用“概率论”、“随机过程”、“数理统计”等数学工具→以普遍化的形式提出,以寻求普遍化的答案和结论,并且理论与工程实践相结合,以雷达系统、通信系统、声纳系统为主要研究对象,主要内容包括:● 随机信号与噪声理论(The Theory of Random Signals and Noise)——分析随机信号与噪声的数学工具● 统计判决(检测)理论(Statistical Decision Theory)——研究在噪声干扰背景中,所关心的信号是属于哪种状态的最佳判决问题(Detection of Signals in Noise)● 参量估计理论(Estimation Theory of Signal Parameters)——研究在噪声干扰背景中,通过对信号的观测,如何构造待估计参数的最佳估计量问题(Estimation of Signal Parameters)● 滤波理论(Filtering Theory)——为了改善信号质量,研究在噪声干扰中所感兴趣信号波形的最佳恢复问题,或离散状态下表征信号在各离散时刻状态的最佳动态估计问题(Estimation of Signal Waveform) 复习重点:信号检测与参量估计 ● 信号检测:根据有限观测,“最佳”区分一个物理系统不同状态的理论 ● 参量估计:根据有限观测,“最佳”找出一个物理系统不同参数的理论如何选择一个估计量&估计量选择的决策过程信号处理否估计量LSE经典方法贝叶斯方法如何选择一个检测器-二元信号检测如何选择一个检测器-多元信号检测*注:ARMA:自回归滑动平均BLUE:最佳线性无偏估计CFAR:恒虚警率CRLB :Cramer-Rao下限EM:数学期望最大化GLRT:广义似然比检验IID:独立同分布LLR:对数似然比LMMSE:线性最小均方误差LMP:局部最大势LRT:似然比检验LSE:最小二乘估计LSI:线性时不变MAP:最大后验概率MLE:最大似然估计MMSE:最小均方误差估计MVU:最小方差无偏NP:Neyman-Pearson准则PRN:伪随机噪声RBLS:Rao-Blackwell-Lehmann-Scheffe定理ROC:接收机工作特性UMP:一致最大势WGN:白色高斯噪声WSS:广义平稳2011《信号检测与估计》复习参考题参数估计部分:1.基本概念理解:最小方差无偏估计,最佳线性无偏估计,最大似然估计,最小二乘估计,矩方法估计,最小均方误差估计,最大似然估计,线性最小均方误差估计,一般(经典)线性模型和贝叶斯线性模型。
一、概念:1. 匹配滤波器。
概念:所谓匹配滤波器是指输出判决时刻信噪比最大的最佳线性滤波器。
应用:在数字信号检测和雷达信号的检测中具有特别重要的意义。
在输出信噪比最大准则下设计一个线性滤波器是具有实际意义的。
2. 卡尔曼滤波工作原理及其基本公式(百度百科)首先,我们先要引入一个离散控制过程的系统。
该系统可用一个线性随机微分方程(Linear StochasticDifference equation)来描述:X(k)=A X(k-1)+BU(k)+W(k)再加上系统的测量值:Z(k)=HX(k)+V(k)上两式子中,X(k)是k时刻的系统状态,U(k)是k时刻对系统的控制量。
A和B是系统参数,对于多模型系统,他们为矩阵。
Z(k)是k时刻的测量值,H是测量系统的参数,对于多测量系统,H为矩阵。
W(k)和V(k)分别表示过程和测量的噪声。
他们被假设成高斯白噪声(White GaussianNoise),他们的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)的covarian ce :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)表示,功率谱密度 为常数,具有这种特性的噪声称为白噪声。
说明:这种称呼来源于光学。
因为光学中将包括全部可见波长的光称为白光,所以我们也将包括了全部频率成分的噪声称为白噪声。
说明:实际上,完全理想的白噪声是不存在的,但只要噪声功率谱均匀分布的范围超过电子系统工作的频率范围很多时,就可以近似认为是白噪声。
例如,热噪声功率谱密度均匀分布的部分高达1013H z,因此可将它看成白噪声。
4. 二元信号检测理论模型 第三章ppt 第4页开始5. 谱、功率谱、频谱的联系与区别。
区别:1、 一个信号的频谱,只是这个信号从时域表示转变为频域表示,只是同一种信号的不同的表示方式而已, 而功率谱是从能量的观点对信号进行的研究,其实频谱和功率谱的关系归根揭底还是信号和功率,能量等之间的关系。
2、 频谱是个很不严格的东西,常常指信号的Fourier 变换,是一个时间平均(tim e averag e)概念;功率谱的概念是针对功率有限信号的(能量有限信号可用能量谱分析),所表现的是单位频带内信号功率随频率的变换情况。
保留频谱的幅度信息,但是丢掉了相位信息,所以频谱不同的信号其功率谱是可能相同的。
3、功率谱是随机过程的统计平均概念,平稳随机过程的功率谱是一个确定函数;而频谱是随机过程样本的Fo uri er变换,对于一个随机过程而言,频谱也是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)。
已知噪声样本nk是均值为零、方差为σ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()2NNkN k k k n NN k Nk k k n xp x H p x H x A 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 ησσ=≥-<∑进一步整理得121ln 12H Nn kk H A x NNA ση=≥+<∑ (3) 检验统计量11Nkk xN=∑是N 个信号的平均值,它是xk (k=1,2,…,N)的函数,是个随机变量。
121ln 1(),()2H Nnk k H A l x x l x NNAσηγγ=≥==+∴<∑令,说明:由于N 次采样的样本x k 之间是独立同分布(iid)的,因此 l(x ) 在两种假设情况下均服从高斯分布,均值和方差计算过程如下。
假设H0情况下,均值和方差分别为:[][][]00112220001111()||011()||()|N N k k k k N N nk 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 N Nσ====⎡⎤⎡⎤===⎢⎥⎢⎥⎣⎦⎣⎦⎡⎤⎡⎤⎛⎫⎛⎫=-==⎢⎥⎢⎥ ⎪ ⎪⎝⎭⎝⎭⎢⎥⎢⎥⎣⎦⎣⎦∑∑∑∑ 假设H 1情况下,同样的方法计算均值和方差为:[]()[][]()11112221111111()||11()||()|N N k k k k N N nk 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=C 11=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:在二元数字通信系统中,假设为H 1时信源输出为1,假设为H0时信源输出为0,信号在通信信道上传输时叠加了均值为零、方差为1的高斯噪声。