用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
相关文档：
如果一个变量你需要几种可能存在的值，那么就可以被定义成为枚举类型。之所以叫枚举就是说将变量或者叫对象可能存在的情况也可以说是可能的值一一例举出来。
举个例子来说明一吧，为了让大家更明白一点，比如一个铅笔盒中有一支笔，但在没有打开之前你并不知道它是什么笔，可能是铅笔也可能是钢笔，这里有两种可能，� ......
一、概述及示例代码
Delphi中包括许多已经封装好的类及控件，其中的非可视化控件库以功能方式划分可处理诸多应用需求。若使用C＋＋实现系统时对某些功能简单调用delphi中现成的库时即可。因此将delphi中的库以DLL形式封装好之后如何将方法导出可供C＋＋调用是本文记录的重点。C＋＋调用的方式有多种，在这里只讨论一种静 ......
#include <stdio.h>
#define MAXHIST 15
#define MAXWORD 11
#define IN 1
#define OUT 0
//打印输入中单词长度的水平方向直方图，
void main()
{
int c, i, nc, state;
int len; //length of each bar
int maxvalue; &n ......
在C语言中，所有传递给函数的参数都是按值传递的。
#include <iostream>
using namespace std;
void Out(int* p)
{
int j = 11;
p = &j;
*p = 12;
cout<<*p<<endl;
}
int main(int argc,char*argv[])
{
&n ......
typedef 声明，简称 typedef，为现有类型创建一个新的名字。比如人们常常使用 typedef 来编写更美观和可读的代码。所谓美观，意指 typedef 能隐藏笨拙的语法构造以及平台相关的数据类型，从而增强可移植性和以及未来的可维护性。
第一、四个用途
用途一：
定义一种类型的别名，而不只是简单的宏替换。可以用作同时声明指 ......
