【正文】 但期間 得到各位老師的幫助和指導，尤其是王惠琴老師，在此表示真誠的感謝！并對同組成員對我的幫助一并表示感謝！ 。離開了實踐的理論是沒有任何意義可言的，理論是需要伴隨著實踐才能完善。 36 致謝 在本次課程設計過程中， 我深深體會到了從書本學習與實際應用中的 不同，這種感同身受必將對我們今后的學習與生活帶來很大的幫助。誤碼率 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)。 end rddd=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)。 34 end end ddd22222=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)。 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 32 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)。 [f5,y_16qamf]=T2F(t,y_16qam)。)。 title(39。幅度 39。)。 xlabel(39。 plot(t,y_16QAM)。 figure(2)。 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])。 28 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)。)。 title(39。幅度 39。)。 27 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。)。 24 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)。 %功率譜圖縱