【正文】 ure(6) semilogy(SNRindB2,theo_err_awgn_prb,39。實際值（接收端）39。理論值（發(fā)射端）39。hold on。*39。axis([2 2 2 2])。)。),title(39。figure(5)for i=1:nb/2plot(idata(i),qdata(i),39。調制信號功率譜密度(Ray)39。fs39。subplot(2,2,2)。)。end%通過高斯信道再通過瑞利衰落信道 figure(4) subplot(2,2,1)plot(Ray_s),title(39。實際值（接收端）39。理論值（發(fā)射端）39。hold on。*39。axis([2 2 2 2])。)。),title(39。figure(3)for i=1:nb/2plot(idata(i),qdata(i),39。)。plot(s111),title(39。)。,fs),title(39。psd(h,s1,39。axis([0 500 5 5])。調制信號(Awgn)39。% 通過高斯信道figure(2)subplot(2,2,1)。解調輸出功率譜密度39。fs39。subplot(3,2,6)。)。plot(demodata1),title(39。)。,fs),title(39。psd(h,s,39。axis([0 500 3 3])。調制信號39。subplot(3,2,3)?；鶐盘柟β首V密度39。fs39。subplot(3,2,2)。)。plot(data0),title(39。h = 。 echo off。 theo_err_awgn_prb(i)=*erfc(sqrt(SNR))。% 理論曲線echo on。 echo off。 smld_symbol_awgn_err_prb(i)=ps1。for i=1:length(SNRindB1), [pb1,ps1]=cm_sm33(SNRindB1(i))。 end。 disp([ps,pb])。 % 比特誤碼率 smld_bit_ray_err_prb(i)=pb。SNRindB2=0::6。 end % 累計誤碼數% abs(demodatadata)求接收端和發(fā)射端% 數據差的絕對值，累計之后就是誤碼個數Ray_num_BER=sum(abs(Ray_demodatadata)) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 誤碼率計算%% 調用了cm_sm32()。 %為了顯示，將它變成波形信號（即傳輸一個1代表單位寬度的高電平）Ray_demodata1=zeros(1,nb/delta_T)。Ray_demodata(1:ml:(nb1))=Ray_idata4。 else Ray_qdata4(n)=0。 end Ray_qchsum(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Ray_ichsum(n)=0 Ray_idata4(n)=1。Ray_qdata4=zeros(1,nb/2)。 Ray_qdata2=Ray_s.*b。 for q=1:nb demodata1((q1)/delta_T+1:q/delta_T)=demodata(q)。 demodata(2:ml:nb)=qdata3。 endend% 將判決之后的數據存放進數組demodata=zeros(1,nb)。 if Awgn_qchsum(n)=0 qdata4(n)=1。 else idata4(n)=0。for n=1:nb/2 Awgn_ichsum(n)=sum(idata2((n1)/delta_T+1:n/delta_T))*delta_T。 endend % 為了顯示星座圖,將信號進行處理idata4=zeros(1,nb/2)。 end if sum(qdata2((n1)/delta_T+1:n/delta_T))=0 qdata3(n)=1。% 抽樣判決的過程，與0作比較，data=0,則置1，否則置0for n=1:nb/2 if sum(idata2((n1)/delta_T+1:n/delta_T))=0 idata3(n)=1。 idata3=zeros(1,nb/2)。 % 通過高斯信道再通過瑞利衰落信道% QPSK 解調部分% 解調部分（高斯信道）idata2=s1.*a。 s111=s1s。% 高斯信道 s1=awgn(s,SNR)。Ray_qdata=qdata1.*ray_qch。ray_qch=raylrnd(,1,nb/2/delta_T)。 ss=abs(fft(s))。endqdata1=qch.*b。for j1=1:nb/2 qch((j11)/delta_T+1:j1/delta_T)=qdata(j1)。 endidata1=ich.*a。 for i=1:nb/2 ich((i1)/delta_T+1:i/delta_T)=idata(i)。 qdata=datanrz(2:ml:nb)。 for q=1:nb data1((q1)/delta_T+1:q/delta_T)=datanrz(q)。 datanrz=data.*21。 % 限定t的取值范圍N=length(t)。 % 采樣頻率SNR=0。 % 傳輸的比特數delta_T=T/200。 % 載波頻率ml=2。八、 參考文獻[1] 李人厚、西安交通大學 [2] 肖明波，通信系統(tǒng)仿真原理與無線應用機械工業(yè)出版社[3] 萬永革編著，通信系統(tǒng)仿真原理與無線應用，北京科學出版社[4] 郭文彬、桑林，通信原理—基于MATLAB的計算機仿真，北京郵電大學出版社，[5]網上資料附件：主程序：T=1。我把程序發(fā)給老師后老師狠狠的教育了我，不應該依靠高版本的優(yōu)化功能把在低版本不能實現的功能強制執(zhí)行，所以我還有很多地方需要學習，包括對誤碼率進行計算的兩個函數我也是借鑒的網上的程序。本文設計出一個QPSK仿真模型，以衡量QPSK在理想信道，高斯白噪聲信道，以及先通過瑞利再通過高斯信道的性能，通過仿真，可以更好地了解QPSK系統(tǒng)的工作原理，而且為硬件的研制提供一定的參考作用。通過圖8可知瑞利信道下的誤碼率和誤比特率明顯高于高斯信道，并與理論值基本符合。六、 總結由此次仿真實驗可知，由于高斯信道和瑞利信道的影響，波形發(fā)生了明顯的變化，功率譜密度也發(fā)生變化，星座圖可反映在噪聲影響下產生的誤碼。由圖可知瑞利衰落信道下的誤碼率比高斯信道下的誤碼率高。實驗結果如圖8所示：圖6圖7圖8實驗結果分析： 由圖7可以得到瑞利衰落信道前后的星座圖，調制信號的曲線圖及其功率譜密度。星座圖反映可接收信號早高斯噪聲的影響下發(fā)生了誤碼，但是大部分還是保持了原來的特性。、高斯信道下的仿真，結果如圖5所示：圖4圖5實驗結果分析： 由圖5可以得到高斯信道下的調制信號，高斯噪聲，調制輸出功率譜密度曲線和QPSK信號的星座圖。程序見附錄五、 設計結果及分析根據圖1和圖2的流程框圖設計仿真程序，得出結果并且分析如下：、理想信道下的仿真，實驗結果如圖3所示：圖3實驗結果分析： 如圖上結果顯示，完成了QPSK信號在理想信道上的調制，傳輸，解調的過程，由于調制過程中加進了載波，因此調制信號的功率譜密度會發(fā)生變化。 MA