數(shù)字信號的最佳接收性能的研究

1、數(shù)字信號的最佳接收性能的研究 科研訓練論文實驗總成績： 裝 訂 線報告份數(shù)： 西安郵電大學 通信與信息工程學院 科研訓練論文專業(yè)班級： 學生姓名： 學號(班內序號): 2013 年 9月 22日數(shù)字信號的最佳接收性能的研究The research about the optimum reception of digital signals摘要數(shù)字信號的最佳接收是通信理論中一個重要的問題。本文主要從接收機性能的角度，分析了數(shù)字信號如何在同樣信道噪聲的條件下，使得正確接收信號的概率最大，而錯誤接收信號的概率減到最小，即最佳接收的問題。應用Matlab進行仿真、分析，通過對系統(tǒng)信號誤碼率的仿真，更加

2、直觀的了解數(shù)字信號的最佳接收性能，以便于比較，評論和改進。AbstractThe optimum reception of digital signals is an important problem in the theory of communication. This article from the receiver's perspective, analyses how the digital signal on the same channel noise conditions, makes the highest probability of correctly rec

3、eived signal data, and to minimize the probability of receiving signals, namely best reception problems. Application of Matlab for simulation, analysis, through to system signal error rate simulation, a more intuitive understanding of the best receiving performance of digital signals in order to com

4、pare, review and improvement.關鍵詞：數(shù)字信號, 最佳接收, 誤碼率。Key words：digital siginal, the optimum reception, probability of error.引言 在數(shù)字通信系統(tǒng)中，接收端收到的的發(fā)送信號和信道噪聲之和。噪聲對數(shù)字信號的影響表現(xiàn)在使接收碼元發(fā)生錯誤。由于數(shù)字通信系統(tǒng)傳輸質量的主要指標是錯誤概率，因此，將錯誤概率最小作為“最佳”的準則是恰當?shù)?。在此次研究中暫不考慮失真的影響，主要討論在二進制數(shù)字通信系統(tǒng)中如何使噪聲引起的錯誤概率最小，從而達到最佳接收的效果。其分析基本原理是將一個接收信號碼元的全部抽

5、樣當作為k維接收矢量空間中的一個矢量，并將接收矢量空間劃分為兩個區(qū)域。按照接收矢量落入哪個區(qū)域來判決是否發(fā)生錯誤。由判決準則可以得出最佳接收機的原理方框圖和計算誤碼率。這個誤碼率在理論上是最佳的，即理論上最小可能達到的。根據(jù)控制載波參量方式的不同，數(shù)字調制主要有調幅（ASK），調頻(FSK)，調相(PSK) 三種基本形式。我們本次的研究主要針對于二進制的2ASK、2FSK、2PSK進行討論，通過對系統(tǒng)信號誤碼率的仿真，了解數(shù)字信號的最佳接收性能。數(shù)字信號的傳輸方式分為基帶傳輸和帶通傳輸，在實際應用中，大多數(shù)信道具有帶通特性而不能直接傳輸基帶信號。為了使數(shù)字信號在帶通信道中傳輸，必須使用數(shù)字基帶

6、信號對載波進行調制，以使信號與信道的特性相匹配。這種用數(shù)字基帶信號控制載波，把數(shù)字基帶信號變換為數(shù)字帶通信號的過程稱為數(shù)字調制。通常使用鍵控法來實現(xiàn)數(shù)字調制，比如對載波的振幅、頻率和相位進行鍵控。2ASK: 2ASK信號的產(chǎn)生方法通常有兩種：模擬調制和鍵控法。解調有相干解調和非相干解調。P=1時f(t)=Acoswt;p=0時f(t)=0;其功率譜密度是基帶信號功率譜的線性搬移2FSK: 一個2FSK信號可以看成是兩個不同載波的2ASK信號的疊加。其解調和解調方法和ASK差不多。2FSK信號的頻譜可以看成是f1和f2的兩個2ASK頻譜的組合。2PSK: 2PSK以載波的相位變化作為參考基準的，

