IT实验--熵计算编程实验.docVIP

云南大学数学与与统计学院 上机实践报告 课程名称：信息论基础实验 年级：2011级 上机实践成绩： 指导教师：陆正福 姓名：凌建新 上机实践名称：熵的计算编程实验 学号：20111910101 上机实践日期： 2014.3.12 上机实践编号：No.2 组号：第七组 上机实践时间： 一、实验目的 二、内容1. 给定二维分布函数，计算联合熵、条件熵、互信息、各变量的熵。 ??? 可选择课本例题2.2.1作为程序测试用例。 2. 自行设定原始数据（如一段文本、一幅图像、一个数据表等），按照频率计算符号的分布，??? 进而计算有关的熵。 三、环境 （注意记录实验中遇到的问题） x y 1 2 3 4 　 1 1/8 1/16 1/32　 1/32　 　 2 1/16　 1/8 1/32　 1/32　 　 3 1/16　 1/16　 1/16　 1/16 　 4 1/4 0　 0　 　0 　 　 　 　 　 　 　 2.实验代码 public class Information { public static void main(String[] args) { double Hx = 0,Hy=0,Hxy=0,Hx_y=0,Hy_x=0; double Ixy=0; double [] num2=new double[4]; double [] num3=new double[4]; double [][] num = {{(double)1/8 , (double)1/16 , (double)1/32, (double)1/32},{(double)1/16, (double)1/8 , (double)1/32 , (double)1/32 },{(double)1/16 , (double)1/16 , (double)1/16, (double)1/16},{(double)1/4,(double)0.0,(double)0.0, (double)0.0 }}; for(int i=0;inum.length;i++){ for(int j=0;jnum[i].length;j++) { System.out.print(num[i][j] + ); } System.out.print(\n); } for(int i=0;inum.length;i++){ for(int j=0;jnum[i].length;j++){ num2[i]= num[j][0]+num[j][1]+ num[j][2]+num[j][3]; Hy =Hy -(num2[i] * Math.log(num2[i]) / Math.log(2)); // H = -∑Pi*log2(Pi) } System.out.println(y的熵为：H(y)= + Hy); break; } for(int i=0;i4;i++){ for(int j=0;j4;j++){ num3[i]= num[0][j] +num[1][j]+ num[2][j]+num[3][j]; Hx =Hx -(num3[i] * Math.log(num3[i]) / Math.log(2)); } System.out.println(y的熵为：H(x)= + Hx); break; } for(int i=0;i4;i++){ for(int j=0;j4;j++){ if(num[i][j]!=0){ Hxy=Hxy-(num[i][j] )* (Math.log(num[i][j]) / Math.log(2));// H = -∑Pi*log2(Pi) } } } Hy_x=Hxy-Hx; Ixy=Hy-Hy_x;