pba=[Pba(:,:,:)];

RS=[RS R(n)];

DS=[DS D(n)];

m=m+1;%外层循环次数加1

end %外层循环结束

%打印输出结果

s0='很好！输入正确，选代结果如下：';

s1='最佳转移概率分布Pba：';

s2='最小信息率Rmin：';

s3='最大失真度Dmax：';

s4='最大拉氏乘于Smax：';

disp(s0);

disp(s1),disp(Pba);

disp(s2),disp(Rmin);

disp(s3),disp(Dmax);

disp(s4),disp(Smax);

%画出信息率失真函数R(D)

plot(DS,RS)

xlabel('允许的失真度D')

ylabel('信息率失真函数R(D)')

title('信息率失真函数R(D)的曲线图')

（四）建立文件exercise6_2.m通过调用ratedf.m计算教材第161页练习6.2（%注意拉式因子s的取值范围为负无穷大到零（小于零），可从一个绝对值充分大的负数取起，如s=-100。） 文件exercise6_2.m的代码和运行结果：

clc;

Pa=[1/2,1/2];

d=[0,1,1/4;

1,0,1/4];

S=-1000000;

ratedf(Pa,d,S);

运行结果:

通过过本次试验我掌握了利用matlab语言编程实现霍夫曼编码、计算率失真函数，这次试验还让我对书本上的本部分知识有了跟深刻的了解，更加熟悉了MatLab软件的使用。通过老师的提问和讲解，还让我明白了很多关于霍夫曼编码的实质性，这次试验使我懂得了不少知识点。