噴泉碼生成矩陣草稿_01

clc;clear;
% k =  ;
x = rand(1,k)<0.5;  %隨機生成碼元序列x

%***************初始度分佈函數開始***************
InitD = 2：K;
if find(d==InitD)
	P = 1/d/(d-1);
elseif d==1
	P = 1/K;
end
PP(d) = P;
%***************初始度分佈函數結束****************
%***************修正度分佈函數開始****************
%K =  ,o = ;
S = c*sqrt(K)*log(K/o);
Initd = 1:round(K/S)-1; %round(K/S)
if find(d==Initd)
	t = S/K/d;
elseif d==round(K/S)      %round(K/S)
	t = (S/K)*log(S/o);
else
	t = 0;
end
tt(d) = t;
Z = 0;
for d = 1:K
	Z = Z + PP(d) + tt(d);
end
Pa = (PP(d)+tt(d))/Z;
%***************修正度分佈函數結束_***************