【文章內容簡介】 Awgn_qchsum(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Awgn_qchsum(n)=0 qdata4(n)=1。 else qdata4(n)=0。 endend% 將判決之后的數據存放進數組demodata=zeros(1,nb)。demodata(1:ml:(nb1))=idata3。 % 存放奇數位demodata(2:ml:nb)=qdata3。 % 存放偶數位%為了顯示，將它變成波形信號（即傳輸一個1代表單位寬度的高電平）demodata1=zeros(1,nb/delta_T)。 % 創(chuàng)建一個1*nb/delta_T的零矩陣for q=1:nb demodata1((q1)/delta_T+1:q/delta_T)=demodata(q)。 % 將極性碼變成對應的波形信號end % 累計誤碼數% abs(demodatadata)求接收端和發(fā)射端% 數據差的絕對值，累計之后就是誤碼個數Awgn_num_BER=sum(abs(demodatadata)) %%%%%%%%%%%%%%%%%%%% 解調部分（瑞利+高斯）Ray_idata2=Ray_s.*a。 % 這里面其實隱藏了一個串并轉換的過程Ray_qdata2=Ray_s.*b。 % 對應的信號與正余弦信號相乘% Ray_idata3=zeros(1,nb/2)。 % 建立1*nb數組，以存放解調之后的信號% Ray_qdata3=zeros(1,nb/2)。% 抽樣判決的過程，與0作比較，data=0,則置1，否則置0% for n=1:nb/2% if Ray_sum(Ray_idata2((n1)/delta_T+1:n/delta_T))=0% Ray_idata3(n)=1。% else Ray_idata3(n)=0。 % end% if Ray_sum(Ray_qdata2((n1)/delta_T+1:n/delta_T))=0% Ray_qdata3(n)=1。% else Ray_qdata3(n)=0。% end% end % 為了顯示星座圖,將信號進行處理Ray_idata4=zeros(1,nb/2)。Ray_qdata4=zeros(1,nb/2)。for n=1:nb/2 Ray_ichsum(n)=sum(idata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Ray_ichsum(n)=0 Ray_idata4(n)=1。 else Ray_idata4(n)=0。 end Ray_qchsum(n)=sum(qdata2((n1)/delta_T+1:n/delta_T))*delta_T。 if Ray_qchsum(n)=0 Ray_qdata4(n)=1。 else Ray_qdata4(n)=0。 endend % 將判決之后的數據存放進數組Ray_demodata=zeros(1,nb)。Ray_demodata(1:ml:(nb1))=Ray_idata4。 % 存放奇數位Ray_demodata(2:ml:nb)=Ray_qdata4。 % 存放偶數位%為了顯示，將它變成波形信號（即傳輸一個1代表單位寬度的高電平）Ray_demodata1=zeros(1,nb/delta_T)。 % 創(chuàng)建一個1*nb/delta_T的零矩陣for q=1:nb Ray_demodata1((q1)/delta_T+1:q/delta_T)=Ray_demodata(q)。 % 將極性碼變成對應的波形信號end % 累計誤碼數% abs(demodatadata)求接收端和發(fā)射端% 數據差的絕對值，累計之后就是誤碼個數Ray_num_BER=sum(abs(Ray_demodatadata)) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 誤碼率計算%% 調用了cm_sm32()。和cm_sm33(）函數%%聲明： 函數聲明在另外倆個M文件中%%作用： cm_sm32()用于瑞利信道誤碼率的計算%% cm_sm33()用于高斯信道誤碼率的計算%% ecoh on/off 作用在于決定是否顯示指令內容%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%SNRindB1=0:1:6。SNRindB2=0::6。% 瑞利衰落信道 for i=1:length(SNRindB1), [pb,ps]=cm_sm32(SNRindB1(i))。 % 比特誤碼率 smld_bit_ray_err_prb(i)=pb。 smld_symbol_ray_err_prb(i)=ps。 disp([ps,pb])。 echo off。 end。% 高斯信道 echo on。for i=1:length(SNRindB1), [pb1,ps1]=cm_sm33(SNRindB1(i))。 smld_bit_awgn_err_prb(i)=pb1。 smld_symbol_awgn_err_prb(i)=ps1。 disp([ps1,pb1])。 echo off。end。% 理論曲線echo on。for i=1:length(SNRindB2), SNR=exp(SNRindB2(i)*log(10)/10)。 % 信噪比 theo_err_awgn_prb(i)=*erfc(sqrt(SNR))。 % 高斯噪聲理論誤碼率 theo_err_ray_prb(i)=*(11/sqrt(1+1/SNR))。 % 瑞利衰落信道理論誤碼率 echo off。end。%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%h = 。 % 類似于C語言的宏定義，方便以下的調用 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 輸出顯示部分%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 第一部分（理想）figure(1)subplot(3,2,1)。plot(data0),title(39?；鶐盘?9。)。axis([0 20000 2 2])。subplot(3,2,2)。psd(h,data1,39。fs39。,fs),title(39。基帶信號功率譜密度39。)。subplot(3,2,3)。plot(s),title(39。調制信號39。)。axis([0 500 3 3])。sub