找回密碼
 立即注冊(cè)

QQ登錄

只需一步，快速開(kāi)始

搜索
查看: 2245|回復(fù): 0
打印 上一主題 下一主題
收起左側(cè)

統(tǒng)計(jì)信源熵C程序

[復(fù)制鏈接]
跳轉(zhuǎn)到指定樓層
樓主
ID:206798 發(fā)表于 2017-6-1 17:37 | 只看該作者 回帖獎(jiǎng)勵(lì) |倒序?yàn)g覽 |閱讀模式
  1. 統(tǒng)計(jì)信源熵程序：
  2. #include<math.h>
  3. #include<stdio.h>
  4. void main()
  5. {
  6. int ch[50]={0};
  7. int sum=0,n=0;
  8. int i;
  9. double p[50]={0};
  10. double H=0;
  11. char c;
  12. printf("please input an essay:\n");

  14. while((c=getchar())!='\n')
  15. {
  16. for(i=65;i<=90;i++)
  17.      {
  18.        if(i==c)
  19.        {
  20.    ch[i-65]++;
  21.    sum++;
  22.        }

  24.     }
  25.     for(i=97;i<=112;i++)
  26.      {
  27.              if(c==i)
  28.              {
  29.               ch[i-65]++;
  30.               sum++;
  31.              }
  32.      }
  33. }
  34. for(i=0;i<50;i++)
  35. {
  36.    p[ i]=(double)ch[ i]/(double)sum;
  37.    if(p[ i]!=0)
  38.    {
  39.            printf("p(%c)=%1.2f",i+65,p[ i]);
  40.            n++;
  41.    }
  42.    if(n==3)
  43.    {
  44.            printf("\n");

  46.            n=0;
  47.    }
  48. }
  49. for(i=0;i<=25;i++)
  50. {
  51. if(p[ i]!=0)
  52. H=H+p[ i]*(log(p[ i]));

  54. }
  55. printf("\n");
  56. printf("信息熵=%f",-H);
  57. printf("\n");
  58. }
復(fù)制代碼


分享到:  QQ好友和群QQ好友和群 QQ空間QQ空間 騰訊微博騰訊微博 騰訊朋友騰訊朋友
收藏收藏 分享淘帖 頂頂 踩踩
回復(fù)

使用道具 舉報(bào)

本版積分規(guī)則

手機(jī)版|小黑屋|51黑電子論壇 |51黑電子論壇6群 QQ 管理員QQ:125739409;技術(shù)交流QQ群281945664

Powered by 單片機(jī)教程網(wǎng)

快速回復(fù) 返回頂部 返回列表