通信原理随机信号分析
通信原理第二章 隨機(jī)信號(hào)分析
一 隨機(jī)過(guò)程
定義
測(cè)試n臺(tái)性能相同的接收機(jī),在同樣條件下,不加信號(hào)測(cè)試其輸出噪聲,波形如圖
(1)每一條曲線 ξi(t)\xi_i(t)ξi?(t) 都是一個(gè)隨機(jī)起伏的時(shí)間函數(shù)——樣本函數(shù)(確知信號(hào)),ξi(t)\xi_i(t)ξi?(t) 稱之為隨機(jī)過(guò)程ξ(t)\xi(t)ξ(t)的一個(gè)實(shí)現(xiàn)/樣本。這是對(duì)于其中一臺(tái)接收機(jī)觀察。
(2)全體樣本函數(shù)的集合稱作隨機(jī)過(guò)程 ξ(t)\xi(t)ξ(t)。
ξ(t)={ξ1(t),ξ2(t),...,ξn(t)}\xi(t)=\{ \xi_1(t),\xi_2(t),...,\xi_n(t)\}ξ(t)={ξ1?(t),ξ2?(t),...,ξn?(t)}
(3)在某一特定時(shí)刻 t1t_1t1?觀察各臺(tái)接收機(jī)的輸出噪聲值 ξ(t1)\xi(t_1)ξ(t1?) ,此時(shí)所有的輸出噪聲值是隨機(jī)過(guò)程ξ(t)\xi(t)ξ(t)一個(gè)隨機(jī)量(隨機(jī)變量):
ξ(t1)={ξ1(t1),ξ2(t1),...,ξn(t1)}\xi(t_1)=\{\xi_1(t_1),\xi_2(t_1),...,\xi_n(t_1)\}ξ(t1?)={ξ1?(t1?),ξ2?(t1?),...,ξn?(t1?)}
因此,隨機(jī)過(guò)程 ξ( t) 是由無(wú)窮多個(gè)隨機(jī)變量構(gòu)成的。
ξ(t)={ξ(t1)+ξ(t2)+...+ξ(tn)+...}\xi(t)=\{\xi(t_1)+\xi(t_2)+...+\xi(t_n)+...\}ξ(t)={ξ(t1?)+ξ(t2?)+...+ξ(tn?)+...}
數(shù)字特征
| 均值 | a(t)=E[ξ(t)]=∫?∞+∞xf1(x,t)dxa(t)=E[\xi(t)]=\int_{-\infty}^{+\infty}xf_1(x,t)\,dxa(t)=E[ξ(t)]=∫?∞+∞?xf1?(x,t)dx | 隨機(jī)過(guò)程的擺動(dòng)中心; 均值的平方是直流功率 |
| 均方值 | E[ξ2(t)]=∫?∞+∞x2f1(x,t)dxE[{\xi^2(t)}]=\int_{-\infty}^{+\infty}x^2f_1(x,t)\,dxE[ξ2(t)]=∫?∞+∞?x2f1?(x,t)dx | 隨機(jī)過(guò)程的平均功率 |
| 方差 | D[ξ(t)]=E{[ξ(t)?a(t)]2}=E[ξ2(t)]?a2(t)=σ2(t)D[\xi(t)]=E\{[\xi(t)-a(t)]^2\}=E[\xi^2(t)]-a^2(t)=\sigma^2(t)D[ξ(t)]=E{[ξ(t)?a(t)]2}=E[ξ2(t)]?a2(t)=σ2(t) | 隨機(jī)過(guò)程的交流功率,相對(duì)于均值的振動(dòng)程度 |
| 自相關(guān)函數(shù) | R(t1,t2)=E[ξ(t1)ξ(t2)]=∫?∞+∞∫?∞+∞x1x2f2(x1,x2;t1,t2)dt1dt2=R(t1,t1+τ)R(t_1,t_2)=E[\xi(t_1)\xi(t_2)]=\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}x_1x_2f_2(x_1,x_2;t_1,t_2)\,dt_1dt_2=R(t_1,t_1+\tau)R(t1?,t2?)=E[ξ(t1?)ξ(t2?)]=∫?∞+∞?∫?∞+∞?x1?x2?f2?(x1?,x2?;t1?,t2?)dt1?dt2?=R(t1?,t1?+τ) | 同一隨機(jī)過(guò)程在兩個(gè)不同時(shí)刻的隨機(jī)變量之間的關(guān)聯(lián)程度 |
| 互相關(guān)函數(shù) | Rξη(t1,t2)=E[ξ(t1)η(t2)]=Rξη[t1,t1+τ]R_{\xi\eta}(t_1,t_2)=E[\xi(t_1)\eta(t_2)]=R_{\xi\eta}[t_1,t_1+\tau]Rξη?(t1?,t2?)=E[ξ(t1?)η(t2?)]=Rξη?[t1?,t1?+τ] | 兩個(gè)隨機(jī)過(guò)程在兩個(gè)不同時(shí)刻的隨機(jī)變量之間的關(guān)聯(lián)程度 |
平均功率=直流功率+交流功率
上表中ξ(t)\xi(t)ξ(t)是隨機(jī)變量,其中的ttt取任意時(shí)刻.
二 平穩(wěn)隨機(jī)過(guò)程
平穩(wěn)隨機(jī)過(guò)程
| 嚴(yán)平穩(wěn) | 隨機(jī)過(guò)程 ξ(t)\xi(t)ξ(t) 的任意 n維分布與時(shí)間起點(diǎn)無(wú)關(guān) | 一維分布與時(shí)間t無(wú)關(guān);二維分布只與時(shí)間間隔τ\tauτ有關(guān) |
| 廣義平穩(wěn) | a(t)=a,R(t1,t1+τ)=R(τ)a(t)=a,R(t_1,t_1+\tau)=R(\tau)a(t)=a,R(t1?,t1?+τ)=R(τ) | 數(shù)學(xué)期望是個(gè)常數(shù),和時(shí)間t無(wú)關(guān);自相關(guān)函數(shù)只與時(shí)間間隔τ\tauτ有關(guān) |
| 兩者關(guān)系 | 嚴(yán)平穩(wěn)一定是廣義平穩(wěn) | 廣義平穩(wěn)不一定是嚴(yán)平穩(wěn) |
各態(tài)歷經(jīng)性
任取平穩(wěn)隨機(jī)過(guò)程 ξ(t)\xi(t)ξ(t)的任一樣本函數(shù) x(t)x(t)x(t),其時(shí)間均值和時(shí)間自相關(guān)滿足
a ̄=x(t) ̄=lim?T→∞1T∫?T/2T/2x(t)dt=a\overline{a}=\overline{x(t)}=\lim_{T\rightarrow\infty}\frac{1}{T}\int_{-T/2}^{T/2}x(t)\,dt=aa=x(t)?=T→∞lim?T1?∫?T/2T/2?x(t)dt=a
R(τ) ̄=x(t)x(t+τ) ̄=lim?T→∞1T∫?T/2T/2x(t)x(t+τ)dt=R(τ)\overline{R(\tau)}=\overline{x(t)x(t+\tau)}=\lim_{T\rightarrow\infty}\frac{1}{T}\int_{-T/2}^{T/2}x(t)x(t+\tau)\,dt=R(\tau)R(τ)?=x(t)x(t+τ)?=T→∞lim?T1?∫?T/2T/2?x(t)x(t+τ)dt=R(τ)
則稱平穩(wěn)隨機(jī)過(guò)程 ξ(t)\xi(t)ξ(t)具有各態(tài)歷經(jīng)性。
意義:可用任意一次實(shí)現(xiàn)的“樣本平均”來(lái)取代隨機(jī)過(guò)程的“統(tǒng)計(jì)平均”,可用任意一次實(shí)現(xiàn)的功率譜密度來(lái)取代隨機(jī)過(guò)程的功率譜密度,簡(jiǎn)化測(cè)量和計(jì)算問(wèn)題;具有各態(tài)歷經(jīng)性的隨機(jī)過(guò)程一定是平穩(wěn)隨機(jī)過(guò)程,反之不一定成立。
三 高斯過(guò)程
(1)定義:任意 n維概率密度都服從正態(tài)分布的隨機(jī)過(guò)程。
(2)重要性質(zhì):
高斯過(guò)程若廣義平穩(wěn),則必狹義平穩(wěn);
高斯過(guò)程中的隨機(jī)變量之間若不相關(guān),則它們統(tǒng)計(jì)獨(dú)立;
若干個(gè)高斯過(guò)程之和仍是高斯過(guò)程;
高斯過(guò)程經(jīng)線性變換后,仍是高斯過(guò)程。
均值為a,方差為σ2\sigma^2σ2的高斯過(guò)程的一維概率密度函數(shù)
當(dāng)均值a=0,標(biāo)準(zhǔn)擦σ=1\sigma=1σ=1時(shí),稱f(x)f(x)f(x)為標(biāo)準(zhǔn)正態(tài)分布函數(shù)
下面給出一些常用的求通信系統(tǒng)的誤碼率的函數(shù)
Q函數(shù)、誤差函數(shù)、互補(bǔ)誤差函數(shù),定義如下
Q(α)Q(\alpha)Q(α)的幾何意義如下:表示大于等于α\alphaα的概率
四.窄帶隨機(jī)過(guò)程
1.定義和表達(dá)式
窄帶隨機(jī)過(guò)程是指其頻帶寬度Δf\Delta fΔf遠(yuǎn)遠(yuǎn)小于中心頻率fcf_cfc?的過(guò)程。用示波器觀察它的波形,是一個(gè)頻率近似fcf_cfc?,包絡(luò)aξ(t)a_{\xi}(t)aξ?(t)和相位?c(t)\phi_c(t)?c?(t) 隨機(jī)緩變的正弦波。
窄帶隨機(jī)過(guò)程的一般表示式:
ξ(t)=aξ(t)cos[wct+?ξ(t)]\xi(t)=a_\xi (t)cos[w_ct+\phi_\xi(t)]ξ(t)=aξ?(t)cos[wc?t+?ξ?(t)]
等價(jià)式
ξ(t)=ξc(t)coswct?ξs(t)sinwct\xi(t)=\xi_c (t)cosw_ct-\xi_s(t)sinw_ctξ(t)=ξc?(t)coswc?t?ξs?(t)sinwc?t
式中ξc(t)\xi_c (t)ξc?(t)和ξs(t)\xi_s (t)ξs?(t)意思如下:
同相分量ξc(t)=aξ(t)cos?ξ(t)\xi_c (t)=a_\xi (t)cos\phi_\xi(t)ξc?(t)=aξ?(t)cos?ξ?(t)
正交分量
ξs(t)=aξ(t)sin?ξ(t)\xi_s (t)=a_\xi (t)sin\phi_\xi(t)ξs?(t)=aξ?(t)sin?ξ?(t)
2.統(tǒng)計(jì)特性
窄帶隨機(jī)過(guò)程ξ(t)\xi(t)ξ(t)的統(tǒng)計(jì)特性可有aξ(t)a_\xi (t)aξ?(t),?ξ(t)\phi_\xi (t)?ξ?(t)來(lái)決定,或者由ξc(t)\xi_c (t)ξc?(t)和ξs(t)\xi_s (t)ξs?(t)的統(tǒng)計(jì)特性決定。反過(guò)來(lái),可以由ξ(t)\xi(t)ξ(t)的統(tǒng)計(jì)特性可確定aξ(t)a_\xi (t)aξ?(t),?ξ(t)\phi_\xi (t)?ξ?(t)的統(tǒng)計(jì)特性,或者ξc(t)\xi_c (t)ξc?(t)和ξs(t)\xi_s (t)ξs?(t)的統(tǒng)計(jì)特性。
下面給出兩個(gè)重要的結(jié)論
-
一個(gè)均值為0,方差為σ?2\sigma_\epsilon^2σ?2?的窄帶平穩(wěn)高斯過(guò)程ξ(t)\xi(t)ξ(t),它的同相分量ξc(t)\xi_c (t)ξc?(t),正交分量ξs(t)\xi_s (t)ξs?(t)同樣是平穩(wěn)隨機(jī)過(guò)程,且均值都為零,方差也相同。另外,在同一時(shí)刻得到的同相分量ξc(t)\xi_c (t)ξc?(t)和正交分量ξs(t)\xi_s (t)ξs?(t)是統(tǒng)計(jì)獨(dú)立的。
數(shù)學(xué)公式為
-
一個(gè)均值為0,方差為σ?2\sigma_\epsilon^2σ?2?的窄帶平穩(wěn)高斯過(guò)程ξ(t)\xi(t)ξ(t),其包絡(luò)aξ(t)a_\xi (t)aξ?(t)的一維分布是瑞利分布,相位?ξ(t)\phi_\xi (t)?ξ?(t)的一維分布是均勻分布,且就一維分布而言,aξ(t)a_\xi (t)aξ?(t)和?ξ(t)\phi_\xi (t)?ξ?(t)是相互獨(dú)立的。
數(shù)學(xué)表述為
3.白噪聲
凡是功率譜密度在整個(gè)頻域內(nèi)都是均勻分布的噪聲,稱為白噪聲。
白噪聲的功率譜密度為
Pξ(w)=n02P_\xi(w)=\frac{n_0}{2}Pξ?(w)=2n0??
由于R(τ)?Pξ(w)R(\tau) \Leftrightarrow P_\xi(w)R(τ)?Pξ?(w),且δ(τ)?1\delta(\tau)\Leftrightarrow1δ(τ)?1
白噪聲的自相關(guān)函數(shù)為
R(τ)=n02δ(τ)R(\tau)=\frac{n_0}{2}\delta(\tau)R(τ)=2n0??δ(τ)
由上面自相關(guān)表達(dá)式,我們知道,白噪聲只有在τ=0\tau=0τ=0的時(shí)候才相關(guān),它在任意兩個(gè)時(shí)刻上的隨機(jī)變量都是不相關(guān)的。
這里補(bǔ)充一點(diǎn),以后討論的熱噪聲和散彈噪聲都視為白噪聲。
4.帶限白噪聲
白噪聲被限制在(?f0,f0)(-f_0,f_0)(?f0?,f0?)內(nèi),則這樣的白噪聲被稱為帶限白噪聲。
功率譜密度
自相關(guān)函數(shù)
五.正弦波加窄帶高斯過(guò)程
通信系統(tǒng)中經(jīng)常會(huì)遇到正弦波加窄帶高斯噪聲的情況。
考慮窄帶高斯噪聲的一般表達(dá)式
ξ(t)=ξc(t)coswct?ξs(t)sinwct\xi(t)=\xi_c (t)cosw_ct-\xi_s(t)sinw_ctξ(t)=ξc?(t)coswc?t?ξs?(t)sinwc?t
得到正弦波+窄帶高斯噪聲的時(shí)域表達(dá)式
其中,s(t)表示正弦波,n(t)表示窄帶高斯噪聲
六.平穩(wěn)隨機(jī)過(guò)程通過(guò)線性系統(tǒng)
設(shè)線性系統(tǒng)的沖激響應(yīng)是h(t)h(t)h(t),輸入為ξi(t)\xi_i(t)ξi?(t),則輸出為
ξo(t)=ξi(t)?h(t)=∫?∞∞h(τ)ξi(t?τ)dτ=∫?∞∞ξi(τ)h(t?τ)dτ\xi_o(t)=\xi_i(t)* h(t)=\int_{-\infin}^{\infin}h(\tau)\xi_i(t-\tau)d\tau=\int_{-\infin}^{\infin}\xi_i(\tau)h(t-\tau)d\tauξo?(t)=ξi?(t)?h(t)=∫?∞∞?h(τ)ξi?(t?τ)dτ=∫?∞∞?ξi?(τ)h(t?τ)dτ
因?yàn)橄到y(tǒng)是物理可實(shí)現(xiàn)的,有
這是要保證系統(tǒng)沖激響應(yīng)h(t)為正值。
輸出信號(hào)的均值
輸出信號(hào)的功率譜密度
參考資料:樊昌信《通信原理》考點(diǎn)精講
總結(jié)
以上是生活随笔為你收集整理的通信原理随机信号分析的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 位运算求a+b
- 下一篇: vectorint[]和vectorin