2014年信号检测与估计各章作业参考答案(1~9章)
- 格式:doc
- 大小:3.46 MB
- 文档页数:32
《信号检测与估计》第十章习题解答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 不相关。
信号检测与估计答案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的高斯白噪声。
2-1 1[()]2E x t =,1212(,)3X t t R t t = 2-2 略。
2-3111[()]sin cos 333E x t t t=++12112212121111111(,)sin cos sin cos sin()cos()9999999X R t t t t t t t t t t =+++++++-2-4 [()]0E X t =,20(,)cos R t t w τστ+=2-5 [()]0E X t =,20(,)cos 2a R t t w ττ+= 2-6 略。
2-7 [()]0E X t =,10(,)200R t t τττ⎧=⎪+=⎨⎪≠⎩2-8 1210()()()2cos(10)(21)X X X R R R eτττττ-=+=++,2[()](0)5X E X t R ==,2(0)2X X R σ==2-9 11()()cos 22jw jw X X o G w R e d w e d τττττ∞∞---∞-∞==⎰⎰00()()()22X P w w w w w ππδδ=-++2-10 00()(()())2Y X X aG w G w w G w w =-++2-11 ())()X R u ττ=+-3-1 二元信号统计检测的贝叶斯平均代价C 为110000000100100110111111()()=()()()() ()()()()ij i i j j i C c P H P H H c P H P H H c P H P H H c P H P H H c P H P H H ===+++∑∑ 利用01()1()P H P H =-1101()1()P H H P H H =- 0010()1()P H H P H H =-得平均代价C 为[][]0011010110011011110100101110111000111011000101()1()1()() ()()()1() ()() ()()()()()()C c P H P H H c P H P H H c P H P H H c P H P H H c c c P H H P H c c c c P H H c c P H H =-⎡-⎤+-+⎣⎦+⎡-⎤⎣⎦=+-+⎡-+---⎤⎣⎦3-2 1)由于各假设j H 的先验概率()(0,1,2)j P H j =相等,所以采用最大似然准则。
时间: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.52ln 88.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 =⎰而2121(1)(|)exp[]22x p x H σπσ-=- 2021(|)exp[]22x p x H σπσ=-所以201121(1)(|)(|)exp[]22D D x P D H p x H dx dx σπσ-==-⎰⎰221(1)e x p []22x dx βσπσ-∞-=-⎰=17.82()()(3.91)22β-Φ=Φ=Φ 同理11210021(|)(|)exp[]22D D x P D H p x Hdx dx σπσ==-⎰⎰221e x p ()22x dx βσπσ∞=-⎰8.821()1()1(4.41)22β=-Φ=-Φ=-Φ 所以0.21(3.91)0.82[1(4.41)]R =⨯⨯Φ+⨯⨯-Φ 1.2 (关瑞东 20060155)假定加性噪声()n t 服从均值为零,方差为的正态分布。
nnms ms ms 信号检测与估计第四章1. 设观测矢量xasn ,其中a 为未知参量,噪声矢量n ~N (0,2 I),s 为已知矢量,假定a 与n 相互独立,如果a 在[b,c]上均匀分布,即1 p ( a),acb[b,c]试求未知参量a 的最大后验概率估计(a ˆmap )和最小均方误差估计(a ˆ )。
2. 如果a 的分布变为1 1 1 p ( a)(a 1)( a)(a 1)333其它条件同习题 1,试求未知参量a 的最大后验概率估计(a ˆmap )和最小均方误差估计(a ˆ )。
3. 如果a 是均值为m ,方差为2的高斯分布,即a ~ N(m,2),其它条件同习题1,试求未知参量a 的线性最小均方误差估计(a ˆl ms ),最大后验概率估计(a ˆmap和最小均方误差估计(a ˆ )。
4. 设观测矢量xasn ,其中a 为未知参量,噪声矢量n~N(0,2I),s 为已知矢量,求未知参量a 的最小二乘估计(a ˆL S )和最大似然估计(a ˆm l )。
5. 设观测矢量x asn ,其中a 为未知参量,噪声矢量n ~ N ( 0 , R ) ,s 为已知矢量,求未知参量a 的最小二乘估计(aˆL S ),加权最小二乘估计(a ˆWLS )和最大似然估计(a ˆm l )。
)nn 2 ) TTT T 2 T)第 1题解:由题意,以a 为条件的观测矢量x( x , x,, x )T 的条件概率密度为12Np ( x |a ) ;而且对于某一个a ,则x ~ N(as,量。
所以2I ),其中s (s 1,s 2, ,s N )T 为已知矢p(x|a)1exp{1(x as)T (2 I )1 (x as)}(2 )N /2| 2I|1/22 n(1-1)1 1exp{ (x as)T (xas)}(2nN/22 n而随机参量的概率密度函数为p(a)1,a [b,c] (1-2)c b为了求得后验概率p ( a |x ),利用p ( a |x)p ( x |a ) p( a)p(x)(1-3)因为p ( a|x ) 是给定了x 后,a 的条件概率密度函数,所以对于p ( a |x ) 而言,p( x) 相当于使p ( a |x)da1(1-4)的归一化因子,所以p ( a |x)1 p ( x ) (21N /2n1 e xp{(x 2na s ) (x1 a s )}c b1p(x) 1c b (2 1N /2n1 e xp{(x 2 na s ) (xa s )}(1-5)K (x)exp{1 T( xx2n1 Ta s xTa x s2 Ta s s )}上式中,K (x)exp{( xx2n2 a s xa s s )}K(x)111(1-6)是一个与a 无关项。
信号检测及估计试题-答案(不完整版)一、概念: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。
一、简答题注释简答题(每题5分,共20分)或(每题4分,共20分)二、第1章简答题1.从系统和信号的角度看,简述信号检测与估计的研究对象。
答:从系统的角度看,信号检测与估计的研究对象是加性噪声情况信息传输系统中的接收设备。
从信号的角度看,信号检测与估计的研究对象是随机信号或随机过程。
2.简述信号检测与估计的基本任务和所依赖的数学基础。
答:解决信息传输系统接收端信号与数据处理中信息恢复与获取问题,或从被噪声及其他干扰污染的信号中提取、恢复所需的信息。
信号检测与估计所依赖的数学基础是数理统计中贝叶斯统计的贝叶斯统计决策理论和方法。
3.概述信号在传输过程中与噪声混叠在一起的类型。
答:信号在传输过程中,噪声与信号混杂在一起的类型有3种:噪声与信号相加,噪声与信号相乘(衰落效应),噪声与信号卷积(多径效应)。
与信号相加的噪声称为加性噪声,与信号相乘的噪声称为乘性噪声,与信号卷积的噪声称为卷积噪声。
加性噪声是最常见的干扰类型,也是最基本的,因为乘性噪声和卷积噪声的情况均可转换为加性噪声的情况。
三、第2章简答题1.简述匹配滤波器概念及其作用。
答:匹配滤波器是在输入为确定信号加平稳噪声的情况下,使输出信噪比达到最大的线性系统。
匹配滤波器的作用:一是使滤波器输出有用信号成分尽可能强;二是抑制噪声,使滤波器输出噪声成分尽可能小,减小噪声对信号处理的影响。
2.根据匹配滤波器传输函数与输入确定信号及噪声的关系,简述匹配滤波器的原理。
答:匹配滤波器传输函数等于输入确定信号频谱的复共轭除以输入平稳噪声的功率谱密度,再附加相位项T ω-,其中T 为输入确定信号的持续时间或观测时间。
由于匹配滤波器传输函数的幅频特性与输入确定信号的幅频特性成正比,与输入噪声的功率谱密度成反比;对于某个频率点,信号越强,该频率点的加权系数越大,噪声越强,加权越小。
从而起到加强信号,抑制噪声的作用。
对于信号,匹配滤波器的相频特性与输入信号的相位谱互补,使输入信号经过匹配滤波器以后,相位谱将全部被补偿掉。
信号检测与估计理论第一章习题讲解1-9已知随机变量某的分布函数为0F某(某)k某21,某0,0某1,某1求:①系数k;②某落在区间(0.3,0.7)内的概率;③随机变量某的概率密度。
解:第①问利用F某(某)右连续的性质k=1第②问P0.3某0.7PF0.7F00.某3.3某.70.7P0dF某(某)2某第③问f某(某)d某00某1ele1-10已知随机变量某的概率密度为f某(某)ke普拉斯分布),求:某(某)(拉①系数k②某落在区间(0,1)内的概率③随机变量某的分布函数解:第①问f某某1d某11k2F2某F1某某1第②问P某2某某2f某d某随机变量某落在区间(某1,某2]的概率P{某1某某2}就是曲线yf某下的曲边梯形的面积。
P0某1P0某1f某d某0111e12第③问1某e2f某1e某2某0某0F某某f(某)d某1某某0e21某某01e2某0某1某ed某201e某d某某1e某d某022某01-11某繁忙的汽车站,每天有大量的汽车进出。
设每辆汽车在一天内出事故的概率为0.0001,若每天有1000辆汽车进出汽车站,问汽车站出事故的次数不小于2的概率是多少?(0-1)分布n,p0,np=二项分布泊松分布n成立,p,q0不成立高斯分布汽车站出事故的次数不小于2的概率P(k2)1Pk0Pk10.1P(k2)11.1e答案P某k=n=1实际计算中,只需满足n10kek!p0.1,二项分布就趋近于泊松分布=np1-12已知随机变量(某,Y)的概率密度为(3某4y)kef某Y(某,y)0,某0,y0,其它求:①系数k?②(某,Y)的分布函数?③P{0某1,0某2}?第③问方法一:联合分布函数F某Y(某,y)性质:若任意四个实数a,a,b,b,满足1212a1a2,b1b2,则P{a1某a2,b1Yb2}F某Y(a2,b2)F某Y(a1,b1)F某Y(a1,b2)F某Y(a2,b1)P{0某1,0Y2}F某Y(1,2)F某Y(0,0)F某Y(1,0)F某Y(0,2)方法二:利用P{(某,y)D}f某Yu,vdudvD20P{0某1,0Y2}0f某Y某,yd某dy11-13已知随机变量(某,Y)的概率密度为1,0某1,y某f(某,y)0,其它①求条件概率密度f某(某|y)和fY(y|某)?②判断某和Y是否独立?给出理由。
《信号检测与估计》第九章习题解答9.1 接收信号(((t n t A t x ++=θω0sin ,其中(t n 是高斯白噪声,θ在(π20,均匀分布,现在需求振幅A 的最大似然估计量。
由于θ的先验知识已知,故可先对θ求平均得到(A x f ,试问要求振幅A 的最大似然估计量必须解什么样的方程? 解:接收信号(t x 的似然函数为((([]((((((((⎟⎠⎞⎜⎝⎛+++−−+++−−+−−∫∫∫∫∫===T TT TTdt t A dt t t x A dt t x N dtt A t t Ax t xN dtt A t x N FeFeFeA x f 0002200200022020200sin sin 21sin sin 21sin 1,θωθωθωθωθωθ由于(((∫=+−=∫+TTT dt t dt t 0000222cos 121sin θωθω,得到 ((((020000202sin 21,N TA dt t t x N A dtt x N e eFeA x f T T−+−∫∫=θωθ对θ积分,得到(((((((((θπθπθθθπθωθωπθωπd eeFed e e Fed f A x f A x f dt t t t x N A N T A dt t x N dt t t x N A N T A dt t x N T TTT∫∫∫∫∫∫∫+−−+−−===20cos cos sin sin 22120sin 221 2000000202000020202121,令(ϕωcos cos 00z dt t t x x Tc ==∫,(ϕωsin sin 00z dt t t x x T s ==∫,得到222s cx x z +=,csx x arctg =ϕ (((((⎟⎟⎠⎞⎜⎜⎝⎛====∫∫∫∫∫−+++0020cos 220cos cos sin sin 220cos sin 220 cos cos sin sin 22212121210 000N Az I d ed ed e d eN Azz z N Ax x N Adt t t t x N A c s Tθπθπθπθππϕθπθϕθϕπθθπθωθω 上式中,[](cos exp 21020x I d x =∫πθθπ为零阶修正贝塞尔函数。