基于matlab的oqpsk調(diào)制解調(diào)實(shí)現(xiàn)課程設(shè)計(jì)-文庫(kù)吧
2025-05-14 08:56 本頁(yè)面
【正文】 圍是 0~3?；鶐盘?hào)產(chǎn)生后可以利用 fft(x,1024)對(duì)其進(jìn)行傅里葉變換，得到基帶信號(hào)的頻譜 圖。再利用函數(shù) oqpskmod(x)對(duì)基帶信號(hào)進(jìn)行調(diào)制，由于MATLAB中的 oqpskmod(x)函數(shù)功能有限，此時(shí)的調(diào)制信號(hào)還不能用波形圖表示出來(lái)，需要自己編寫部分代碼（見附錄），所以使用調(diào)制信號(hào)的星座圖代替調(diào)制信號(hào)圖，使用函數(shù) scatter(y)可生成星座圖。代碼編寫完成后即可得到相應(yīng)的 OQPSK調(diào)制波形圖，同樣使用 fft(n,1024)函數(shù)來(lái)產(chǎn)生波形信號(hào)的頻譜圖，其中 n為調(diào)制信號(hào)，即在調(diào)制信號(hào)中 《 基于 MATLAB的 OQPSK調(diào)制解調(diào)實(shí)現(xiàn) 》 第 5頁(yè) 共 22頁(yè) 取 1024個(gè)點(diǎn)進(jìn)行傅里葉變換。調(diào)制波形圖出來(lái)以后進(jìn)行加噪聲處理，此處利用的是awgn(n,snr)函數(shù)， snr 為信噪比的值，先設(shè)定 snr 為 1。隨后產(chǎn)生加入噪聲后的調(diào)制信號(hào)頻譜圖。下面先介紹幾個(gè)關(guān)鍵變量及重要函數(shù)，然后將展示出調(diào)制過(guò)程中出現(xiàn)的比較重要的圖形。在這里為了更清晰地看出調(diào)制后信號(hào)的波形圖、頻譜圖，及方便解調(diào)時(shí)對(duì)比加入噪聲后出現(xiàn)的錯(cuò)碼位置，只選取 8個(gè)基帶信號(hào)。 a=8。 %基帶信號(hào)碼元個(gè)數(shù) x = randint(a,1,[0 3])。 % 隨機(jī)產(chǎn)生 a個(gè)四進(jìn)制數(shù) y =oqpskmod(x)。% 進(jìn)行 oqpsk調(diào)制 scatterplot(y) %解調(diào)信號(hào)的星座圖 fc=100。 %載波頻率 fs=1600。 %抽樣頻率 b=0。 %設(shè)置初始相位為 0 s=0。 %調(diào)制信號(hào)的波形圖縱坐標(biāo)初始化 n=0。 %調(diào)制信號(hào)波形圖二維初始化 基帶信號(hào)（如圖 31）： 圖 31 基帶信號(hào) 基帶信號(hào)的頻譜圖如下（圖 32 ）： r=fft(x,1024)。 %取 1024個(gè)點(diǎn)對(duì)基帶信號(hào)進(jìn)行傅里葉變換 fr=(0:length(r)1)*fs/length(r)fs/2。 %橫坐標(biāo)為頻域 plot(fr,abs(r))。 % OQPSK基帶信號(hào)頻譜圖 《 基于 MATLAB的 OQPSK調(diào)制解調(diào)實(shí)現(xiàn) 》 第 6頁(yè) 共 22頁(yè) 圖 32基帶信號(hào)頻譜圖 函數(shù)調(diào)制信號(hào)的星座圖如下（圖 33）所示 ： scatterplot(y) %產(chǎn)生調(diào)制信號(hào)的星座圖 圖 33調(diào)制信號(hào)星座圖 函數(shù)調(diào)制信號(hào)的星座圖會(huì)隨輸入基帶信號(hào)初值的變化而改變，即不同的第一個(gè)信號(hào)碼元對(duì)應(yīng)不同的星座圖。因?yàn)?OQPSK信號(hào)是由兩路原本相互正交、但其中一路比另一路早半個(gè)周期的信號(hào)組成，這就使得第一路信號(hào)剛進(jìn) 去時(shí)和第二路信號(hào)最后進(jìn)去時(shí)系統(tǒng)中都只有一路信號(hào)，此時(shí)會(huì)出現(xiàn)純實(shí)數(shù)或純虛數(shù)，即會(huì)出現(xiàn)上面調(diào)制信號(hào)星座圖中非角落的點(diǎn)。 調(diào)制信號(hào)波形圖如下（圖 34）所示 ： m=(2*pi*fc*k)/fs+b。 %合相位 s=sin(m)。 %生成波形圖 n=[n sin(m)]。 %將一維轉(zhuǎn)化為二維 g=1:length(n)。 % g為調(diào)制波形圖的橫坐標(biāo) plot(g,n) %調(diào)制信號(hào)波形圖 《 基于 MATLAB的 OQPSK調(diào)制解調(diào)實(shí)現(xiàn) 》 第 7頁(yè) 共 22頁(yè) 圖 34 調(diào)制信號(hào)波形圖 加了噪聲之后的調(diào)制信號(hào)波形圖如下（圖 35）所示 ： snr=。 %信噪比 noi=n +awgn(n,snr)。%加入噪聲后的調(diào)制信號(hào) 圖 35 加噪聲后的調(diào)制信號(hào) 調(diào)制信號(hào)頻譜圖如下（圖 36）所示 ： s1=fft(n,1024)。 %去 1024個(gè)點(diǎn)對(duì)調(diào)制信號(hào)進(jìn)行傅里葉變化 fs1=(0:length(s1)1)*fs/length(s1)fs/2。 %橫坐標(biāo)為頻域 plot(fs1,abs(s1))。 %OQPSK調(diào)制信號(hào)頻譜圖 《 基于 MATLAB的 OQPSK調(diào)制解調(diào)實(shí)現(xiàn) 》 第 8頁(yè) 共 22頁(yè) 圖 36調(diào)制信號(hào)頻譜圖 調(diào)制信號(hào)加噪聲后的頻譜圖如下（圖 37）所示 ： s2=fft(noi,1024)。 %取 1024個(gè)點(diǎn)對(duì)加噪聲后的調(diào)制信號(hào)進(jìn)行傅里葉變換 fs2=(0:length(s2)1)*fs/length(s2)fs/2。 %橫坐標(biāo)為頻域 plot(fs2,abs(s2)) %加入噪聲后 OQPSK調(diào)制信號(hào)的頻譜圖 圖 37調(diào)制信號(hào)加入噪聲后的頻譜圖 OQPSK 解調(diào)實(shí)現(xiàn) 此步建立在已經(jīng)進(jìn)行完 OQPSK調(diào)制的基礎(chǔ)之上。在 OQPSK的解調(diào)過(guò)程中，我們使用的是 oqpskdemod(y)這個(gè)函數(shù)，與函數(shù) oqpskmod(x)一樣，這個(gè)函數(shù)的功能也非常有 《 基于 MATLAB的 OQPSK調(diào)制解調(diào)實(shí)現(xiàn) 》 第 9頁(yè) 共 22頁(yè) 限，僅是 oqpskmod(x)的逆過(guò)程，即 oqpskdemod(y)的輸入信號(hào)必須是 oqpskmod(x)的輸出信號(hào)，否則系統(tǒng)就會(huì)報(bào)錯(cuò)。在沒有加噪聲之前，僅依靠簡(jiǎn)單的 oqpskdemod(y)函數(shù)就可以將調(diào)制后的信號(hào)解調(diào)出來(lái)。 k=oqpskdemod(y)。 %解調(diào)信號(hào) 解調(diào)信號(hào)如下（圖 38）所示： 圖 38解調(diào)信號(hào) 解調(diào)信號(hào)的頻譜圖如下（圖 39）所示： jtpp=fft(k,1024)。 %解調(diào)信號(hào)的傅里葉變換 fj=(0:length(jtpp)1
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1