標題: 統(tǒng)計信源熵C程序 [打印本頁]

---

作者: 胡頹葉子    時間: 2017-6-1 17:37
標題: 統(tǒng)計信源熵C程序

1. 統(tǒng)計信源熵程序：
   
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");
    
13. 
    
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.        }
    
23. 
    
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");
    
45. 
    
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]));
    
53. 
    
54. }
    
55. printf("\n");
    
56. printf("信息熵=%f",-H);
    
57. printf("\n");
    
58. }
    

復(fù)制代碼

---

歡迎光臨 (http://www.torrancerestoration.com/bbs/)

Powered by Discuz! X3.1