费诺编码的matlab实现Word下载.docx
- 文档编号:19327579
- 上传时间:2023-01-05
- 格式:DOCX
- 页数:7
- 大小:103.09KB
费诺编码的matlab实现Word下载.docx
《费诺编码的matlab实现Word下载.docx》由会员分享,可在线阅读,更多相关《费诺编码的matlab实现Word下载.docx(7页珍藏版)》请在冰豆网上搜索。
实验一:
1.实验目的
1)掌握费诺编码的思想和具体方法。
2)用MATLAB语言实现费诺编码。
2.实验原理及编码思想:
费诺编码属于概率匹配编码,但不是最佳的编码方法。
在编N进制码时首先将信源消息符号按其出现的概率依次由大到小排列开来,并将排列好的信源符号按概率值分N大组,使N组的概率之和近似相同,并对各组赋予一个N进制码元0、1……N-1。
之后再针对每一大组内的信源符号做如上的处理,即再分为概率和相同的N组,赋予N进制码元。
如此重复,直至每组只剩下一个信源符号为止。
此时每个信源符号所对应的码字即为费诺码。
具体过程如下:
[1]将信源消息符号按其出现的概率大小依次排列:
P1>
=P2>
=…>
=Pn。
[2]依次排列的信源符号按概率值分为两大组,使两个组的概率之和近似相同,并对各组赋予一个二进制码元“0”和“1”。
[3]使划分后的两个组的概率之和近似相同,并对各组赋予一个二进制符号“0”和“1”。
[4]如此重复,直至每个组只剩下一个信源
4.程序代码
clc;
clear;
A=[0.19,0.18,0.17,0.16,0.13,0.10,0.06,0.01];
A=fliplr(sort(A));
%降序排列
[m,n]=size(A);
fori=1:
n
B(i,1)=A(i);
%生成B的第1列
end
%生成B第2列的元素
a=sum(B(:
1))/2;
fork=1:
n-1
ifabs(sum(B(1:
k,1))-a)<
=abs(sum(B(1:
k+1,1))-a)
break;
end
n%生成B第2列的元素
ifi<
=k
B(i,2)=0;
else
B(i,2)=1;
%生成第一次编码的结果
END=B(:
2)'
;
END=sym(END);
%生成第3列及以后几列的各元素
j=3;
while(j~=0)
p=1;
while(p<
=n)
x=B(p,j-1);
forq=p:
ifx==-1
ifB(q,j-1)==x
y=1;
continue;
y=0;
ify==1
q=q+1;
ifq==p|q-p==1
B(p,j)=-1;
ifq-p==2
B(p,j)=0;
END(p)=[char(END(p)),'
0'
];
B(q-1,j)=1;
END(q-1)=[char(END(q-1)),'
1'
a=sum(B(p:
q-1,1))/2;
fork=p:
q-2
ifabs(sum(B(p:
=abs(sum(B(p:
k+1,1))-a);
fori=p:
q-1
B(i,j)=0;
END(i)=[char(END(i)),'
B(i,j)=1;
p=q;
end
C=B(:
j);
D=find(C==-1);
[e,f]=size(D);
ife==n
j=0;
j=j+1;
B
A
END
[u,v]=size(char(END(i)));
L(i)=v;
avlen=sum(L.*A)
H=-sum(A.*log2(A));
disp(['
ENtropyofthesource:
'
num2str(H),'
bits'
]);
5.实验结果
①实验运行结果(以[0.19,0.18,0.17,0.16,0.13,0.10,0.06,0.01]为例)
②计算:
信息熵是多少?
上述编码方法的平均码长是多少?
信息熵是2.7829;
平均码长是2.88;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 编码 matlab 实现