LDPC在AWGN性能仿真

clear
clc
clear

%--------------------------------------------------------------------------
%-----Set-Simulation-Parameters--------------------------------------------
%--------------------------------------------------------------------------

snrs = [10 13 15];      %SNR values
codeRate = 9/10; %Possible values for codeRate are 1/4 1/3 2/5 1/2 3/5 2/3 3/4 4/5 5/6 8/9 and 9/10. The block length of the code is 64800
mod_order = 4;  %PSK Modulation Order
frames = 2500;  %Number of frames （fame size is 64800 bits） to be simulated

%--------------------------------------------------------------------------

rounds = size（snrs2）;

messageLength = round（64800*codeRate）;

for run = 1:1:rounds

framepattern = [];    
    
snrvalue = snrs（run）;

H = dvbs2ldpc（codeRate）;

%spy（H）;   % Visualize the location of nonzero elements in H.

errors = 0;

hEnc = comm.LDPCEncoder（H）;
hMod = comm.PSKModulator（mod_order ‘BitInput‘true）;
hChan = comm.MIMOChannel（‘MaximumDopplerShift‘ 0 ‘NumTransmitAntennas‘1‘NumReceiveAntennas‘1 ‘TransmitCorrelationMatrix‘ 1 ‘ReceiveCorrelationMatrix‘ 1 ‘PathGainsOutputPort‘ true）;
hAWGN = comm.AWGNChannel（‘NoiseMethod‘‘Signal to noise ratio （SNR）‘‘SNR‘snrvalue）;
% hDemod = comm.PSKDemodulator（4 ‘BitOutput‘true‘DecisionMethod‘‘Approximate log-likelihood ratio‘...
%                              ‘Variance‘ 1/10^（hChan.SNR/10））;
hDemod = comm.PSKDemodulator（4 ‘BitOutput‘true‘DecisionMethod‘‘Approximate log-likelihood ratio‘）;                         
hDec = comm.LDPCDecoder（H‘DecisionMethod‘ ‘Soft decision‘）;
%hError = comm