7、當基帶信號為0時相位相對于初始相位為0，當基帶信號為1時相對于初始相位為180°。一、最佳接收的概念及常見的接收信號方法最佳接收理論研究從噪聲中如何最好地提取有用信號?！白罴选钡臏蕜t：錯誤概率最小。產(chǎn)生錯誤的原因：暫不考慮失真的影響，主要討論在二進制數(shù)字通信系統(tǒng)中如何使噪聲引起的錯誤概率最小。判決規(guī)則：設在一個二進制通信系統(tǒng)中發(fā)送碼元“1”的概率為P(1)，發(fā)送碼元“0”的概率為P(0)，則總誤碼率Pe等于式中Pe1 = P(0/1) 發(fā)送“1”時，收到“0”的條件概率； Pe0 = P(1/0) 發(fā)送“0”時，收到“1”的條件概率；上面這兩個條件概率稱為錯誤轉移概率。裝 訂 線A0

8、A1r由接收矢量決定的兩個聯(lián)合概率密度函數(shù)f0(r)和f1(r)的曲線畫在下圖中當先驗概率相等時，即P(1) = P(0)時，f0(r0) = f1(r0)，所以最佳分界點位于圖中兩條曲線交點處的r 值上。在判決邊界確定之后，按照接收矢量r 落在區(qū)域A0應判為收到的是“0”的判決準則，這時有：若 ， 則判為“0” ； 反之，若 ， 則判為“1” 。在發(fā)送“0”和發(fā)送“1”的先驗概率相等時，上兩式的條件簡化為： 若f0(r) > f1(r)，則判為“0” 若f0(r) < f1(r)，則判為“1”這個判決準則常稱為最大似然準則。按照這個準則判決就可以得到理論上最佳的誤碼率，即達到理論

9、上的誤碼率最小值最佳接收機的原理方框圖： 確知信號接收信號 隨參信號 隨機相位信號 起伏信號二、研究討論出適應于數(shù)字信號的最佳接收實際接收機和最佳接收機的性能比較：由此可以得知，實際接收機與最佳接收機的性能在公式的形式上是一樣的。 普通接收系統(tǒng)的r(r=S/N)與最佳接收系統(tǒng)的Eb/n0相對應。最佳接收性能特點：誤碼率僅和Eb / n0以及相關系數(shù)r有關，與信號波形及噪聲功率無直接關系。 碼元能量Eb與噪聲功率譜密度n0之比，實際中相當于信號噪聲功率比Ps/Pn。因為若系統(tǒng)帶寬B等于1/Ts，則有 按照能消除碼間串擾的奈奎斯特速率傳輸基帶信號時，所需的最小帶寬為(1/2Ts) Hz。對于已調信

10、號，若采用的是2PSK或2ASK信號，則其占用帶寬應當是基帶信號帶寬的兩倍，即恰好是(1/Ts) Hz。所以，在工程上，通常把(Eb/n0)當作信號噪聲功率比看待。相關系數(shù) r 對于誤碼率的影響。 當兩種碼元的波形相同，相關系數(shù)最大，即r = 1時，誤碼率最大。這時的誤碼率Pe = 1/2。因為這時兩種碼元波形沒有區(qū)別，接收端是在沒有根據(jù)的亂猜。 當兩種碼元的波形相反，相關系數(shù)最小，即r = -1時，誤碼率最小。這時的最小誤碼率等于 例如2PSK信號的相關系數(shù)就等于 -1。當兩種碼元正交，即相關系數(shù) r 等于0時，誤碼率等于 例如2FSK信號的相關系數(shù)就等于或近似等于零若兩種碼元中有一種的能量

