可读性极低,而且其实也没必要用R语言写,图个乐罢了
p=c(0.4,0.2,0.2,0.1,0.1)###输入形如c(0.4,0.2,0.2,0.1,0.1)的概率向量,即每个待编码消息的发生概率
p1=p###将概率向量另存,最后计算编码效率要用
mazijuzhen=matrix(,nrow=length(p),ncol=length(p)-1)###码字矩阵:第i行对应向量p的第i个分量所对应的那个待编码消息的编码后的码字
group=matrix(c(1:length(p),rep(NA,length(p)*(length(p)-1))),nrow=length(p),ncol=length(p))###初始分组:每一行代表一组,每个行向量的所有分量代表此组的所有元素,初始时,有多少个待编码消息就分多少个组,每组只有一个待编码消息,以整数i代表向量p的第i个分量所对应的那个待编码消息
i=1###开始编码
for(i in 1:(length(p)-1))
{
orderp=order(p,decreasing = FALSE)###orderp的分量依次是:p的最小分量的下标,p的第二小分量的下标。。。
mazijuzhen[group[orderp[1],],i]=0###给概率最小的两个消息组编上0和1
mazijuzhen[group[orderp[2],],i]=1
group[min(c(orderp[1],orderp[2])),]=c(na.omit(group[min(c(orderp[1],orderp[2])),]),na.omit(group[max(c(orderp[1],orderp[2])),]),rep(NA,length(p)-length(c(na.omit(group[min(c(orderp[1],orderp[2])),]),na.omit(group[max(c(orde