matlab费诺编码程序

%函数f1存放于
function x=f1(i,j,p,r)
global x;
x=char(x);
if(j<=i)
return;
else
>
q=0;
for t=i:j %对于区间[i,j]自上而下求累加概率值
q=p(t)+q;y(t)=q;
end
for t=i:j%把所有自上而下的累加概率值与该区间总概率值减该累加概率值之差取绝对值存在一数组
v(t)=abs(y(t)-(q-y(t)));
end
for t=i:j
|
if(v(t)==min(v)) %求该数组中最小的一个值来确定分界点位置
for k=i:t %赋值码字
x(k,r)='0';
end
for k=(t+1):j
x(k,r)='1';
end
d=t;
《
f1(i,d,p,r+1); %递归调用及相互调用
f2(d+1,j,p,r+1);
f1(d+1,j,p,r+1);
f2(i,d,p,r+1);
else
end
end
end
?
return;
%函数f2存放于
function x=f2(i,j,p,r)
global x;
x=char(x);
if(j<=i)
return;
else
{
q=0;
for t=i:j %对于区间[i,j]自上而下求累加概率值
q=p(t)+q;y(t-i+1)=q;
end
for t=1:j-(i-1)%把所有自上而下的累加概率值与该区间总概率值减该累加概率值之差取绝对值存在一数组
v(t)=abs(y(t)-(q-y(t)));
end
for t=1:j-(i-1)
…
if(v(t)==min(v)) %求该数组中最小的一个值来确定分界点位置
d=t+i-1;
for k=i:d %赋值码字
x(k,r)='0';
end
for k=(d+1):j
x(k,r)='1';
end
@
f2(d+1,j,p,r+1);%递归调用及相互调用
f1(i,d,p,r+1);
f2(i,d,p,r+1);
f1(d+1,j,p,r+1);
else
end
end
end
?
return;
主函数
clc;clear all;
fprintf('………………费诺编码程序………………\n '); fprintf('小组成员：苏、杨、张\n \n ');
fprintf('请输入信源符号的个数：');
N=input('N=');%输入信源符号的个数
s=0;l=0;H=0;
》
for i=1:N
fprintf('请输入第%d个符号的概率：',i);
p(i)=input('p=');%输入信源符号概率分布矢量,0<p(i)<1 if p(i)<=0|| p(i)>=1
error('请注意P的范围是0<P<1')
end
s=s+p(i)
H=H+(- p(i)*log2(p(i)));%计算信源信息熵
,
end
if (s~=1)
error('信源符号概率和不等1')
end
tic;
for i=1:N-1 %按概率分布大小对信源排序
for j=i+1:N
if p(i)<p(j)
、
m=p(j);p(j)=p(i);p(i)=m;
end
end
end
x=f1(1,N,p,1);
for i=1:N %计算平均码长
L(i)=length(find(x(i,:)));
l=l+p(i)*L(i);
`
end
n=H/l; %计算编码效率
fprintf('编码后所得码字:\n');
disp(x) %显示按概率降序排列的码字
fprintf('平均码长:K=\n');
disp(l)% 显示平均码长
fprintf('信息熵:H(X)=\n');
disp(H)%显示信息熵
fprintf('编码效率: η=\n');
disp(n) %显示编码效率
%fprintf('计算耗时time= %f\n',toc);。