- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
《信息处理与编码》结课大作业
学号:
班级:
姓名:
成绩:
香农—费诺—埃利斯码的MATLAB语言实现
1.编码原理:香农—费诺—埃利斯码不是分组码,它根据信源符号的积累概率分配码字,不是最正确码,但它的编码和译码效率都很高。步骤1、求修正累计概率。2、转换成二进制小数形式。3、求码长。4、截取码字。
2.编码如下:
clc;
clear;
A=[0.4,0.3,0.1,0.09,0.07,0.04];
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
end
fori=1:n%生成B第2列的元素
ifi=k
B(i,2)=0;
else
B(i,2)=1;
end
end
%生成第一次编码的结果
END=B(:,2);
END=sym(END);
%生成第3列及以后几列的各元素
j=3;
while(j~=0)
p=1;
while(p=n)
x=B(p,j-1);
forq=p:n
ifx==-1
break;
else
ifB(q,j-1)==x
y=1;
continue;
else
y=0;
break;
end
end
end
ify==1
q=q+1;
end
ifq==p|q-p==1
B(p,j)=-1;
else
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];
else
a=sum(B(p:q-1,1))/2;
fork=p:q-2
ifabs(sum(B(p:k,1))-a)=abs(sum(B(p:k+1,1))-a);
break;
end
end
fori=p:q-1
ifi=k
B(i,j)=0;
END(i)=[char(END(i)),0];
else
B(i,j)=1;
END(i)=[char(END(i)),1];
end
end
end
end
p=q;
end
C=B(:,j);
D=find(C==-1);
[e,f]=size(D);
ife==n
j=0;
else
j=j+1;
end
end
B
A
END
fori=1:n
[u,v]=size(char(END(i)));
L(i)=v;
end
avlen=sum(L.*A)
运行结果如图1所示:
图1运行结
文档评论(0)