【正文】 入，相信QAM調(diào)制解調(diào)技術(shù)必將有更廣闊的應(yīng)用前景。參考文獻(xiàn)[1 ] 樊昌信、張甫翊、徐炳祥、吳成柯等編著。《通信原理》(第五版)。北京：國防工業(yè)出版社，2001[2 ] 張輝、曹麗娜編著?！冬F(xiàn)代通信原理與技術(shù)》。西安：西安電子科技大學(xué)出版社，2002[3 ] 曹志剛、錢亞生編著。 《現(xiàn)代通信原理。 》北京：清華大學(xué)出版社，1992[4 ] 曾興雯、劉乃安、陳健等編著?！陡哳l電路原理與分析》。西安：西安電子科技大學(xué)出版社，2001[5 ] 通信原理實驗編寫組編著?！锻ㄐ旁韺嶒灐飞虾＃荷虾＝煌ù髮W(xué)出版社，1987[6 ] 趙民建、袁夢濤、李世巨等。全數(shù)字多星座圖、可變符號率QAM調(diào)制器[J]。電路與系統(tǒng)學(xué)報，2001，6[7 ] 樊昌信、宮錦文、劉忠成編著《。通信原理及系統(tǒng)實驗》電子工業(yè)出版社，2007[8 ] 王興亮、達(dá)新宇編?！稊?shù)字通信原理與技術(shù)》(第二版)。 西安電子科技大學(xué)出版社，2004[9 ] 宮豐奎、張力、連冰等《針對QAM信號的解調(diào)系統(tǒng)的設(shè)計》[J]。電視技術(shù)，2003，12[10] 傅海陽、趙品勇《SDH微波通信系統(tǒng)》北京：人民郵電出版社，2003[11] 王立寧，MATLAB與通信仿真[M]北京：人民郵電出版社，2000[12] 張葛祥、李娜。Matlab仿真技術(shù)與應(yīng)用[M]，北京：清華大學(xué)出版社，2003[13] 鄧華。Matlab通信仿真及應(yīng)用實例[M]，北京：人民郵電大學(xué)出版社，2003[14] 于風(fēng)云，張平。QAM調(diào)制與解調(diào)的全數(shù)字實現(xiàn)[J]，現(xiàn)代電子技術(shù)，2005，（3）=5355[15] 陶為戈，張娟，朱日失華，基于FPGA的QAM調(diào)制器的設(shè)計[J]，江蘇技術(shù)師范學(xué)院學(xué)報，2007,13（2）；3539[16] Bygan David A。QAM for terrestrial and cable transmission。IEEE Trans。On Consumer Elecetronics。1995,41(3)。383391[17]李白萍。現(xiàn)代通信理論[M]。西安電子科技大學(xué)出版社，2006；143144[18] Peebles PZ Communication System Principles,AddisonWesley Publishing Company,1976,3642[19] 張森、張正亮等，MATLAB仿真技術(shù)與實力應(yīng)用教程。機(jī)械工業(yè)出版社，2004[20] 樊昌信等，通信原理；國防工業(yè)出版社，2004畢業(yè)設(shè)計小結(jié)經(jīng)過這一次的畢業(yè)設(shè)計，我學(xué)習(xí)并學(xué)會了使用應(yīng)用軟件Matlab，能夠在該環(huán)境下進(jìn)行簡單的電路仿真。基本上懂得了數(shù)字調(diào)制解調(diào)技術(shù)的發(fā)展?fàn)顩r，并對振幅調(diào)制解調(diào)原理有了更深一步的認(rèn)識理解了MQAM的原理，知道了眼圖、星座圖、各個頻譜的產(chǎn)生原理、誤碼性能測試和相關(guān)參數(shù)的設(shè)定。畢業(yè)設(shè)計帶給我了不少收獲，從學(xué)習(xí)MATLAB的使用，到仿真程序，程序剛開始的時候有不少錯誤，一些錯誤是由于粗心造成的，一些錯誤是由于對于原理理解上的偏差造成的。這也讓我明白了，做什么事情都要嚴(yán)謹(jǐn)，不可以有半點馬虎。同時此次畢業(yè)設(shè)計也暴露出自己專業(yè)基礎(chǔ)知識的很多不足。比如缺乏綜合應(yīng)用和理論聯(lián)系實際的能力，等等。這次實踐是對自己大學(xué)四年所學(xué)只是的一次檢閱，使我明白自己的知識還很淺薄，雖然馬上畢業(yè)了，但是自己的求學(xué)之路還很漫長，以后更應(yīng)該在工作中學(xué)習(xí)，不明白的地方一定要不恥下問，努力使自己成為一個對社會有貢獻(xiàn)的人。感謝老師對我的諄諄教導(dǎo)，讓我學(xué)到了很多的知識。附 錄附源程序代碼： main_plot。mclear。clc。echo off。close all。N=10000。 %設(shè)定碼元數(shù)量fb=1。 %基帶信號頻率fs=32。 %抽樣頻率fc=4。 %載波頻率，為便于觀察已調(diào)信號，我們把載波頻率設(shè)的較低Kbase=2。 % Kbase=1，不經(jīng)基帶成形濾波，直接調(diào)制。 % Kbase=2，基帶經(jīng)成形濾波器濾波后，再進(jìn)行調(diào)制info=random_binary(N)。 %產(chǎn)生二進(jìn)制信號序列[y，I，Q]=qam(info，Kbase，fs，fb，fc)。 %對基帶信號進(jìn)行16QAM調(diào)制y1=y。 y2=y。 %備份信號，供后續(xù)仿真用T=length(info)/fb。 m=fs/fb。 nn=length(info)。dt=1/fs。 t=0:dt:Tdt。 subplot(211)。 %便于觀察，這里顯示的已調(diào)信號及其頻譜均為無噪聲干擾的理想情況%由于測試信號碼元數(shù)量為10000個，在這里我們只顯示其總數(shù)的1/10plot(t(1:1000)，y(1:1000)，t(1:1000)，I(1:1000)，t(1:1000)，Q(1:1000)，[0 35]，[0 0]，39。b:39。)。title(39。已調(diào)信號(In:red，Qn:green)39。)。%傅里葉變換，求出已調(diào)信號的頻譜n=length(y)。 y=fft(y)/n。 y=abs(y(1:fix(n/2)))*2。q=find(y1e04)。 y(q)=1e04。 y=20*log10(y)。f1=m/n。 f=0:f1:(length(y)1)*f1。subplot(223)。 plot(f，y，39。r39。)。 grid on。 title(39。已調(diào)信號頻譜39。)。 xlabel(39。f/fb39。)。 %畫出16QAM調(diào)制方式對應(yīng)的星座圖subplot(224)。 constel(y1，fs，fb，fc)。 title(39。星座圖39。)。SNR_in_dB=8:2:24。 %AWGN信道信噪比for j=1:length(SNR_in_dB) y_add_noise=awgn(y2，SNR_in_dB(j))。 %加入不同強(qiáng)度的高斯白噪聲 y_output=qamdet(y_add_noise，fs，fb，fc)。 %對已調(diào)信號進(jìn)行解調(diào) numoferr=0。 for i=1:N if (y_output(i)~=info(i))， numoferr=numoferr+1。 end。 end。 Pe(j)=numoferr/N。 %統(tǒng)計誤碼率end。figure。semilogy(SNR_in_dB，Pe，39。red*39。)。 grid on。xlabel(39。SNR in dB39。)。 ylabel(39。Pe39。)。title(39。16QAM調(diào)制在不同信道噪聲強(qiáng)度下的誤碼率39。)。random_binary。m%產(chǎn)生二進(jìn)制信源隨機(jī)序列function [info]=random_binary(N) if nargin == 0， %如果沒有輸入?yún)?shù)，則指定信息序列為10000個碼元 N=10000。end。for i=1:N， temp=rand。 if (temp0。5)， info(i)=0。 % 1/2的概率輸出為0 else info(i)=1。 % 1/2的概率輸出為1 endend。qam。mfunction [y，I，Q]=qam(x，Kbase，fs，fb，fc)。%T=length(x)/fb。 m=fs/fb。 nn=length(x)。dt=1/fs。 t=0:dt:Tdt。%串/并變換分離出I分量、Q分量，然后再分別進(jìn)行電平映射I=x(1:2:nn1)。 [I，In]=two2four(I，4*m)。 Q=x(2:2:nn)。 [Q，Qn]=two2four(Q，4*m)。 if Kbase==2。 %基帶成形濾波 I=bshape(I，fs，fb/4)。 Q=bshape(Q，fs，fb/4)。 end。 y=I。*cos(2*pi*fc*t)Q。*sin(2*pi*fc*t)。 %調(diào)制qamdet。m%QAM信號解調(diào)function [xn，x]=qamdet(y，fs，fb，fc)。dt=1/fs。 t=0:dt:(length(y)1)*dt。I=y。*cos(2*pi*fc*t)。 Q=y。*sin(2*pi*fc*t)。 [b，a]=butter(2，2*fb/fs)。 %設(shè)計巴特沃斯濾波器I=filtfilt(b，a，I)。 Q=filtfilt(b，a，Q)。m=4*fs/fb。 N=length(y)/m。 n=(。6:1:N)*m。 n=fix(n)。In=I(n)。 Qn=Q(n)。 xn=four2two([In Qn])。 %I分量Q分量并/串轉(zhuǎn)換，最終恢復(fù)成碼元序列xnnn=length(xn)。 xn=[xn(1:nn/2)。xn(nn/2+1:nn)]。 xn=xn(:)。 xn=xn39。bshape。m%基帶升余弦成形濾波器function y=bshape(x，fs，fb，N，alfa，delay)。%設(shè)置默認(rèn)參數(shù)if nargin6。 delay=8。 end。if nargin5。 alfa=0。5。 end。if nargin4。 N=16。 end。b=firrcos(N，fb，2*alfa*fb，fs)。y=filter(b，1，x)。two2four。m%二進(jìn)制轉(zhuǎn)換成四進(jìn)制function [y，yn]=two2four(x，m)。T=[0 1。3 2]。 n=length(x)。 ii=1。for i=1:2:n1。 xi=x(i:i+1)+1。 yn(ii)=T(xi(1)，xi(2))。 ii=ii+1。end。yn=yn1。5。 y=yn。 for i=1:m1。 y=[y。yn]。end。 y=y(:)39。 %映射電平分別為1。5；0。5；0。5；1。5four2two。m%四進(jìn)制轉(zhuǎn)換成二進(jìn)制function xn=four2two(yn)。y=yn。 ymin=min(y)。 ymax=max(y)。 ymax=max([ymax abs(ymin)])。ymin=abs(ymax)。 yn=(yymin)*3/(ymaxymin)。 %設(shè)置門限電平，判決I0=find(yn 0。5)。 yn(I0)=zeros(size(I0))。I1=find(yn=0。5 amp。 yn1。5)。 yn(I1)=ones(size(I1))。I2=find(yn=1。5 amp。 yn2。5)。 yn(I2)=ones(size(I2))*2。I3=find(yn=2。5)。 yn(I3)=ones(size(I3))*3。%一位四進(jìn)制碼元轉(zhuǎn)換為兩位二進(jìn)制碼元T=[0 0。0 1。1 1。1 0]。 n=length(yn)。 for i=1:n。 xn(i，:)=T(yn(i)+1，:)。end。 xn=xn39。 xn=xn(:)。 xn=xn39。constel。m%畫出星座圖function c=constel(x，fs，fb，fc)。N=length(x)。 m=2*fs/fb。 n=fs/fc。 i1=mn。 i=1。 ph0=(i11)*2*pi/n。 while i = N/m。 xi=x(i1:i1+n1)。 y=2*fft(xi)/n。 c(i)=y(2)。 i=i+1。 i1=i1+m。end。 %如果無輸出，則作圖if nargout1。 cmax=max(abs(c))。 ph=(0:5:360)*pi/180。 plot(1。414*cos(ph)，1。414*sin(ph)，39。c39。)。 hold on。 for i=1:length(c)。 ph=ph0angle(c(i))。 a=abs(c(i))/cmax*1。414。 plot(a*cos(ph)，a*sin(ph)，39。r*39。)。 end。 plot([1。5 1。5]，[0 0]，39。k:39。，[0 0]，[1。5 1。5]，39。k:39。)。 hold off。 axis equal。 axis([1。5 1。5 1。5 1。5])。end。