鲁棒孤子分布的仿真程序

% Raptor codes中对RSD进行仿真
% 参考文献：《喷泉码的编译码算法研究与仿真实现》 西安电子科技大学 高宽

clc
clear all;
close all;


% 设置输入符号数k
k = 132;
% 设置允许的译码失败概率
dlt = 0.2;
% 设置常数c
% 给出c的取值范围
c_min = sqrt（k）/log（k/dlt）/（k-1）;
c_max = sqrt（k）/log（k/dlt）/2;
c = c_min;
% ripple大小的期望值
R = c * log（k/dlt） * sqrt（k）;

for i = 1:k
    if i < k/R - 1
        Tau（i） = R/（i*k）;
    elseif i == k/R
        Tau（i） = R*log（R/dlt）/k;
    else
        Tau（i） = 0;
    end
end

for j = 1:k
   if j == 1
       rho（j） =