多符号差分球形译码算法的matlab仿真.包含仿真操作录像

clc close all clear all addpath('functions'); R=1e4;% sphere radios M=2; % MPSK symbols N=10; % multiple symbol window size Ns=(N-1)*round(1E5/(N-1));% number of symbols % option for known statistics or unknown stats stats_opt='known'; % number of basis functions in unknown stats I=3; EbN0_dB=0:2:40;% SNR range in dB EbN0=10.^(EbN0_dB/10); EsN0=log2(M)*EbN0; %fD=100;% doppler frequency fdTs=.03;% normalized doppler frequency ch_dis=1; % this loop scans the SNR range for ind=1:length(EbN0) ind nbits=0;%bit counter err_ms=0;% error counter err_cdd=0;% error counter err_coh=0;% error counter err_dcoh=0;% error counter clc EbN0(ind) % AWGN noise variance N0=1/EsN0(ind); % create and factorize C^-1 to find L Hermitian, Eq.2,3,4 Lh=computeL(N,fdTs,N0/2,stats_opt,I); % this loop keeps going to get a certain amount of errors ERR_TH=100; while err_ms< ERR_TH %|| err_dcoh<100 % binary data,0/1 xb=bits(log2(M)*Ns); % QPSK symbols v=bin2mpsk(xb,M); % DPSK modulation s=dpsk(v); Nd=length(s); % correlated channel coefficients, Jakes model h=flat_cos(Nd,fdTs,ch_dis); % AWGN noise CN(0,N0) n=cxn(Nd,N0); % received signals r=h.*s+n;% diff symbols %r_coh=h(1:Ns).*v.'+n(1:Ns);%coherent % multiple symbol sphere detection for k=1:N-1:Nd-N+1 rN=r(k:k+N-1).'; U=conj(Lh*diag(rN)); shat=MSDSD(U,M,R); vh_ms(k:k+N-2)=diff_detector(shat); end % conventional non-coherent detection vh_cdd=diff_detector(r); % binary detection bh_cdd=mpsk2bin(vh_cdd,M); bh_ms=mpsk2bin(vh_ms,M); err_ms=err_ms+sum(abs(xb-bh_ms)); err_cdd=err_cdd+sum(abs(xb-bh_cdd)); nbits=log2(M)*Ns+nbits;%Ns is number of symbols end % compute practical BER ber_ms(ind)=err_ms/nbits;% multiple symbol ber_cdd(ind)=err_cdd/nbits;% conventional end % BER theory value from book "digital comm on fading channels" page 293&259 Pb_coh=.5*(1-sqrt(EbN0./(1+EbN0))); roh=besselj(0,2*pi*fdTs); if M==2 Pb_cdd=0.5*(1+EbN0*(1-roh))./(1+EbN0); elseif M==4 Pb_cdd=.5*(1+EbN0*(1-roh))./(1+EbN0);% this is not exact! end for k=1:length(EsN0) N0=1/EsN0(k); Pb_ms(k)=ber_msd(M,N,N0,fdTs); end %% plot % marker mk=['r-s';'k-*';'g--';'b-o';'c--';'m--']; semilogy(EbN0_dB,ber_ms,mk(1,:),'LineWidth',2,'MarkerSize',6); hold on semilogy(EbN0_dB,ber_cdd,mk(2,:),'LineWidth',2,'MarkerSize',6); semilogy(EbN0_dB,Pb_cdd,mk(3,:),'LineWidth',2,'MarkerSize',6); semilogy(EbN0_dB,Pb_ms,'b-','LineWidth',2,'MarkerSize',6); grid on xlabel('E_b/N_0,[dB]') ylabel('Bit Error Rate') leg1=['MSD ',stats_opt,', N=',num2str(N)]; legend(leg1,'Simulation CDD','Theory CDD','Theory MSDD'); set(gca,'XTick',EbN0_dB(1):5:EbN0_dB(end))