shannon编码的matlab代码

clc
clear all
P=[0.5540.2250.1120.109];

p=-sort（-P）;%降序排列
n=length（p）;
pp=cumsum（p）;
pp=[0pp];%概率累加
bi=ceil（log2（1./p））;%计算编码位数
L=bi*p‘
H_x=-p*log2（p‘）
R=H_x/L

 %将十进制转化为二进制
for i=1:n
    N=bi（i）;
    count=0;
    tempnum=pp（i）;
    
    
    while （N）
        count=count+1;
        if（count>N）
            b