ÓÃc±àдÇóÐÅÏ¢ÂÛÖи÷ÖÖìغÍÐÅÏ¢Á¿µÄ³ÌÐò

#include<stdio.h>
#include<math.h>
#define u 20
int i,j,n,m;
float Ix[u],I[u][u],H_X,H_Y,H_XY,I_XY,H_XpY,Pypx[u][u],Px[u],H_YpX,Py[u],Pxpy[u][u],Pxy[u][u];
/*¶¨ÒåµÄIx[]=I(x[i])×ÔÐÅÏ¢£»I[][]=I(x[i];y[j])»¥ÐÅÏ¢£»H_X=H(X)ƽ¾ù×ÔÐÅÏ¢£»H_XY=H(XY)ÁªºÏìØ;
I_XY=I(X;Y)ƽ¾ù»¥ÐÅÏ¢£»H_XpY=H(X|Y)¡¢H_YpX=H(Y|X)Ìõ¼þìØ;Pypx[i][j]=P(y[j]|x[i])Ìõ¼þ¸ÅÂÊ;
Px[i]=P(x[i])·¢x[i]µÄ¸ÅÂÊ£»H_XpY=H(Y/X)Ìõ¼þìØ;Py[j]=P(y[j])ÊÕµ½y[j]µÄ¸ÅÂÊ£»Pxpy[i][j]=P(x[i]/y[j])Ìõ¼þ¸ÅÂÊ£»
Pxy[i][j]=P(x[i]y[j])ÁªºÏ¸ÅÂÊ*/
/*¶¨ÒåÒÔ2Ϊµ×µÄ¶ÔÊýº¯Êý*/
float log2(float x)
{
 float z;
 z=(float)(log(x)/log(2));
 return z;
}
/*Çó×ÔÐÅÏ¢µÄº¯Êý*/
float self_inf(float x)
{
  float z;
  z=-log2(x);
  return z;
}
/*Çó»¥ÐÅÏ¢µÄº¯Êý*/
float mut_inf(int x,int y)
{
 float z;
 z=log2(Pxpy[x][y]/Px[x]);
 return z;
}
/*Çóƽ¾ù×ÔÐÅÏ¢µÄº¯Êý*/
float entropy(float *x,int n)
{
  float z=0;
  for(i=1;i<=n;i++)
  {
   z+=(*(x+i))*log2(1/(*(x+i)));
  }
  return z;
}
/*ÇóÁªºÏìصĺ¯Êý*/
float joint_entropy(float (*p)[u])
{
 float z=0;
 for(i=1;i<=n;i++)
  for(j=1;j<=m;j++)
   {
    z+=*(p[i]+j)*log2(1/(*(p[i]+j)));
   }
 return z; 
}
main()
{
 float s=0;
 printf("\nplease in put the dimension of 'X' and 'Y'\n");
 scanf("%d %d",&n,&m);
 printf("\nThe dimension of X is n=%d\nThe dimension of Y is m=%d\nPlease input the condition probability:P(y[j]/x[i]),",n,m);
 printf("(after you input one number please click the 'enter')\n");
/*Ìõ¼þ¸ÅÂÊP(y[j]/x[i])¸³Öµ*/
 for(i=1;i<=n;i++)
  {for(j=1;j<=m;j++)
   {printf("P(y[%d]/x[%d])=",j,i);
          scanf("%f",&Pypx[i][j]);
         }
  }
 printf("please input P(x[i]):(after you input one number please click the 'e