资源描述：
《信息论作业全》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、信息论作业信号与信息处理专业姚艳红(200511547)要求：处理图像kma.bmp,并画出各图像的灰度层次的分布图，并对256位和16位信源数据进行霍夫曼编码1.程序如下：clear;%提取原图像的数据x=imread('lena.bmpr);%一..显不原图像figure;subplot(2,2,1)imshow(x);xl=double(x);%计算图像的爛y=zeros(1,256);forp=1:256*256;y(x1(p)+1)=y(x1(p)+1)+1;endyl=y;y=y/(256*256);h=0;forq=1
2、:256ify(q)>0h=h-y(q)*log2(y(q));endend%二...用四位来表示图像，即将后四位屏蔽掉，(变为零)y2=zeros(1,256);forp=1:256;forq=1:256;x2(p,q)=x1(p,q)-rem(x1(p,q),16);endendforp=1:256*256;y2(x2(p)+1)=y2(x2(p)+1)+1;%分布endsubplot(2,2,2);imshow(uinl8(x2));%三….也是四位表示y3=zeros(1,256);forp=1:256;forq=1:256
3、;x3(p,q)二x2(p,q)/16;endendforp=l:256*256;y3(x3(p)+1)=y3(x3(p)+1)+1;%分布endsubplot(2,2,3)imshow(uint8(x3));%四….四位表示y4=zeros(1,256);forp=l:256;forq=l:256;x4(p,q)=255・x3(p,q);endendforp=1:256*256;y4(x4(p)+1)=y4(x4(p)+1)+1;%分布endsubplot(2,2,4)imshow(uint8(x4));%画出四个图象屮灰度层次的
4、分布图figure;subplot(2,2,l)stemCyl//);subplot(2,2,2)stem(y2,T);subplot(2,2,3)stem(y3,7);subplot(2,2,4)slcm(y4,'.');%霍夫曼编码%对256个灰度层次的信源数据进行编码mapl=huffman(double(yl))%对16个灰度层次的信源数据进行编码y5=zeros(l,16);j"；fori=1:256ify2(i)>0;y5(j)二y2(i);j=j+l；endendmap2=huffman(double(y5))1.运行
5、结果：(1)计算原图像的嫡为:h=7.3482(2)所显示的图像：(3)各图像中各个灰度层次的概率分布图（4）16位（后四位被屏蔽掉的数据）霍夫曼编码:0uooioooor1•looior2*000*3•mr4•oo10*5‘110’6•ioi*7'OlO*8•oir9•ooir10nooo*11•mo112•looir13'1001001'14uooiooor151100100000*（5）256位霍夫曼编码:0111110100011001100000000000000000000000000000000000000000000
6、000000,1*1111010001100110000000000000000000000000000000000000000000000000r2*111101000110011000000000000000000000000000000000000000000000000r3UlllOlOOOllOOll00000000000000000000000000000000000000000000000r4'1111010001100110000000000000000000000000000000000000000000000r5
7、UlllOlOOOllOOll000000000000000000000000000000000000000000000r6U1110100011001100000000000000000000000000000000000000000000r7*1111010001100110000000000000000000000000000000000000000000r811111010001100110000000000000000000000000000000000000000001*9*11110100011001100000000
8、000000000000000000000000000000000r1041110100011001100000()0()0()0()0()0()0()0()0()0()0000000000000000r11*111101000110
