【正文】 [f3,s_16qam2f]=T2F(t,s_16qam2)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,s_16qam2)。 figure(2)。)。 title(39。功率譜密度 (dB/Hz)39。)。 xlabel(39。 plot(f2,10*log10(abs(s_16qam1f).^2/T))。 figure(2)。 grid。上支路已調(diào)信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,4,1)。 %生成上支路頻帶信號 s_16qam2=pd_NRZ22 .* h2t。 pd_NRZ22((w1)*64+k)=phxh2(k)。 phxh2=d22(w1)+(d22(w)d22(w1))**(1+sin(pi*tm/8*tao*Ts))。 tm=0:dt:8*Ts*taodt。 pd_NRZ22=d_NRZ22(1:Lt)。 grid。24轉(zhuǎn)換后下支路信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,2,3)。 grid。24轉(zhuǎn)換后下支路信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,3,6)。 % 下支路 d_NRZ22=conv(dd22,gt2)。 grid。24轉(zhuǎn)換后上支路信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,2,1)。 grid。24電平轉(zhuǎn)換后上支路信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,3,3)。 d_NRZ11=conv(dd11,gt2)。 end dd11=sigexpand(d11,4*fc*N_sample)。 d11(m)=2*d1(2*m1)+d1(2*m)。 % 下支路 % 上下支路 24電平轉(zhuǎn)換 d11=[]。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,h1t)。 figure(1)。 % 載波 h1t=cos(2*pi*fc*t)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,d_NRZ2(1:Lt))。 figure(1)。 % 下支路 dd2=sigexpand(d2,2*fc*N_sample)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,d_NRZ1(1:Lt))。 figure(1)。 gt1=ones(1,2*fc*N_sample)。 d2(i)=d(2*(i1)+2)。 d2=[]。 grid?；鶐盘柟β首V圖 39。)。 ylabel(39。頻率 (Hz)39。 axis([fx1,fx2,fy1,fy2])。 subplot(2,2,2)。 [f1,d_NRZf]=T2F(t,d_NRZ(1:Lt))。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,d_NRZ(1:Lt))。 figure(5)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,d_NRZ(1:Lt))。 % 基帶信號 figure(1)。 % 雙極性 gt=ones(1,fc*N_sample)。 %功率譜圖縱坐標(biāo)起點(diǎn) fy2=25。 %功率譜圖橫坐標(biāo)起點(diǎn) fx2=10。 %時域波形圖縱坐標(biāo)起點(diǎn) ty2=。 % 時域波形圖橫坐標(biāo)起點(diǎn) tx2=8。 % 時間向量 Lt=length(t)。 % 抽樣時間間隔 T=N*Ts。 % 碼元數(shù) Ts=1。 % 載波頻率 N_sample=8。 sf=T/N*fftshift(sf)。 f=N/2*df:df:N/2*dfdf。 df=1/T。 dt=t(2)t(1)。 out=reshape(out,1,M*N)。 out=zeros(M,N)。 另外，我還要感謝在論文撰寫過程中給我?guī)椭耐瑢W(xué)們，特別是和我同宿舍的她們，她們給予了我很大的幫助和支持，是她們無私的付出才成就了我今天的成功，謝謝！ 最后，我要感謝我的家人，是他們始終如一地支持著我，他們是我學(xué)習(xí)和生活中的摯友，他們是我人生中最寶貴的財富，這份彌足珍貴的情誼將會永遠(yuǎn)深藏在我的心中。 致 謝 論文寫到這里已接近尾聲，在本文完成之際，首先，我要特別感謝我的指導(dǎo)老師王玲老師，在她的諄諄教導(dǎo)和督促之下，我才順利的完成了我的畢業(yè)設(shè)計(jì)工作，在這段時間里，我又學(xué)到了新的知識。 5. 將 16QAM調(diào)制解調(diào)系統(tǒng)與 2DPSK系統(tǒng)的抗噪聲性能進(jìn)行對比，獲得了他們兩者的誤碼率曲線，進(jìn)行比較，發(fā)現(xiàn) 16QAM的抗噪聲性能不如 2DPSK，這是與理論相符的，也即，當(dāng)信噪比一定的情況下， M越大，誤碼率 Pe也越大。 3. 利用 MATLAB實(shí)現(xiàn)了 16QAM調(diào)制與解調(diào)系統(tǒng)的設(shè)計(jì)，實(shí)現(xiàn)與仿真， 并得到相應(yīng)的調(diào)制解調(diào)波形，發(fā)現(xiàn)解調(diào)信號波形與輸入信號波形存在一定時延，所以該系統(tǒng)的實(shí)時性有不足，但并不影響對誤碼率的檢測，以及系統(tǒng)能夠的抗噪聲性能。 AM d2 AM d1 第五章 結(jié)論與總結(jié) 本文研究的重點(diǎn)是對基于 MATLAB的 16QAM 調(diào)制解調(diào)系統(tǒng)進(jìn)行設(shè)計(jì)與仿真，得到以下的結(jié)論 1. 對 16QAM調(diào)制解調(diào)系統(tǒng)基本原理進(jìn)行了較為深入地理解與分析，并且根據(jù)其原理編制了仿真程序。這樣在高斯白噪聲信道處調(diào)整信噪比，得到如圖 43 所示的誤碼率圖。 這里我們采用的是差分相干解調(diào)的方法，所以并不需要做載波恢復(fù)。 （ 2） 16QAM 與 2DPSK 和 64QAM 這里 我們設(shè)計(jì)一個 2DPSK 調(diào)制和解調(diào)系統(tǒng)和一個 64QAM 調(diào)制和解調(diào)系統(tǒng)的仿真模型，以觀察其與 16QAM 信號的抗噪聲性能，并對它們的誤碼率進(jìn)行比較。而 16QAM 信號，在等概率出現(xiàn)條件下，可以計(jì)算出其最大功率和平均功率之比等于 倍，即 dB。但是，這時是在最大功率（振幅）相等的條件下比較的，沒有考慮這兩種體制的平均功率差別。設(shè)其最大振幅為 AM，則 16PSK 信號的相鄰矢量端點(diǎn)的歐氏距離等于 AAd in21 ?? ? 式 41 而 16QAM 信號的相鄰點(diǎn)歐氏距離等于 AMd 2122 ????? 式 42 d2 和 d1 的比值就代表這兩種體制的噪聲容限之比。 通過調(diào)整 高斯白噪聲信道的信噪比 snr（ Eb/No），可以得到如圖 52 所示的誤碼率圖： 1 0 . 5 0 0 . 5 1 1 . 5 2 2 . 5103102101100Q A M 信號誤碼率分析信噪比誤碼率 圖 49 QAM信號誤碼率分析 可見 16QAM 信號的誤碼率隨著信噪比的增大而逐漸減小，這與理論分析是完全一致的 16QAM 信號與其它調(diào)制信號的性能比較 （ 1） 16QAM 和 16PSK 星座圖中相鄰點(diǎn)歐氏距離直接代表這噪聲容限的大小。 end 將兩路并聯(lián)信號經(jīng)過轉(zhuǎn)換成為一路信號，下圖為基帶 信號和解調(diào)信號的對比及它們的頻率譜密度圖： 0 2 4 6 842024時間 ( S )幅度基帶信號時域波形圖 1 0 5 0 5 10 4 0 2 0020頻率 ( H z )功率譜密度(dB/Hz)基帶信號功率譜圖0 2 4 6 842024時間 ( S )幅度解調(diào)信號時域波形圖 P n = 1 0 d B 1 0 5 0 5 10 4 0 2 0020頻率 ( H z )功率譜密度(dB/Hz)解調(diào)信號功率譜圖 P n = 1 0 d B 圖 48基帶與解調(diào)信號 可以看出，并沒有錯碼，可見本次仿真是成功的 第 節(jié) QAM 性能分析 16QAM 抗噪聲性能仿真 對于 QAM，可以看成是由兩個相互正交且獨(dú)立的多電平 ASK 信號疊加而成。 映射前數(shù)據(jù) 電平 /V 雙極性 3 00 1 1 1 01 1 1 1 10 1 1 3 11 1 1 0 2 4 6 842024時間 ( S )幅度上支路抽樣判決后信號時域波形圖0 2 4 6 842024時間 ( S )幅度下支路抽樣判決后信號時域波形圖0 2 4 6 842024時間 ( S )幅度上支路 4 2 變換后信號時域波形圖0 2 4 6 842024時間 ( S )幅度下支路 4 2 電平變換后信號時域波形圖 圖 47 上下支路抽樣判決及 42 轉(zhuǎn)換后信號時域波形圖 并串轉(zhuǎn)換 經(jīng)過并串轉(zhuǎn)換即可得到解調(diào)信號，采用循環(huán)語句實(shí)現(xiàn)， for s=1:N/2 ddd(2*s1)=dd111111(s)。 0 2 4 6 8432101234時間 ( S )幅度上支路通過低通濾波器信號時域波形圖0 2 4 6 8432101234時間 ( S )幅度下支路通過低通濾波器信號時域波形圖 圖 46上下支路通過低通濾波器信號時域波形圖 抽樣判決和 42 電平轉(zhuǎn)換 抽樣判決是在每個碼元中間抽樣，并用幾個判決語句進(jìn)行判決，例如當(dāng)碼元幅度大于 1 時，判為 3。在第二行程序中， x 表示輸入序列， y 表示輸出序列。 y=filter(b,a,x)。 濾波器 IIR 濾波器采用的巴特沃斯低通濾波器有現(xiàn)成的模型，我們可以加以利用， 因此 在本文涉及的仿真中濾波器均選擇巴特沃斯低通濾波器。 為使解調(diào)效果較好，采用噪聲的強(qiáng)度較小，設(shè)置 Pn=10dB. 第 節(jié) 16QAM 解調(diào)模塊的建立與仿真 系統(tǒng)先前所得的 16QAM 調(diào)制信號通過高斯白噪聲信道以后便可以解調(diào)了。 對已調(diào)制信號可采用 wgn 函數(shù)添加加性高斯噪聲。 上下支路信號在加載波之前還經(jīng)過平滑處理，以濾除較高頻率的信號，使實(shí)驗(yàn)結(jié)果更加理想。 經(jīng)過串并轉(zhuǎn)換后，并行輸出的每一路碼元傳輸速率降為原來的一半即 Rb/2. 輸入 d： 1 1 1 1 1 1 1 1 上支路 d_NRZ1： 1 1 1 1 下支路 d_NRZ2： 1 1 1 1 0 2 4 6 8432101234時間 ( S )幅度串并轉(zhuǎn)換后上支路信號時域波形圖0 2 4 6 8432101234時間 ( S )幅度串并轉(zhuǎn)換后下支路信號時域波形圖 圖 42 串并轉(zhuǎn)換后上下支路信號時域波形圖 24 電平轉(zhuǎn)換 24 電平轉(zhuǎn)換就是將輸入信號的 2 電平信號狀態(tài)經(jīng)過轉(zhuǎn)換后變成相應(yīng)的 4 電平信號。 下圖為 16QAM 解調(diào)框圖： 圖 36 16 QAM 解調(diào)框圖 第四章 基于