【正文】 end 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% 將判決之后的數(shù)據(jù)存放進數(shù)組demodata=zeros(1,nb)。demodata(1:ml:(nb1))=idata3。 % 存放奇數(shù)位demodata(2:ml:nb)=qdata3。 % 存放偶數(shù)位%為了顯示，將它變成波形信號（即傳輸一個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 % 累計誤碼數(shù)% abs(demodatadata)求接收端和發(fā)射端% 數(shù)據(jù)差的絕對值，累計之后就是誤碼個數(shù)Awgn_num_BER=sum(abs(demodatadata)) %%%%%%%%%%%%%%%%%%%% 解調(diào)部分（瑞利+高斯）Ray_idata2=Ray_s.*a。 Ray_qdata2=Ray_s.*b。 % Ray_idata3=zeros(1,nb/2)。 % 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 % 將判決之后的數(shù)據(jù)存放進數(shù)組Ray_demodata=zeros(1,nb)。Ray_demodata(1:ml:(nb1))=Ray_idata4。 % 存放奇數(shù)位Ray_demodata(2:ml:nb)=Ray_qdata4。 % 存放偶數(shù)位%為了顯示，將它變成波形信號（即傳輸一個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 % 累計誤碼數(shù)% abs(demodatadata)求接收端和發(fā)射端% 數(shù)據(jù)差的絕對值，累計之后就是誤碼個數(shù)Ray_num_BER=sum(abs(Ray_demodatadata)) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 誤碼率計算%% 調(diào)用了cm_sm32()。和cm_sm33(）函數(shù)%%聲明： 函數(shù)聲明在另外倆個M文件中%%作用： cm_sm32()用于瑞利信道誤碼率的計算%% cm_sm33()用于高斯信道誤碼率的計算%% ecoh on/off 作用在于決定是否顯示指令內(nèi)容%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%SNRindB1=0:1:8。SNRindB2=0::8。% 瑞利衰落信道 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語言的宏定義，方便以下的調(diào)用 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 第一部分（理想）figure(1)subplot(212)。psd(h,data0,39。fs39。,fs),title(39?；鶐盘柟β首V密度39。)。subplot(211)。plot(data0),title(39?；鶐盘?9。)。axis([0 20000 2 2])。figure(2)subplot(212)。psd(h,s,39。fs39。,fs),title(39。調(diào)制信號功率譜密度39。)。subplot(211)。plot(s),title(39。調(diào)制信號39。)。axis([0 500 3 3])。figure(3)subplot(212)。psd(h,demodata1,39。fs39。,fs),title(39。解調(diào)輸出功率譜密度39。)。subplot(211)。plot(demodata1),title(39。解調(diào)輸出39。)。axis([0 20000 2 2])。figure(4)subplot(313)。psd(h,s1,39。fs39。,fs),title(39。調(diào)制信號功率譜密度(Awgn)39。)。subplot(312)。plot(s1),title(39。調(diào)制信號(Awgn)39。)。axis([0 500 5 5])。subplot(311)。plot(s111),title(39。高斯噪聲曲線39。)。axis([0 2000 5 5])。figure(5)for i=1:nb/2plot(idata(i),qdata(i),39。r+39。),title(39。QPSK信號星座圖（Awgn）39。)。hold on。legend(39。理論值（發(fā)射端）39。,39。實際值（接收端）39。)。axis([2 2 2 2])。plot(Awgn_ichsum(i),Awgn_qchsum(i),39。*39。)。hold on。endfigure(6) subplot(111) semilogy(SNRindB2,theo_err_awgn_prb,39。r39。),title(39。誤碼率曲線39。)。hold on。 semilogy(SNRindB1,smld_bit_awgn_err_prb,39。r*39。)。hold on。 semilogy(SNRindB2,theo_err_ray_prb)。hold on。 semilogy(SNRindB1,smld_bit_ray_err_prb,39。*39。)。 xlabel(39。Eb/No39。)。ylabel(39。BER39。)。 legend(39。理論AWGN39。,39。仿真AWGN39。,39。理論Rayleigh39。,39。仿真Rayleigh39。)。二．QPAK和OQPSK對比仿真程序：%QPSK amp。 OQPSKclear all。close all。M=4。Ts=1。fc=10。N_sample=16。N_num=100。dt=1/fc/N_sample。t=0:dt:N_num*Tsdt。T=dt*length(t)。py1f=zeros(1,length(t))。 %功率譜密度1py2f=zeros(1,length(t))。 %功率譜密度2for PL=1:100 %輸入100段N_num個碼字的波形，為了使功率譜密度看起來更加平滑， %可以取這100段信號功率譜密度的平均 d1=sign(randn(1,N_num))。 d2=sign(randn(1,N_num))。 gt=ones(1,fc*N_sample)。 %QPSK調(diào)制 s1=sigexpand(d1,fc*N_sample)。 s2=sigexpand(d2,fc*N_sample)。 b1=conv(s1,gt)。 b2=conv(s2,gt)。 s1=b1(1:length(s1))。 s2=b2(1:length(s2))。 st_qpsk=s1.*cos(2*pi*fc*t)s2.*sin(2*pi*fc*t)。 s2_delay=[ones(1,N_sample*fc/2) s2(1:endN_sample*fc/2)]。 st_oqpsk=s1.*cos(2*pi*fc*t)s2_delay.*sin(2*pi*fc*t)。 %經(jīng)過帶通后，在經(jīng)過非線性電路 [f y1f]=T2F(t,st_qpsk)。 [f y2f]=T2F(t,st_oqpsk)。 [t y1]=bpf(f,y1f,fc1/Ts,fc+1/Ts)。 [t y2]=bpf(f,y2f,fc1/Ts,fc+1/Ts)。 subplot(221)。 plot(t,y1)。xlabel(39。t39。)。ylabel(39。QPSK波形39。)。 axis([5 15 ])。title(39。經(jīng)過帶通后的波形39。)。 subplot(222)。 plot(t,y2)。xlabel(39。t39。)。ylabel(39。OQPSK波形39。)。 axis([5 15 ])。title(39。經(jīng)過帶通后的波形39。)。 %經(jīng)過非線性電路 y1=*tanh(2*y1)。 y2=*tanh(2*y2)。 [f y1f]=T2F(t,y1)。 [f y2f]=T2F(t,y2)。 py1f=py1f+abs(y1f).^2/T。 %QPSK不同段信號功率譜密度相加 py2f=py2f+abs(y2f).^2/T。 %OQPSK不同段信號功率譜密度相加endpy1f=py1f/100 %QPSK100段功率譜密度平均py2f=py2f/100 %OQPSK100段功率譜密度平均subplot(223)。plot(f,10*log10(py1f))。xlabel(39。f39。)。ylabel(39。QPSK功率譜密度(dB/Hz)39。)。title(39。經(jīng)過非線性電路后的功率譜密度39。)。axis([15 15 30 10])。subplot(224)。plot(f,10*log10(py2f))。xlabel(39。f39。)。ylabel(39。OQPSK功率譜密度(dB/Hz)39。)。title(39。經(jīng)過非線性電路后的功率譜密度39。)。axis([15 15 30 10])。fig