【導讀】調制后再傳輸的通信方式，在現代通信中有著十分重要的地位且應用廣泛。擾性、在電路上實現也較為簡單。量節(jié)省發(fā)送功率，設備簡單并易于制造等優(yōu)點。本文主要通過MATLAB軟件針對對QPSK和QDPSK調制系統(tǒng)搭建仿真平臺，從信號波形、功率譜密度與誤碼率三方面對二者展開討論比較。系統(tǒng)進行了仿真，并將仿真結果和相位理想狀態(tài)下的分析數據進行了對比。四進制為例對比MPSK與MDPSK系統(tǒng)的抗噪聲性能。

　　

【正文】 但是對于誤碼率的仿真則較難得到理想結果，這同時也能反映出自身對誤碼率相關知識掌握的欠缺。 本次畢業(yè)設計讓我收獲頗多，不僅加深了對數字通信、 QPSK 系統(tǒng)和 QDPSK系統(tǒng)相關知識的理解，鞏固了自己的專業(yè)知識，還鍛煉了綜合運用所學的知識來解決面臨的實際問題的能力 ；同時也反映出自己在專業(yè)知識的儲備方面還十分的欠缺。希望經過這次設計，自己能夠取長補短，為以后的學習和工作積累寶貴的經驗。 參考文獻 [1] 馬俊杰 .淺談數字通信的優(yōu)點以及應用 .價值工程 .2020 年 09 期： 145 [2] 曹志剛、錢亞生 .現代通信原理 .北京：清華大學出版社， 1992 [3] 姬宗嶺 .數字通信技術的應用與發(fā)展 .兵團教育學院學報 ,1999 年 04 期 [4] 楊震中 .數字通信技術的進展與趨向 .現代有線傳輸 ,1995 年 01 期 [5] 樊昌信、曹麗娜 .通信原理 .6 版 .北京：國防工業(yè)出版， :188~233 [6] 鄧華 .MATLAB 通信仿真及應用實例詳解 .北京：人民郵電出版社， :294~318 [7] , . MultipleSymbol Differential Detection of MPSK[J].IEEETransactions on :300～ 308. [8] , Noncoherent Detection of PSK Signals[J]. :146~155. [9] 侯光華、韋武鋼、楊森 .衛(wèi)星通信中的 OQPSK調制技術分析 .無線互聯(lián)科技 .2020年 01期 :13~15 [10] 求是科技 .MATLAB .人民郵電出版社 .2020 [11] 陳懷琛、吳大正、高西全 .MATLAB 及在電子信息課程中的應用 .3 版 .北京：電子工業(yè)出版社， [12] 約翰 .、馬蘇 德 .薩勒赫 .現代通信系統(tǒng) ——使用 .西安：西安交通大學出版社， :268~278 [13] Andrea Goldsmith、 Stanford University. WIRELESS COMMUNICATIONS. Cambridge University [14] 周炯磐、龐沁華等，通信原理，北京郵電大學出版社， 2020年 6月 [15] 郭文彬、桑林，通信原理－基于 MATLAB的計算機仿真，北京郵電大學出版社， 2020年 6月 [16] 徐明遠、邵玉 斌， MATLAB仿真在通信與電子工程中的應用，西安電子科技大學出版社， 2020年 6月 [17] Proakis John G. Digital Communication, （影印版，第 3版），北京：電子工業(yè)出版社， 1998； 致謝 本論文是在林霏老師的悉心指導下完成的，從選題、前期準備、研究到論文的構思、寫作和修改，無不凝結了老師的心血。在論文完成之際，首先衷心感謝林老師對我的深切關懷和諄諄教導。林老師淵博的學識、嚴謹的治學態(tài)度、謙和的為人、忘我的科研精神深深地感染了我，使我在學術品質與治學態(tài)度上受 益匪淺，在此謹向林老師致以誠摯的謝意和崇高的敬意。 最后衷心感謝電氣學院的所有給予我指導和幫助的領導、老師和同學，衷心的感謝你們！ 附錄 主程序： %QPSK 信號波形 clc。 clear all。 close all。 %全局變量 Ts=。 %時間分辨率 (時點間間隔），抽樣時間間隔 t=0:Ts:5。 %時間坐標 s N=length(t)。 fs=1/Ts。 df=fs/(N1)。 %df 為頻率分辨率（頻率間隔） n=(N1)/2:(N1)/2。 f=n*df。 %頻率坐標 %基帶信號 figure。 p=[1 0 1 1 0 1 0 0 1 1]。 %待傳送的編碼串 d1=[0::t(end)]。 %每個編碼的發(fā)送延遲時間 p1=[p, zeros(length(d1)length(p))]。 %在發(fā)送串后補零 d=[d1。p1]39。 %產生 d 矩陣 m=pulstran(,d,39。rectpuls39。,)。 subplot (4,1,1)。 plot(t,m) axis([t(1) t(end) ]) grid on title(39?；鶐盘?p39。) %抽取 IQ 信號 QPSK 兩個基帶符號為一組 ， 映射到一對 IQ pQI=reshape(p,2,5)。 pQ=pQI(2,:)。 pI=pQI(1,:)。 d2=[0:1:t(end)1]。 %IQ 信號波形 %dQ=[d2。pQ]39。 %dI=[d2。pI]39。 %mI=pulstran(,dI,39。rectpuls39。,1)。 %mQ=pulstran(,dQ,39。rectpuls39。,1)。 %IQ 值映射 ―0‖映射成 1， ―1‖映射成 1 I(pI==1)=1/sqrt(2)。 I(pI==0)=1/sqrt(2)。 Q(pQ==1)=1/sqrt(2)。 Q(pQ==0)=1/sqrt(2)。 dImap=[d2。I]39。 dQmap=[d2。Q]39。 mImap=pulstran(,dImap,39。rectpuls39。,1)。 mQmap=pulstran(,dQmap,39。rectpuls39。,1)。 subplot(4,1,2) plot(t,mImap) axis([t(1) t(end) ]) grid on title(39。QPSK 映射后的 I 路信號 39。) subplot(4,1,3) plot(t,mQmap) axis([t(1) t(end) ]) grid on title(39。QPSK 映射后的 Q 路信號 39。) %載波信號 fc=10 %載波頻率 c1=cos(2*pi*fc*t)。 c2=sin(2*pi*fc*t)。 %已調信號 s1=mImap.*c1。 s2=mQmap.*c2。 s=s1s2。 subplot(4,1,4) plot(t,s) axis([t(1) t(end) ]) title(39。QPSK 進入信道的信號 s=s1s239。) %QDSK 信號波形 clc。 clear all。 close all。 %全局變量 Ts=。 %時間分辨率 (時點間間隔），抽樣時間間隔 t=0:Ts:5。 %時間坐標 s N=length(t)。 fs=1/Ts。 df=fs/(N1)。 %df 為頻率分辨率（頻率間隔） n=(N1)/2:(N1)/2。 f=n*df。 %頻率坐標 %基帶信號 figure。 p=[1 0 1 1 0 1 0 0 1 1]。 %待傳送的編碼串 d1=[0::t(end)]。 %每個編碼的發(fā)送延遲時間 p1=[p, zeros(length(d1)length(p))]。 %在發(fā)送串后補零 d=[d1。p1]39。 %產生 d 矩陣 m=pulstran(,d,39。rectpuls39。,)。 subplot (4,1,1)。 plot(t,m) axis([t(1) t(end) ]) grid on title(39?；鶐盘?p39。) %抽取 IQ 信號 QDPSK 兩個基帶符號為一組 ， 映射到一對 IQ pQI=reshape(p,2,5)。 pQ=pQI(2,:)。 pI=pQI(1,:)。 d2=[0:1:t(end)1]。 %IQ 信號波形 %dQ=[d2。pQ]39。 %dI=[d2。pI]39。 %mI=pulstran(,dI,39。rectpuls39。,1)。 %mQ=pulstran(,dQ,39。rectpuls39。,1)。 %IQ 值映射 “ 0” 映射成 +1，“ 1” 映射成 1 I(pI==1)=1/sqrt(2)。 I(pI==0)=+1/sqrt(2)。 Q(pQ==1)=1/sqrt(2)。 Q(pQ==0)=+1/sqrt(2)。 dImap=[d2。I]39。 dQmap=[d2。Q]39。 mImap=pulstran(,dImap,39。rectpuls39。,1)。 mQmap=pulstran(,dQmap,39。rectpuls39。,1)。 subplot(4,1,2) plot(t,mImap) axis([t(1) t(end) ]) grid on title(39。QDPSK 映射后的 I 路信號 39。) subplot(4,1,3) plot(t,mQmap) axis([t(1) t(end) ]) grid on title(39。QDPSK 映射后的 Q 路信號 39。) %載波信號 fc=10 %載波頻率 c1=cos(2*pi*fc*t+pi/2)。 c2=sin(2*pi*fc*tpi/2)。 %已調信號 s1=mImap.*c1。 s2=mQmap.*c2。 s=s1+s2。 subplot(4,1,4) plot(t,s) axis([t(1) t(end) ]) title(39。QDPSK 進入信道的信號 s=s1+s239。) %QPSK 功率譜密度 m=1。 p1=10。 Fs=1000。 nfft=1024。 cs_n=xcorr(s,39。unbiased39。)。 %計算序列的自相關函數 CSk=fft(cs_n,nfft)。 %傅里葉變換 Pxx=abs(CSk)。 index=0:round(nfft/21)。 k=index*Fs/nfft。 plot_Pxx=10*log10(Pxx(index+1))。 figure plot(k,plot_Pxx)。 title(39。QPSK 功率譜密度 39。) %QDPSK 功率譜密度 m=。 p1=10。 Fs=1000。 nfft=1024。 cs_n=xcorr(s,39。unbiased39。)。 %計算序列的自相關函數 CSk=fft(cs_n,nfft)。 %傅里葉變換 Pxx=abs(CSk)。 index=0:round(nfft/21)。 k=index*Fs/nfft。 plot_Pxx=10*log10(Pxx(index+1))。 figure plot(k,plot_Pxx)。 title(39。QDPSK 功率譜密度 39。) %QPSK 誤碼率 clear all。 clc。 M = 4。 x = randint(100000,1,M)。 % 隨機序列 s= pskmod(x,M)。 % 調制 SNR = 0:2:30。 % 信噪比范圍 for n = 1:length(SNR) r = awgn(s,SNR(n),39。measured39。)。 % 加高斯噪聲 r1 = pskdemod(r,M)。 % 解調 [nErrors, BER(n)] = biterr(x(1:end),r1(1:end))。%計算仿真誤比特率 end for i=1:length(SNR) SNR1(i)=10^(SNR(i)/10)。 % 信噪比 BERtheory(i)=1(*erfc(sqrt(SNR1(i)/2)))^2。 BERtheory(i)=BERtheory(i)/2。 end figure semilogy(SNR,BERtheory,39。b39。)。 %不同顏色畫出誤比特率曲線對比 hold on semilogy(SNR,BER,39。r39。)。 %不同顏色畫出誤比特率曲線對比 legend(39。 理論 BER39。,39。仿真 BER39。)。 %曲線名稱 grid on xlabel(39。SNR (dB)39。)。 ylabel(39。BER39。)。 %橫縱坐標 title(39。QPSK 通過 AWGN 信道 39。)。 %QDPSK 誤碼率 clear all。 close all。 clc。 %