無線實(shí)驗(yàn)指導(dǎo)書(編輯修改稿)
2024-08-31 01:15 本頁面
【文章內(nèi)容簡(jiǎn)介】 K信號(hào)相干解調(diào)器構(gòu)成，其原理框圖如圖15所示。S1(t)低通濾波判決低通濾波判決Cos(w0t)sin(w0t)逆映射并串轉(zhuǎn)換基帶信號(hào)I(t)Q(t)圖15 QPSK解調(diào)原理框圖數(shù)學(xué)推導(dǎo)調(diào)制信號(hào)為：調(diào)制信號(hào)經(jīng)過高斯信道為：I路相干解調(diào)：Q路相干解調(diào)：由公式推導(dǎo)可以知道，只需將)，通過一個(gè)具有低通特性的系統(tǒng)即可獲得只包含I和Q信息的信號(hào)，實(shí)現(xiàn)解調(diào)，這里選用了巴特沃斯型的低通濾波器。附程序代碼%日期 %功能 QPSK的調(diào)制解調(diào)，基帶信號(hào)點(diǎn)數(shù)t（限偶數(shù)），基波頻率w0可設(shè)置clear all。nb=32。 % 傳輸?shù)谋忍財(cái)?shù)T=1。 % 基帶信號(hào)寬度，也就是基波頻率fc=8/T。 % 載波頻率ml=2。 % 調(diào)制信號(hào)類型的一個(gè)標(biāo)志位（選取2的原因見23行）c = 4*nb。 %單周期采樣點(diǎn)數(shù)delta_T=T/c。 % 采樣間隔fs=1/delta_T。 % 采樣頻率t=0:delta_T:nb*Tdelta_T % 限定t的取值范圍 c * nbN=length(t)。 % 采樣數(shù) %%%%%%%%%%%%%%%%%%% 調(diào)制部分% 基帶信號(hào)的產(chǎn)生data=fix(2*rand(1,nb))。 % 調(diào)用一個(gè)隨機(jī)函數(shù)（0 or 1），輸出到一個(gè)1*100的矩陣datanrz=data.*21。 % 變成極性碼for i=1:nb data1((i1)/delta_T+1:i/delta_T)=datanrz(i)。 % 將極性碼變成對(duì)應(yīng)的波形信號(hào)end % 將基帶信號(hào)變換成對(duì)應(yīng)波形信號(hào)for i=1:nb data0((i1)/delta_T+1:i/delta_T)=data(i)。 % 將基帶信號(hào)變成對(duì)應(yīng)的波形信號(hào)end % 串并轉(zhuǎn)換，將奇偶位數(shù)據(jù)分開idata=datanrz(1:ml:(nb1))。 % 將奇偶位分開，因此間隔m1為2 qdata=datanrz(2:ml:nb)。% QPSK信號(hào)的調(diào)制for i=1:nb/2 ich(2*((i1)/delta_T+1):2*(i/delta_T))=idata(i)。endfor ii=1:N/T a(ii)=(1/sqrt(2))*cos(2*pi*fc*t(ii))。 endidata1=ich.*a。 % 奇數(shù)位數(shù)據(jù)與余弦函數(shù)相乘，得到一路的調(diào)制信號(hào)for j=1:nb/2 qch(2*((j1)/delta_T+1):2*(j/delta_T))=qdata(j)。endfor jj=1:N/T b(jj)=(1/sqrt(2))*sin(2*pi*fc*t(jj))。endqdata1=qch.*b。% 偶數(shù)位數(shù)據(jù)與余弦函數(shù)相乘，得到另一路的調(diào)制信號(hào)st = idata1 qdata1。%%%%%%%%%%%%%%%%%%%%%%%%信道中SNR=0。 % 信噪比stn = awgn(st,SNR)。%%%%%%%%%%%%%%%%%%%%%%%%解調(diào)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%設(shè)計(jì)濾波器%%%%%%%%[B,A] = butter(3,39。low39。)。[h1,w] = freqz(B,A)。%%%%%%%%%%%%%%%%%%%%%%%%%%%%相干解調(diào)ist = stn .* a。p =length(ist)qst = stn .* (b)。%%%%%%%%%%%%%%%%%%%%濾波istl = filter(B,A,ist)。qstl = filter(B,A,qst)。%%%%%%%%%%%%%%%%%%%%%%%%%抽樣判決%%%%%%%%%%for i = 1 : nb/2 if istl(2*(p/nb)*(i1)+(1*(p/nb))) = 0 in(i) = 1。 else in(i) = 0。 end if qstl(2*(p/nb)*(i1)+(1*(p/nb))) = 0 qn(i) = 1。 else qn(i) = 0。 endend%%%%%%%%%%%%%%%%%并串轉(zhuǎn)換for i = 1 : nb/2 y(2*i1) = in(i)。 y(2*i) = qn(i)。endfor i = 1 : nb yy((i1)/delta_T+1:i/delta_T) = y(i)。enddatayNfigure。subplot(4,1,1)plot(data0*),title(39。基帶信號(hào)，4096維二進(jìn)制序列，對(duì)應(yīng)向量是data039。)。subplot(4,1,2)plot(data1*),title(39。雙極性信號(hào)，4096維雙極性序列，對(duì)應(yīng)向量是data139。)。subplot(4,1,3)plot(ich*),title(39。I路數(shù)據(jù)，4096維雙極性序列，對(duì)應(yīng)向量是ich39。)。subplot(4,1,4)plot(qch*),title(39。Q路數(shù)據(jù)，4096維雙極性序列，對(duì)應(yīng)向量是qch 39。)。figure。subplot(4,1,1)plot(ist),title(39。相干解調(diào)I路信號(hào)，4096維且值為(1,1)的序列，對(duì)應(yīng)向量是ist 39。)。subplot(4,1,2)plot(qst),title(39。相干解調(diào)Q路信號(hào)，4096維且值為(1,1)的序列，對(duì)應(yīng)向量是qst 39。)。subplot(4,1,3)plot(istl),title(39。I路解調(diào)波形，4096維且值為(1,1)的序列，對(duì)應(yīng)向量是istl 39。)。subplot(4,1,4)plot(qstl),title(39。Q路解調(diào)波形，4096維且值為(1,1)的序列，對(duì)應(yīng)向量是qstl 39。)。%%%%%%%%%%%%%%%%%畫圖%%%%%%%%%%%%%%%%%%%%%%%%%%%figure。subplot(4,2,1)。plot(data0*),title(39?；鶐盘?hào)39。)。subplot(4,2,2)。psd(abs(fft(data0))),title(39?；鶐盘?hào)頻譜39。)。subplot(4,2,3)。plot(st),title(39。調(diào)制信號(hào)39。)。subplot(4,2,4)。psd(abs(fft(st))),title(39。調(diào)制信號(hào)頻譜39。)。subplot(4,2,5)。pl
試題試卷相關(guān)推薦
鄂ICP備17016276號(hào)-1