無碼間串擾基帶傳輸系統(tǒng)的仿真設(shè)計

湯巧治,陳雙燕

（閩南理工學院電子與電氣工程學院,福建泉州 362700）

由于各種噪聲和隨機因素的影響，進行通信系統(tǒng)實際電路及性能的試驗與研究已經(jīng)變得不可行。若系統(tǒng)若干參數(shù)不滿足工程要求，意味著整個通信系統(tǒng)需要重新建設(shè)，如此代價實在太大。因此，Matlab作為強大的計算機輔助分析與設(shè)計工具，結(jié)合其系統(tǒng)仿真方法，人們可以迅速構(gòu)建一個通信系統(tǒng)模型，提供一個便捷、高效和精確的評估平臺。此法將有助于提高新技術(shù)理論成果轉(zhuǎn)化為實際產(chǎn)品的效率，進一步降低成本，因而越來越受到業(yè)界的青睞[1]。數(shù)字基帶傳輸是數(shù)字頻帶通信系統(tǒng)的基礎(chǔ)，其研究內(nèi)容涉及碼間串擾（InterSymbol Interference，ISI）的成因以及消除技術(shù)、功率譜（power spectral density，PSD）分析和抗噪聲性能分析等[2]，因此研究并改善基帶傳輸系統(tǒng)的性能意義重大。

1 基帶傳輸系統(tǒng)原理分析

基帶傳輸系統(tǒng)框圖如圖1所示。發(fā)送濾波器至接收濾波器總的傳輸特性為H(ω)，則H(ω)=GT(ω)C(ω)GR(ω)。其中，GT(ω)、C(ω)、GR(ω)分別對應發(fā)送濾波器、信道和接收濾波器的頻譜[2-3]。

圖1 基帶傳輸系統(tǒng)框圖

所謂基帶傳輸系統(tǒng)的碼間串擾ISI是由于系統(tǒng)傳輸總特性不理想，導致前后碼元的波形畸變、展寬，并使前面波形出現(xiàn)很長的拖尾，蔓延到當前碼元的采樣時刻上，從而對當前碼元的判決造成干擾[4]。相關(guān)理論分析表明：奈奎斯特第一準則是判斷一個基帶系統(tǒng)能否實現(xiàn)無ISI的理論依據(jù)。工程上，具有滾升余弦頻率特性的傳輸信道Hrcos(ω)可滿足無碼間串擾傳輸要求[1-2]，其頻譜為：

其中，Ts為一個碼元的持續(xù)時間，0 ≤α≤1稱為滾降系數(shù)。其沖激響應為：

對于物理信道是加性高斯白噪聲信道的情況，可以證明，當發(fā)送濾波器與接收濾波器相互匹配的時候，即GT(ω)=GR?(ω)時，通信性能（誤碼率最?。┻_到最佳[1]。對于理想的物理信道(C(ω)=1)，收發(fā)濾波器相互匹配時有：

此時可求得發(fā)送濾波器和接收濾波器的傳遞函數(shù)的實數(shù)解為：

無ISI條件下，信道傳遞函數(shù)是滾升余弦的，匹配的收發(fā)濾波器稱為平方根滾升余弦濾波器(Square Root Raised Cosine Filter)[1]，其傳遞函數(shù)的實數(shù)解為：

匹配的收發(fā)濾波器的沖激響應是：

工程上，滾升余弦濾波器和平方根滾升余弦濾波器通常用FIR濾波器來近似實現(xiàn)。Matlab通信工具箱中提供了設(shè)計升余弦濾波器的函數(shù)“rcosine”，用于計算升余弦或平方根升余弦濾波器的寫法為[1]：num=rcosine(Fd,Fs,'fir/normal',r,delay)或者num=rco?sine(Fd,Fs,'fir/sqrt',r,delay)。

其中，'fir/normal'用于FIR滾升余弦濾波器設(shè)計，'fir/sqrt'用于FIR平方根滾升余弦濾波器設(shè)計；r是滾降系數(shù)，在0到1之間取值；Fd為輸入數(shù)字序列的碼元速率；Fs為濾波器采樣率，其值必須是Fd的正整數(shù)倍；delay是濾波器輸入到響應峰值之間的時延（單位是碼元周期個數(shù)）。

2 基帶系統(tǒng)仿真設(shè)計

