【正文】 end endend%%零內(nèi)插zero=5。 else I_diff(i)=1。Q_diff(i1)==1 I_diff(i)=1。 elseif I_diff(i1)==1amp。Q_diff(i1)==1 I_diff(i)=1。 elseif I_diff(i1)==1amp。Q_diff(i1)==1 I_diff(i)=1。 end else if I_diff(i1)==1amp。 else I_diff(i)=1。Q_diff(i1)==1 I_diff(i)=1。 elseif I_diff(i1)==1amp。Q_diff(i1)==1 I_diff(i)=1。 elseif I_diff(i1)==1amp。Q_diff(i1)==1 I_diff(i)=1。Q(i)==1 if I_diff(i1)==1amp。 end elseif I(i)==1amp。 else I_diff(i)=1。Q_diff(i1)==1 I_diff(i)=1。 elseif I_diff(i1)==1amp。Q_diff(i1)==1 I_diff(i)=1。 elseif I_diff(i1)==1amp。Q_diff(i1)==1 I_diff(i)=1。Q(i)==1 if I_diff(i1)==1amp。 end elseif I(i)==1amp。 else I_diff(i)=1。Q_diff(i1)==1 I_diff(i)=1。 elseif I_diff(i1)==1amp。Q_diff(i1)==1 I_diff(i)=1。 elseif I_diff(i1)==1amp。Q_diff(i1)==1 I_diff(i)=1。Q(i)==1 % 0+1 11 if I_diff(i1)==1amp。for i=2:data_len if I(i)==1amp。I_diff(1)=I(1)。 end endend%% 碼變換I_diff=zeros(1,data_len)。 else Q(i1)=1。 end else if input(i)==1 Q(i1)=1。 else I(i)=1。 endend%%串并轉(zhuǎn)換for i=1:data_len if rem(i,2)==1 if input(i)==1 I(i)=1。for i=1:data_len if rand_data_len(i)= input(i)=1。data_len=10000。clc。%QDPSK 誤碼率clear all。QPSK 通過 AWGN 信道39。)。 ylabel(39。SNR (dB)39。)。,39。 %不同顏色畫出誤比特率曲線對比legend(39。r39。)。endfiguresemilogy(SNR,BERtheory,39。 % 信噪比 BERtheory(i)=1(*erfc(sqrt(SNR1(i)/2)))^2。 % 解調(diào)[nErrors, BER(n)] = biterr(x(1:end),r1(1:end))。)。 % 信噪比范圍for n = 1:length(SNR)r = awgn(s,SNR(n),39。 % 隨機(jī)序列s= pskmod(x,M)。M = 4。)%QPSK 誤碼率clear all。 title(39。 plot_Pxx=10*log10(Pxx(index+1))。 index=0:round(nfft/21)。 %計算序列的自相關(guān)函數(shù) CSk=fft(cs_n,nfft)。unbiased39。nfft=1024。p1=10。QPSK 功率譜密度39。 figureplot(k,plot_Pxx)。 k=index*Fs/nfft。 %傅里葉變換Pxx=abs(CSk)。)。 cs_n=xcorr(s,39。Fs=1000。)%QPSK 功率譜密度m=1。subplot(4,1,4)plot(t,s)axis([t(1) t(end) ])title(39。s2=mQmap.*c2。c2=sin(2*pi*fc*tpi/2)。QDPSK 映射后的 Q 路信號39。QDPSK 映射后的 I 路信號39。,1)。mQmap=pulstran(,dQmap,39。rectpuls39。Q]39。I]39。Q(pQ==0)=+1/sqrt(2)。I(pI==0)=+1/sqrt(2)。,1)。%mQ=pulstran(,dQ,39。rectpuls39。pI]39。pQ]39。d2=[0:1:t(end)1]。pQ=pQI(2,:)。基帶信號 p39。subplot (4,1,1)。rectpuls39。p1]39。 %每個編碼的發(fā)送延遲時間p1=[p, zeros(length(d1)length(p))]。p=[1 0 1 1 0 1 0 0 1 1]。f=n*df。df=fs/(N1)。 %時間坐標(biāo) sN=length(t)。%全局變量Ts=。clear all。QPSK 進(jìn)入信道的信號 s=s1s239。s=s1s2。%已調(diào)信號s1=mImap.*c1。)%載波信號fc=10 %載波頻率c1=cos(2*pi*fc*t)。)subplot(4,1,3)plot(t,mQmap)axis([t(1) t(end) ])grid ontitle(39。 subplot(4,1,2)plot(t,mImap)axis([t(1) t(end) ])grid ontitle(39。rectpuls39。,1)。mImap=pulstran(,dImap,39。dQmap=[d2。dImap=[d2。Q(pQ==1)=1/sqrt(2)。 %IQ 值映射 “0”映射成1， “1”映射成 1I(pI==1)=1/sqrt(2)。rectpuls39。,1)。%mI=pulstran(,dI,39。%dI=[d2。%IQ 信號波形%dQ=[d2。pI=pQI(1,:)。)%抽取 IQ 信號 QPSK 兩個基帶符號為一組，映射到一對 IQpQI=reshape(p,2,5)。plot(t,m)axis([t(1) t(end) ])grid ontitle(39。,)。 %產(chǎn)生 d 矩陣m=pulstran(,d,39。 %在發(fā)送串后補(bǔ)零d=[d1。 %待傳送的編碼串d1=[0::t(end)]。 %頻率坐標(biāo)%基帶信號figure。 %df 為頻率分辨率（頻率間隔）n=(N1)/2:(N1)/2。fs=1/Ts。 %時間分辨率(時點(diǎn)間間隔） ，抽樣時間間隔t=0:Ts:5。close all。 最后衷心感謝電氣學(xué)院的所有給予我指導(dǎo)和幫助的領(lǐng)導(dǎo)、老師和同學(xué)，衷心的感謝你們！附錄主程序：%QPSK 信號波形clc。在論文完成之際，首先衷心感謝林老師對我的深切關(guān)懷和諄諄教導(dǎo)。希望經(jīng)過這次設(shè)計，自己能夠取長補(bǔ)短，為以后的學(xué)習(xí)和工作積累寶貴的經(jīng)驗(yàn)。QPSK 和 QDPSK 信號波形和功率譜密度的仿真相對順利，但是對于誤碼率的仿真則較難得到理想結(jié)果，這同時也能反映出自身對誤碼率相關(guān)知識掌握的欠缺。仿真的過程大致上分為設(shè)計仿真思路、編程調(diào)試、得到仿真結(jié)果并分析三大步。本次設(shè)計課題內(nèi)容涉及的 QPSK 和 QDPSK 調(diào)制都是在專業(yè)課學(xué)習(xí)當(dāng)中基礎(chǔ)且重要的知識內(nèi)容，要想對其進(jìn)行進(jìn)一步深入的模擬仿真和性能比較，就不僅僅局限于課本上的簡單介紹的理論知識了。第五章 結(jié)束語QPSK 調(diào)制解調(diào)和 QDPSK 調(diào)制解調(diào)因其相對較高的性能而廣泛應(yīng)用于數(shù)字通信的調(diào)制解調(diào)當(dāng)中。并且對于 QPSK 系統(tǒng)隨著信噪比的增大仿真值越接近理論值。0 2 4 6 8 10 12 14 16 18 2022151010105100SNR (dB)BERQPSK仿QDPSK仿仿仿仿仿 QPSK BERQDPSK BER圖 49 QPSK 與 QDPSK 誤碼率曲線比較QPSK 與 QDPSK 誤碼率比較仿真條件輸入 100000 個隨機(jī)數(shù)，信噪比范圍在 020dB 下，QPSK 系統(tǒng)與 QDPSK 系統(tǒng)的誤碼率比較QPSK 與 QDPSK 誤碼率比較從圖 49 可以看出，對于 QPSK 與 QDPSK 系統(tǒng)這兩個系統(tǒng)的誤碼率都是隨著信噪比的增大而減小，但是同一信噪比下，QPSK的誤碼率比 QDPSK 的誤碼率要小，要使 QDPSK 的誤碼率與 QPSK 的誤碼率相同，要通過增大 QDPSK 系統(tǒng)的信噪比才可以達(dá)到，因此 QPSK 的抗噪聲性能要優(yōu)于 QDPSK 系統(tǒng)。在得出的理論值與仿真值做比較上還是有差別的。0 2 4 6 8 10 12 14 16 18 2022151010105100SNR (dB)BERQDPSK仿仿AWGN仿仿 仿仿BER仿仿 BER 圖 48 QDPSK 誤碼率曲線QDPSK 系統(tǒng)仿真條件：即對于進(jìn)制數(shù) M=4 時，信噪比范圍在 020dB 下，輸入 40000 個隨機(jī)數(shù)，通過蒙特卡洛仿真方法，得出 QDPSK 系統(tǒng)的誤碼率曲線。0 50 100 150 200 250 300 350 400 450 500403020100102030 QPSK仿仿仿仿仿圖 44QPSK 功率譜密度 圖 45QDPSK 功率譜密度分析 QPSK 系統(tǒng)和 QDPSK 系統(tǒng)的誤碼率,仿真實(shí)現(xiàn)采用蒙特卡洛方法，仿真思路如圖 46 所示 r c 2 比特符號 n s均勻隨機(jī)數(shù)發(fā)生器高斯隨機(jī)數(shù)發(fā)生器QPSK 或QDPSK 映射檢測器高斯隨機(jī)數(shù)發(fā)生器比較比特差錯計數(shù)器符號差錯計數(shù)器圖 46 誤碼率仿真框圖0 2 4 6 8 10 12 14 16 1810151010105100SNR (dB)BERQPSK仿仿AWGN仿仿 仿仿BER仿仿 BER圖 47 QPSK 誤碼率曲線QPSK 系統(tǒng)仿真條件：即對于進(jìn)制數(shù) M=4 時，信噪比范圍在 030dB 下，輸入 100000 個隨機(jī)數(shù)，通過蒙特卡洛仿真方法，得出 QPSK 系統(tǒng)的誤碼率曲線。除內(nèi)部函數(shù)以外，所有MATLAB 的核心文件和工具箱文件都是可讀可改的源文件。 MATLAB 包含兩個部分：核心部分和各種可選工具箱。由于MATLAB 的程序不用編譯等預(yù)處理，也不生成可執(zhí)行文件，程序?yàn)榻忉寛?zhí)行，所以速度較慢。MATLAB 還具有較強(qiáng)的編輯圖形界面的能力。MATLAB 的圖形功能強(qiáng)大。由于 MATLAB 是用 C 語言編寫的， MATLAB 提供了和 C 語言幾乎一樣多的運(yùn)算符，靈活使用 MATLAB 的運(yùn)算符將使程序變得極為簡短。MATLAB 程序書寫形式自由，利用其豐富的庫函數(shù)避開繁雜的子程序編程任務(wù)，壓縮了一切不必要的編程工作。 以 下 簡 單 介紹 一 下 MATLAB 的 主 要 特 點(diǎn) 。 MATLAB 語 言 是 一 種 高 級 的 基 于 矩 陣 /數(shù) 組 的 語 言 ， 其 最 突 出 的 特 點(diǎn) 就是 簡 潔