標(biāo)題: 機(jī)器學(xué)習(xí) 花卉識(shí)別分類(lèi)C++例程 [打印本頁(yè)]
作者: ffffaaaa    時(shí)間: 2020-6-1 19:21
標(biāo)題: 機(jī)器學(xué)習(xí) 花卉識(shí)別分類(lèi)C++例程


機(jī)器學(xué)習(xí) 花卉識(shí)別分類(lèi)源程序如下:
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. const int maxn=155;
  4. vector<int> v[maxn];
  5. map<string, int> mp;
  6. struct Iris
  7. {
  8.     double prop[4];
  9.     //0_Sepal_Length 1_Sepal_Width
  10.     //2_Petal_Length 3_Petal_Width
  11.     int sum;//原本分類(lèi);
  12.     int pv;//預(yù)測(cè)分類(lèi);
  13. }all[maxn];
  14. void print(Iris a)
  15. {
  16.     printf("%12.1f %11.1f %12.1f %11.1f %12d %15d ",a.prop[0],a.prop[1],a.prop[2],a.prop[3],a.sum,a.pv);
  17.     if(a.sum==a.pv) printf("\n");
  18.     else printf(" 預(yù)測(cè)錯(cuò)誤\n");
  19. }
  20. void Iris_Division(int t)
  21. //分層隨機(jī)分成t份，用于t折交叉驗(yàn)證；
  22. {
  23.     srand(time(NULL));//重置隨機(jī)數(shù)種子；
  24.     int book[maxn], cnt;
  25.     memset(book, 0, sizeof(book));
  26.     for(int i=1; i<=t; i++)
  27.     {
  28.         for(int j=1; j<=3; j++)
  29.         {
  30.             cnt=0;
  31.             while(cnt<150/(t*3))
  32.             {
  33.                 int x=rand()%50+1+50*(j-1);
  34.                 if(book[x]==1) continue;
  35.                 book[x]=1;
  36.                 v[i].push_back(x);
  37.                 cnt++;
  38.             }
  39.         }
  40.     }
  41. }
  42. struct Tree_Node
  43. {
  44.     int sum;//-1表示非葉節(jié)點(diǎn)， 0、1、2表示該節(jié)點(diǎn)劃分種類(lèi)
  45.     int basic;//依據(jù)什么劃分
  46.     double div;//分界值
  47.     struct Tree_Node* left;//小于等于分界值的兒子節(jié)點(diǎn)
  48.     struct Tree_Node* right;//大于分界值的兒子節(jié)點(diǎn)
  49. };
  50. double log2(double x) { return log(x)/log((double)(2.0)); }
  51. double Cal_Ent_D(int cnt[], int n, int sum)
  52. {
  53.     double Ent_D=0;
  54.     for(int i=0; i<n; i++)
  55.         if(cnt[i]!=0)
  56.         {
  57.             double x;
  58.             x=cnt[i]*1.0/sum;
  59.             Ent_D=Ent_D-x*log2(x);
  60.         }
  61.     return Ent_D;
  62. }
  63. void Cal_Max_Gain(set<int>train, int* basic, double* div)
  64. //計(jì)算最大的信息增益，選擇劃分依據(jù);
  65. {
  66.     int cnt[3]={0};//3種分類(lèi)結(jié)果的數(shù)量；
  67.     int sum=train.size();
  68.     //cout<<sum<<endl;
  69.     for(auto x : train) cnt[all[x].sum]++;

  71.     double Ent_D=Cal_Ent_D(cnt, 3, sum);

  73.     vector<double> t;
  74.     double maxx=-1e5;
  75.     for(int i=0; i<4; i++)
  76.     {
  77.         t.clear();
  78.         for(auto x : train) t.push_back(all[x].prop[i]);
  79.         sort(t.begin(), t.end());//連續(xù)值排序，用于二分離散化
  80.         for(int j=0; j<(int)(t.size()-1); j++)
  81.         {
  82.             double y=(t[j]+t[j+1])/2.0;
  83.             double ans=Ent_D;
  84.             memset(cnt, 0, sizeof(cnt));
  85.             sum=0;
  86.             for(auto x : train) if(all[x].prop[i]<=y) cnt[all[x].sum]++, sum++;
  87.             if(sum!=0) ans=ans-(double)sum/(double)(t.size())*Cal_Ent_D(cnt, 3, sum);

  89.             memset(cnt, 0, sizeof(cnt));
  90.             sum=0;
  91.             for(auto x : train) if(all[x].prop[i]>y) cnt[all[x].sum]++, sum++;
  92.             if(sum!=0) ans=ans-(double)sum/(double)(t.size())*Cal_Ent_D(cnt, 3, sum);

  94.             if(ans>maxx)
  95.             {
  96.                 maxx=ans;
  97.                 *basic=i;
  98.                 *div=y;
  99.             }
  100.         }
  101.     }
  102. }
  103. void Build_Tree(Tree_Node** T, set<int>train)//建立決策樹(shù)
  104. {
  105.     *T=(Tree_Node*)malloc(sizeof(Tree_Node));
  106.     ((*T)->sum)=-1; ((*T)->basic)=0; ((*T)->div)=0;
  107.     set<int> s;
  108.     for(auto x : train) s.insert(all[x].sum);
  109.     if(s.size()==1) { ((*T)->sum)=(*s.begin()); s.clear(); return; }

  111.     Cal_Max_Gain(train, &((*T)->basic), &((*T)->div));

  113.     set<int> left_train, right_train;
  114.     for(auto x : train)
  115.     {
  116.         if(all[x].prop[(*T)->basic] <= ((*T)->div)) left_train.insert(x);
  117.         else right_train.insert(x);
  118.     }
  119.     Build_Tree(&(*T)->left, left_train);
  120.     Build_Tree(&(*T)->right, right_train);
  121. }

  123. int Test_Data(Tree_Node* T, int x)
  124. {
  125.     if((T->sum)!=-1) return T->sum;
  126.     if(all[x].prop[T->basic]<=(T->div))
  127.         return Test_Data(T->left, x);
  128.     else return Test_Data(T->right, x);
  129. }
  130. int main()
  131. {
  132.     freopen("iris.txt","r",stdin);
  133.     int cnt=0;
  134.     for(int i=1; i<=150; i++)
  135.     {
  136.         string s;
  137.         scanf("%lf,%lf,%lf,%lf,",&all[i].prop[0],&all[i].prop[1],&all[i].prop[2],&all[i].prop[3]);
  138.         cin>>s;
  139.         if(!mp.count(s)) { mp[s]=cnt; cnt++;}
  140.         all[i].sum=mp[s];
  141.     }
  142.     int t=10;
  143.     Iris_Division(t);//分層隨機(jī)分成t份，用于t折交叉驗(yàn)證；
  144.     set<int> train, test;
  145.     double Ans=0.0;
  146.     for(int i=1; i<=10; i++)
  147.     {
  148.         train.clear();  test.clear();
  149.         for(int j=1; j<=t; j++)
  150.         {
  151.             if(j==i) for(int k=0; k<(int)v[j].size(); k++) test.insert(v[j][k]);
  152.             else for(int k=0; k<(int)v[j].size(); k++) train.insert(v[j][k]);
  153.         }
  154.         Tree_Node* tree;
  155.         Build_Tree(&tree, train);
  156.         printf("Sepal_Length Sepal_Width Petal_Length Petal_Width Actual_Value Predicted_Value\n");
  157.         int Cor_Num=0, All_Num=0;
  158.         for(auto x: test)
  159.         {
  160.             all[x].pv=Test_Data(tree, x);
  161.             if(all[x].sum==all[x].pv) Cor_Num++;
  162.             All_Num++;
  163.             print(all[x]);
  164.         }
  165.         printf("第%d次 預(yù)測(cè)正確率：%.3f\n\n", i, (double)Cor_Num/(double)All_Num);
  166.         Ans+=(double)Cor_Num/(double)All_Num;
  167.     }
  168.     printf("平均預(yù)測(cè)正確率：%.3f\n\n", Ans/(double)t);
  169.     return 0;
  170. }
復(fù)制代碼

所有資料51hei提供下載:
機(jī)器學(xué)習(xí)程序（c++）.7z (251.07 KB, 下載次數(shù): 17)


作者: tieq1952    時(shí)間: 2020-6-7 08:18
謝謝分享!!!



歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1