【正文】 。在大學階段，理論的學習和實踐是密不可分的。同時，和班里同學一起討論，共同解決問題，查到了有用的資料，就大家一起分享，我們合理的分工，才能高效的完成這次課設(shè)，而且促進了同學之間的相互了解，加深了我們的友誼，提高了我們的團隊協(xié)作精神。致謝課設(shè)到這里已接近尾聲，在本文完成之際，首先，我要特別感謝我們同組的同學，是他們陪著我一起發(fā)現(xiàn)問題，解決問題，才幫助一起度過了一個有一個難關(guān)，最終獲得現(xiàn)在的成果。誤碼率39。)。) xlabel(39。) title(39。 figure(7) semilogy(log(p1),ber,39。)p1=0::10。) ylabel(39。 xlabel(39。16QAM星座圖39。 end end end axis([4 4 4 4])。grid。解調(diào)信號功率譜圖39。)。ylabel(39。頻率(Hz)39。axis([fx1,fx2,fy1,fy2])。subplot(2,2,4)。 [f8,r_16qamf]=T2F(t,r_16qam(1:Lt))。)。title(39。幅度39。)。xlabel(39。plot(t,r_16qam(1:Lt))。figure(5)。endrddd=sigexpand(ddd,fc*N_sample)。for s=1:N/2 ddd(2*s1)=dd111111(s)。grid。下支路42電平變換后信號時域波形圖39。)。ylabel(39。時間(S)39。plot(t,d_NRZ24(1:Lt))axis([tx1,tx2,ty1,ty2])。 figure(4)。 endendddd22222=sigexpand(dd222222,2*fc*N_sample)。 elseif dd2222(q)==3 dd222222(2*q1)=1。 elseif dd2222(q)==1 dd222222(2*q1)=1。 elseif dd2222(q)==1 dd222222(2*q1)=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)。 endendddd11111=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)。 endenddd22222=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)。%通過乘法器2r_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)。[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)。 endend% 生成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)。)。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。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])。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?；鶐盘枙r域波形圖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