LDPC的matlab 仿真 BP算法 BPSK调制

LDPC仿真，matlab BP算法 LDPC的matlab 仿真 BP算法 BPSK调制

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个文件信息