费诺编码
功能简介:首先输入信源符号的个数,然后输入各个信源符号的概率。(概率之和等于1)
主程序: 李松林
#include<stdio.h> 湖北师范学院
#include<math.h>
#define G 20
int function1(int c,int M);
void function2(int c,int num1,int M,int cs);
void function3(int c,int M);
float a[G]={0},tmp=0,m[G]={0},k=0,H=0,num[G]={0},sum1=0,sum2=0;
int K[G]={0},i,j,N,s[G][G],c,cs=0,num1=0,ks[G],zj[G],zh[G];
void main()
{
loop:printf("请输入信源符号个数N:");
scanf("%d",&N);
for(i=0;i<N;i++)
{
scanf("%f",&a[i]);
tmp+=a[i];
}
if(tmp<0.9999||tmp>1.0001)
{
printf("输入的数据不符合要求,请重新输入\n");
tmp=0;
goto loop;
}
else
{
for(i=0;i<N-1;i++)
for(j=i+1;j<N;j++)
{
if(a[i]<a[j]) //从大到小排序
{
tmp=a[i];
a[i]=a[j];
a[j]=tmp;
}
}
for(i=0;i<N;i++)
{
m[i]=-log10(a[i])/log10(2.0);//求出-log p(ai)并保存在数组m里
}
/
function3(0,N);
for(i=0;i<N;i++)
{
H+=a[i]*m[i];//求信源熵
k+=a[i]*K[i];//平均码长
}
printf("信源消息符号ai 符号概率p(ai) 码长Ki 二元码字\n");
for(i=0;i<N;i++)
{
printf(" a%d %-4.2f %d ",i+1,a[i],K[i]);
for(j=0;j<K[i];j++)
printf("%d",s[i][j]);
printf("\n");
}
printf("信源熵H=%5.3f bit/符号 费诺码的平均码长k=%5.3f 码元/符号 编码效率n=%4.1f%% \n",H,k,100*(H/k));
}
}
int function1(int c,int M)//求出分界点
{
int d=c;
for(i=c;i<M-1;i++)
{
for(j=c;j<i+1;j++)
sum1+=a[j];
for(j=i+1;j<M;j++)
sum2+=a[j];
num[i]=fabs(sum1-sum2);//取绝对值
sum1=0;
sum2=0;
}
for(i=c;i<M-1;i++)
{
if(num[d]>num[i])
d=i;
}
d+=1;
return d;
}
void function2(int c,int num1,int M,int cs)//进行分组并计算码长、码字
{
for(i=c;i<num1;i++)
{
s[i][cs]=0;
K[i]+=1;
}
for(j=num1;j<M;j++)
{
s[j][cs]=1;
K[j]+=1;
}
}
void function3(int c,int M)
{
num1=function1(c,M);
function2(c,num1,M,cs);
ks[cs]=c;
zj[cs]=num1;
zh[cs]=M;
cs++;
if((zj[cs-1]-ks[cs-1])>1)
{
function3(ks[cs-1],zj[cs-1]);
}
if((zh[cs-1]-zj[cs-1])>1)
{
function3(zj[cs-1],zh[cs-1]);
}
cs--;
}
第一次输入:6 0.32 0.22 0.18 0.16 0.08 0.04
结果:
第二次输入:7 0.20 0.19 0.18 0.17 0.15 0.10 0.01
结果:
第三次输入:8 0.1 0.18 0.4 0.05 0.06 0.1 0.07 0.04
结果: