基于MATLAB的多进制LDPC译码算法的仿真

最新推荐文章于 2024-06-24 10:20:22 发布

fpga和matlab

最新推荐文章于 2024-06-24 10:20:22 发布

阅读量3.9k

点赞数 4

分类专栏： MATLAB 板块4:编码译码文章标签：算法多进制LDPC

本文链接：https://blog.csdn.net/ccsss22/article/details/123540847

版权

MATLAB 同时被 2 个专栏收录

1606 篇文章 1630 订阅

订阅专栏

板块4:编码译码

69 篇文章 43 订阅

订阅专栏

本文介绍了一种使用伽罗达域的多进制准循环LDPC编码在不同信噪比下的仿真过程，通过详细步骤展示了编码、信道处理、解码及误码率计算。结果以 EbN0 与 BER 曲线的形式展示，适用于不同 q 值的性能评估。

摘要由CSDN通过智能技术生成

1.仿真预览

2.部分核心代码

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 < 40)%仿真越长，效果越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

C33