基于Matlab/Simulink對數(shù)字基帶無ISI傳輸系統(tǒng)進行系統(tǒng)建模，具體模型如圖2所示[5]。

圖2 無ISI基帶傳輸系統(tǒng)的Simulink模型（高斯噪聲Var=0.01 W）

基帶傳輸系統(tǒng)模型的發(fā)送部分各模塊說明如下，選用貝努利二進制序列發(fā)生器（Bernoulli Bina?ry Generator）模塊作為系統(tǒng)信號源，來產(chǎn)生單極性非歸零矩形脈沖（NRZ），設(shè)置信源傳碼率為1000 Baud（數(shù)值上等于比特率1000 bit/s），其輸出為沖激脈沖形式的數(shù)據(jù)序列。單雙極性轉(zhuǎn)換模塊將NRZ轉(zhuǎn)換成BRZ。系統(tǒng)仿真時，系統(tǒng)的仿真速率設(shè)置為10000（即1e4）Hz，對應的固定仿真步長為1e-4 s。因此BRZ數(shù)據(jù)進入發(fā)送濾波器之前需要進行采樣率（速率）轉(zhuǎn)換，由上采樣和速率轉(zhuǎn)換兩個模塊實現(xiàn)。發(fā)送濾波器為平方根升余弦濾波器，以DSP系統(tǒng)工具箱的“Discrete Filter”模塊實現(xiàn)，參數(shù)設(shè)置如圖3。

圖3 平方根滾升余弦濾波器參數(shù)設(shè)置

模塊關(guān)鍵參數(shù)說明如下：分子系數(shù)通過“rco?sine”函數(shù)計算，滾降系數(shù)為0.5，濾波器延時為10個碼元時間，故設(shè)置為：rcosine(1,10,'fir/sqrt',0.5,10)；其分母系數(shù)設(shè)置為1；采樣周期與系統(tǒng)仿真步長一致。信道為加性高斯信道，采用隨機數(shù)發(fā)生器和加法器來實現(xiàn)。隨機數(shù)發(fā)生器產(chǎn)生的高斯噪聲均值（Mean）設(shè)為0，Var設(shè)為固定值0.01 W。由隨機過程分析可知，噪聲的方差即是噪聲平均功率[2]。

接收部分，接收濾波器與發(fā)送濾波器相匹配，參數(shù)設(shè)置與發(fā)送濾波器相同。將接收信號進行眼圖觀察。眼圖是一系列數(shù)字信號在示波器上累積而顯示的圖形，因形狀像人眼睛而得名。在工程上，為了便于觀察接收波形中的碼間串擾和噪聲的情況，可在采樣判決設(shè)備的輸入端口處以恢復的采樣時鐘作為同步，用示波器觀察該端口的接收波形。Simulink眼圖模型參數(shù)設(shè)置如圖4所示。

圖4 眼圖模塊參數(shù)設(shè)置

眼圖關(guān)鍵參數(shù)設(shè)置如下：每次掃描顯示（Sym?bols per trace）的碼元個數(shù)若設(shè)置為3，則顯示3只眼睛；每個碼元的采樣點數(shù)（Samples per symbol）設(shè)置為10個；并調(diào)整采樣延時偏移量（Sample offset），使圖位于觀察范圍的正中央。

接收端采樣判決時需要提供位同步時鐘脈沖，由時鐘提取子系統(tǒng)來恢復，其內(nèi)部結(jié)構(gòu)如圖5所示。

圖5 位同步時鐘提取子系統(tǒng)

其中，Abs求模模塊的作用是倍頻。由于信源的傳碼率為1000 Baud，故IIR帶通濾波器的中心頻率應為2000 Hz，其通頻帶設(shè)置為1900～2100 Hz。采用鎖相環(huán)來鎖定定時脈沖的二次諧波分量，再通過二分頻得出位定時脈沖。故壓控振蕩器VCO的靜態(tài)頻率可設(shè)置為1995 Hz左右，靈敏度8 Hz/V，輸出信號幅度為1 V，采樣時間同系統(tǒng)仿真步長。計數(shù)器最大值設(shè)置為1，用于二分頻。延時模塊用于調(diào)整時鐘起始時刻對準眼圖的最佳采樣時刻。

