第一章 离散随机信号 笔记
第一章 离散随机信号
1.1 引言
定义 :设已给概率空间 (Ω,ζ,P)(\Omega,\zeta,P)(Ω,ζ,P) ,ZZZ 为整数集,若对每一整数 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,ξ)(ξ∈Ω) 与之对应,则称依赖与参数 nnn 的一列随机变量 x(n,ξ)x(n,\xi)x(n,ξ) 为一离散时间随机过程或随机序列,记为 {x(n,ξ),ξ∈Ω,n∈Z}\{x(n,\xi),\xi \in\Omega,n\in Z\}{x(n,ξ),ξ∈Ω,n∈Z} ,简记为 {xn}\{x_n\}{xn} 。
随机序列有一下特点:
(1)随机序列中的任何一个点上的取值都是不能先验确定的随机变量。
(2)随机序列可以用它的统计平均特性来表征。随机序列在个时间点上的随机变量取值是服从某种确定的概率分布。
(3)平稳随机信号的能量化表示。一随机信号的能量谱称为功率谱密度。一个平稳随机信号的功率谱是确定的,因此功率谱可以统计表征一个随机过程的谱特性。
综上所述,对于离散随机信号的概念和表征问题,有:
(1)一个随机信号在各个时间点上的取值以及不同点上取值之间的相互关联性只能用概率特性或统计平均特性来表征,它的确定值是无法先验表达的。
(2)一个平稳随机信号在各频率点上的能量的取值可以用功率谱,密度函数与自相关函数统计描述。
1.2 离散时间随机信号的时域(统计)表示
一般情况下,我们大多只对离散随机信号的一维概率分布和二维概率分布等有限维度概率分布进行研究。
1.2.1 离散时间随机过程的概率分布
(宽)平稳随机过程:若一个随机过程的数学期望及方差与时间无关,相关函数仅与时间间隔有关,则称这个随机过程为广义平稳随机过程。
(宽)平稳随机过程的三个特征:
(1) E(xn)=CξE(x_n) = C_\xiE(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[x2(n)]<+∞E[x^2(n)] < +\inftyE[x2(n)]<+∞ 均方值有界。
严平稳随机过程:所谓随机过程严格平稳,是指它的任何n维分布函数或概率密度函数与时间起点无关。
严平稳随机过程特点:其数学期望、方差与时间无关,自相关函数仅与时间间隔有关。
- (宽)平稳随机过程 和 严平稳随机过程 的区别和联系
(1)一个宽平稳过程不一定是严平稳过程,一个严平稳过程也不一定宽平稳过程。
(2)宽平稳过程定只涉及与一维、二维分布有关的数字特征,所以一个严平稳过程只要二阶矩存在,则必定是宽平稳过程。但反过来,一般是不成立的。
(3)正态过程是一个重要特例,一个宽平稳的正态过程必定是严平稳的。
各态历经随机过程:如果一个平稳随机过程,只要满足一些较宽的条件,则一个样本函数在整个时间轴上的平均值可以用来代替其集平均(统计平均值和自相关函数等),这就是各态历经随机过程。
一般来说,在一个随机过程中,不同样本函数的时间平均值是不一定相同的,而集平均则是一定的。因此,一般的随机过程的时间平均≠集平均,只有平稳随机过程才有可能是具有各态历经性的。即各态历经的随机过程一定是平稳的,而平稳的随机过程则需要满足一定条件才是各态历经的。
1.2.2 离散时间随机过程的数字特征
- 数学期望(均值) 在n点上的
随机变量 x(n)x(n)x(n) 的均值(用 mx(n)m_x(n)mx(n) 表示)定义为
mx(n)=E[x(n)]≜∫−∞∞xp1(x)dxm_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∣2p1(x)dxE[|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) 的方差(用 σx2(n)\sigma_x^2(n)σx2(n) 表示)定义为 [x(n)−E(x(n))][x(n)-E(x(n))][x(n)−E(x(n))] 的均方值
σx2(n)≜E{∣x(n)−E(x(n))∣2}≜E{∣x(n)−mx(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}
通过以上公式可以推导出:
平均功率 = 交流功率 + 直流功率
σx2(n)=E[∣x(n)∣2]−∣mx2(n)∣E[∣x(n)∣2]=σx2(n)+∣mx2(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)∣}
以上三个特征量仅与一维概率密度 p1(x)p_1(x)p1(x) 有关。对于平稳随机过程,其一维概率密度与时间无关,故一个平稳随机序列的 mx(n)、E(∣x(n)∣2)、σx2(n)m_x(n)、E(|x(n)|^2)、\sigma_x^2(n)mx(n)、E(∣x(n)∣2)、σx2(n) 均是与时间无关的常数,可以将变量 nnn 去掉。
与二维概率分布有关的统计特性主要有自相关函数和自协方差函数。
- 自相关函数
一个平稳随机过程的自相关函数用 ϕxx(m)\phi_{xx}(m)ϕxx(m) 表示
ϕxx(m)≜E[x(n1)x∗(n2)]=E[x(n)x∗(n−m)]=∫−∞∞∫−∞∞x1x2∗p2(x1,x2;m)dx1dx2\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
- 自协方差函数
一个平稳随机序列的协方差函数用 γxx(m)\gamma_{xx}(m)γxx(m) 表示
γxx(m)≜E[(x(n)−mx)(x(n+m)−mx)∗]\gamma_{xx}(m)\triangleq E[(x(n)-m_x)(x(n+m)-m_x)^*] γxx(m)≜E[(x(n)−mx)(x(n+m)−mx)∗]
自相关函数 与 自协方差函数 的关系
γxx(m)=ϕxx(m)−∣mx∣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)} 的自相关、自协方差及两个序列的互相关和互协方差序列定义为
ϕxx(m)=E[x(n)x(n−m)]γxx(m)=E[(x(n)−mx)(x(n−m)−mx)]\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)]
ϕxy(m)=E[x(n)y(n−m)]γxy(m)=E[(x(n)−mx)yx(n−m)−my)]\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
γxx(m)=ϕxx(m)−mx2γxy(m)=ϕxy(m)−mxmy\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
ϕxx(0)=E[x2(n)]均方值γxx(0)=ϕxx(0)−mx2方差\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 ϕxx(m)\phi_{xx}(m)ϕxx(m) 和 γxx(m)\gamma_{xx}(m)γxx(m) 都是偶函数
ϕxx(m)=ϕxx(−m)γxx(m)=γxx(−m)\phi_{xx}(m)=\phi_{xx}(-m)\\ \gamma_{xx}(m)=\gamma_{xx}(-m) ϕxx(m)=ϕxx(−m)γxx(m)=γxx(−m)
性质 4
ϕxy(m)=ϕyx(−m)γxy(m)=γyx(−m)\phi_{xy}(m)=\phi_{yx}(-m)\\ \gamma_{xy}(m)=\gamma_{yx}(-m) ϕxy(m)=ϕyx(−m)γxy(m)=γyx(−m)
性质 5
∣ϕxx(m)∣≤ϕxx(0)∣γxx(m)∣≤γxx(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) ,则
ϕyy(m)=ϕxx(m)γyy(m)=γxx(m)\phi_{yy}(m)=\phi_{xx}(m)\\ \gamma_{yy}(m)=\gamma_{xx}(m) ϕyy(m)=ϕxx(m)γyy(m)=γxx(m)
性质 7 对于大多随机序列,mmm 愈大则相关性愈小,当 mmm 趋近无穷大时,可以认为不相关。若平稳随机序列 {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) γxx(m)\gamma_{xx}(m)γxx(m) 和 ϕxx(m)\phi_{xx}(m)ϕxx(m) 除相差一个常量 mx2m_x^2mx2 外,其它特性相同。
(2) ϕxx(m)\phi_{xx}(m)ϕxx(m) 是一个随机过程的最主要的统计特征,它不仅说明了相关性还隐含了 mx,σx2,E[x2(n)]m_x,\sigma_x^2,E[x^2(n)]mx,σx2,E[x2(n)] 等主要特征量,例:
E[x2(n)]=ϕxx(0)mx2=ϕxx(∞)σx2=E[x2(n)]−mx2=ϕxx(0)−ϕxx(∞)}\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=[x1x2⋮xn]\pmb{X}= \left[ \begin{matrix} x_1\\ x_2\\ \vdots\\ x_n \end{matrix} \right] XXX=⎣⎢⎢⎢⎡x1x2⋮xn⎦⎥⎥⎥⎤
其均值向量为
E[X]=[E[x1]E[x2]⋮E[xn]]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]⎦⎥⎥⎥⎤
其自协方差阵定义为
cov[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}
其互协方差阵定义为
cov(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…
在实践中需要多次大量的实验以获得足够多的( NNN 足够大)序列来保证计算的可靠性,显然这是十分困难的。但是我们假设一个平稳过程,集合平均等于一个样本函数在整个时间轴上的时间平均,这样就能使用时间平均对随机序列的均值和自相关函数做出合理近似估计。
定义 设 {x(n)}\{x(n)\}{x(n)} 是一个平稳随机序列,如果均值极限
<x(n)>≜limN→∞12N+1∑n=−NNx(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 有
mx=<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)>≜limN→∞12N+1∑n=−NN[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 有
ϕxx(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)成立,因为均方收敛必依概率收敛,那么当 NNN 足够大时,就可以用一次抽样所得的一个现实 {xn(ξ),n∈Z}(ξ∈Ω)\{x_n(\xi),n\in Z\}(\xi\in\Omega){xn(ξ),n∈Z}(ξ∈Ω) 所作的时间平均
12N+1∑n=−NNx(n,ξ)\frac{1}{2N+1}\sum_{n=-N}^{N}x(n,\xi) 2N+11n=−N∑Nx(n,ξ)
与
12N+1∑n=−NN[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,ξ)]
分别作为 mxm_xmx 和 ϕxx(m)\phi_{xx}(m)ϕxx(m) 的近似估计。
思考 :证明,∣ϕxy∣≤ϕxx(0)ϕyy(0)|\phi_{xy}|\leq \sqrt{\phi_{xx}(0)\phi_{yy}(0)}∣ϕxy∣≤ϕxx(0)ϕyy(0) .
构造:
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)\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 离散时间随机信号的 zzz 域及频域(统计)表示
平稳随机信号无始无终,能量无限,且不满足绝对可和(绝对值的和是一个有限常数)。
故不存在傅里叶变换,也不存在 zzz 变换,但其功率有限,存在功率谱。
傅里叶变换条件(狄利克雷条件):
1.函数在任意有限区间内连续,或只有有限个第一类间断点(当t从左或右趋于这个间断点时,函数有有限的左极限和右极限);
2.在一个周期内,函数有有限个极大值或极小值。
3.x(t)在单个周期内绝对可积。
zzz 变换条件(充要条件):
级数绝对可和。使级数绝对可和的成立的所有 zzz 值称为收敛域。
如果 ϕxx(m)\phi_{xx}(m)ϕxx(m) 满足绝对可和,∑m=−∞∞∣ϕxx(m)∣≤∞\sum_{m=-\infty}^\infty |\phi_{xx}(m)|\leq \infty∑m=−∞∞∣ϕxx(m)∣≤∞ ,则 ϕxx(m)\phi_{xx}(m)ϕxx(m) 存在离散傅里叶变换,且均值 mx=0m_x=0mx=0 。
1.3.1 自相关函数的 zzz 变换及其收敛域
设矩形窗口函数:
wR(n)={1,0≤n≤N0,其它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其它
则,
xN(n)=x(n)wR(n)x_N(n)=x(n)w_R(n) xN(n)=x(n)wR(n)
傅里叶及其共轭的傅里叶变换:
XN(ω)=∑n=0N−1x(n)e−jωnXN∗(ω)=∑k=0N−1x∗(k)e−jωkX_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[∣XN(ω)∣2]=∑n=0N−1∑k=0N−1E[x(n)x∗(k)]e−jω(n−k)=∑m=−(N−1)N−1(N−∣m∣)ϕxx(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[∣1NXN(ω)∣2]=∑m=−(N−1)N−1(1−∣m∣N)ϕxx(m)e−jωmlimN→∞E[∣1NXN(ω)∣2]=∑m=−∞∞ϕxx(m)e−jωmE[|\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
这里 ∣1NXN(ω)∣2|\frac{1}{N}X_N(\omega)|^2∣N1XN(ω)∣2 叫做加窗序列周期图, E[∣1NXN(ω)∣2]E[|\frac{1}{N}X_N(\omega)|^2]E[∣N1XN(ω)∣2] 平均功率谱密度。
即,功率谱密度和自相关函数互为傅里叶变换和逆变换
Pxx(ω)=∑m=−∞∞ϕxx(m)e−jωmϕxx(m)=12π∫−ππPxx(ω)ejωmdω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ω
对于 zzz 变换,我们有
Γxx(z)=∑m=−∞∞γxx(m)z−1,Rz−≤∣z∣≤Rz+γxx(m)=12πj∮CΓxx(z)zm−1dz,C∈(Rz−,Rz+)\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+)
Φxx(z)=∑m=−∞∞ϕxx(m)z−1,Rz−≤∣z∣≤Rz+ϕxx(m)=12πj∮CΦxx(z)zm−1dz,C∈(Rz−,Rz+)\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+)
由于要使得随机信号存在 zzz 变换,故均值 mx=0m_x=0mx=0 此时 γxx(m)=ϕxx(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∗)
因为 ϕxx(m)\phi_{xx}(m)ϕxx(m) 是偶函数,所以 Φxx(z)=Φxx(z−1)\Phi_{xx}(z)=\Phi_{xx}(z^{-1})Φxx(z)=Φxx(z−1) 。如果 z=z1z=z_1z=z1 是 Φxx(z)\Phi_{xx}(z)Φxx(z) 的一个极点,则 z=1/z1z=1/z_1z=1/z1 也是 Φxx(z)\Phi_{xx}(z)Φxx(z) 的一个极点。由此可见 Φxx(z)\Phi_{xx}(z)Φxx(z) 和 Γxx(z)\Gamma_{xx}(z)Γxx(z) 的收敛域一定是包含单位圆的圆环。
1.3.2 功率谱分析
既然 Γxx(z)\Gamma_{xx}(z)Γxx(z) 的收敛域包含单位元,那么就可以取单位圆作积分围线,并将 z=ejωz=e^{j\omega}z=ejω 代入得:
γxx(m)=12πj∮CΓxx(z)zm−1dz=12π∫−ππΓxx(ejω)ejωmdω\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ω
令 Pxx(ω)=Γxx(ejω)P_{xx}(\omega) = \Gamma_{xx}(e^{j\omega})Pxx(ω)=Γxx(ejω) ,则有
γxx(m)=12π∫−ππPxx(ω)ejωmdωγxx(0)=E[x2(n)]=12π∫−ππPxx(ω)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)} 的平均功率密度,而且因为它是相对于频率的平均功率密度,所以又被称为功率谱密度。功率谱密度和自相关函数互为傅里叶变换和逆变换
Pxx(ω)=∑m=−∞∞ϕxx(m)e−jωmϕxx(m)=12π∫−ππPxx(ω)ejωmdω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ω 处的功率谱密度为
Pxx(ω)=limN→∞12N+1∣∑n=−NNx(n)e−jωm∣2P_{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
对于平稳随机序列,则有以下定理
设 ϕxx(m)\phi_{xx}(m)ϕxx(m) 是平稳随机序列 {x(n),n∈Z}\{x(n),n\in Z\}{x(n),n∈Z} 的自相关函数,如果 ϕxx(m)\phi_{xx}(m)ϕxx(m) 绝对可和,则有
Pxx(ω)=limN→∞12N+1E[∣∑n=−NNx(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]
证明:
12N+1E[∣∑n=−NNx(n)e−jωm∣2]=12N+1E[∑k=−NNx∗(k)ejωk∑n=−NNx(n)e−jωn]=12N+1E[∑k=−NN∑n=−NNx∗(k)x(n)e−jω(n−k)]=12N+1∑k=−NN∑n=−NNϕxx(n−k)e−jω(n−k)=12(2N+1)∑l−m=−2N2N∑l+m=−2N2Nϕxx(m)e−jωm,m=n−k,l=n+k=12(2N+1)[∑m=−2N0∑l=−2N−m2N+mϕxx(m)e−jωm+∑m=02N∑l=−2N+m2N−mϕxx(m)e−jωm]=∑m=−2N0(1+m2N+1)ϕxx(m)e−jωm+∑m=02N(1−m2N+1)ϕxx(m)e−jωm=∑m=−2N0(1−∣m∣2N+1)ϕxx(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
故
Pxx(ω)=limN→∞12N+1E[∣∑n=−NNx(n)e−jωm∣2]=∑m=−∞∞ϕxx(m)e−jωmP_{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\pi2π ,功率谱密度在一个周期内的积分与随机序列的平均功率成正比。
4)随机序列的功率谱密度是周期信号,因此可以展开为傅里叶级数,ϕxx(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∈Zy(n)=\sum_{k=-\infty}^{\infty}h(k)x(n-k),\quad n\in Z y(n)=k=−∞∑∞h(k)x(n−k),n∈Z
则称 LLL 为线性时不变系统,称 {h(k),k∈Z}\{h(k),k \in Z\}{h(k),k∈Z} 为线性系统的冲激响应,它是线性系统的时域表征,如果 k<0k<0k<0 时 h(k)=0h(k)=0h(k)=0 ,则 LLL 是因果的。此时
y(n)=∑k=0∞h(k)x(n−k),n∈Zy(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)∣<∞
就称 LLL 是稳定的,稳定意味着有界的输入导致了有界的输出,在实践中,常见的系统都是稳定和因果的线性系统。
1.4.1 线性时不变系统对随机输入的响应
LLL 是线性时不变系统,其冲激响应为 h(n)h(n)h(n) ,频率响应为 H(ω)H(\omega)H(ω) ,且满足
1)∑k=−∞∞∣h(k)∣<∞2)∑k=−∞∞∑n=−∞∞h∗(k)h(n)ϕxx(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)} 的均值、自相关函数和功率谱密度
my=mx∑k=−∞∞h(k)ϕyy(n,n+m)=∑k=−∞∞h(k)∑r=−∞∞h(r)ϕxx(m+k−r)=ϕyy(m)Pyy(ω)=∣H(ω)∣2Pxx(ω)\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 系统输入、输出的互相关函数和互谱密度
现在让我们来讨论关于线性非时变系统的输入和输出之间的互相关函数 ϕxy(m)\phi_{xy}(m)ϕxy(m) 及互相关谱密度函数 Pxy(ω)P_{xy}(\omega)Pxy(ω) 。
设 线性时不变系统 LLL 的输入和输出分别为平稳随机序列 {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)} 存在谱密度 Pxx(ω)P_{xx}(\omega)Pxx(ω) ,则系统的输入 {x(n)}\{x(n)\}{x(n)} 和输出 {y(n)}\{y(n)\}{y(n)} 平稳相关,且它们的互相关谱密度函数为
Pxy(ω)=H(ω)Pxx(ω),Pyx=H∗(ω)Pxx(ω)P_{xy}(\omega)=H(\omega)P_{xx}(\omega),P_{yx}=H^*(\omega)P_{xx}(\omega) Pxy(ω)=H(ω)Pxx(ω),Pyx=H∗(ω)Pxx(ω)
ϕxy(m)=ϕxx(m)∗h(m),ϕyy(m)=ϕxy(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)
第一章 离散随机信号 笔记相关推荐
- 【食品加工技术】第一章 食品加工技术概述 笔记
[食品加工技术]第一章 食品加工技术概述 笔记 1.1 食品行业常用术语 食品 Food 传统食品 Traditional Food 天然食品 Natural Food 模拟食品 Imitation ...
- 数字信号处理第一章:离散时间信号与离散时间系统
参考书:<数字信号处理-理论. 算法与实现>第二版 胡广书 清华大学出版社 1. 离散时间信号的基本概念 1) 离散信号概述 a)连续时间信号x(t) -->离散时间信号 x(nTs ...
- 工程经济学第一章工程经济学概述笔记
工程经济学第一章工程经济学概述 1.1.1 工程与经济的概念 1. 工程的概念 工程不同于科学,也不同于技术. 科学本质上属于认识世界的范畴.科学是人类探索自然和社会现象并取得认识的过程和结果.这 里 ...
- 电磁场第一章——矢量分析工具 复习笔记
笔记不是用来初学,而是用来过知识点的,比较简洁,提纲挈领,细节不太多,更不涉及计算.后面章节会陆续更. 场:指物体在空间中的分布情况.在空间区域中某一点有确定物理量与之对应,称在该区域定义了该物理量的 ...
- 深度学习入门基于Python的理论与实现_第一章_Python入门(原创笔记)
前言 此书使用Python作为编程语言,尽可能地少使用外部库,从零开始实现深度学习的程序. 此书从简单的机器学习问题开始,最终实现一个能高精度地识别图像的系统. 此书以图像识别为主题,主要学习使用深度 ...
- 【转】数据库软考笔记(一)第一章 计算机硬件基础知识笔记总结
转自:https://blog.csdn.net/qq_36411874/article/details/115057569 目录 运算器:ALU.AC.DR.PSW 控制器:IR.PC.AR.ID ...
- 第一章 C++初步知识笔记(上)
1.1 从C到C++ C++对C 的增强: 在C语言的基础上对C语言做了扩充 在面向对象的基础上增加了面向对象的支持 1.2 最简单的C++程序 #include <iostream> / ...
- 《机器学习》 —— 第一章:绪论 学习笔记
1.1 引言 机器学习所研究的主要内容,是关于在计算机上从数据中产生"模型"(model)的算法,即"学习算法"(learning algorithm). 1. ...
- 《DirectX 9.0 3D游戏开发编程基础》 第一章 初始化Direct3D 读书笔记
REF设备 参考光栅设备,他能以软件计算方式完全支持Direct3D Api.借助Ref设备,可以在代码中使用那些不为当前硬件所支持的特性,并对这此特性进行测试. D3DDEVTYPE 在程序代码中, ...
最新文章
- consul agent的一个故障
- hdu 4533(树状数组区间更新+单点查询)
- IHttpModule 与IHttpHandler的区别
- C#创建桌面快捷方式
- python正态分布函数_python3-正态分布
- hystrix基于request cache请求缓存技术优化批量数据查询接口
- win10找回win7的windows照片查看器
- 电脑每次开机都要重新设置时间解决方法
- oracle的创建视图,oracle 创建视图
- MATLAB显示图像为什么一片白
- 图及其应用实验报告 c语言,图的应用 实验报告
- html5自动换行,html5文字自动换行怎么设置
- 系统分析与设计小组项目总结报告
- vivo系统升级服务器无响应,vivo系统升级教程
- 游戏时弹出内存不能为read的提示框
- MySQL之库表设计篇:一到五范式、BC范式与反范式详解
- 屏幕的背光是什么呢?
- 用google突破图书馆入口IP限制小技巧
- 爲什麽一个标准的反相器中PMOS管的寬長比要比N管大
- 联邦学习实战-1:用python从零开始实现横向联邦学习