11、等于零，例如2ASK信號，則誤碼率為比較以上3種情況，它們之間的性能差3dB，即2ASK信號的性能比2FSK信號的性能差3dB，而2FSK信號的性能又比2PSK信號的性能差3dB。三、誤碼率曲線的繪制方法： 一般蒙特卡洛仿真是針對最佳接收機的，對于來自同一信道的接收信號，信道高斯噪聲的功率譜密度是不變的，最佳接收機性能之所以比一般接收機好，是由于其輸入信噪比是一般接收機的2倍（一般接收機帶通濾波器帶寬按2/Ts計算）；故信道高斯白噪聲的單邊帶功率譜密度即可按最佳接收機設計，也可按普通接收機設計。蒙特卡羅方法也稱為隨機模擬方法,有時也稱為隨機抽樣技術或統(tǒng)計實驗方法。它的基本思想是:為了求解數(shù)學、

12、物理、工程技術以及生產(chǎn)管理等方面的問題,首先建立一個概率模型或隨機過程,使它的參數(shù)等于問題的解;然后通過對模型或過程的觀察或抽樣試驗來計算所求參數(shù)的統(tǒng)計特征,最后給出所求解的近似值。而解得精確度可用估計值的標準誤差來表示。蒙特卡羅方法可以解決各種類型的問題,但總的來說,視其是否涉及隨機過程的性態(tài)和結果,該方法處理的問題可以分為兩類:第一類是確定性的數(shù)學問題,首先建立一個與所求解有關的概率模型,使所求的解就是我們所建立模型的概率分布或數(shù)學期望;然后對其進行隨機抽樣觀察,即產(chǎn)生隨機變量;最后用其算術平均值作為所求解的近似估計值。第二類是隨機性問題,被考察的元素更多的受到隨機性的影響,一般情況下采用

13、直接模擬方法,即根據(jù)實際物理情況的概率法則,用電子計算機進行抽樣試驗。在應用蒙特卡羅方法解決實際問題的過程中,大體有如下幾個內容:(1)對求解的問題建立簡單而又便于實現(xiàn)的概率統(tǒng)計模型,使所求的解恰好是所建立模型的概率分布或數(shù)學期望。(2)根據(jù)概率統(tǒng)計模型的特點和計算實踐的需要,盡量改進模型,以便減小方差和費用,提高計算效率。(3)建立對隨機變量的抽樣方法,其中包括建立產(chǎn)生偽隨機數(shù)的方法和建立對所遇到的分布產(chǎn)生隨機變量的隨機抽樣方法。(4)給出獲得所求解的統(tǒng)計估計值及其方差或標準誤差的方法。四、仿真實現(xiàn)MATLAB是一種功能強大的科學計算和工程仿真軟件，它的交互式集成界面能夠幫助用戶快速的完成數(shù)

14、值分析、數(shù)字信號處理、仿真建模、和優(yōu)化等功能。本課程設計需要運用MATLAB編程實現(xiàn)2ASK,2FSK,2PSK調制解調過程，并且輸出其調制后的波形，畫出頻譜、功率譜密度圖、誤碼率曲線，并重點比較各種調制信號的誤碼率情況，討論最佳接收性能。通信系統(tǒng)的抗噪聲性能是指系統(tǒng)克服加性噪聲影響的能力。在數(shù)字通信系統(tǒng)中，信道的加性噪聲能使傳輸碼元產(chǎn)生錯誤，錯誤程度通常用誤碼率來衡量。與數(shù)字基帶系統(tǒng)一樣，分析二進制數(shù)字調制系統(tǒng)的抗噪聲性能，也就是要計算系統(tǒng)由加性噪聲產(chǎn)生的總誤碼率。五、結論：由matlab仿真繪制出實際接收機和最佳接收機的誤碼率曲線對比圖，如下。從橫向和縱向的角度分別對比分析數(shù)字信號的最佳接

