W=reshape(W',1,m*n);

s0='很好！输入正确，Huffman编码1结果如下：';

s1='Huffman编码所得码字W：';

s2='Huffman编码平均码字长度L:';

s3='Huffman编码所得码字W的方差V：';

s4='Huffman编码的编码效率q：';

disp(s0);

disp(s1),disp(B),disp(W);

disp(s2),disp(L);

disp(s3),disp(V);

disp(s4),disp(q);

（二）建立文件exercise4_8_1.m通过调用huffman.m计算教材第117页练习4.8（1）信源Y

文件exercise4_8_1.m的代码和运行结果：

clc;

P=[0.2,0.19,0.18,0.17,0.15,0.10,0.01];

Y=huffman(P);

运行结果:

（三）参照教材第232页信息率失真函数的迭代计算程序建立文件ratedf.m

文件ratedf.m的代码：

%信息率失真函数的迭代计算，选代精度取为10^(-7) %

%在信源的输入概率分布Pa和失真矩阵d已知的条件下求出信息率失真函数

%