第一章 离散随机信号
1.1 引言
定义 :设已给概率空间 ( Ω , ζ , P ) (\Omega,\zeta,P) (Ω,ζ,P) , Z Z Z 为整数集,若对每一整数 n ( n ∈ Z ) n(n\in Z) n(n∈Z) ,均有定义在 ( Ω , ζ , P ) (\Omega,\zeta,P) (Ω,ζ,P) 上的一个随机变量 x ( n , ξ ) ( ξ ∈ Ω ) x(n,\xi)\quad(\xi \in \Omega) x(n,ξ)(ξ∈Ω) 与之对应,则称依赖与参数 n n n 的一列随机变量 x ( n , ξ ) x(n,\xi) x(n,ξ) 为一离散时间随机过程或随机序列,记为 { x ( n , ξ ) , ξ ∈ Ω , n ∈ Z } \{x(n,\xi),\xi \in\Omega,n\in Z\} {x(n,ξ),ξ∈Ω,n∈Z} ,简记为 { x n } \{x_n\} {xn} 。
随机序列有一下特点:
(1)随机序列中的任何一个点上的取值都是不能先验确定的随机变量。
(2)随机序列可以用它的统计平均特性来表征。随机序列在个时间点上的随机变量取值是服从某种确定的概率分布。
(3)平稳随机信号的能量化表示。一随机信号的能量谱称为功率谱密度。一个平稳随机信号的功率谱是确定的,因此功率谱可以统计表征一个随机过程的谱特性。
综上所述,对于离散随机信号的概念和表征问题,有:
(1)一个随机信号在各个时间点上的取值以及不同点上取值之间的相互关联性只能用概率特性或统计平均特性来表征,它的确定值是无法先验表达的。
(2)一个平稳随机信号在各频率点上的能量的取值可以用功率谱,密度函数与自相关函数统计描述。
1.2 离散时间随机信号的时域(统计)表示
一般情况下,我们大多只对离散随机信号的一维概率分布和二维概率分布等有限维度概率分布进行研究。
1.2.1 离散时间随机过程的概率分布
(宽)平稳随机过程:若一个随机过程的数学期望及方差与时间无关,相关函数仅与时间间隔有关,则称这个随机过程为广义平稳随机过程。
(宽)平稳随机过程的三个特征:
(1) E ( x n ) = C ξ E(x_n) = C_\xi E(xn)=Cξ 均值是一个常数(随机变量)。
(2) E [ x ( n ) x ∗ ( n − m ) ] = R ( m ) E[x(n)x^*(n-m)] = R(m) E[x(n)x∗(n−m)]=R(m) 自相关函数只与两点的差有关。
(3) E [ x 2 ( n ) ] < + ∞ E[x^2(n)] < +\infty E[x2(n)]<+∞ 均方值有界。
严平稳随机过程:所谓随机过程严格平稳,是指它的任何n维分布函数或概率密度函数与时间起点无关。
严平稳随机过程特点:其数学期望、方差与时间无关,自相关函数仅与时间间隔有关。
- (宽)平稳随机过程 和 严平稳随机过程 的区别和联系
(1)一个宽平稳过程不一定是严平稳过程,一个严平稳过程也不一定宽平稳过程。
(2)宽平稳过程定只涉及与一维、二维分布有关的数字特征,所以一个严平稳过程只要二阶矩存在,则必定是宽平稳过程。但反过来,一般是不成立的。
(3)正态过程是一个重要特例,一个宽平稳的正态过程必定是严平稳的。
各态历经随机过程:如果一个平稳随机过程,只要满足一些较宽的条件,则一个样本函数在整个时间轴上的平均值可以用来代替其集平均(统计平均值和自相关函数等),这就是各态历经随机过程。
一般来说,在一个随机过程中,不同样本函数的时间平均值是不一定相同的,而集平均则是一定的。因此,一般的随机过程的时间平均≠集平均,只有平稳随机过程才有可能是具有各态历经性的。即各态历经的随机过程一定是平稳的,而平稳的随机过程则需要满足一定条件才是各态历经的。
1.2.2 离散时间随机过程的数字特征
- 数学期望(均值) 在n点上的
随机变量 x ( n ) x(n) x(n) 的均值(用 m x ( n ) m_x(n) mx(n) 表示)定义为
m x ( n ) = E [ x ( n ) ] ≜ ∫ − ∞ ∞ x p 1 ( x ) d x m_x(n)=E[x(n)]\triangleq \int_{-\infty}^{\infty}xp_1(x)dx mx(n)=E[x(n)]≜∫−∞∞xp1(x)dx
- 均方值(功率)
随机变量 x ( n ) x(n) x(n) 的均方值定义为 ∣ x ( n ) ∣ 2 |x(n)|^2 ∣x(n)∣2 的期望
E [ ∣ x ( n ) ∣ 2 ] ≜ ∫ − ∞ ∞ ∣ x ∣ 2 p 1 ( x ) d x E[|x(n)|^2]\triangleq\int_{-\infty}^{\infty}|x|^2p_1(x)dx E[∣x(n)∣2]≜∫−∞∞∣x∣2p1(x)dx
- 方差
随机变量 x ( n ) x(n) x(n) 的方差(用 σ x 2 ( n ) \sigma_x^2(n) σx2(n) 表示)定义为 [ x ( n ) − E ( x ( n ) ) ] [x(n)-E(x(n))] [x(n)−E(x(n))] 的均方值
σ x 2 ( n ) ≜ E { ∣ x ( n ) − E ( x ( n ) ) ∣ 2 } ≜ E { ∣ x ( n ) − m x ( n ) ∣ 2 } \sigma_x^2(n)\triangleq E\{|x(n)-E(x(n))|^2\}\triangleq E\{|x(n)-m_x(n)|^2\} σx2(n)≜E{∣x(n)−E(x(n))∣2}≜E{∣x(n)−mx(n)∣2}
通过以上公式可以推导出:
平均功率 = 交流功率 + 直流功率
σ x 2 ( n ) = E [ ∣ x ( n ) ∣ 2 ] − ∣ m x 2 ( n ) ∣ E [ ∣ x ( n ) ∣ 2 ] = σ x 2 ( n ) + ∣ m x 2 ( n ) ∣ } \left. \begin{matrix} \sigma_x^2(n)=E[|x(n)|^2]-|m_x^2(n)|\\ E[|x(n)|^2]=\sigma_x^2(n)+|m_x^2(n)|\\ \end{matrix} \right\} σx2(n)=E[∣x(n)∣2]−∣mx2(n)∣E[∣x(n)∣2]=σx2(n)+∣mx2(n)∣}
以上三个特征量仅与一维概率密度 p 1 ( x ) p_1(x) p1(x) 有关。对于平稳随机过程,其一维概率密度与时间无关,故一个平稳随机序列的 m x ( n ) 、 E ( ∣ x ( n ) ∣ 2 ) 、 σ x 2 ( n ) m_x(n)、E(|x(n)|^2)、\sigma_x^2(n) mx(n)、E(∣x(n)∣2)、σx2(n) 均是与时间无关的常数,可以将变量 n n n 去掉。
与二维概率分布有关的统计特性主要有自相关函数和自协方差函数。
- 自相关函数
一个平稳随机过程的自相关函数用 ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 表示
ϕ x x ( m ) ≜ E [ x ( n 1 ) x ∗ ( n 2 ) ] = E [ x ( n ) x ∗ ( n − m ) ] = ∫ − ∞ ∞ ∫ − ∞ ∞ x 1 x 2 ∗ p 2 ( x 1 , x 2 ; m ) d x 1 d x 2 \phi_{xx}(m)\triangleq E[x(n_1)x^*(n_2)]=E[x(n)x^*(n-m)]\\ =\int_{-\infty}^\infty\int_{-\infty}^\infty x_1x_2^*p_2(x_1,x_2;m)dx_1dx_2 ϕxx(m)≜E[x(n1)x∗(n2)]=E[x(n)x∗(n−m)]=∫−∞∞∫−∞∞x1x2∗p2(x1,x2;m)dx1dx2
- 自协方差函数
一个平稳随机序列的协方差函数用 γ x x ( m ) \gamma_{xx}(m) γxx(m) 表示
γ x x ( m ) ≜ E [ ( x ( n ) − m x ) ( x ( n + m ) − m x ) ∗ ] \gamma_{xx}(m)\triangleq E[(x(n)-m_x)(x(n+m)-m_x)^*] γxx(m)≜E[(x(n)−mx)(x(n+m)−mx)∗]
自相关函数 与 自协方差函数 的关系
γ x x ( m ) = ϕ x x ( m ) − ∣ m x ∣ 2 \gamma_{xx}(m)=\phi_{xx}(m)-|m_x|^2 γxx(m)=ϕxx(m)−∣mx∣2
1.2.3 离散时间平稳过程相关序列与协方差性质
两个实的、平稳相关的平稳随机序列
{
x
(
n
)
}
\{x(n)\}
{x(n)} 和
{
y
(
n
)
}
\{y(n)\}
{y(n)} 的
{
x
(
n
)
}
\{x(n)\}
{x(n)} 的自相关、自协方差及两个序列的互相关和互协方差序列定义为
ϕ
x
x
(
m
)
=
E
[
x
(
n
)
x
(
n
−
m
)
]
γ
x
x
(
m
)
=
E
[
(
x
(
n
)
−
m
x
)
(
x
(
n
−
m
)
−
m
x
)
]
\phi_{xx}(m)=E[x(n)x(n-m)]\\ \gamma_{xx}(m)=E[(x(n)-m_x)(x(n-m)-m_x)]\\
ϕxx(m)=E[x(n)x(n−m)]γxx(m)=E[(x(n)−mx)(x(n−m)−mx)]
ϕ x y ( m ) = E [ x ( n ) y ( n − m ) ] γ x y ( m ) = E [ ( x ( n ) − m x ) y x ( n − m ) − m y ) ] \phi_{xy}(m)=E[x(n)y(n-m)]\\ \gamma_{xy}(m)=E[(x(n)-m_x)yx(n-m)-m_y)]\\ ϕxy(m)=E[x(n)y(n−m)]γxy(m)=E[(x(n)−mx)yx(n−m)−my)]
性质 1
γ
x
x
(
m
)
=
ϕ
x
x
(
m
)
−
m
x
2
γ
x
y
(
m
)
=
ϕ
x
y
(
m
)
−
m
x
m
y
\gamma_{xx}(m)=\phi_{xx}(m)-m_x^2\\ \gamma_{xy}(m)=\phi_{xy}(m)-m_xm_y
γxx(m)=ϕxx(m)−mx2γxy(m)=ϕxy(m)−mxmy
性质 2
ϕ
x
x
(
0
)
=
E
[
x
2
(
n
)
]
均
方
值
γ
x
x
(
0
)
=
ϕ
x
x
(
0
)
−
m
x
2
方
差
\phi_{xx}(0)=E[x^2(n)]\qquad {均方值}\\ \gamma_{xx}(0)=\phi_{xx}(0)-m_x^2\qquad{方差}
ϕxx(0)=E[x2(n)]均方值γxx(0)=ϕxx(0)−mx2方差
性质 3
ϕ
x
x
(
m
)
\phi_{xx}(m)
ϕxx(m) 和
γ
x
x
(
m
)
\gamma_{xx}(m)
γxx(m) 都是偶函数
ϕ
x
x
(
m
)
=
ϕ
x
x
(
−
m
)
γ
x
x
(
m
)
=
γ
x
x
(
−
m
)
\phi_{xx}(m)=\phi_{xx}(-m)\\ \gamma_{xx}(m)=\gamma_{xx}(-m)
ϕxx(m)=ϕxx(−m)γxx(m)=γxx(−m)
性质 4
ϕ
x
y
(
m
)
=
ϕ
y
x
(
−
m
)
γ
x
y
(
m
)
=
γ
y
x
(
−
m
)
\phi_{xy}(m)=\phi_{yx}(-m)\\ \gamma_{xy}(m)=\gamma_{yx}(-m)
ϕxy(m)=ϕyx(−m)γxy(m)=γyx(−m)
性质 5
∣
ϕ
x
x
(
m
)
∣
≤
ϕ
x
x
(
0
)
∣
γ
x
x
(
m
)
∣
≤
γ
x
x
(
0
)
|\phi_{xx}(m)|\leq\phi_{xx}(0)\\ |\gamma_{xx}(m)|\leq\gamma_{xx}(0)
∣ϕxx(m)∣≤ϕxx(0)∣γxx(m)∣≤γxx(0)
性质 6 如果
y
(
n
)
=
x
(
n
−
m
)
y(n) = x(n-m)
y(n)=x(n−m) ,则
ϕ
y
y
(
m
)
=
ϕ
x
x
(
m
)
γ
y
y
(
m
)
=
γ
x
x
(
m
)
\phi_{yy}(m)=\phi_{xx}(m)\\ \gamma_{yy}(m)=\gamma_{xx}(m)
ϕyy(m)=ϕxx(m)γyy(m)=γxx(m)
性质 7 对于大多随机序列,
m
m
m 愈大则相关性愈小,当
m
m
m 趋近无穷大时,可以认为不相关。若平稳随机序列
{
x
(
n
)
}
\{x(n)\}
{x(n)} 中不含有任何周期分量,则有
KaTeX parse error: No such environment: align at position 8: \begin{̲a̲l̲i̲g̲n̲}̲ \lim\limits_{m…
由以上性质可以得出一下几点:
(1) γ x x ( m ) \gamma_{xx}(m) γxx(m) 和 ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 除相差一个常量 m x 2 m_x^2 mx2 外,其它特性相同。
(2) ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 是一个随机过程的最主要的统计特征,它不仅说明了相关性还隐含了 m x , σ x 2 , E [ x 2 ( n ) ] m_x,\sigma_x^2,E[x^2(n)] mx,σx2,E[x2(n)] 等主要特征量,例:
E [ x 2 ( n ) ] = ϕ x x ( 0 ) m x 2 = ϕ x x ( ∞ ) σ x 2 = E [ x 2 ( n ) ] − m x 2 = ϕ x x ( 0 ) − ϕ x x ( ∞ ) } \left. \begin{aligned} &E[x^2(n)]=\phi_{xx}(0)\\ &m_x^2=\phi_{xx}(\infty)\\ &\sigma_x^2=E[x^2(n)]-m_x^2=\phi_{xx}(0)-\phi_{xx}(\infty) \end{aligned} \right\} E[x2(n)]=ϕxx(0)mx2=ϕxx(∞)σx2=E[x2(n)]−mx2=ϕxx(0)−ϕxx(∞)⎭⎪⎬⎪⎫
对于多为随机向量
X
=
[
x
1
x
2
⋮
x
n
]
\pmb{X}= \left[ \begin{matrix} x_1\\ x_2\\ \vdots\\ x_n \end{matrix} \right]
XXX=⎣⎢⎢⎢⎡x1x2⋮xn⎦⎥⎥⎥⎤
其均值向量为
E
[
X
]
=
[
E
[
x
1
]
E
[
x
2
]
⋮
E
[
x
n
]
]
E[\pmb{X}]= \left[ \begin{matrix} E[x_1]\\ E[x_2]\\ \vdots\\ E[x_n]\\ \end{matrix} \right]
E[XXX]=⎣⎢⎢⎢⎡E[x1]E[x2]⋮E[xn]⎦⎥⎥⎥⎤
其自协方差阵定义为
c
o
v
[
X
]
=
E
{
[
X
−
E
(
X
)
]
[
X
−
E
(
X
)
]
T
}
cov[\pmb{X}]=E\{[\pmb{X}-E(\pmb{X})][\pmb{X}-E(\pmb{X})]^T\}
cov[XXX]=E{[XXX−E(XXX)][XXX−E(XXX)]T}
其互协方差阵定义为
c
o
v
(
X
,
Y
)
=
E
{
[
X
−
E
(
X
)
]
[
Y
−
E
(
Y
)
]
T
}
cov(\pmb{X},\pmb{Y})=E\{[\pmb{X}-E(\pmb{X})][\pmb{Y}-E(\pmb{Y})]^T\}
cov(XXX,YYY)=E{[XXX−E(XXX)][YYY−E(YYY)]T}
如果
X
,
Y
\pmb{X},\pmb{Y}
XXX,YYY 不是向量,而是复数,则公式(25)(26)中的转置符号应该为共轭符号,如果
X
,
Y
\pmb{X},\pmb{Y}
XXX,YYY 既是向量,又是复数,则转置符号应该为共轭转置符号。
1.2.4 平稳序列的时间平均与遍历性
随机序列某一点处的均值和自相关函数,可以用平均集来计算,即
KaTeX parse error: No such environment: align at position 8: \begin{̲a̲l̲i̲g̲n̲}̲ &m_x=\lim\limi…
在实践中需要多次大量的实验以获得足够多的(
N
N
N 足够大)序列来保证计算的可靠性,显然这是十分困难的。但是我们假设一个平稳过程,集合平均等于一个样本函数在整个时间轴上的时间平均,这样就能使用时间平均对随机序列的均值和自相关函数做出合理近似估计。
定义 设
{
x
(
n
)
}
\{x(n)\}
{x(n)} 是一个平稳随机序列,如果均值极限
<
x
(
n
)
>
≜
lim
N
→
∞
1
2
N
+
1
∑
n
=
−
N
N
x
(
n
)
<x(n)> \triangleq \lim\limits_{N\rightarrow\infty}\frac{1}{2N+1}\sum_{n=-N}^{N}x(n)
<x(n)>≜N→∞lim2N+11n=−N∑Nx(n)
存在,称它为
{
x
(
n
)
}
\{x(n)\}
{x(n)} 的时间均值。进而,如果以概率为 1 有
m
x
=
<
x
(
n
)
>
m_x=<x(n)>
mx=<x(n)>
则称
{
x
(
n
)
}
\{x(n)\}
{x(n)} 的均值具有遍历性。
对于平稳随机序列
{
x
(
n
)
}
\{x(n)\}
{x(n)} ,如果对
∀
m
∈
Z
\forall m \in Z
∀m∈Z ,均方极限
<
x
(
n
)
x
∗
(
n
−
m
)
>
≜
lim
N
→
∞
1
2
N
+
1
∑
n
=
−
N
N
[
x
(
n
)
x
∗
(
n
−
m
)
]
<x(n)x^*(n-m)> \triangleq \lim\limits_{N\rightarrow\infty}\frac{1}{2N+1}\sum_{n=-N}^{N}[x(n)x^*(n-m)]
<x(n)x∗(n−m)>≜N→∞lim2N+11n=−N∑N[x(n)x∗(n−m)]
存在,称它为
{
x
(
n
)
}
\{x(n)\}
{x(n)} 的时间自相关函数。进而,如果以概率为 1 有
ϕ
x
x
(
m
)
=
<
x
(
n
)
x
∗
(
n
−
m
)
>
\phi_{xx}(m)=<x(n)x^*(n-m)>
ϕxx(m)=<x(n)x∗(n−m)>
则称
{
x
(
n
)
}
\{x(n)\}
{x(n)} 的自相关函数具有遍历性。
这里注意,当(30)(32)成立,因为均方收敛必依概率收敛,那么当 N N N 足够大时,就可以用一次抽样所得的一个现实 { x n ( ξ ) , n ∈ Z } ( ξ ∈ Ω ) \{x_n(\xi),n\in Z\}(\xi\in\Omega) {xn(ξ),n∈Z}(ξ∈Ω) 所作的时间平均
1 2 N + 1 ∑ n = − N N x ( n , ξ ) \frac{1}{2N+1}\sum_{n=-N}^{N}x(n,\xi) 2N+11n=−N∑Nx(n,ξ)
与
1 2 N + 1 ∑ n = − N N [ x ( n , ξ ) x ∗ ( n − m , ξ ) ] \frac{1}{2N+1}\sum_{n=-N}^{N}[x(n,\xi)x^*(n-m,\xi)] 2N+11n=−N∑N[x(n,ξ)x∗(n−m,ξ)]
分别作为 m x m_x mx 和 ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 的近似估计。
思考 :证明, ∣ ϕ x y ∣ ≤ ϕ x x ( 0 ) ϕ y y ( 0 ) |\phi_{xy}|\leq \sqrt{\phi_{xx}(0)\phi_{yy}(0)} ∣ϕxy∣≤ϕxx(0)ϕyy(0) .
构造:
E [ λ x ( n ) + y ( n + m ) ] 2 ≥ 0 ⇒ λ 2 E [ x 2 ( n ) ] + 2 λ E [ x ( n ) y ( n + m ) ] + E [ y 2 ( n + m ) ] ≥ 0 ⇒ ϕ x x ( 0 ) λ 2 + 2 ϕ x y ( m ) λ + ϕ y y ( 0 ) ≥ 0 , n = 0 因 为 ϕ x x ( 0 ) > 0 , 故 Δ = 4 ∣ ϕ x y ( m ) ∣ 2 − 4 ϕ x x ( 0 ) ϕ y y ( 0 ) > 0 即 ∣ ϕ x y ∣ ≤ ϕ x x ( 0 ) ϕ y y ( 0 ) \begin{aligned} &E[\lambda x(n)+y(n+m)]^2\geq 0\\ \Rightarrow &\lambda^2E[x^2(n)]+2\lambda E[x(n)y(n+m)]+E[y^2(n+m)] \geq 0\\ \Rightarrow &\phi_{xx}(0) \lambda^2 +2\phi_{xy}(m)\lambda+\phi_{yy}(0) \geq 0 ,\quad n=0 \\ 因为 \\ &\phi_{xx}(0) > 0,故 \Delta = 4|\phi_{xy}(m)|^2-4\phi_{xx}(0) \phi_{yy}(0) > 0\\ 即 \\ &|\phi_{xy}|\leq \sqrt{\phi_{xx}(0)\phi_{yy}(0)} \end{aligned} ⇒⇒因为即E[λx(n)+y(n+m)]2≥0λ2E[x2(n)]+2λE[x(n)y(n+m)]+E[y2(n+m)]≥0ϕxx(0)λ2+2ϕxy(m)λ+ϕyy(0)≥0,n=0ϕxx(0)>0,故Δ=4∣ϕxy(m)∣2−4ϕxx(0)ϕyy(0)>0∣ϕxy∣≤ϕxx(0)ϕyy(0)
1.3 离散时间随机信号的 z z z 域及频域(统计)表示
平稳随机信号无始无终,能量无限,且不满足绝对可和(绝对值的和是一个有限常数)。
故不存在傅里叶变换,也不存在 z z z 变换,但其功率有限,存在功率谱。
傅里叶变换条件(狄利克雷条件):
1.函数在任意有限区间内连续,或只有有限个第一类间断点(当t从左或右趋于这个间断点时,函数有有限的左极限和右极限);
2.在一个周期内,函数有有限个极大值或极小值。
3.x(t)在单个周期内绝对可积。
z z z 变换条件(充要条件):
级数绝对可和。使级数绝对可和的成立的所有 z z z 值称为收敛域。
如果 ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 满足绝对可和, ∑ m = − ∞ ∞ ∣ ϕ x x ( m ) ∣ ≤ ∞ \sum_{m=-\infty}^\infty |\phi_{xx}(m)|\leq \infty ∑m=−∞∞∣ϕxx(m)∣≤∞ ,则 ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 存在离散傅里叶变换,且均值 m x = 0 m_x=0 mx=0 。
1.3.1 自相关函数的 z z z 变换及其收敛域
设矩形窗口函数:
w
R
(
n
)
=
{
1
,
0
≤
n
≤
N
0
,
其
它
w_R(n)= \left\{ \begin{aligned} 1,\qquad &0\leq n\leq N\\ 0,\qquad&其它 \end{aligned} \right.
wR(n)={1,0,0≤n≤N其它
则,
x
N
(
n
)
=
x
(
n
)
w
R
(
n
)
x_N(n)=x(n)w_R(n)
xN(n)=x(n)wR(n)
傅里叶及其共轭的傅里叶变换:
X
N
(
ω
)
=
∑
n
=
0
N
−
1
x
(
n
)
e
−
j
ω
n
X
N
∗
(
ω
)
=
∑
k
=
0
N
−
1
x
∗
(
k
)
e
−
j
ω
k
X_N(\omega)=\sum_{n=0}^{N-1}x(n)e^{-j\omega n} \\ X_N^*(\omega)=\sum_{k=0}^{N-1}x^*(k)e^{-j\omega k}
XN(ω)=n=0∑N−1x(n)e−jωnXN∗(ω)=k=0∑N−1x∗(k)e−jωk
故
E
[
∣
X
N
(
ω
)
∣
2
]
=
∑
n
=
0
N
−
1
∑
k
=
0
N
−
1
E
[
x
(
n
)
x
∗
(
k
)
]
e
−
j
ω
(
n
−
k
)
=
∑
m
=
−
(
N
−
1
)
N
−
1
(
N
−
∣
m
∣
)
ϕ
x
x
(
m
)
e
−
j
ω
m
\begin{aligned} E[|X_N(\omega)|^2]&=\sum_{n=0}^{N-1}\sum_{k=0}^{N-1}E[x(n)x^*(k)]e^{-j\omega(n-k)}\\ &=\sum_{m=-(N-1)}^{N-1}(N-|m|)\phi_{xx}(m)e^{-j\omega m} \end{aligned}
E[∣XN(ω)∣2]=n=0∑N−1k=0∑N−1E[x(n)x∗(k)]e−jω(n−k)=m=−(N−1)∑N−1(N−∣m∣)ϕxx(m)e−jωm
由此
E
[
∣
1
N
X
N
(
ω
)
∣
2
]
=
∑
m
=
−
(
N
−
1
)
N
−
1
(
1
−
∣
m
∣
N
)
ϕ
x
x
(
m
)
e
−
j
ω
m
lim
N
→
∞
E
[
∣
1
N
X
N
(
ω
)
∣
2
]
=
∑
m
=
−
∞
∞
ϕ
x
x
(
m
)
e
−
j
ω
m
E[|\frac{1}{N}X_N(\omega)|^2]=\sum_{m=-(N-1)}^{N-1}(1-\frac{|m|}{N})\phi_{xx}(m)e^{-j\omega m}\\ \lim\limits_{N\rightarrow\infty}E[|\frac{1}{N}X_N(\omega)|^2]=\sum_{m=-\infty}^{\infty}\phi_{xx}(m)e^{-j\omega m}
E[∣N1XN(ω)∣2]=m=−(N−1)∑N−1(1−N∣m∣)ϕxx(m)e−jωmN→∞limE[∣N1XN(ω)∣2]=m=−∞∑∞ϕxx(m)e−jωm
这里
∣
1
N
X
N
(
ω
)
∣
2
|\frac{1}{N}X_N(\omega)|^2
∣N1XN(ω)∣2 叫做加窗序列周期图,
E
[
∣
1
N
X
N
(
ω
)
∣
2
]
E[|\frac{1}{N}X_N(\omega)|^2]
E[∣N1XN(ω)∣2] 平均功率谱密度。
即,功率谱密度和自相关函数互为傅里叶变换和逆变换
P
x
x
(
ω
)
=
∑
m
=
−
∞
∞
ϕ
x
x
(
m
)
e
−
j
ω
m
ϕ
x
x
(
m
)
=
1
2
π
∫
−
π
π
P
x
x
(
ω
)
e
j
ω
m
d
ω
P_{xx}(\omega)=\sum_{m=-\infty}^{\infty}\phi_{xx}(m)e^{-j\omega m}\\ \phi_{xx}(m)=\frac{1}{2\pi}\int_{-\pi}^{\pi}P_{xx}(\omega)e^{j\omega m}d\omega
Pxx(ω)=m=−∞∑∞ϕxx(m)e−jωmϕxx(m)=2π1∫−ππPxx(ω)ejωmdω
对于
z
z
z 变换,我们有
Γ
x
x
(
z
)
=
∑
m
=
−
∞
∞
γ
x
x
(
m
)
z
−
1
,
R
z
−
≤
∣
z
∣
≤
R
z
+
γ
x
x
(
m
)
=
1
2
π
j
∮
C
Γ
x
x
(
z
)
z
m
−
1
d
z
,
C
∈
(
R
z
−
,
R
z
+
)
\Gamma_{xx}(z)=\sum_{m=-\infty}^{\infty}\gamma_{xx}(m)z^{-1}, \quad R_{z^-}\leq |z| \leq R_{z^+}\\ \gamma_{xx}(m)=\frac{1}{2\pi j}\oint_{C}\Gamma_{xx}(z)z^{m-1}dz, \quad C\in (R_{z^-},R_{z^+})
Γxx(z)=m=−∞∑∞γxx(m)z−1,Rz−≤∣z∣≤Rz+γxx(m)=2πj1∮CΓxx(z)zm−1dz,C∈(Rz−,Rz+)
Φ x x ( z ) = ∑ m = − ∞ ∞ ϕ x x ( m ) z − 1 , R z − ≤ ∣ z ∣ ≤ R z + ϕ x x ( m ) = 1 2 π j ∮ C Φ x x ( z ) z m − 1 d z , C ∈ ( R z − , R z + ) \Phi_{xx}(z)=\sum_{m=-\infty}^{\infty}\phi_{xx}(m)z^{-1}, \quad R_{z^-}\leq |z| \leq R_{z^+}\\ \phi_{xx}(m)=\frac{1}{2\pi j}\oint_{C}\Phi_{xx}(z)z^{m-1}dz, \quad C\in (R_{z^-},R_{z^+}) Φxx(z)=m=−∞∑∞ϕxx(m)z−1,Rz−≤∣z∣≤Rz+ϕxx(m)=2πj1∮CΦxx(z)zm−1dz,C∈(Rz−,Rz+)
由于要使得随机信号存在
z
z
z 变换,故均值
m
x
=
0
m_x=0
mx=0 此时
γ
x
x
(
m
)
=
ϕ
x
x
(
m
)
\gamma_{xx}(m)=\phi_{xx}(m)
γxx(m)=ϕxx(m) .
x
(
n
)
⟷
X
(
z
)
x
(
−
n
)
⟷
X
(
z
−
1
)
x
∗
(
n
)
⟷
X
∗
(
z
∗
)
x(n)\longleftrightarrow X(z)\\ x(-n)\longleftrightarrow X(z^{-1})\\ x^*(n)\longleftrightarrow X^*(z^*)\\
x(n)⟷X(z)x(−n)⟷X(z−1)x∗(n)⟷X∗(z∗)
因为
ϕ
x
x
(
m
)
\phi_{xx}(m)
ϕxx(m) 是偶函数,所以
Φ
x
x
(
z
)
=
Φ
x
x
(
z
−
1
)
\Phi_{xx}(z)=\Phi_{xx}(z^{-1})
Φxx(z)=Φxx(z−1) 。如果
z
=
z
1
z=z_1
z=z1 是
Φ
x
x
(
z
)
\Phi_{xx}(z)
Φxx(z) 的一个极点,则
z
=
1
/
z
1
z=1/z_1
z=1/z1 也是
Φ
x
x
(
z
)
\Phi_{xx}(z)
Φxx(z) 的一个极点。由此可见
Φ
x
x
(
z
)
\Phi_{xx}(z)
Φxx(z) 和
Γ
x
x
(
z
)
\Gamma_{xx}(z)
Γxx(z) 的收敛域一定是包含单位圆的圆环。
1.3.2 功率谱分析
既然
Γ
x
x
(
z
)
\Gamma_{xx}(z)
Γxx(z) 的收敛域包含单位元,那么就可以取单位圆作积分围线,并将
z
=
e
j
ω
z=e^{j\omega}
z=ejω 代入得:
γ
x
x
(
m
)
=
1
2
π
j
∮
C
Γ
x
x
(
z
)
z
m
−
1
d
z
=
1
2
π
∫
−
π
π
Γ
x
x
(
e
j
ω
)
e
j
ω
m
d
ω
\begin{aligned} \gamma_{xx}(m)&=\frac{1}{2\pi j}\oint_{C}\Gamma_{xx}(z)z^{m-1}dz\\ &=\frac{1}{2\pi}\int_{-\pi}^{\pi}\Gamma_{xx}(e^{j\omega})e^{j\omega m}d\omega \end{aligned}
γxx(m)=2πj1∮CΓxx(z)zm−1dz=2π1∫−ππΓxx(ejω)ejωmdω
令
P
x
x
(
ω
)
=
Γ
x
x
(
e
j
ω
)
P_{xx}(\omega) = \Gamma_{xx}(e^{j\omega})
Pxx(ω)=Γxx(ejω) ,则有
γ
x
x
(
m
)
=
1
2
π
∫
−
π
π
P
x
x
(
ω
)
e
j
ω
m
d
ω
γ
x
x
(
0
)
=
E
[
x
2
(
n
)
]
=
1
2
π
∫
−
π
π
P
x
x
(
ω
)
d
ω
\gamma_{xx}(m)=\frac{1}{2\pi}\int_{-\pi}^{\pi}P_{xx}(\omega)e^{j\omega m}d\omega\\ \gamma_{xx}(0)=E[x^2(n)]=\frac{1}{2\pi}\int_{-\pi}^{\pi}P_{xx}(\omega)d\omega
γxx(m)=2π1∫−ππPxx(ω)ejωmdωγxx(0)=E[x2(n)]=2π1∫−ππPxx(ω)dω
此时正好反应的是
{
x
(
n
)
}
\{x(n)\}
{x(n)} 的平均功率密度,而且因为它是相对于频率的平均功率密度,所以又被称为功率谱密度。功率谱密度和自相关函数互为傅里叶变换和逆变换
P
x
x
(
ω
)
=
∑
m
=
−
∞
∞
ϕ
x
x
(
m
)
e
−
j
ω
m
ϕ
x
x
(
m
)
=
1
2
π
∫
−
π
π
P
x
x
(
ω
)
e
j
ω
m
d
ω
P_{xx}(\omega)=\sum_{m=-\infty}^{\infty}\phi_{xx}(m)e^{-j\omega m}\\ \phi_{xx}(m)=\frac{1}{2\pi}\int_{-\pi}^{\pi}P_{xx}(\omega)e^{j\omega m}d\omega
Pxx(ω)=m=−∞∑∞ϕxx(m)e−jωmϕxx(m)=2π1∫−ππPxx(ω)ejωmdω
1.3.3 谱密度的物理意义
在信号与系统中,若设
x
(
n
)
x(n)
x(n) 为一确定性的功率信号,则
x
(
n
)
x(n)
x(n) 在频率
ω
\omega
ω 处的功率谱密度为
P
x
x
(
ω
)
=
lim
N
→
∞
1
2
N
+
1
∣
∑
n
=
−
N
N
x
(
n
)
e
−
j
ω
m
∣
2
P_{xx}(\omega)=\lim\limits_{N\rightarrow\infty}\frac{1}{2N+1}|\sum_{n=-N}^N x(n)e^{-j\omega m}|^2
Pxx(ω)=N→∞lim2N+11∣n=−N∑Nx(n)e−jωm∣2
对于平稳随机序列,则有以下定理
设 ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 是平稳随机序列 { x ( n ) , n ∈ Z } \{x(n),n\in Z\} {x(n),n∈Z} 的自相关函数,如果 ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 绝对可和,则有
P x x ( ω ) = lim N → ∞ 1 2 N + 1 E [ ∣ ∑ n = − N N x ( n ) e − j ω m ∣ 2 ] P_{xx}(\omega)=\lim\limits_{N\rightarrow\infty}\frac{1}{2N+1}E[|\sum_{n=-N}^N x(n)e^{-j\omega m}|^2] Pxx(ω)=N→∞lim2N+11E[∣n=−N∑Nx(n)e−jωm∣2]
证明:
1 2 N + 1 E [ ∣ ∑ n = − N N x ( n ) e − j ω m ∣ 2 ] = 1 2 N + 1 E [ ∑ k = − N N x ∗ ( k ) e j ω k ∑ n = − N N x ( n ) e − j ω n ] = 1 2 N + 1 E [ ∑ k = − N N ∑ n = − N N x ∗ ( k ) x ( n ) e − j ω ( n − k ) ] = 1 2 N + 1 ∑ k = − N N ∑ n = − N N ϕ x x ( n − k ) e − j ω ( n − k ) = 1 2 ( 2 N + 1 ) ∑ l − m = − 2 N 2 N ∑ l + m = − 2 N 2 N ϕ x x ( m ) e − j ω m , m = n − k , l = n + k = 1 2 ( 2 N + 1 ) [ ∑ m = − 2 N 0 ∑ l = − 2 N − m 2 N + m ϕ x x ( m ) e − j ω m + ∑ m = 0 2 N ∑ l = − 2 N + m 2 N − m ϕ x x ( m ) e − j ω m ] = ∑ m = − 2 N 0 ( 1 + m 2 N + 1 ) ϕ x x ( m ) e − j ω m + ∑ m = 0 2 N ( 1 − m 2 N + 1 ) ϕ x x ( m ) e − j ω m = ∑ m = − 2 N 0 ( 1 − ∣ m ∣ 2 N + 1 ) ϕ x x ( m ) e − j ω m \begin{aligned} \frac{1}{2N+1}E[|\sum_{n=-N}^N x(n)e^{-j\omega m}|^2] &=\frac{1}{2N+1}E[\sum_{k=-N}^N x^*(k)e^{j\omega k}\sum_{n=-N}^N x(n)e^{-j\omega n}]\\ &=\frac{1}{2N+1}E[\sum_{k=-N}^N \sum_{n=-N}^N x^*(k)x(n)e^{-j\omega (n-k)}]\\ &=\frac{1}{2N+1}\sum_{k=-N}^N \sum_{n=-N}^N \phi_{xx}(n-k)e^{-j\omega (n-k)}\\ &=\frac{1}{2(2N+1)}\sum_{l-m=-2N}^{2N} \sum_{l+m=-2N}^{2N}\phi_{xx}(m)e^{-j\omega m},\quad m=n-k,l=n+k\\ &=\frac{1}{2(2N+1)}[\sum_{m=-2N}^{0} \sum_{l=-2N-m}^{2N+m}\phi_{xx}(m)e^{-j\omega m}\\ &+\sum_{m=0}^{2N} \sum_{l=-2N+m}^{2N-m}\phi_{xx}(m)e^{-j\omega m}]\\ &=\sum_{m=-2N}^{0}(1+\frac{m}{2N+1})\phi_{xx}(m)e^{-j\omega m}\\ &+\sum_{m=0}^{2N}(1-\frac{m}{2N+1})\phi_{xx}(m)e^{-j\omega m}\\ &=\sum_{m=-2N}^{0}(1-\frac{|m|}{2N+1})\phi_{xx}(m)e^{-j\omega m}\\ \end{aligned} 2N+11E[∣n=−N∑Nx(n)e−jωm∣2]=2N+11E[k=−N∑Nx∗(k)ejωkn=−N∑Nx(n)e−jωn]=2N+11E[k=−N∑Nn=−N∑Nx∗(k)x(n)e−jω(n−k)]=2N+11k=−N∑Nn=−N∑Nϕxx(n−k)e−jω(n−k)=2(2N+1)1l−m=−2N∑2Nl+m=−2N∑2Nϕxx(m)e−jωm,m=n−k,l=n+k=2(2N+1)1[m=−2N∑0l=−2N−m∑2N+mϕxx(m)e−jωm+m=0∑2Nl=−2N+m∑2N−mϕxx(m)e−jωm]=m=−2N∑0(1+2N+1m)ϕxx(m)e−jωm+m=0∑2N(1−2N+1m)ϕxx(m)e−jωm=m=−2N∑0(1−2N+1∣m∣)ϕxx(m)e−jωm
故
P x x ( ω ) = lim N → ∞ 1 2 N + 1 E [ ∣ ∑ n = − N N x ( n ) e − j ω m ∣ 2 ] = ∑ m = − ∞ ∞ ϕ x x ( m ) e − j ω m P_{xx}(\omega)=\lim\limits_{N\rightarrow\infty}\frac{1}{2N+1}E[|\sum_{n=-N}^N x(n)e^{-j\omega m}|^2]=\sum_{m=-\infty}^{\infty}\phi_{xx}(m)e^{-j\omega m} Pxx(ω)=N→∞lim2N+11E[∣n=−N∑Nx(n)e−jωm∣2]=m=−∞∑∞ϕxx(m)e−jωm
综合上面的讨论,可知随机序列的功率谱密度具有以下几个特点:
1)随机序列的自相关函数与它的功率谱密度构成一对离散时间傅里叶变换。
2)随机序列的自相关函数具有离散时间的形式,也称自相关序列,但它的功率谱密度却是变量 ω \omega ω 的连续函数。
3)对于非周期离散时间自相关序列,其功率谱密度为一周期信号,周期为 2 π 2\pi 2π ,功率谱密度在一个周期内的积分与随机序列的平均功率成正比。
4)随机序列的功率谱密度是周期信号,因此可以展开为傅里叶级数, ϕ x x ( m ) \phi_{xx}(m) ϕxx(m) 正是傅里叶级数的各谐波分量系数。
5)随机序列的自相关函数与功率谱密度分别从时域、频域反映了随机序列的二阶统计特性,虽然形式不同,但两者包含的信息是一致的。
1.4 线性系统对随机信号的影响
所谓系统,用数学的语言描述,就是从输入序列
{
x
(
n
)
,
n
∈
Z
}
\{x(n),n \in Z\}
{x(n),n∈Z} 到输出序列
{
y
(
n
)
,
n
∈
Z
}
\{y(n),n \in Z\}
{y(n),n∈Z} 的映射,记为
{
y
(
n
)
,
n
∈
Z
}
=
L
{
x
(
n
)
,
n
∈
Z
}
\{y(n),n \in Z\}=L\{x(n),n \in Z\}
{y(n),n∈Z}=L{x(n),n∈Z} 。如果输入序列与输出序列满足:
y
(
n
)
=
∑
k
=
−
∞
∞
h
(
k
)
x
(
n
−
k
)
,
n
∈
Z
y(n)=\sum_{k=-\infty}^{\infty}h(k)x(n-k),\quad n\in Z
y(n)=k=−∞∑∞h(k)x(n−k),n∈Z
则称
L
L
L 为线性时不变系统,称
{
h
(
k
)
,
k
∈
Z
}
\{h(k),k \in Z\}
{h(k),k∈Z} 为线性系统的冲激响应,它是线性系统的时域表征,如果
k
<
0
k<0
k<0 时
h
(
k
)
=
0
h(k)=0
h(k)=0 ,则
L
L
L 是因果的。此时
y
(
n
)
=
∑
k
=
0
∞
h
(
k
)
x
(
n
−
k
)
,
n
∈
Z
y(n)=\sum_{k=0}^{\infty}h(k)x(n-k),\quad n\in Z
y(n)=k=0∑∞h(k)x(n−k),n∈Z
若冲激响应
{
h
(
k
)
,
k
∈
Z
}
\{h(k),k \in Z\}
{h(k),k∈Z} 满足条件
∑
k
=
−
∞
∞
∣
h
(
k
)
∣
<
∞
\sum_{k=-\infty}^{\infty}|h(k)|<\infty
k=−∞∑∞∣h(k)∣<∞
就称
L
L
L 是稳定的,稳定意味着有界的输入导致了有界的输出,在实践中,常见的系统都是稳定和因果的线性系统。
1.4.1 线性时不变系统对随机输入的响应
L L L 是线性时不变系统,其冲激响应为 h ( n ) h(n) h(n) ,频率响应为 H ( ω ) H(\omega) H(ω) ,且满足
1 ) ∑ k = − ∞ ∞ ∣ h ( k ) ∣ < ∞ 2 ) ∑ k = − ∞ ∞ ∑ n = − ∞ ∞ h ∗ ( k ) h ( n ) ϕ x x ( n − k ) < ∞ \begin{aligned} &1)\sum_{k=-\infty}^{\infty}|h(k)|<\infty \\ &2)\sum_{k=-\infty}^{\infty}\sum_{n=-\infty}^{\infty}h^*(k)h(n)\phi_{xx}(n-k)<\infty \end{aligned} 1)k=−∞∑∞∣h(k)∣<∞2)k=−∞∑∞n=−∞∑∞h∗(k)h(n)ϕxx(n−k)<∞
则有:
y
(
n
)
=
∑
k
=
−
∞
∞
h
(
k
)
x
(
n
−
k
)
=
h
(
n
)
∗
x
(
n
)
y(n)=\sum_{k=-\infty}^{\infty}h(k)x(n-k)=h(n)*x(n)
y(n)=k=−∞∑∞h(k)x(n−k)=h(n)∗x(n)
也是平稳过程。
{
y
(
n
)
}
\{y(n)\}
{y(n)} 的均值、自相关函数和功率谱密度
m
y
=
m
x
∑
k
=
−
∞
∞
h
(
k
)
ϕ
y
y
(
n
,
n
+
m
)
=
∑
k
=
−
∞
∞
h
(
k
)
∑
r
=
−
∞
∞
h
(
r
)
ϕ
x
x
(
m
+
k
−
r
)
=
ϕ
y
y
(
m
)
P
y
y
(
ω
)
=
∣
H
(
ω
)
∣
2
P
x
x
(
ω
)
\begin{aligned} &m_y=m_x\sum_{k=-\infty}^{\infty}h(k)\\ &\phi_{yy}(n, n+m)=\sum_{k=-\infty}^{\infty}h(k)\sum_{r=-\infty}^{\infty}h(r)\phi_{xx}(m+k-r)=\phi_{yy}(m)\\ &P_{yy}(\omega)=|H(\omega)|^2P_{xx}(\omega)\\ \end{aligned}
my=mxk=−∞∑∞h(k)ϕyy(n,n+m)=k=−∞∑∞h(k)r=−∞∑∞h(r)ϕxx(m+k−r)=ϕyy(m)Pyy(ω)=∣H(ω)∣2Pxx(ω)
1.4.2 系统输入、输出的互相关函数和互谱密度
现在让我们来讨论关于线性非时变系统的输入和输出之间的互相关函数 ϕ x y ( m ) \phi_{xy}(m) ϕxy(m) 及互相关谱密度函数 P x y ( ω ) P_{xy}(\omega) Pxy(ω) 。
设 线性时不变系统
L
L
L 的输入和输出分别为平稳随机序列
{
x
(
n
)
,
n
∈
Z
}
\{x(n),n \in Z\}
{x(n),n∈Z} 和
{
y
(
n
)
,
n
∈
Z
}
\{y(n),n \in Z\}
{y(n),n∈Z} 且
{
x
(
n
)
}
\{x(n)\}
{x(n)} 存在谱密度
P
x
x
(
ω
)
P_{xx}(\omega)
Pxx(ω) ,则系统的输入
{
x
(
n
)
}
\{x(n)\}
{x(n)} 和输出
{
y
(
n
)
}
\{y(n)\}
{y(n)} 平稳相关,且它们的互相关谱密度函数为
P
x
y
(
ω
)
=
H
(
ω
)
P
x
x
(
ω
)
,
P
y
x
=
H
∗
(
ω
)
P
x
x
(
ω
)
P_{xy}(\omega)=H(\omega)P_{xx}(\omega),P_{yx}=H^*(\omega)P_{xx}(\omega)
Pxy(ω)=H(ω)Pxx(ω),Pyx=H∗(ω)Pxx(ω)
ϕ x y ( m ) = ϕ x x ( m ) ∗ h ( m ) , ϕ y y ( m ) = ϕ x y ( m ) ∗ h ( − m ) \phi_{xy}(m)=\phi_{xx}(m)*h(m),\phi_{yy}(m)=\phi_{xy}(m)*h(-m) ϕxy(m)=ϕxx(m)∗h(m),ϕyy(m)=ϕxy(m)∗h(−m)