接收端數(shù)據(jù)恢復部分由乘法器、滯環(huán)比較器（Relay）和觸發(fā)子系統(tǒng)（Triggered Subsystem）構(gòu)成，分別實現(xiàn)采樣、判決以及數(shù)據(jù)脈沖的再生。由于基帶系統(tǒng)傳輸?shù)碾p極性信號，所以Relay模塊的判決電平設(shè)置為0，大于0判為1，否則判為0。觸發(fā)子系統(tǒng)設(shè)置為上升沿觸發(fā)，用來恢復NRZ矩形脈沖數(shù)據(jù)。

3 基帶傳輸系統(tǒng)結(jié)果測試

3.1 仿真環(huán)境說明

本次設(shè)計在Matlab 2014b版本的仿真環(huán)境中進行。仿真求解器采用固定步長方式，步長為1e-4 s。由于發(fā)送濾波器、接收濾波器以及采樣判決過程均存在一定的延時，為使示波器各個窗口的信號波形能精確對齊，需要精心設(shè)置各部分觀測點信號的延時量。各點延時量具體值見圖2中5個Delay模塊所示。

3.2 基帶傳輸系統(tǒng)測試結(jié)果

圖2為無ISI基帶傳輸系統(tǒng)模型某次仿真所得結(jié)果?，F(xiàn)將該次仿真結(jié)果補充說明如下：隨機數(shù)發(fā)生器中高斯噪聲Var=0.01 W，Matlab仿真時間設(shè)置為20 s（200000個仿真步長），所得誤碼率Pe=0.007005。

（1）功能分析

基帶系統(tǒng)模型仿真結(jié)果的各觀測點信號均已送至示波器進行顯示，示波器各信號波形如圖6所示。

圖6 基帶傳輸系統(tǒng)各點信號波形

現(xiàn)對基帶傳輸系統(tǒng)各點信號波形進行分析。示波器設(shè)置信號波形顯示的時間范圍是0～0.01 s。示波器共設(shè)置7個顯示窗口。第1個波形為系統(tǒng)信號源產(chǎn)生的單極性非歸零矩形脈沖（NRZ）。由于碼率為1000 Baud，因此一個碼元持續(xù)時間（碼元周期）為0.001 s，所以圖6第一個窗口顯示的NRZ個數(shù)剛好10個，分別表示二進制的1010011001。第2個波形為對應的雙極性歸零矩形脈沖（BRZ）。第3個波形為經(jīng)發(fā)送濾波器變換后的升余弦信號，滾降系數(shù)為0.5。第4個波形為疊加上高斯白噪聲后的升余弦信號，即經(jīng)信道傳輸?shù)男盘?。?個波形為接收濾波器輸出信號。第6個波形為位同步時鐘信號，對比發(fā)送端的BRZ波形，同步保持時間約為6個碼元持續(xù)時間，有4個時鐘同步脈沖出現(xiàn)微小偏差（見圖6中紅色箭頭所示）。當時鐘同步出現(xiàn)偏差的值大于一個閾值后，同步提取子系統(tǒng)中的鎖相環(huán)開始工作以獲取同步時鐘。第7個波形為經(jīng)采樣判決后恢復的NRZ矩形脈沖。在這一階段雖然出現(xiàn)微小的時鐘偏差，但接收端仍能夠正確無誤地恢復信源發(fā)送的數(shù)據(jù)。

（2）眼圖分析

將接收端經(jīng)接收濾波器濾波之后的升余弦信號進行眼圖觀察和分析，所得眼圖波形如圖7所示。

圖7 接收信號眼圖波形分析

從圖7可知，當方差Var=0.01 W時，接收信號的眼圖規(guī)整，信號線跡清晰。圖7(a)中，中央的紅色箭頭表示該碼元的最佳采樣判決時刻。當噪聲方差增大至1 W時，眼圖線跡模糊，此時將導致接收端恢復數(shù)據(jù)時誤碼率急劇增大。

（3）誤碼率分析

誤碼率是數(shù)字通信系統(tǒng)一個極為重要的性能指標。仿真過程中發(fā)現(xiàn)，單次啟動模型仿真時，不同的仿真時間會得出差別較大的誤碼率值，這是由Matlab仿真環(huán)境本身的影響導致的。每一次仿真結(jié)果需要滿足一定的數(shù)據(jù)點數(shù)才能得出較為精確、可靠的結(jié)果，故需研究系統(tǒng)不同的單次仿真時間對此基帶系統(tǒng)抗噪性能的影響。故設(shè)置噪聲方差Var=0.01 W，改變不同的單次仿真時間t，執(zhí)行一次系統(tǒng)仿真所得誤碼率Pe的結(jié)果如圖8。

