%?Raptor?codes中對RSD進行仿真
%?參考文獻：《噴泉碼的編譯碼算法研究與仿真實現(xiàn)》?西安電子科技大學?高寬

clc
clear?all;
close?all;


%?設置輸入符號數(shù)k
k?=?132;
%?設置允許的譯碼失敗概率
dlt?=?0.2;
%?設置常數(shù)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?????????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）?=