第二章高斯噪声背景下谐波恢复数学模型.docx
- 文档编号:9896511
- 上传时间:2023-02-07
- 格式:DOCX
- 页数:23
- 大小:137.61KB
第二章高斯噪声背景下谐波恢复数学模型.docx
《第二章高斯噪声背景下谐波恢复数学模型.docx》由会员分享,可在线阅读,更多相关《第二章高斯噪声背景下谐波恢复数学模型.docx(23页珍藏版)》请在冰豆网上搜索。
第二章高斯噪声背景下谐波恢复数学模型
第二章高斯噪声背景下谐波恢复数学模型
2.1高斯过程
式中
|B|---归一化协方差矩阵的行列式,即
Bjk---行列式|B|中元素bjk的代数余因子;
bjk---归一化协方差函数:
bjk
E{[(tj)-aj][(tQ-ak]}
,方差和协方差决
由上式可以看出,正态随机过程的n维分布仅由各随机变量的数学期望
疋.
一维高斯正态分布的概率密度函数可写为:
f(X)=
2eXP;
(x_a)2
2厂
式中,a及二是两个常量(均值和方差).
高斯噪声一般分为白噪声和有色噪声。
功率谱密度在整个频域内都是均匀分布的噪声,被称之为白噪声,即
显见,白噪声的自相关函数仅在.=0时才不为零;这说明,白噪声只有在零点才相关,而它在任意两个时刻上的随机变量都是不相关的。
有色噪声与白噪声不同,它的功率谱在整个频段上不是均匀分布的。
2.2谐波恢复的数学模型
高斯噪声背景下的谐波恢复,主要是利用特征子空间分析的方法,对观测值进行处理从而估计出原始信号的频率等参量特征,即完成了在噪声背景下对信号的恢复。
我们首先对特征子空间进行分析。
从几何意义上说
协方差几何空间=信号子空间+噪声子空间我们所要做的就是从大的空间抽取低秩子空间,对信号进行分解处理。
下面将要介绍的Pisarenko,MUSIC,Prony,ESPRIT等方法其核心思想都是由此而来的。
设观测信号模型:
y(n)=x(n)+w(n)
j(wz)
x(n)八:
ieJ(ii)
iA
其中,:
i为信号幅度
w为谐波信号频率
<为相位在I」,二均匀分布的参量
p为谐波个数
w(n)为零均值方差为二2的高斯白噪声
构造mm协方差矩阵
_Ry(0)Ry(T)…RyLM杓)[
Ry
(1)Ry(0)…Ry(—M+2)
aa+a
.Ry(M-1)Ry(M)…Ry(0)_
设R=S+W:
S表示信号协方差矩阵,W表示噪声协方差矩阵。
p2
S=送aisisiH
i=1
其中:
.exp(j(M
「°
F面我们对协方差矩阵R进行分析:
R=S+W
p
s=l:
|«i
i4
H
SS
首先我们来说明R为什么等于S+W.y(n)的自相关函数
ry(m)
=E\
Jn)y(n+m)】
-p*
=E*
aiexp^join-
■-p*
=E*
匡ctiexp(—阿门-
p
p
*
=E
S》c(|exp[jn(t\
i壬
p
p
*
=工
cti送ct|exp[jn(们1
1仝
p
2
=工
«iexp(joim)
i)7"*(n)".「*iexp(j■l(n-m^l)-w(n-m)」l_y」
i)':
lexp(j」(nm)jEW*丄yJ.l
一■」j・|m]E[exp(|-[)];「2、(m)
)j(I-i)匚(m)
2、(m)
(n)w(nm)
«i
2$SH
c2|
p
=11
i4
«i
exp(j-i)
exp(j■i)
exp(j(M-1)「)I2I
p
i=1
ai
exp(j(M-1岡)
-
exp(j-i)
exp(j2i)
exp(-j•i)
1
exp(j-i)
exp(-j(M-1)「)
exp(—jM国i)
exp(—j(M+1)①j)
对应矩
exp(j(M-1)⑷J
exp(j(M-2)i)
匚2
阵各项可知现在对矩阵
「°
R=S+W.
R进行分析:
•、二2
(1)设-i,Vi分别是矩阵S的特征值和特征向量,则
i=1~M
讨论:
因为
值
2
(2)W=;「I
讨论:
(1)
M>p,信号阵的秩必为p。
所以
(i=p+1~M)。
M
由于1=7所以W二二
iA
S阵有
p个非零特征值(i=1~p),M-p个零特征
H
ViVi
M
2亠-■-MM
i=1
M
2H
p
R八■iviviH
i=1
p
八('i二「VjVjH亠二■-2vivil
idi»1
R阵与S阵具有相同的特征向量vi,且
(2)Ci,vJi=1〜p
是信号特征对
(打,vj…i=p1~M是噪声特征对
分别记:
Vs=[Vi,Vp]
G=[Vpi,,vm]
由信号向量张成的子空间叫做信号子空间,而有噪声向量张成的子空间叫做噪声
子空间。
(3)矩阵R的特征值
厲丄2..
[為+仃…i=1〜p為={2
戶i二p1〜M
由以上结论我们可以看出,信号向量和噪声子空间中的所有向量(包括它们的线性组合)
正交的。
即
M
SH(瓦GkVk)=0
k-p1
由于此结论十分重要我们在此进行证明。
2H
SiSi
?
p
Si
exp(ji)
jexp(jp^i)一
S二EBEh
0pn
[S1,,Sp]
S2H
2.H
ap丄Sp
2SiSiH
SVpi=EBEHVpi=0v;iEBEHVpi=0(EHVpd)HB(EHVp.J=0-B是正定的
EHVpi=0
展开上式即有证毕;
利用信号向量和噪声空间的正交性进行信号恢复的方法称之为噪声子空间方法。
第三章算法概述与分析
定义观测信号的空间协方差矩阵为
(3.1)
(3.2)
H
R=E[y(n)y(n)]
设噪声方差为二2,由假设得
H
E{[As(n)+w(n)][As(n)+w(n);};
E{[As(n)+w(n)][[*s(n)];+w;(n)]}E{[As(n)+w(n)][s(n)A+;wE[As(n)s(n)A+w(^)s(n)
As(n)w;H(n)+w(n)w(n)];;
=E[As(n)s;(n)A]+E[w(n)(n)A]+
HH
E[As(n)w(n)]+E[w(n)w(n)]
=ASA+;「Im
其中,S=E[s(n)sH(n)],Im为单位阵。
注意到(3.1)式表达的是观测信号的“统计平均”。
我们认为通信中的随机过程是平稳随
机过程,而平稳随机过程具有各态历经性,所以在这里,我们可以用观测信号的“时间平均”
来代替“统计平均”,即
R=E[y(n)yH(n)]
yi(n)yi(n)%(n)y2(n)…yi(n)yM(n)
_yM(n)yi(n)yM(n皿(n)…yM(n)yM(n)
「E[yi(n)%(n)]
_|E[y2(n)yi(n)]
=9
_E[yM(n)yi(n)]
■Rll(0)尺2(0)
R21(0)R22(0)
aa
_Rm1(0)Rm2(0)
R中各元素的估计为:
E[yi(n)y2(n)]
E[y2(n)y2(n)]
a
E[yM(n)y2(n)]
R1M(0)
R2M(0)
Rmm(0)
E[yi(n)yM(n)]
E[y2(n)yM(n)]
a
E[yM(n)yM(n)]
1n
Rj(0)y(n)y「(n)
Nn」
故可得R的估计为
A1NH
Ry(n)y(n)(3.3)
Nn4
3.1Pisarenko谐波分解法
在Pisarenko谐波分解法中,考虑的是由p个实正弦组成的确定性过程
p
x(n)=送Asin(2兀佃+d)(3.1.1)
我们假定,初始相位“是在(-二,二)均匀分布的独立随机变量,它在一次实现中为常量。
先来推导过程{x(n)}满足的差分方程。
为此,先来考虑单个正弦波的情况,即
x(n)=sin(2二fn二).
将三角恒等式
sin(2二fn旳sin[2f(n—2)討=2cos(2二f)sin[2f(n-1)寸]
中的正弦函数换乘x(n)后,得到二阶差分方程式x(n)-2cos(2二f)x(n-1)x(n-2)=0
上式两边去z变换,得
[1-2cos(2二f)z,z°]X(z)=0
这样,我们就得到特征多项式
1-2cos(2二f)z4=0
它有一对共轭负根即
z=cos(2二f)一jsin(2f)=ej2f
共轭根的模为1,即\z\=\z2\=1,由它们可决定正弦频率£二arctan[Im(zi)/Re(zi)]/2■:
通常我们取正的频率。
显然,如果p个实正弦波没有重复频率的话,征多项式
则着p个频率也应由特
P2p
I丨(Z-Zi)(z-Z*)akZ2p上=0
i4k=0
2p
或vakz2p上=0
k卫
(3.1.2)
的根决定。
易知,\z\=1,且系数ak是对称的,即
ak=a2p±(0-k一p)
与式(3.1.2)对应的差分方程为
2p
x(n)亠二aix(n_i)=0
i£
进一步的,我们来考虑白噪声中的正弦波过程:
y(n)=x(n)+w(n)
其中{w(n)}是一零均值、方差为cW的高斯白噪声,即
2
{w(n)}~WN(0,二w)。
将式(3.1.5)代入式(3.1.4),立即可得
2p2p
y(n)…工aiy(n-i)二w(n)'、'aiw(n-i)
这表明白噪声中的正弦波过程{y(n)}是一个特殊的ARMA(2p,2p)
数相同。
为了推导AR参数满足的方程,令
T、
y=[y(n),y(n-1),…,y(n-2p)]
a=[1,a1,…,a2p]T1
w=[w(n),w(n—1),…,w(n-2p)]T
J
于是式(3.1.6)可以写成一下矩阵形式:
TT
ya二wa
用向量y左乘式(3.1.8)两边,并取数学期望得到
E{yyT}a二E{ywT}a
(3.1.3)
(3.1.4)
(3.1.5)
(3.1.6)
过程,其AR参数与MA参
(3.1.7)
(3.1.8)
(3.1.9)
若Ry(l)=E{y(n)y(nI)},则显知,
Ry(O)Ry(_1)
Ry(-2p)
E{yyT}二
Ry
(1)
Ry(O)
Ry(-2p1)
=Ry
.Ry(2p)Ry(2p-1)
Ry(O)
TTT2
E{yw}二E{(xw)w}二E{ww}-.■wI
将上述关系代入式(3.1.9),得特征方程
(3.1.10)
Ryada
式(3.1.10)表明,二W是自相关矩阵Ry的特征值,而a是对应特征值二W的特征向量。
式(3.1.10)
组成了由Pisarenko发展的谐波分解方法的基础。
这样一来,谐波恢复问题就转化成了自相
关矩阵Ry的特征分解。
(1)计算矩阵
(2)求矩阵R最小特征值,即为
即取前m+1八.行和列组成新矩阵
的特征分解。
%,并用表示%的多重度。
然后将R将结为m+1-i,
R1,重复这一步直至从某一阶降至下一阶,最小特征值不
变。
这样,即可定出p,与最小特征值对应的特征向量各分量用a0,a1/',a2p表示。
(3)求多项式的根,有a0■dz亠■亠a2pz2p=0
***
(4)将上式的根记为乙,乙,Z2,Z2,…,Zp,Zp,其中Zi=exp(j2二f「t),由此定出fi(取正值)无噪声的正弦波过程(2.1—1)式的自相关函数为
p
Rx()八Pcos(2£)
i4
其中,R=Aj2/2。
因此,白噪声中的正弦波y(n)=x(n)+w(n)的自相关函数
Ry()二Rx()Rw(),故有
Ry(0)八Pi二
i=1
Ry(.)-vPicos(2二fj.)(.^0)
i二1
P=[R,…,pp]t、
r=[Ry
(1),…,Ry(p)]T卜
{Dj}={cos(2Tfij)}
DP二rD二{Djj}P二D」r
写成一式为
p
&()-;•、()'Rcos(2t)»W、()Rx()
i=1
因此接收信号y(n)的功率谱由上式的傅里叶变换可得
2pP
Py(f)=%+2:
寸[B(f+fi)W(f-fi)]
12
下面进行仿真实验:
信号:
x(n)冷©5®5)•e"5"25)
噪声:
高斯白噪声,方差为1,均值(E)为0
4
2
0
-2
-4
020406080100
frequencyestimation
211[r
1--
0I1LLk1
020406080100
-4
0
1.5
0.5
2
0
-2
100
20406080
4pisarenkoobservalsignalSNR=0
frequencyestimation
20
40
60
80
100
pisarenkoobservalsignalSNR=-2
机*由卄K卄H+Hth址社rUHHHi4吐H卄H+KI•住冉屮卄根卜碎H~H比宙郴啊卄HH齢屮1乂屮忖朴昌
1.5
0.5
20
40
60
80
100
0
20
100
E=0.4715
1.4349
Pisarenko谐波分解法在理论上首次揭示出了相关矩阵对应于最小特征值的特征向量的
系数即为正弦波过程的ARMA模型的AR参数。
但是作为一种算法,pisarenko算法并不有
效,实用。
3.3MUSIC法
从原理上讲,Pisarenko分解时将自相关矩阵分解成信号子空间(对应与大特征值的特
征向量)和噪声子空间(对应于噪声方差的特征向量)。
由于特征多项式
12d
1aiz-•亠a?
pZ=0的根决定正弦频率,而多项式的系数ai又是由对应于噪声方差的
特征向量决定的,所以Pisarenko分解实质是一种噪声子空间方法。
Pisarenko方法不是一种
十分有效的实际算法,但MUSIC和下面将要介绍的ESPRIT方法却十分有效。
MUSIC是多重信号分类(MultipleSignalClassificationCharacterization)的英文缩写。
信号处理中的几个重要问题都可以归结为估计下列模型的参数:
y(t)=A("x(t)e(t)(t=1,…,N)(331)
式中,y(t)・Cmn是带噪声的数据向量,x(t)・Cn1是信号幅值的向量,e(t)・Cm1代表
加性噪声,而矩阵Cmn具有下面的特殊结构:
A(R二[a(J,,a(n)](3.3.2)
这里,{J是实的参数,aCikCm1是第I个信号与y(t)之间的传输向量,且
二二[J,「冷]
对以上模型,我们做以下假定。
A1:
m>n,且对应于不同「值的向量a(■)是线性独的。
A2:
E{e(t)}=0,E{e(t)eH(t)}=;:
2l,E{e(t)eT(t)}=0
这里H代表向量的共轭转置,T代表转置。
A3:
矩阵P=E{x(t)xH(t)}是非奇异的,且N>m.当满足假设条件时,观测向量y(t)的协方
差矩阵由
R?
E{y(t)yH(t)}二A⑺PAHQ)‘1(3.3.3)
给定。
为方便计,我们将A(^)简记为A。
注意到R是一对称阵,其特征值分解具有下列形式:
R?
U32UH,其中,
——2222
Z=diag[;_,…,;「m]对角线元素、二二i叫做R的特征值。
在A1条件下,矩阵A显然是非奇异的,即rank(A)=n从而,在A3条件下APAH的秩
也为n。
因此R的特征值必须满足下列关系:
2
i匚(i二1,,n)
i-;「(i=n1「,m)
不妨令’1,…,’n分别对应的特征向量为S,…,Sn,而’n1,i,'m,对应的特征向量为
5,…,gm_n,定义
S=[3,,Sn]
G=[g1,,gm_n]
U分为两个子矩阵U=[S|G]。
S和G分别叫做信号子空间和噪声子空间。
于是,特征矩阵
我们的兴趣是研究噪声子空间G和矩阵A之间的关系。
一方面,由于匚2和G分别是R
的特征值和对应的特征向量,故有特征方程
(334)
另一方面,用G右乘式(3.3.3),又有
综合上面两式给出
从而有
由于P是非奇异的,所以
上式又可写作
由于U是酉阵,故
于是式(3.3.7)也可写作
ghAPAHG=(AHG)HP(AHG)=0。
_SH1
UUH=[S|G].|歸
-G-
SSHGGH二I
MUSIC算法的基本思想是对真实相关阵R使用性质(337)或(338)。
在实际应用中,R是未知的,但它可从观测数据一致估计,即用下式计算
1N
R=-y(t)yH(t)
N14
类似于R的特征分解,令{$,…,Sn,gir,gm』}表示R的归一化特征向量,且按特征值降序排列,定义
Hh
f(,)二a()GGa()(3.3.9)
=aH()(I-SSH)a()(3.3.10)
一般地说,当n>m-n时使用(3.3.9);反之,使用(3.3.10)。
使f(,)最小的n个•.值即为所求。
下面叙述MUSIC算法的步骤:
(1)利用观测值作出样本自相关矩阵R。
(2)对R进行特征值分解,求出n以及,1,…,n,;「,3,…,sn。
(3)利用式(3.3.9)或(3.3.10)计算出n个值。
下面进行仿真实验:
j(0.5t0.15)j(1.5t025)
信号:
x(n)=ee
噪声:
高斯白噪声,方差为1,均值(E)为0
MUSICobservalsignalSNR=0
-4
2
0
-2
4MUSICobservalsignalSNR=-2
0100200300400500600
3000
0
2000
1000
00.511.522.533.544.555.56
00.511.522.533.544.555.56
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第二 章高斯 噪声 背景 谐波 恢复 数学模型