用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
相关文档：
什么是回调函数？
简而言之，回调函数就是一个通过函数指针调用的函数。如果你把函数的指针（地址）作为参数传递给另一个函数，当这个指针被用为调用它所指向的函数时，我们就说这是回调函数。
为什么要使用回调函数？
因为可以把调用者与被调用者分开。调用者不关心谁是被调用者，所有它需知道的，只是存在 ......
在这个例子里，回调函数比较两个值。查找函数向比较函数传递两个指向需要进行比较的值的指针，并检查比较函数的返回值。例如：零表示相等的值，现在查找函
数就与类型无关，因为它本身并不执行实际的比较。确实，调用者必须编写必需的比较函数，但这样做是很容易的，因为调用者知道链表中所包含的值的类型 ......
简单的函数调用，通过简单的函数调用反汇编可以清楚了解如下
1.栈到底是什么，如何操纵栈的？
2.参数和临时变量是以什么形式在哪存放？
3.如何传递返回值？
举例：
#include <stdio.h>
int add(int a,int b)
{
int c=0;
c=a+b;
......
前几天搞一个项目涉及到C/C++混合编程，在Keil下编译时提示：
error: #800:this declaration may not have extern "C" linkage
error: At end of source : #67:expected a "}"
当时立马想到是C和C++编程时需要extern "C"，马上试了试，然后宣告失败。
因为双击后提示出错的地方是template<clase T>,以为是keil不� ......
C/C++招聘的一些感受和经验[转]
前段时间刚参加了n多公司的C/C++软件工程师的面试，有国企，外企，私企（moto，飞思卡尔，港湾，中国卫星XXX（这个牛））等等等等。感受感想颇多，近日终于空闲，在此表述一下。
本人基本条件：3年开发经验，2year+ Windows development experence，1 ......
