• 大小: 360KB
    文件类型: .rar
    金币: 2
    下载: 1 次
    发布日期: 2021-06-10
  • 语言: Matlab
  • 标签: LDPC  

资源简介

多进制纠错码LDPC编译码matlab算法仿真程序,有注释说明

资源截图

代码片段和文件信息

clc;
clear all;
m=250;n=[3005001000];%校验矩阵的n是码长
frame=500;
Eb_N0=[0:0.5:5];

for ii=1:3
R(ii)=(n(ii)-m)/n(ii)
H(ii)=getH(mn(ii))       %load H; H-->G
[G(ii)valid]=H2G(H(ii)); 
while valid==0                   %valid作为校验矩阵是否为满秩的标志若为非满秩(valid=0)则返回重新利用1A随机构造
H(ii)=getH(mn(ii));                     %构造校验矩阵H
[G(ii)valid]=H2G(H(ii));                %将校验矩阵H转化为生成矩阵G
end
end




    for i=1:length(Eb_N0)
     EN=10^(Eb_N0(i)/10)
    sigma=sqrt(1./(2*EN*R));
    ber(i)=0;

      for num=1:frame  
        num
        x = (sign(randn(1size(G1)))+1)/2; % random bits
        y = mod(x*G2);                     % coding 
        bpskmod =bpsk(y);                          %BPSK modulation
        
        z=bpskmod + sigma*randn(1length(bpskmod));   % AWGN transmission
        f1=1./(1+exp(-2*z/(sigma^2)));         % likelihoods
        f0=1-f1;
        [z_hat success k] = ldpc_decode(zf0f1H);
        x_hat = z_hat(size(G2)+1-size(G1):size(G2));%译码
        x_hat = x_hat‘;                           %转置,保证x_hat与x的矩阵维数一致 
        %erro=find(x_hat~=x);                     %寻找错误信息位
        err(num)=length(find(x_hat~=x))          %求出错误信息位位数
        
        end   %for num
        ber(i)=sum(err)/(length(x)*frame)        %计算比特误码率BER
        
    end %for i
   

       
semilogy(Eb_N0ber‘b-o‘);
xlabel(‘Eb/No(dB)‘‘fontweight‘‘bold‘);
ylabel(‘BER‘‘fontweight‘‘bold‘);
%title(‘ldpc在awgn下性能仿真‘‘fontsize‘12‘fontweight‘‘bold‘‘fontname‘‘黑体‘)
%c=legend(‘r=3/4L=512‘1);
hold on
grid on;
%hold off;
 %for j

        
        
        
        
        
        
        
        

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1731  2011-05-22 19:52  LDPC_程序1\1.asv

     文件       1767  2011-05-22 21:21  LDPC_程序1\1.m

     文件        347  2011-05-20 00:22  LDPC_程序1\bpsk.m

     文件      64037  2011-05-21 21:32  LDPC_程序1\duibi.jpg

     文件       3302  2011-05-22 20:05  LDPC_程序1\frame.asv

     文件      12010  2011-05-20 21:04  LDPC_程序1\G.mat

     文件        333  2010-05-09 11:31  LDPC_程序1\getG.asv

     文件       4255  2011-05-22 15:54  LDPC_程序1\getH.asv

     文件       4219  2010-05-08 10:14  LDPC_程序1\getH.m

     文件       8440  2011-05-22 22:09  LDPC_程序1\H.mat

     文件        847  2010-05-12 19:53  LDPC_程序1\H2G.m

     文件       3523  2011-05-20 20:46  LDPC_程序1\ldpc_decode.asv

     文件       3523  2011-05-21 22:03  LDPC_程序1\ldpc_decode.m

     文件       1457  2011-05-22 21:57  LDPC_程序1\ldpc_decodeber.asv

     文件       1465  2011-05-22 22:01  LDPC_程序1\ldpc_decodeber.m

     文件          0  2011-05-21 17:38  LDPC_程序1\Untitled.asv

     文件      63282  2011-05-21 17:13  LDPC_程序1\untitled.jpg

     文件      65922  2011-05-20 22:30  LDPC_程序1\图.jpg

     文件     236278  2010-05-09 08:59  LDPC_程序1\图r=0.5\252_504_1000.bmp

     文件     252518  2010-05-09 11:27  LDPC_程序1\图r=0.5\252_504_10000(补充).bmp

     文件     236278  2010-05-09 11:00  LDPC_程序1\图r=0.5\252_504_10000.bmp

     文件     236278  2010-05-11 12:48  LDPC_程序1\图r=0.5\256_512_20000(补充).bmp

     文件       6784  2010-05-15 11:11  LDPC_程序1\图r=0.5\256_512_3000.fig

     文件     236278  2010-05-13 12:15  LDPC_程序1\图r=0.5\256_512_5000.bmp

     文件       8963  2010-05-14 18:24  LDPC_程序1\图r=0.5\256_512_5000.fig

     文件     252518  2010-05-09 15:26  LDPC_程序1\图r=0.5\96_192_10000.bmp

    ..A.SH.     27136  2010-05-13 15:03  LDPC_程序1\图r=0.5\Thumbs.db

     文件       6776  2010-05-15 16:49  LDPC_程序1\图r=0.5\untitled.fig

     文件       6740  2010-05-16 18:37  LDPC_程序1\图r=0.667\1.fig

     文件       6790  2010-05-14 20:52  LDPC_程序1\图r=0.667\2.fig

............此处省略38个文件信息

评论

共有 条评论