北郵通原gmsk實(shí)驗(yàn)(編輯修改稿)
2025-05-14 00:00 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 化循環(huán)次數(shù)，象限數(shù)，及路徑等，設(shè)立兩級(jí)循環(huán)計(jì)算每次輸入情況下的相位值（注意此處連續(xù)相位），進(jìn)行采樣八次畫(huà)出一條正弦值的圖，然后移位bn數(shù)組及改變其他輸入量，重新畫(huà)出下一次的正弦值的圖，循環(huán)多次在屏幕上畫(huà)出一個(gè)眼圖，其代碼實(shí)現(xiàn)如下：①a=randint(1,1024)。 %等概率生成 0/1 b=2*a1。screen=40。%刷屏 screengmsk=zeros(1,80)。 %初始化 pathtemp=[b(1),b(2),b(3),b(4),b(5 )]。%賦初值%LL 表示絕對(duì)的 LL=0。 %模4后的L for i=1:screenfor k=1:10for j=1/8:1/8:1 %采樣 8 次gmsk((k1)*8+j*8)=sin(Ph(j,temp,L))。%不用加上因?yàn)槟?4 而帶來(lái)的 相位上的損失endL=mod(L+temp(1),4)。 %因?yàn)閷?xiě)的函數(shù) Phi是要輸入模 4的 L 的,所以也要同時(shí)求出for kk=1:4 %移位,得到新的 temptemp(kk)=temp(kk+1)。temp(5)=b((i1)*10+k+5)。 %給 temp(5)賦值,得到完整的 temp endendplot(gmsk)。title(39。GMSK 基帶正弦信號(hào)眼圖39。)。 hold on。endhold off。繪圖如下：圖10 量化前GMSK基帶余弦，正弦信號(hào)眼圖（余）弦表的均勻量化及Bin文件的生成 對(duì)正（余）弦表的量化范圍為[1,1],且每個(gè)樣值量化為 10bit,故量化為1024 個(gè)電平。若采用均勻量化,則量化間隔為 ,其誤差較小,不足以對(duì)整個(gè)編解碼產(chǎn)生不可恢復(fù)的影響，且考慮到其實(shí)現(xiàn)簡(jiǎn)單的因素，所以采用均勻量化,既能滿足實(shí)驗(yàn)對(duì)精確度的要求,又將運(yùn)算復(fù)雜度降低到了能夠接受的程度。具體步驟如下:1)對(duì)正 (余) 弦表中1024個(gè)點(diǎn)的每個(gè)樣值進(jìn)行10bit 量化。2)實(shí)驗(yàn)板上只有插接 8KB 的 ROM,所以 ROM 的數(shù)據(jù)位只有 8 位,共有 3 塊 ROM。存儲(chǔ)方式是: ROM1存儲(chǔ) 量化值的低 8 位；ROM3存儲(chǔ) 量化值的低8位；ROM2存儲(chǔ)和的高兩位；由于ROM每次讀8位，所以存儲(chǔ)也要按此方式，其代碼如下①M(fèi)=1024。 %量化delta=2/M。 quancos=zeros(1,8*1024)。quansin= zeros(1,8*1024) 。 quansc=zeros(1,8*1024)。for n=1:1024qcos=floor((cos(phi(n))+1)/delta)。%量化后的 cosqsin=floor((sin(phi(n))+1)/delta)。%量化后的 sincosh=floor(qcos/256)。%cos 高兩位cosl=floor(mod(qcos,256))。%cos 低 8 位quancos(n*87)=cosl。 sinh=floor(qsin/256 )。sinl= floor(mod(qsin,256))。 quansin(n*87)=sinl。end②cosfile= fopen(39。39。,39。w39。)。 %寫(xiě)入bin文件fwrite(cosfile,quancos)。fclose(cosfile)。sinfile= fopen(39。39。,39。w39。)。fwrite(sinfile,quansin)。fclose(sinfile)。scfile= fopen(39。39。,39。w39。)。 fwrite(scfile,quansc)。 fclose(scfile )。 圖11 SIN低八位bin表 圖12 COS低八位bin表 圖13 SIN和COS的高位繪制原理與量化前眼圖原理相同，區(qū)別在于繪制時(shí)正弦和余弦值要通過(guò)查表得到，其代碼實(shí)現(xiàn)和繪制眼圖如下：①M(fèi)=1024。delta=2/M。 cos_table_lh(1024)=0。sin_table_lh(1024)=0。for n=1:1024qcos=floor((cos(phi(n))+1)/delta)。%量化后的 cos cos_table_lh(n)=qcos*delta1。qsin=floor((sin(phi(n))+1)/delta)。%量化后的sin sin_table_lh(n)=qsin*delta1。enda=randint(1,1024)。 %等概率生成 0/1b=2*a1。screen=40。%刷屏 screengmsk=zeros(1,80)。 %初始化 pathtemp=[b(1),b(2),b(3),b(4),b(5)]。%賦初值L=0。 %模4后的Lc(5)=0。 for i=1:screen for k=1:10temp(5)=b((i1)*10+k+5)。 %給 temp(5)賦值,得到完整的 tempfor j=0:1:7 %采樣 8 次for x=1:5if(temp(x)==1) c(x)=0。else c(x)=1。endendn=16*c(1)+8*c(2)+4*c(3)+2*c(4)+c(5)+1。%第 n個(gè)狀態(tài)gmsk((k1)*8+j+1)=cos_table_lh((n1)*32+L*8+j+1)。 % endL=mod(L+temp(1),4)。%因?yàn)閷?xiě)的函數(shù) Phi是要輸入模 4 的L 的,所以也要同時(shí)求出 for kk=1:4 %移位,得到新的 temptemp(kk)=temp(kk+1)。 end endplot(gmsk)。title(39。量化后 GMSK 基帶余弦信號(hào)眼圖39。)。 hold on。 end hold off。圖14 量化后 GMSK 基帶余弦信號(hào)、正弦信號(hào)圖 mif文件的生成for i=1:8192 data(i)=quansc(i)。endfid=fopen(39。/Users/zhouliang/Desktop/GMSK/39。,39。w39。)。fprintf(fid,39。WIDTH=
教學(xué)教案相關(guān)推薦
鄂ICP備17016276號(hào)-1