15、受機性能： 由圖可以得出結論：三種信號之間的性能差3dB，即2ASK信號的性能比2FSK信號的性能差3dB，而2FSK信號的性能又比2PSK信號的性能差3dB。將實際接收機和最佳接收機的誤碼率曲線作比較可以看出，若實際接收機中的信號噪聲功率比r等于最佳接收機中的碼元能量和噪聲功率譜密度之比Eb/n0，則兩者的誤碼率性能一樣。但是，由于實際接收機總不可能達到這一點。所以，實際接收機的性能總是比不上最佳接收機的性能。參考文獻1張志涌 徐彥琴等編著.MATLAB教程. 北京：北京航空航天大學出版社，2001年2郭文彬 桑林.通信原理基于Matlab的計算機仿真.北京：北京郵電大學出版社2009年3樊

16、昌信,曹麗娜. 通信原理. 第6版.北京：國防工業(yè)出版社，2009年：194-205，298-311.4葛哲學. 精通MATLAB. 電子工業(yè)出版社 .2008年5 孫學軍等. 通信原理. 電子工業(yè)出版社. 2001年6 曹志剛等. 現(xiàn)代通信原理. 北京：清華大學出版社 . 1992年8月7 曾一凡. 擴頻通信原理. 機械工業(yè)出版社 . 2005年9月8 宮錦文. 通信原理實驗指導書. 電子工業(yè)出版社 . 2007年9 Cahn C R. Performance of Digital Phase-Modulation Communication Systems.IRE Trans.On Com

17、mun.Systems 1957.7源程序主程序1：clear all;close all;echo on%-系統(tǒng)仿真參數(shù)A=1; %載波振幅fc=2; %載波頻率（Hz）snr=1; %信噪比dBN_sample=8;% 基帶信號中每個碼元的的采樣點數(shù)N=10000; % 碼元數(shù)Ts=1; % 碼元寬度df=0.01%頻率分辨率B=1/Ts;f_start=fc-B;f_cutoff=fc+B;fs=fc*N_sample%系統(tǒng)采樣頻率，即考慮載波后，一個碼元內的采樣點數(shù)ts=Ts/fs; % 系統(tǒng)采樣間隔t=0:ts:N*Ts-ts;Lt=length(t);%-畫出調制信號波形及功率譜%

