Simulink Simulation of DPSK (Differential PSK) passband signaling to measure the SER

[@WonYYang]

Watch kzbin.info/www/bejne/hXiqo4FnitqXhMU for Simulink Simulation of OFDM communication system to get its BER Watch kzbin.info/www/bejne/Z4vEXphpd8qhn9U for Simulink Simulation of coherent 4-FSK passband signaling to measure the SER Watch kzbin.info/www/bejne/a568gKakaMx9i9k for MATLAB Simulation of DPSK (Differential PSK) passband signaling to measure the BER Watch kzbin.info/www/bejne/eIeseZyhataioa8 for Simulink Simulation of DPSK (Differential PSK) passband signaling to measure the SER Watch kzbin.info/www/bejne/aoW3eoufjbiIb8U for Simulink Simulation of noncoherent FSK passband signaling to measure the SER Watch kzbin.info/www/bejne/h3a9nHqwg9yWfMk for MATLAB Simulation of noncoherent FSK passband signaling to measure the SER Watch kzbin.info/www/bejne/bXXbfGlpZrCqq80 for MATLAB Simulation of coherent FSK passband signaling to measure the SER Watch kzbin.info/www/bejne/apS0dJufadNki7s for MATLAB Simulation of QAM passband signaling to measure the BER Watch kzbin.info/www/bejne/gZWvoICfYst8i5Y for Simulink Simulation of QAM passband signaling to measure the SER Watch kzbin.info/www/bejne/nJu3nJWEnd51nZo for MATLAB Simulation of PSK passband signaling to measure the BER Watch kzbin.info/www/bejne/h5u0hZJmqtFqb80 for Simulink Simulation of PSK passband signaling to measure the BER The Simulink model file is available at kr.mathworks.com/matlabcentral/fileexchange/173550-simulink-simulation-of-dpsk-passband-signaling. %do_DPSK_passband_sim.m clear, clf b=2; M=2^b; % Number of bits per symbol and Modulation order Ns=160; % This is crucial to getting a realistic BER performance Ts=1e-5; T=Ts/Ns; % Symbol duration and Sampling interval A=sqrt(2/Ts); wc=2*pi*10/Ts; % Carrier frequency[rad/s] d=5; % Some delay to see how well the DPSK combats a delay Target_no_of_error=100; Transmitted_Signal_Power=1/Ts; Simulink_mdl='DPSK_passband_sim'; SNRdBs=[1:2:11]; Received_Signal_Power=Transmitted_Signal_Power*(1+10^(-SNRdBs(end)/10)*Ts/b/T) for i=1:length(SNRdBs) SNRdB=SNRdBs(i); EbN0dB=SNRdB-3; %10*log10(b); out=sim(Simulink_mdl,2e5*Ts); SERs(i)=ser(end,1); fprintf('%d-ary %s with SNRdB=%5.2f[dB] yields SER=%10.4e ',M,Simulink_mdl,SNRdB,SERs(i)) end SERs SNRdBt=[0:0.1:SNRdBs(end)]; EbN0dBt=SNRdBt-3; poset=prob_error(SNRdBt,'DPSK',b,'sym'); [BERt,SERt]=berawgn(EbN0dBt,'DPSK',M); %SERt=b*BERt; semilogy(SNRdBt,poset,'k', SNRdBs,SERs,'*', SNRdBt,SERt,'m:') title(sprintf('SER for %2d-DPSK passband signaling',M)) xlabel('SNR[dB]'), shg function mo=detector_DPSK(dth,M) phases=2*pi/M*[0:M-1]; [dmin,mmin]=min(abs(exp(j*dth)-exp(j*phases))); mo=mmin-1; function p=prob_error(SNRbdB,signaling,b,opt1,opt2) % Finds the symbol/bit error probability for given SNRbdB=Eb/(N0/2)[dB](Table 7.1) % Note that EbN0dB=SNRbdB-3. % Copyleft: Won Y. Yang, [e-mail], CAU for academic use only if nargin1))*Q(sqrt(b*SNRb)*sin(pi/M)); %Eq.(7.3.7) if lower(opt1(1))=='b'&b>1, p=p/b; end elseif signaling(1:3)=='DPS' %DPSK if b==1 p=exp(-SNRb/2)/2; %Eq.(7.4.9) else for i=1:NSNR EbN0=SNRb(i)/2; %b; % From the MATLAB built-in function berawgn() tol=1e-6; % tol=max(min(1e-4/EbN0.^5,1e-4),eps); int_dpsk=@(y,EbN0,M,b)exp(-b*EbN0*(sin(pi/M))^2./(1+sqrt(1-(sin(pi/M)).^2).*cos(y))); %?????? p(i)=1/pi*quad(int_dpsk,1e-5,pi*(1-1/M),tol,[],EbN0,M,b); end end if lower(opt1(1))=='b', p=p/b; end elseif signaling(1:3)=='QAM' L = 2^(ceil(b/2)); N = M/L; tmpL = 1-2*(L-1)/L*Q(sqrt(3*b/2/(L^2-1)*SNRb)); tmpN = 1-2*(N-1)/N*Q(sqrt(3*b/2/(N^2-1)*SNRb)); p = 1-tmpL.*tmpN; %Eq.(7.5.6) if lower(opt1(1))=='b'&b>1, p = p/b; end end end function I=Gauss_Hermite(f,N,varargin) [t,w]=Gausshp(N); ft=feval(f,t,varargin{:}); I=w*ft'; end