永利皇宫官网网址

网站最新上线:http://www.hslogic.com/

matlab代做仿真案例-多进制LDPC编译码仿真

上一篇 / 下一篇  2017-04-19 17:49:38

多进制LDPC编译码仿真
演示地址:

clc;
clear;
close all;
warning off;
addpath 'functions\'

q         = 4;                      %定义伽罗达域,支持4,8,16
nbits     = log2(q);                %多进制符号值
EbNo      = [0  0.5  1  1.5 2];   %信噪比

%以下两个值设置小点,这样仿真可以快点,值较大的时候,仿真速度非常慢,但性能较好
M         = 24;              %校验数
N         = 32;             %矩阵长度
t         = 3;               %权值                 
h         = ldpc_generate(M,N,t,q,1); %产生矩阵H
[H,G]     = ldpc_h2g(h,q);                    %产生矩阵G

for kk=1:length(EbNo)
    kk
    totalNumErr = 0;
    count       = 0;
    SNR         = 10^(EbNo(kk)/10);
    sigma       = 1/(sqrt(SNR));    

    while (totalNumErr < 100000)%仿真越长,效果越OK
        kk
        totalNumErr
        
        %产生一组随机数
        x  = floor(rand(1,size(G,1))*q); 
        
        %进行多进制准循环LDPC编码
        y  = func_Ldpc_enc(x,G,q);   
        
        %多进制信道处理
        yb = (fliplr(de2bi(y,nbits)))';      
        yb = yb(:);                          
        zb = 2*yb-1;                      
        zb = zb + sigma*randn(size(zb));      

        f1=1./(1+exp(-2*zb/sigma^2));        
        f1 = f1(:);                         
        f1 = reshape(f1,nbits,length(y));                    
        f0=1-f1;
        junk = ones(q,length(y));           
        [v0, v1, pp] = func_message_passing(f0,f1,junk);    
        
        [z_hat,success, k] = func_Ldpc_dec_log(pp,H,q);
        x_hat = z_hat(size(G,2)+1-size(G,1):size(G,2));
        x_hat = x_hat';    
        
        %统计误码率
        totalNumErr = totalNumErr+(sum (x_hat~=x'));
       
        count = count + 1;
    end
    
    BERs(kk)=totalNumErr/(count*length(x));
end

figure;
semilogy(EbNo,BERs,'r-o');
xlabel('EbN0');
ylabel('BER');
grid on;

if q == 4
   save ERR4_log.mat EbNo BERs
end
if q == 8
   save ERR8_log.mat EbNo BERs
end
if q == 16
   save ERR16_log.mat EbNo BERs
end


TAG: Matlab matlab MATLAB

引用 删除 邓依   /   2018-11-29 10:53:58
请问您有完整的多进制LDPC码编译码matlab程序吗?我正在找这方面的程序
 

评分:0

我来说两句

显示全部

:loveliness: :handshake :victory: :funk: :time: :kiss: :call: :hug: :lol :'( :Q :L ;P :$ :P :o :@ :D :( :)

Open Toolbar