【正文】 sample:end)。 %在每個(gè)碼元中間抽樣 dd1111=[]。for n=1:N/4 if dd111(n)1 dd1111(n)=3。 elseif 0dd111(n) amp。 dd111(n)1 dd1111(n)=1。 elseif 1dd111(n) amp。 dd111(n)0 dd1111(n)=1。 else dd1111(n)=3。 endend dd11111=sigexpand(dd1111,4*fc*N_sample)。 d_NRZ13=conv(dd11111,gt2)。figure(4)。subplot(2,2,1)。plot(t,d_NRZ13(1:Lt))axis([tx1,tx2,ty1,ty2])。xlabel(39。時(shí)間(S)39。)。ylabel(39。幅度39。)。title(39。上支路抽樣判決后信號(hào)時(shí)域波形圖39。)。grid。%通過(guò)乘法器2r_16QAM21 = y_16qam .* h2t。%通過(guò)低通濾波器[b,a]=butter(3,)。r_16QAM21=filter(b,a,r_16QAM21)。figure(3)。subplot(2,2,4)。plot(t,r_16QAM21)axis([tx1,tx2,ty1,ty2])。xlabel(39。時(shí)間(S)39。)。ylabel(39。幅度39。)。title(39。下支路通過(guò)低通濾波器信號(hào)時(shí)域波形圖39。)。grid。%抽樣判決 dd222=r_16QAM21(2*fc*N_sample:4*fc*N_sample:end)。 %在每個(gè)碼元中間抽樣dd2222=[]。for n=1:N/4 if dd222(n)1 dd2222(n)=3。 elseif 0dd222(n) amp。 dd222(n)1 dd2222(n)=1。 elseif 1dd222(n) amp。 dd222(n)0 dd2222(n)=1。 else dd2222(n)=3。 endenddd22222=sigexpand(dd2222,4*fc*N_sample)。 d_NRZ23=conv(dd22222,gt2)。figure(4)。subplot(2,2,3)。plot(t,d_NRZ23(1:Lt))axis([tx1,tx2,ty1,ty2])。xlabel(39。時(shí)間(S)39。)。ylabel(39。幅度39。)。title(39。下支路抽樣判決后信號(hào)時(shí)域波形圖39。)。grid。% 42電平轉(zhuǎn)換dd111111=[]。for p=1:N/4 if dd1111(p)==3 dd111111(2*p1)=1。 dd111111(2*p)=1。 elseif dd1111(p)==1 dd111111(2*p1)=1。 dd111111(2*p)=1。 elseif dd1111(p)==1 dd111111(2*p1)=1。 dd111111(2*p)=1。 elseif dd1111(p)==3 dd111111(2*p1)=1。 dd111111(2*p)=1。 endendddd11111=sigexpand(dd111111,2*fc*N_sample)。 d_NRZ14=conv(ddd11111,gt1)。 figure(4)。subplot(2,2,2)。plot(t,d_NRZ14(1:Lt))axis([tx1,tx2,ty1,ty2])。xlabel(39。時(shí)間(S)39。)。ylabel(39。幅度39。)。title(39。上支路42變換后信號(hào)時(shí)域波形圖39。)。grid。dd222222=[]。for q=1:N/4 if dd2222(q)==3 dd222222(2*q1)=1。 dd222222(2*q)=1。 elseif dd2222(q)==1 dd222222(2*q1)=1。 dd222222(2*q)=1。 elseif dd2222(q)==1 dd222222(2*q1)=1。 dd222222(2*q)=1。 elseif dd2222(q)==3 dd222222(2*q1)=1。 dd222222(2*q)=1。 endendddd22222=sigexpand(dd222222,2*fc*N_sample)。 d_NRZ24=conv(ddd22222,gt1)。 figure(4)。subplot(2,2,4)。plot(t,d_NRZ24(1:Lt))axis([tx1,tx2,ty1,ty2])。xlabel(39。時(shí)間(S)39。)。ylabel(39。幅度39。)。title(39。下支路42電平變換后信號(hào)時(shí)域波形圖39。)。grid。% 并串轉(zhuǎn)換ddd=[]。for s=1:N/2 ddd(2*s1)=dd111111(s)。 ddd(2*s)=dd222222(s)。endrddd=sigexpand(ddd,fc*N_sample)。 %解調(diào)信號(hào)輸出r_16QAM=conv(rddd,gt)。figure(5)。subplot(2,2,3)。plot(t,r_16qam(1:Lt))。axis([tx1,tx2,ty1,ty2])。xlabel(39。時(shí)間(S)39。)。ylabel(39。幅度39。)。title(39。解調(diào)信號(hào)時(shí)域波形圖39。)。grid。 [f8,r_16qamf]=T2F(t,r_16qam(1:Lt))。figure(5)。subplot(2,2,4)。plot(f8,10*log10(abs(r_16QAMf).^2/T))。axis([fx1,fx2,fy1,fy2])。xlabel(39。頻率(Hz)39。)。ylabel(39。功率譜密度(dB/Hz)39。)。title(39。解調(diào)信號(hào)功率譜圖39。)。grid。figure(6)hold on for n1=3:3 for n2=3:3 if ~((n1==0)|(n2==0)|(n1==2)|(n1==2)|(n2==2)|(n2==2)) scatter(n1,n2)。 end end end axis([4 4 4 4])。 hold off title(39。16QAM星座圖39。)。 xlabel(39。同相支路39。) ylabel(39。正交支路39。)p1=0::10。 ber=3/4*erfc(sqrt(*p1))。 figure(7) semilogy(log(p1),ber,39。b*39。) title(39。QAM信號(hào)誤碼率分析39。) xlabel(39。信噪比39。)。 ylabel(39。誤碼率39。)。致謝 這次課程設(shè)計(jì)，給我留下了很深的印象。在這期間，讓我受益匪淺 。 通過(guò)這次課程設(shè)計(jì)，使我對(duì)正交幅度調(diào)制有了全面的認(rèn)識(shí)，對(duì)通信系統(tǒng)的知識(shí)又有了深刻的理解，在之前通信原理的學(xué)習(xí)以及完成課后的作業(yè)的過(guò)程中，已經(jīng)使用過(guò)MATLAB，對(duì)其有了一些基礎(chǔ)的了解和認(rèn)識(shí)。本次經(jīng)過(guò)一周的課程設(shè)計(jì)，讓我們更熟悉了通信原理的基本知識(shí)和MATLAB的語(yǔ)言，讓我們把課上的理論知識(shí)運(yùn)用到實(shí)際中去，讓我們更近一步地鞏固了課堂上所學(xué)的理論知識(shí)，并能很好地理解與掌握通信系統(tǒng)仿真中的基本概念、基本原理、基本分析方法。 在此次課程設(shè)計(jì)過(guò)程中，要特別感謝我的指導(dǎo)老師的指導(dǎo)與督促，同時(shí)感謝她的諒解與包容。沒(méi)有老師的幫助也就沒(méi)有今天的這篇文章。求學(xué)歷程是艱苦的，但又是快樂(lè)的。 總的來(lái)說(shuō)，通過(guò)這次設(shè)計(jì)使我了解了MATLAB的使用方法，學(xué)會(huì)分析通信系統(tǒng)性能和誤碼率，提高了分析和動(dòng)手實(shí)踐能力。同時(shí)我相信，進(jìn)一步加強(qiáng)對(duì)MATLAB的學(xué)習(xí)與研究對(duì)我今后的學(xué)習(xí)將會(huì)起到很大的幫助！本文參考了大量的文獻(xiàn)資料，在此，向各學(xué)術(shù)界的前輩們致敬！32