18、 產(chǎn)生二進制信源d=sign(randn(1,N);dd=sigexpand(d+1)/2,fc*N_sample);gt=ones(1,fc*N_sample); % NRZ波形d_NRZ=conv(dd,gt);d_NRZ1=d_NRZ(1:Lt);pause%畫出單極性NRZ波形及其功率譜figure(1)subplot(221); plot(t,d_NRZ1);% 畫出單極性NRZ信號波形axis(0 50 0 1.2);xlabel('t');ylabel('單極性信號');subplot(222); d_NRZ1f,d_NRZ1,df1,f=T2F(

19、d_NRZ1,ts,df,fs);%求出單極性NRZ信號功率譜plot(f,10*log10(abs(fftshift(d_NRZ1f).2/length(f);% 畫出單極性NRZ信號功率譜 axis(-3*B 3*B -50 0);xlabel('f');ylabel('單極性信號PDF');pause%畫出雙極性NRZ波形及其功率譜d_sjx=2*d_NRZ-1;%生成雙極性NRZ信號d_sjx1=d_sjx(1:Lt);subplot(223); plot(t,d_sjx1);% 畫出雙極性NRZ信號波形axis(0 50 0 1.2);xlabel(&

20、#39;t');ylabel('雙極性信號');subplot(224); d_sjx1f,d_sjx1,df1,f=T2F(d_sjx1,ts,df,fs);%求出雙極性NRZ信號功率譜plot(f,10*log10(abs(fftshift(d_sjx1f).2/length(f);% 畫出雙極性NRZ信號功率譜axis(-3*B 3*B -50 0);xlabel('f');ylabel('雙極性信號PDF');%-畫出數(shù)字頻帶信號及其功率譜% 對數(shù)字基帶信號進行2ASK調制 ht=A*sin(2*pi*fc*t);%載波s_2as

21、k=d_NRZ(1:Lt).*ht;%生成已調信號2ASKpause%畫出已調信號2ASK及其功率譜figure(2)subplot(211);plot(t,s_2ask);%畫出2ASK信號axis(0 50 -1.2 1.2);xlabel('t');ylabel('2ASK');%求2ASK信號功率譜s_2askf,s_2ask,df1,f=T2F(s_2ask,ts,df,fs);%求出單極性NRZ信號功率譜subplot(212);plot(f,10*log10(abs(fftshift(s_2askf).2/length(f);% 畫出單極性NRZ信

22、號功率譜axis(-fc-3*B fc+3*B -50 0);xlabel('f');ylabel('2ASK信號PDF');% 對數(shù)字基帶信號進行2FSK調制 s_2fsk=A*cos(2*pi*fc*t+2*pi*d_sjx(1:Lt).*t);%生成2FSK信號 pause%畫出已調信號2FSK及其功率譜subplot(223)plot(t,s_2fsk);%畫出2FSK波形axis(0 10 -1.2 1.2);xlabel('t');ylabel('2FSK');subplot(224)s_2fskf,s_2fsk,df

23、1,f=T2F(s_2fsk,ts,df,fs);%求2FSK信號功率譜plot(f,10*log10(abs(fftshift(s_2fskf).2/length(f);% 畫出2FSK功率譜axis(-fc-3*B fc+3*B -50 0);xlabel('f');ylabel('2FSK信號PDF');% 對數(shù)字基帶信號進行2PSK調制 s_2psk=d_sjx(1:Lt).*ht;%生成2PSK信號pause%畫出已調信號2PSK及其功率譜figure(3)subplot(321)plot(t,s_2psk);%畫出2PSK波形axis(0 10 -1

24、.2 1.2);xlabel('t');ylabel('2PSK'); subplot(322)s_2pskf,s_2psk,df1,f=T2F(s_2psk,ts,df,fs);%求2PSK信號功率譜plot(f,10*log10(abs(fftshift(s_2pskf).2/length(f);% 畫出2PSK信號功率譜axis(-fc-3*B fc+3*B -50 0);xlabel('f');ylabel('2PSK信號PDF');%-將2PSK信號送入信道進行傳輸，先生成信道加性高斯白噪聲噪聲snr_lin=10(sn

25、r/10); %換算成倍數(shù)signal_energy=0.5*A2*Ts;%求出接收信號平均能量noise_power=( signal_energy *fs)/(snr_lin*4);%求出噪聲方差（噪聲均值為0）noise_std=sqrt(noise_power);%求出噪聲均方差noise=noise_std.*randn(1,Lt);%以噪聲均方差作為幅度產(chǎn)生高斯白噪聲%-將已調信號送入信道pause%畫出信道中的高斯白噪聲及其功率譜subplot(323)plot(t,noise(1:Lt);%畫出噪聲xlabel('t');ylabel('信道噪聲'

26、;);axis(0 10 -3 3);noisef,noise,df1,f=T2F(noise,ts,df,fs);%求信道噪聲功率譜subplot(324)plot(f,10*log10(abs(fftshift(noisef).2/length(f);% 畫出信道噪聲功率譜axis(-fs/2-2 fs/2+2 -50 0);xlabel('f');ylabel('信道噪聲PDF');r=s_2ask(1:Lt)+noise(1:Lt);%疊加了噪聲的已調信號，相當于將已調信號送入理想信道pause%畫出加噪后的已調信號2PSK及其功率譜subplot(32

27、5)plot(t,r);%畫出加噪后的已調信號2PSKxlabel('t');ylabel('加噪2ASK信號');axis(0 10 -3 3);rf,r,df1,f=T2F(r,ts,df,fs);%求加噪后的已調信號2ASK功率譜subplot(326)%畫出加噪后已調信號的功率譜plot(f,10*log10(abs(fftshift(rf).2/length(f);% 畫出已調信號2ASK功率譜axis(-fc-3*B fc+3*B -50 0);xlabel('f');ylabel('加噪2ASK信號PDF');%-在

28、接收端準備進行解調，先通過帶通濾波器pause%畫出帶通濾波器H,f=bp_f(length(s_2askf),f_start,f_cutoff,df1,fs,1);%經(jīng)過理想帶通濾波器figure(4)subplot(322)plot(f,abs(fftshift(H);% 畫出理想帶通濾波器axis(-fc-3*B fc+3*B -0.2 1.2);xlabel('f');ylabel('理想BPF');DEM = H.*s_2askf; %濾波器輸出的頻譜dem=F2T(DEM,fs);%濾波器的輸出波形 dem1=dem(1:Lt)pause%經(jīng)過理想帶

29、通濾波器后的信號波形及功率譜subplot(323)%經(jīng)過理想帶通濾波器后的信號波形plot(t,dem1)%畫出經(jīng)過理想帶通濾波器后的信號波形axis(0 50 -4 4);xlabel('t');ylabel('理想BPF輸出信號');demf1,dem1,df1,f=T2F(dem1,ts,df,fs);%求經(jīng)過理想帶通濾波器后信號功率譜subplot(324)plot(f,10*log10(abs(fftshift(demf1).2/length(f);% 畫出經(jīng)過理想帶通濾波器后信號功率譜axis(-fc-3*B fc+3*B -50 0);xlabe

30、l('f');ylabel('理想BPF輸出信號PDF');%-進行相干解調，先和本地載波相乘,即混頻subplot(325)%畫出同頻同相的本地載波plot(t,ht);axis(0 50 -1.2 1.2);xlabel('t');ylabel('本地載波');subplot(326)% 本地載波頻譜htf,ht,df1,f=T2F(ht,ts,df,fs);plot(f,fftshift(abs(htf)% 畫出載波頻譜axis(-fc-3*B fc+3*B 0 15);xlabel('f');ylabel(

31、'本地載波頻譜');pause%畫出混頻后的信號及功率譜figure(5)der=dem1(1:Lt).*ht(1:Lt);%和本地載波相乘，即混頻subplot(221)%畫出混頻后的波形plot(t,der);axis(0 50 -2 2);xlabel('t');ylabel('混頻后的信號');derf,der,df1,f=T2F(der,ts,df,fs);%求混頻后信號的功率譜subplot(222)plot(f,10*log10(abs(fftshift(derf).2/length(f);%畫出混頻后的功率譜axis(-fc-3*

32、B fc+3*B -50 0);xlabel('f');ylabel('混頻后信號的PDF');%-再經(jīng)過低通濾波器pause%畫出理想低通濾波器LPF,f=lp_f(length(derf),B,df1,fs,1);%求低通濾波器subplot(224) % 畫出理想低通濾波器plot(f,fftshift(abs(LPF);axis(-fc-3*B fc+3*B -0.2 1.2);xlabel('f');ylabel('理想LPF');pause%混頻信號經(jīng)理想低通濾波器后的波形及功率譜DM = LPF.*derf; %理想

33、低通濾波器輸出信號頻譜dm=F2T(DM,fs); %理想低通濾波器的輸出波形figure(6)subplot(221)plot(t,dm(1:Lt);%畫出經(jīng)過低通濾波器后的解調出的波形axis(0 50 -1.2 1.2);xlabel('t');ylabel('LPF輸出信號');subplot(222)dmf,dm,df1,f=T2F(dm,ts,df,fs);%求LPF輸出信號的功率譜plot(f,10*log10(abs(fftshift(dmf).2/length(f);%畫出LPF輸出信號的功率譜axis(-fc-3*B fc+3*B -50 0

34、);xlabel('f');ylabel('LPF輸出信號功率譜');%-最后對LPF輸出信號抽樣判決dm= dm(1:Lt);panjue=zeros(1,N);%建立存儲判決值的矩陣%抽樣判決，規(guī)則：大于等于0判1，小于0判-1for i=1:N; if dm(fc*N_sample*(i-1)+fc*N_sample/2+1)>=0;%抽樣判決時刻 panjue(i)=1; else panjue(i)=-1; endend%-畫出判決出的基帶信號波形，并和調制信號比較rr=sigexpand(panjue,fc*N_sample);rrt=ones

35、(1,fc*N_sample); % NRZ波形huifu_NRZ=conv(rr,rrt);pause%觀察調制信號和解調信號波形subplot(224)plot(t,d_sjx(1:Lt);%調制信號波形axis(0 50 -1.2 1.2);xlabel('t');ylabel('調制信號');subplot(223)plot(t,huifu_NRZ(1:Lt);%解調信號波形axis(0 50 -1.2 1.2);xlabel('t');ylabel('解調信號');%-統(tǒng)計誤碼數(shù)numoferr=sum(abs(panj

36、ue-d)/2)/N%計算出錯誤碼元數(shù)主程序2 誤碼率：clear allclose allA=1; %載波振幅fc=2; %載波頻率（Hz）SNRindB1=-5:1:12; % 信噪比取值向量，dB為單位SNRindB2=-5:0.2:12; % 信噪比取值向量，dB為單位N_sample=100;%每個碼元的的采樣點數(shù)N=10000; % 碼元數(shù)Ts=1; % 碼元寬度d=sign(rand(1,N)-0.5+eps);%產(chǎn)生雙極性二進制碼元df=0.01;simu_err_prb=zeros(1,length(SNRindB1);%理論誤碼率for i=1:length(SNRindB

37、2) % 計算信噪比值 SNR=10(SNRindB2(i)/10); % 計算普通接收機的理論誤碼率。theo_err_prba(i)=0.5*erfc(sqrt(0.25*SNR);%一般蒙特卡洛仿真是針對最佳接收機的，對于來自同一信道的接收信號，信道高斯噪聲的功率譜密度是不變的，最佳接收機性能之所以比一般接收機好，是由于其輸入信噪比是一般接收機的2倍（一般接收機帶通濾波器帶寬按2/Ts計算）；故信道高斯白噪聲的單邊帶功率譜密度即可按最佳接收機設計，也可按普通接收機設計。%互補誤差函數(shù)theo_err_prba1(i)=0.5*erfc(sqrt(0.5*SNR);%最佳接收機理論誤碼率曲

38、線 end%計算普通接收機實際誤碼率for i=1:length(SNRindB2) % 計算信噪比值 SNR=10(SNRindB2(i)/10); % 計算普通接收機的理論誤碼率。theo_err_prbp(i)=0.5*erfc(sqrt(SNR);%一般蒙特卡洛仿真是針對最佳接收機的，對于來自同一信道的接收信號，信道高斯噪聲的功率譜密度是不變的，最佳接收機性能之所以比一般接收機好，是由于其輸入信噪比是一般接收機的2倍（一般接收機帶通濾波器帶寬按2/Ts計算）；故信道高斯白噪聲的單邊帶功率譜密度即可按最佳接收機設計，也可按普通接收機設計。%互補誤差函數(shù)theo_err_prbp1(i)=

39、0.5*erfc(sqrt(2*SNR);%最佳接收機理論誤碼率曲線 end% 誤碼率曲線圖：估計值和理論值曲線對比圖for i=1:length(SNRindB2) % 計算信噪比值 SNR=10(SNRindB2(i)/10); % 計算普通接收機的理論誤碼率。theo_err_prbf(i)=0.5*erfc(sqrt(0.5*SNR);%一般蒙特卡洛仿真是針對最佳接收機的，對于來自同一信道的接收信號，信道高斯噪聲的功率譜密度是不變的，最佳接收機性能之所以比一般接收機好，是由于其輸入信噪比是一般接收機的2倍（一般接收機帶通濾波器帶寬按2/Ts計算）；故信道高斯白噪聲的單邊帶功率譜密度即可按最佳接收機設計，也可按