通信原理gmsk實(shí)驗(yàn)報(bào)告-文庫(kù)吧
2025-07-08 17:18 本頁(yè)面
【正文】 編寫計(jì)算cos?t及sin?t的程序，并設(shè)計(jì)余弦及正弦ROM表。 將余弦和正弦碼表中的每個(gè)樣值的10bit碼字，按照地址邏輯進(jìn)行存放， 并用matlab程序?qū)崿F(xiàn)。 得到ROM存儲(chǔ)的基本波形表，利用信號(hào)源產(chǎn)生的偽隨機(jī)序列｛an｝經(jīng)預(yù)編碼后得到｛bn｝，經(jīng)過地址邏輯運(yùn)算由ROM中順序取出cos?t及sin?t的離散值，然后利用計(jì)算機(jī)繪圖程序（功能相當(dāng)于數(shù)模變換DAC），得到基帶波形的輸出，觀察仿真眼圖。 把得到的正余弦表進(jìn)行數(shù)字量化，寫入BIN文件，下載到硬件系統(tǒng)中，通過示波器觀察實(shí)際硬件實(shí)現(xiàn)的GMSK信號(hào)眼圖。五、 系統(tǒng)設(shè)計(jì)總體設(shè)計(jì)軟件部分 系統(tǒng)軟件流程圖如下圖所示。（1）gt函數(shù)的產(chǎn)生思路：了解到Matlab內(nèi)置有erfc函數(shù)，可以通過erfc(x)的表達(dá)式，得到Q(t)的表達(dá)似式，從而根據(jù)設(shè)定的時(shí)間等參數(shù)得到gt的表達(dá)式，然后在[,]內(nèi)用linspace取1000個(gè)點(diǎn)作圖可以得到gt的圖像。具體編程實(shí)現(xiàn)如下（）：代碼：function gt = gt_m( t )%計(jì)算g(t)clc。clear。T=1/270833。B=。a=1/B*sqrt(log(2)/2)。t=linspace(*T,*T,1000)。gt=1/(2*T)*1/2*(erfc(pi/a*(tT/2))erfc(pi/a*(t+T/2)))。plot(t,gt)。title(39。 BT=，高斯濾波器矩形脈沖響應(yīng)g(t)39。)。xlabel(39。t(限定在5T范圍內(nèi))39。)。ylabel(39。g(t)39。)。end可以得到圖像如下圖所示：（2）相位路徑?t與cos?t及sin?t的碼表的產(chǎn)生 思路：首先由于每一個(gè)相位路徑是由5個(gè)碼元所組成的，先要編寫四個(gè)象限的｛bn｝表，將所有情況全部羅列出來。然后根據(jù)相位路徑產(chǎn)生的公式，借助matlab中的積分quad函數(shù)，計(jì)算相位路徑?t，最后得到cos?t及sin?t，分別用一個(gè)1*1024的矩陣來表示。具體編程實(shí)現(xiàn)如下（）。 流程圖如下所示：代碼：clc。clear。T=1/270833。B=。a=sqrt(log(2)/2)/B。t=linspace(*T,*T,1024)。g=@(t)(*erfc(pi*(tT/2)/a)*erfc(pi*(t+T/2)/a))/(2*T)。 %g(t)b=zeros(32,5)。for n=0:31 %將n轉(zhuǎn)化為二進(jìn)制 m=n。 for j=4:1:1 if floor(m/2^j)==1 b(n+1,5j)=1,m=m2^j,end。 end。 b(n+1,5)=m。 for i=1:5 if b(n+1,i)==0 b(n+1,i)=1。 end。 end。end。pa=[]。for n=0:31 for L=0:3 fai=0。 for j=2:2 fai=fai+pi*b(n+1,j+3)*quad(g,*T,j*TT/2)。 end fai=fai+L*pi/2。 for i=0:7 delta=0。 for j=2:2 delta=delta+pi*b(n+1,j+3)*quad(g,j*TT/2,i*T/8j*TT/2)。 end pb=fai+delta。 pa=[pa,pb]。 end endend%plot(t,pa)s=sin(pa)。 %正弦表c=cos(pa)。 %余弦表save(39。39。,39。s39。)。save(39。39。,39。c39。)。(3)眼圖的繪制思路：得到ROM存儲(chǔ)的基本波形表后，利用信號(hào)源產(chǎn)生的偽隨機(jī)序列｛an｝經(jīng)預(yù)編碼后得到｛bn｝，經(jīng)過地址邏輯運(yùn)算由ROM中順序取出cos?t及sin?t的離散值，然后利用計(jì)算機(jī)繪圖程序（功能相當(dāng)于數(shù)模變換DAC），得到基帶波形的輸出，觀察仿真眼圖。流程圖如下：說明：A、 開辟空間存入ROM表以及偽隨機(jī)序列。B、根據(jù)所用的ROM表確定采樣頻率fs=8fb(fb=), 即一個(gè)碼元時(shí)間內(nèi)有8個(gè)抽樣值，每樣值的量化電平是10，并設(shè)初始相位是0，即起始地址是0000000000。C、確定讀取信號(hào)的數(shù)目為10*TIMES，即每十個(gè)碼元時(shí)間的波形顯示于同一屏幕。共顯示TIMES次。D、由bk2,bk1,bk,bk+1,bk+25個(gè)碼元及象限L形成地址邏輯獲得ROM表中的7位地址；E、再取三位地址碼，順序取出i=0,1,2,3,4,5,6,7個(gè)抽樣量化值，由10位地址邏輯找到ROM中基帶波形的位置，將其轉(zhuǎn)換成實(shí)際值，存入一個(gè)數(shù)據(jù)文件用于眼圖仿真，共有10*TIMES*n個(gè)值（n=8）；F、根據(jù)產(chǎn)生眼圖的原理，將每次掃描結(jié)果疊加而成，刪除程序中清屏幕命令，即可看到眼圖；代碼：close all。%仿真眼圖p。 %將ROM表導(dǎo)入m=2^231。 %偽隨機(jī)序列個(gè)數(shù)TIMS=500。 %顯示次數(shù)k=5。t=1:80。bcos=zeros(1,80)。bsin=zeros(1,80)。seq=round(rand(1,m))。 %產(chǎn)生偽隨機(jī)序列an=[0,seq(1:4)]。 %初始4位隨機(jī)碼figure(1)。for i=1:TIMS for j=0:9 an=[an(2:5),seq(k)]。 %取下一
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1