圖8 單次仿真時間與基帶系統(tǒng)誤碼率的關(guān)系

從圖8可知，系統(tǒng)單次仿真時間需要設(shè)置大于10 s，誤碼率才能小于7e-3。另外，當單次仿真時間取15 s、20 s和25 s這三個值時，系統(tǒng)的誤碼率相當接近；并且當單次仿真時間設(shè)為25 s時的誤碼率反而比單次仿真時間為20 s的誤碼率更大了?？紤]到系統(tǒng)仿真分析的效率，故選定此基帶系統(tǒng)單次仿真時間為20 s，并在此基礎(chǔ)上進一步研究此基帶系統(tǒng)的誤碼率與噪聲方差的關(guān)系。

在Matlab中可通過錯誤率統(tǒng)計模塊（Error Rate Calculation）繪制此無ISI基帶傳輸系統(tǒng)的誤碼率與噪聲方差的關(guān)系曲線。為了獲得此曲線圖，需要編寫一個M文件，通過運行此M文件多次調(diào)用系統(tǒng)仿真模型，從仿真結(jié)果中得到不同噪聲方差值時的誤碼率，從而繪制出誤碼率曲線[6-7]。相應的M文件代碼如下：

程序執(zhí)行之前，需要將隨機數(shù)發(fā)生器中噪聲的方差Var設(shè)為變量sigma2。并將錯誤率統(tǒng)計模塊輸出的誤碼率值以變量Error形式傳輸?shù)組atlab的工作空間。

M文件程序執(zhí)行完畢，獲得不同噪聲方差值時的誤碼率曲線，相關(guān)曲線圖如9。

圖9 不同噪聲方差與誤碼率曲線（模型單次仿真時間t=20 s）

從噪聲方差與誤碼率曲線來看，通過程序調(diào)用基帶傳輸系統(tǒng)的Simulink模型進行仿真所得的誤碼率要比單次執(zhí)行基帶傳輸系統(tǒng)的Simulink模型進行仿真所得的誤碼率大，例如當噪聲方差Var=0.01 W時，誤碼率曲線圖中的誤碼率Pe≈0.011；而單次執(zhí)行Simulink模型所得的誤碼率Pe=0.007005。對比分析原因有二：其一是程序調(diào)用方式得到的誤碼率是平均值，而單次執(zhí)行Simulink模型所得的誤碼率是最后較穩(wěn)定較精確的值。其二，通過程序調(diào)用基帶傳輸系統(tǒng)的Simulink模型進行仿真所得的誤碼率要通過變量返回到Matlab工作空間，再調(diào)用畫圖語句畫圖所得，這過程也導致誤碼率增加。但是通過程序調(diào)用基帶傳輸系統(tǒng)的Simulink模型進行仿真獲得誤碼率方式具有明顯的工程指導意義，從圖9中的（a）和（b）來看，噪聲方差Var<0.1 W時，誤碼率較小，且變化平緩。從（c）圖中可知，噪聲方差Var<0.7 W時，誤碼率Pe<0.02，而且變化也較平緩；當Var>0.7 W時，誤碼率Pe隨噪聲增大而明顯地增大。

4 結(jié)論

首先闡述了數(shù)字基帶傳輸系統(tǒng)的組成以及相關(guān)原理。其次利用Matlab/Simulink仿真平臺，對無ISI基帶傳輸系統(tǒng)進行系統(tǒng)建模、功能仿真和抗高斯噪聲的性能分析。系統(tǒng)仿真步長為1e-4 s，信源傳碼率為1000 Baud。

仿真結(jié)果表明：在高斯噪聲方差Var=0.01 W條件下，Matlab單次仿真的最佳時間為20 s，此時仿真數(shù)據(jù)點數(shù)充足，所得誤碼率Pe=0.007005。并且在Var=0.01 W條件下，接收信號的線跡清晰，眼圖規(guī)整。從系統(tǒng)的噪聲方差與誤碼率曲線來看，當Var<0.1 W時，誤碼率小，且變化平緩；當Var<0.7 W時，誤碼率Pe<0.02，而且變化也較平緩；當Var>0.7 W時，誤碼率Pe隨噪聲增大而明顯地增大；此結(jié)果在工程領(lǐng)域的參數(shù)計算及性能優(yōu)化方面極具指導價值。