单径瑞利信道中的BPSK相干解调的(理论)误码率性能
目錄
- 引言
- 理論差錯(cuò)率推導(dǎo)
- 理論差錯(cuò)率編程繪圖
- 時(shí)分的導(dǎo)引輔助的信道估計(jì)方法
引言
最近完成老師給的作業(yè),題目如下:
(1)請(qǐng)推導(dǎo)出單徑瑞利信道中的BPSK相干解調(diào)的理論誤碼率性能,并畫(huà)出比特信噪比與誤碼率的關(guān)系曲線。
(2)在單徑瑞利信道中,請(qǐng)?jiān)O(shè)計(jì)一種時(shí)分的導(dǎo)引輔助的信道估計(jì)方法,用Simulink進(jìn)行仿真,測(cè)量BPSK的誤碼率性能,畫(huà)出比特信噪比與誤碼率的關(guān)系曲線,并與(1)的曲線進(jìn)行對(duì)比。
對(duì)于BPSK的調(diào)制解調(diào),可以參考之前的文章:
BPSK調(diào)制解調(diào)
理論差錯(cuò)率推導(dǎo)
BPSK一般是輸入±1進(jìn)行調(diào)制。在此處,我們假設(shè)輸入信號(hào)為:s1(t)=g(t)s_1(t)=g(t)s1?(t)=g(t)和s2(t)=?g(t)s_2(t)=-g(t)s2?(t)=?g(t)。其中g(t)g(t)g(t)是在區(qū)間[0,Tb][0,T_b][0,Tb?]非零而在其他處為零的任意脈沖,其能量為ξg\xi_gξg?。
則輸入信號(hào)可表示為s1(t)=ξbs_1(t)=\sqrt{\xi_b}s1?(t)=ξb??和s2(t)=?ξbs_2(t)=-\sqrt{\xi_b}s2?(t)=?ξb??。假設(shè)兩個(gè)信號(hào)是等概發(fā)送的,再假設(shè)此時(shí)發(fā)送信號(hào)s1s_1s1?。
在小尺度衰落信道中,發(fā)送信號(hào)s1(t)s_1(t)s1?(t)將發(fā)送乘性失真,單徑瑞利衰落信道意味著至少在一個(gè)信號(hào)傳輸間隔內(nèi),乘法過(guò)程可以看做是乘一個(gè)常數(shù)。因此,當(dāng)發(fā)送信號(hào)s1(t)s_1(t)s1?(t)時(shí),在一個(gè)信號(hào)傳輸間隔內(nèi)的等效低通接收信號(hào)為:
r1(t)=ae?j?s1(t)+z(t)(0≤t≤T)(1)r_1(t)=ae^{-j\phi}s_1(t)+z(t) (0\le t \le T) \tag{1}r1?(t)=ae?j?s1?(t)+z(t)(0≤t≤T)(1)
式(1)中,z(t)z(t)z(t)表示惡化信號(hào)的復(fù)高斯白噪聲過(guò)程。
假設(shè)信道衰落足夠慢,以致于相移?\phi?能夠從接收信號(hào)中無(wú)誤差地估計(jì)出來(lái)。在這種情況下,能夠?qū)崿F(xiàn)接收信號(hào)的理想相干檢測(cè)。于是,接收信號(hào)可用一個(gè)匹配濾波器來(lái)處理。
那么接收端經(jīng)過(guò)匹配濾波器得到的解調(diào)信號(hào)應(yīng)該為:r=as1+n=aξb+n(2)r=as_1+n=a\sqrt{\xi_b}+n \tag{2}r=as1?+n=aξb??+n(2)
式(1)中,n表示均值為0,方差為σn2=12N0\sigma_n^2=\frac{1}{2}N_0σn2?=21?N0?的加性高斯分量。a表示衰減。此時(shí),要判斷接收到的rrr究竟是s1s_1s1?還是s2s_2s2?,就要進(jìn)行判決。在先驗(yàn)概率相同的情況下,最大后驗(yàn)概率準(zhǔn)則和最大似然準(zhǔn)則的效果相同。
在本實(shí)驗(yàn)中,將r與閾值0進(jìn)行比較,如果r>0,則判為s1s_1s1?,否則判為s2s_2s2?
顯然,r被判決為s1和s2概率分布函數(shù)(PDF)分別是:
p(r∣s1)=?1πN0e?(r?aξb)2/N0(3)p(r|s_1)=-\frac{1}{\sqrt{\pi N_0}} e^{-(r-a\sqrt{\xi_b})^2/N_0} \tag{3}p(r∣s1?)=?πN0??1?e?(r?aξb??)2/N0?(3)
p(r∣s2)=?1πN0e?(r+aξb)2/N0(4)p(r|s_2)=-\frac{1}{\sqrt{\pi N_0}} e^{-(r+a\sqrt{\xi_b})^2/N_0} \tag{4}p(r∣s2?)=?πN0??1?e?(r+aξb??)2/N0?(4)
在給定了發(fā)送符號(hào)s1(t)s_1(t)s1?(t)的情況下,錯(cuò)誤概率就是r<0的概率。即
p(e∣s1)=∫?∞0p(r∣s1)dr=Q(2a2ξbN0)(5)p(e|s_1)=\int_{-\infty}^{0}p(r|s_1)dr \\ =Q\begin{pmatrix} \sqrt{\frac{2a^2\xi_b}{N_0}} \end{pmatrix} \tag{5}p(e∣s1?)=∫?∞0?p(r∣s1?)dr=Q(N0?2a2ξb????)(5)
式(5)是一個(gè)Q函數(shù),定義為Q(x)=12π∫x∞e?t2/2dt,x≥0Q(x)=\frac{1}{\sqrt{2 \pi}} \int_{x}^{\infty}e^{-t^2/2}dt,x\ge 0Q(x)=2π?1?∫x∞?e?t2/2dt,x≥0。同理,當(dāng)發(fā)射信號(hào)s2s2s2時(shí),錯(cuò)誤概率也是p(e∣s2)=Q(2a2ξbN0)p(e|s_2)=Q\begin{pmatrix}\sqrt{\frac{2a^2\xi_b}{N_0}}\end{pmatrix}p(e∣s2?)=Q(N0?2a2ξb????)。因?yàn)閟1和s2是等概率發(fā)送的,所以平均錯(cuò)誤概率仍然是
pb=12p(e∣s1)+12p(e∣s2)=Q(2a2ξbN0)(6)p_b=\frac{1}{2}p(e|s_1)+\frac{1}{2}p(e|s_2) \\ =Q\begin{pmatrix}\sqrt{\frac{2a^2\xi_b}{N_0}}\end{pmatrix} \tag{6}pb?=21?p(e∣s1?)+21?p(e∣s2?)=Q(N0?2a2ξb????)(6)
從式(6)我們可以發(fā)現(xiàn),錯(cuò)誤概率只跟比特信噪比:ξbN0\frac{\xi_b}{N_0}N0?ξb??和衰減a有關(guān)。我們定義接收信噪比γb=a2ξbN0\gamma_b=\frac{a^2\xi_b}{N_0}γb?=N0?a2ξb??。則BPSK的差錯(cuò)率為:pb(γb)=Q(2γb)(7)p_b(\gamma_b)=Q(\sqrt{2\gamma_b}) \tag{7}pb?(γb?)=Q(2γb??)(7)
我們把式(7)看成是條件差錯(cuò)概率,其條件是衰減系數(shù)a保持不變。但我們?yōu)榱说玫絘變化時(shí)的出錯(cuò)率,所以我們需要將(7)中的pb(γb)p_b(\gamma_b)pb?(γb?)對(duì)γb\gamma_bγb?求平均,即需要計(jì)算如下積分:
pb=∫0+∞pb(γb)p(γb)dγb(8)p_b=\int_{0}^{+\infty}p_b(\gamma_b)p(\gamma_b)d\gamma_b \tag{8}pb?=∫0+∞?pb?(γb?)p(γb?)dγb?(8)
其中,p(γb)p(\gamma_b)p(γb?)是a為隨機(jī)變量時(shí),γb\gamma_bγb?的概率密度函數(shù)。
此時(shí),我們只需知道p(γb)p(\gamma_b)p(γb?)即可。在單徑瑞利信道中:因?yàn)閍是瑞利分布的,故a2a^2a2為具有兩個(gè)自由度的χ2\chi^2χ2分布。因此, γb\gamma_bγb?也是χ2\chi^2χ2分布的。容易證明
p(γb)=1γbˉe?γbγb(γb≥0)(9)p(\gamma_b)=\frac{1}{\bar{\gamma_b}}e^{-\gamma_b\sqrt{\gamma_b}}(\gamma_b\ge0) \tag{9}p(γb?)=γb?ˉ?1?e?γb?γb??(γb?≥0)(9)
此時(shí)我們定義比特信噪比:BSNR=2ξb/N0BSNR=2\xi_b/N_0BSNR=2ξb?/N0?
式(9)中,γbˉ\bar{\gamma_b}γb?ˉ?是平均信噪比,它定義為:γbˉ=ξbN0E(a2)\bar{\gamma_b}=\frac{\xi_b}{N_0}E(a^2)γb?ˉ?=N0?ξb??E(a2)。由于a2a^2a2為具有兩個(gè)自由度的χ2\chi^2χ2分布,所以E(a2)=2E(a^2)=2E(a2)=2,γbˉ=4BSNR\bar{\gamma_b}=4BSNRγb?ˉ?=4BSNR。
讓我們嘗試把式(7)、(9)帶入式(8),積分結(jié)果為:
pb=12(1?γbˉ1+γbˉ)=12(1?4BSNR1+4BSNR)(10)p_b=\frac{1}{2}\begin{pmatrix}{1-\sqrt{\frac{\bar{\gamma_b}}{1+\bar{\gamma_b}}}}\end{pmatrix}=\frac{1}{2}\begin{pmatrix}{1-\sqrt{\frac{4BSNR}{1+4BSNR}}}\end{pmatrix} \tag{10}pb?=21?(1?1+γb?ˉ?γb?ˉ????)=21?(1?1+4BSNR4BSNR???)(10)
式(10)就是理想情況下的差錯(cuò)率表達(dá)式。
理論差錯(cuò)率編程繪圖
%author: ddatalent %date: 2021/06/28 clear;close;clc; snr=zeros(1,30); snr(1)=1e-3; for i=1:29snr(i+1)=snr(i)*1.5; end Pb=zeros(1,length(snr)); for i=1:length(snr)Pb(i)=1/2*(1-sqrt(4*snr(i)/(1+4*snr(i)))); end;figure() plot(20*log10(snr),Pb); xlabel('SNR(dB)'); ylabel('Pe'); title('理想誤碼率曲線'); grid on時(shí)分的導(dǎo)引輔助的信道估計(jì)方法
單徑瑞利信道,由于只有單徑的信號(hào),不考慮時(shí)延擴(kuò)展,它描述的是信道的時(shí)間選擇性衰落特性,也就是平坦衰落,我們?cè)谝粠?hào)的幀頭插入一個(gè)‘1’比特作為導(dǎo)引估計(jì)信道特性,并且可以認(rèn)為在一幀內(nèi)信道特性不會(huì)變化。
本文參考自:
《MIMO-OFDM Wireless Communications with MATLAB》
《數(shù)字通信》(第四版,第14章)
總結(jié)
以上是生活随笔為你收集整理的单径瑞利信道中的BPSK相干解调的(理论)误码率性能的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: 微型计算机控制赖寿宏pdf,微型计算机控
- 下一篇: 数据库系统的三级模式结构