【正文】 elseif dd2222(q)==1 。 for q=1:N/4 if dd2222(q)==3 dd222222(2*q1)=1。 grid。上支路 42變換后信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 plot(t,d_NRZ14(1:Lt)) axis([tx1,tx2,ty1,ty2])。 figure(4)。 end end ddd11111=sigexpand(dd111111,2*fc*N_sample)。 elseif dd1111(p)==3 dd111111(2*p1)=1。 elseif dd1111(p)==1 dd111111(2*p1)=1。 elseif dd1111(p)==1 dd111111(2*p1)=1。 for p=1:N/4 if dd1111(p)==3 dd111111(2*p1)=1。 grid。下支路抽樣判決后信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 plot(t,d_NRZ23(1:Lt)) axis([tx1,tx2,ty1,ty2])。 figure(4)。 end end dd22222=sigexpand(dd2222,4*fc*N_sample)。 dd222(n)0 dd2222(n)=1。 dd222(n)1 dd2222(n)=1。 for n=1:N/4 if dd222(n)1 dd2222(n)=3。 %抽樣判決 dd222=r_16QAM21(2*fc*N_sample:4*fc*N_sample:end)。)。 title(39。幅度 39。)。 xlabel(39。 subplot(2,2,4)。 r_16QAM21=filter(b,a,r_16QAM21)。 %通過乘法器 2 r_16QAM21 = y_16qam .* h2t。)。 title(39。幅度 39。)。 xlabel(39。 subplot(2,2,1)。 d_NRZ13=conv(dd11111,gt2)。 else dd1111(n)=3。 elseif 1dd111(n) amp。 elseif 0dd111(n) amp。 %在每個碼元中間抽樣 dd1111=[]。 grid。上支路通過低通濾波器信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 plot(t,r_16QAM11) axis([tx1,tx2,ty1,ty2])。 figure(3)。 %通過低通濾波器 [b,a]=butter(3,)。 grid。接收信號功率譜圖 39。)。 ylabel(39。頻率 (Hz)39。 axis([fx1,fx2,fy1,fy2])。 subplot(2,4,8)。 33 [f5,y_16qamf]=T2F(t,y_16qam)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,y_16QAM)。 figure(2)。 %使信號通過瑞利信道 fadedSig = filter(chan,s_16qam)。 noise =wgn(m1,Lt,p1)。 % 信道加入高斯白噪聲進行接受解調(diào) % 產(chǎn)生高斯白噪聲 m1=1。)。 title(39。功率譜密度 (dB/Hz)39。)。 xlabel(39。 plot(f4,10*log10(abs(s_16qamf).^2/T))。 figure(2)。 grid。已調(diào)信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,4,3)。 s_16qam=s_16qam1+s_16qam2。)。 title(39。功率譜密度 (dB/Hz)39。)。 xlabel(39。 plot(f3,10*log10(abs(s_16qam2f).^2/T))。 figure(2)。 grid。下支路已調(diào)信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,4,2)。 grid。上支路已調(diào)信號功率譜圖 39。)。 ylabel(39。頻率 (Hz)39。 axis([fx1,fx2,fy1,fy2])。 subplot(2,4,5)。 [f2,s_16qam1f]=T2F(t,s_16qam1)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,s_16qam1)。 %生成下支路頻帶信號 figure(2)。 end end % 生成 16QAM信號 s_16qam1=pd_NRZ11 .* h1t。 for k=1:8*tao*fc*N_sample pd_NRZ11((w1)*64+k)=phxh1(k)。 for w=2:N/4 phxh1=d11(w1)+(d11(w)d11(w1))**(1+sin(pi*tm/8*tao*Ts))。 tao=3/16。 % 對上下支路信號進行平滑 pd_NRZ11=d_NRZ11(1:Lt)。)。 30 title(39。幅度 39。)。 xlabel(39。 plot(t,d_NRZ22(1:Lt))。 figure(3)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,d_NRZ22(1:Lt))。 figure(1)。 dd22=sigexpand(d22,4*fc*N_sample)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,d_NRZ11(1:Lt))。 figure(3)。)。 title(39。幅度 39。)。 xlabel(39。 29 plot(t,d_NRZ11(1:Lt))。 figure(1)。 % 上支路 gt2=ones(1,4*fc*N_sample)。 d22(m)=2*d2(2*m1)+d2(2*m)。 for m=1:N/4。 % 上支路 d22=[]。 grid。載波信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,3,4)。 h2t=sin(2*pi*fc*t)。 grid。串并轉(zhuǎn)換后下支路信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 28 subplot(2,3,5)。 d_NRZ2=conv(dd2,gt1)。 grid。串并轉(zhuǎn)換后上支路信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,3,2)。 d_NRZ1=conv(dd1,gt1)。 end % 上支路 dd1=sigexpand(d1,2*fc*N_sample)。 for i=1:N/2 d1(i)=d(2*(i1)+1)。 % 串并轉(zhuǎn)換 d1=[]。)。 title(39。功率譜密度 (dB/Hz)39。)。 xlabel(39。 plot(f1,10*log10(abs(d_NRZf).^2/T))。 figure(5)。 grid。基帶信號時域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,2,1)。 grid?；鶐盘枙r域波形圖 39。)。 ylabel(39。時間 (S)39。 axis([tx1,tx2,ty1,ty2])。 subplot(2,3,1)。 % NRZ波形 d_NRZ=conv(dd,gt)。 %功率譜圖縱坐標終點 %產(chǎn)生二進制信源 d=[1,1,1,1,1,1,1,1] dd=sigexpand(d,fc*N_sample)。 %功率譜圖橫坐標終點 fy1=40。 %時域波形圖縱坐標終點 fx1=10。 %時域波形圖橫坐標終點 ty1=。 % 時間向量長度 tx1=0。 % 信號持續(xù)時間長度 t=0:dt:Tdt。 % 碼元寬度 dt=Ts/fc/N_sample。 % 基帶碼元抽樣點數(shù) N=8。 主函數(shù)代碼 fc=2。 sf=fft(st)。 N=length(st)。 T=t(end)。 計算信號的傅里葉變化 function[f,sf]=T2F(t,st)。 out(1,:)=d。： 6578 [6]吳偉陵等著，移動通信原理，北京：電子工業(yè)出版社， 2020： 3057 [7].李建新，現(xiàn)代通信系統(tǒng)分析與仿真 MATLAB 通信工具箱，西安：西安電子科技大學出版社， 2020： 1345 [8]潘子宇， Matlab 通信仿真設計指導書，南京工程學院， 2020： 2436 [9]劉敏， MATLAB 通信仿真與應用，北京：國防工業(yè)出版社： 5660 25 附 錄 源程序代碼 如下 將輸入的序列擴展成間隔為 N1個 0的序列 function [out]=sigexpand(d,M) N=length(d)。 通過本次課程設計使我對 MATLAB 語言有了更進一步的掌握，對無線衰落信道的性能有了更清楚的了解。通過觀察 16QAM 調(diào)制解調(diào)信號的頻譜以及對查閱的資料進行研究發(fā)現(xiàn) 16QAM 的抗噪聲性能相當優(yōu)越。 23 總結 本課程設計主要是在 MATLAB 環(huán)境下，設計仿真了基于 QAM 調(diào)制的無線衰落信道模型，首先對信號源產(chǎn)生的二進制序列進行 QAM調(diào)制，然后使已調(diào)信號分別通過瑞利，高斯信道，在進行 QAM 解調(diào)，對解調(diào)后的信號進行誤碼率分析。因此，利用多電平誤碼率的分析方法，可得到 M 進制 QAM的誤碼率為 [1]： ])(1l o g3[)11(022 nEL Le r f cLP be ??? () 式中， ML? ， Eb 為每碼元能量， n0為噪聲單邊功率譜密度。因此無線通信中的瑞利信道對信號傳輸造成的影響更大。 42 電平轉(zhuǎn)換是 24 電平轉(zhuǎn)換的逆過程，其映射關系如下圖所示 ： 表 42 電平映射關系 映射前數(shù)據(jù) 電平 /V 雙極性 3 00 1 1 1 01 1 1 1 10 1 1 3 11 1 1 圖 47上下支路抽樣判決及 42轉(zhuǎn)換后信號時域波形圖 并串轉(zhuǎn)換 經(jīng)過并串轉(zhuǎn)換即可得到解調(diào)信號 ， 將兩路并聯(lián)信號經(jīng)過轉(zhuǎn)換成 為一路信號，下圖為基帶信號和解調(diào)信號的對比及它們的頻率譜密度圖 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 電平變換后信號時域波形圖 21 圖 48高斯噪聲下的 基帶 信號 與解調(diào)信號 圖 49瑞利信道 下的 基帶 信號 與解調(diào)信號 可以看出， 加高斯白噪聲時 并沒有錯碼， 經(jīng)統(tǒng)計得此時 誤碼率為 0。整個函數(shù)表示信號通過濾波器的