【正文】 end。 axis equal。)。,[0 0],[ ],39。 plot([ ],[0 0],39。)。 plot(a*cos(ph),a*sin(ph),39。 ph=ph0angle(c(i))。 hold on。c39。 ph=(0:5:360)*pi/180。 %如果無輸出 ， 則作圖 if nargout1。 i1=i1+m。 c(i)=y(2)。 xi=x(i1:i1+n1)。 ph0=(i11)*2*pi/n。 i1=mn。 m=2*fs/fb。 通信系統(tǒng)綜合設(shè)計(jì)與實(shí)踐 第 31 頁 %畫出星座圖 function c=constel(x,fs,fb,fc)。 xn=xn(:)。 end。 for i=1:n。1 0]。0 1。 yn(I3)=ones(size(I3))*3。 yn(I2)=ones(size(I2))*2。 I2=find(yn= amp。 yn)。 yn(I0)=zeros(size(I0))。 yn=(yymin)*3/(ymaxymin)。 ymax=max([ymax abs(ymin)])。 ymin=min(y)。 %映射電平分別為 ； ； ； %四進(jìn)制轉(zhuǎn)換成二進(jìn)制 function xn=four2two(yn)。 end。 y=[y。 y=yn。 end。 yn(ii)=T(xi(1),xi(2))。 for i=1:2:n1。 n=length(x)。 T=[0 1。 y=filter(b,1,x)。 end。 if nargin4。 alfa=。 end。 %設(shè)置默認(rèn)參數(shù) if nargin6。 xn=xn39。xn(nn/2+1:nn)]。 通信系統(tǒng)綜合設(shè)計(jì)與實(shí)踐 第 29 頁 %I 分量 Q 分量并 /串轉(zhuǎn)換，最終恢復(fù)成碼元序列 xn nn=length(xn)。 Qn=Q(n)。 n=fix(n)。 N=length(y)/m。 Q=filtfilt(b,a,Q)。 [b,a]=butter(2,2*fb/fs)。 I=y.*cos(2*pi*fc*t)。 dt=1/fs。 y=I.*cos(2*pi*fc*t)Q.*sin(2*pi*fc*t)。 Q=bshape(Q,fs,fb/4)。 if Kbase==2。 Q=x(2:2:nn)。 %串 /并變換分離出 I 分量、 Q 分量，然后再分別進(jìn)行電平映射 I=x(1:2:nn1)。 dt=1/fs。 m=fs/fb。 通信系統(tǒng)綜合設(shè)計(jì)與實(shí)踐 第 28 頁 function [y,I,Q]=qam(x,Kbase,fs,fb,fc)。 % 1/2 的概率輸出為 0 else info(i)=1。 for i=1:N, temp=rand。 %產(chǎn)生二進(jìn)制信源隨機(jī)序列 function [info]=random_binary(N) if nargin == 0, %如果沒有輸入?yún)?shù)，則指定信息序列為 10000 個碼元 N=10000。16QAM 調(diào)制在不同信道噪聲強(qiáng)度下的誤碼率 39。)。 ylabel(39。SNR in dB39。 grid on。red*39。 通信系統(tǒng)綜合設(shè)計(jì)與實(shí)踐 第 27 頁 figure。 Pe(j)=numoferr/N。 end。 %對已調(diào)信號進(jìn)行解調(diào) numoferr=0。 %AWGN 信道信噪比 for j=1:length(SNR_in_dB) y_add_noise=awgn(y2,SNR_in_dB(j))。)。 title(39。 %畫出 16QAM 調(diào)制方式對應(yīng)的星座圖 subplot(224)。f/fb39。)。 title(39。)。 plot(f,y,39。 f=0:f1:(length(y)1)*f1。 y=20*log10(y)。 q=find(y1e04)。 y=fft(y)/n。)。 title(39。b:39。 subplot(211)。 dt=1/fs。 m=fs/fb。 y2=y。 %產(chǎn)生二進(jìn)制信號序列 [y,I,Q]=qam(info,Kbase,fs,fb,fc)。 % Kbase=1,不經(jīng)基帶成形濾波，直接調(diào)制 。 %抽樣頻率 fc=4。 %設(shè)定碼元數(shù)量 fb=1。close all。clc。《詳解 MATLAB/Simulink 通信系統(tǒng)建模與仿真（配視頻教程）》?！锻ㄐ畔到y(tǒng)工程》（第二版）?！锻ㄐ旁怼罚ㄉ希??！锻ㄐ旁怼??！锻ㄐ旁怼罚ǖ?6 版）。隨著研究不斷深入，相信 QAM 調(diào)制解調(diào)技術(shù)必將有更廣闊的應(yīng)用前景。QAM 數(shù)字調(diào)制器作為 DVB 系統(tǒng)的前端設(shè)備，接受來自編碼器、服務(wù)器、 DVB 網(wǎng)關(guān)、視頻服務(wù)器等設(shè)備的 TS 流，進(jìn)行 RS 編碼，卷積編碼和 QAM 數(shù)字調(diào)制，輸出的射頻信號可以直接在有線電視網(wǎng)上傳送，同時也可根據(jù)需要選擇中頻輸出，它以其靈活的配置和優(yōu)越的性能指標(biāo)，廣泛的應(yīng)用于數(shù)字有線電視傳輸域和數(shù)字 MMDS 系統(tǒng)。正交振幅調(diào)制 QAM( Quadrature Amplitude Modulation )就是一種頻譜利用率很高的調(diào)制方式，其在中、大容量數(shù)字微波通信系統(tǒng)、有線 電視網(wǎng)絡(luò)高速數(shù)據(jù)傳輸、衛(wèi)星通信系統(tǒng)等領(lǐng)域得到了廣 泛應(yīng)用。 未來展望 在現(xiàn)代通信中，提高頻譜利用率一直是人們關(guān)注的焦點(diǎn)之一。實(shí)際的信道是很復(fù)雜的，在實(shí)際的應(yīng)用中應(yīng)根據(jù)不同的要求選用不同的調(diào)制方式。 作為一種現(xiàn)代調(diào)制技術(shù)， QAM 相比較傳統(tǒng)調(diào) 制技術(shù)，有著很多優(yōu)于傳統(tǒng)調(diào)制技術(shù)的特性，使得它在中、大容量數(shù)字微波通信系統(tǒng)、有線電視網(wǎng)絡(luò)高速數(shù)據(jù)傳輸、衛(wèi)星通信系統(tǒng)和蜂窩系統(tǒng)等系統(tǒng)中得到大規(guī)模應(yīng)用，大大提高系統(tǒng)通信質(zhì)量和通信效率。 本文的重要貢獻(xiàn) (1) 對 QAM 調(diào)制與解調(diào)技術(shù)的基本知識做了概括性的總結(jié)與分析。近 年來，隨著通信業(yè)務(wù)需求的增長，尋找頻譜利用率高的數(shù)字調(diào)制方式已成為數(shù)字通信系統(tǒng)設(shè)計(jì)，研究的主要目標(biāo)之一。所以，當(dāng)信噪比比較小的時候，理論計(jì)算的誤碼率的值要小于仿真得到的值。 在同樣信噪比時，誤符號率比誤比特率 （誤碼率） 要 大 。理論誤比特率 39。理論誤符號率 39。誤符號率 39。誤比特率 39。誤比特率和誤符號率 39。)。) xlabel(39。 title(39。k.39。k*39。ko39。 %理論誤符號率 通信系統(tǒng)綜合設(shè)計(jì)與實(shí)踐 第 19 頁 ber1=1/log2(M)*ser1。 %誤符號率 end P4=2*(11/sqrt(M))*qfunc(sqrt(3*snr1/(M1)))。 [err,ber(indx)]=biterr(msg,decmsg,log2(M))。 y=qamdemod(rx,M)。 %求每個符號的平均功率 for indx=1:length(EsN0) sigma=sqrt(spow/(2*snr1(indx)))。 %Gray 映射 msgmod=qammod(msg1,M)。 %信噪比轉(zhuǎn)換為線性值 msg=randint(1,nsymbol,M)。 %Gray 編碼規(guī)則 EsN0=5:20。 %每種信噪比下的發(fā)送符號數(shù) M=16。)。 title(39。Pe39。)。 xlabel(39。)。 semilogy(SNR_in_dB,Pe,39。 %統(tǒng)計(jì)誤碼率 end。 end。 for i=1:N if (y_output(i)~=info(i)), numoferr=numoferr+1。 xn=xn39。 xn=xn39。 xn(i,:)=T(yn(i)+1,:)。 n=length(yn)。1 1。 %一位四進(jìn)制碼元轉(zhuǎn)換為兩位二進(jìn)制碼元 T=[0 0。 I3=find(yn=)。 yn)。 yn(I1)=ones(size(I1))。 I1=find(yn= amp。 %設(shè)置門限電平 ， 判決 I0=find(yn )。 ymin=abs(ymax)。 通信系統(tǒng)綜合設(shè)計(jì)與實(shí)踐 第 16 頁 ymax=max(y)。 y=yn。 xn=xn39。 xn(nn/2+1:nn)]。 %I 分量 Q 分量并 /串轉(zhuǎn)換，最終恢復(fù)成碼元序列 xn nn=length(xn)。 Qn=Q(n)。 n=fix(n)。 N=length(y)/m。 Q=filtfilt(b,a,Q)