c语言性别体重身高的程序,1-1 C语言手撕高斯-朴素贝叶斯 - 通过身高和体重推测性别（机器学习）...

目录

目录

目录

项目介绍

项目1-1的准确度

项目可行性：原理简介

高斯分布(正态分布)

朴素贝叶斯

代码流程

获取数据(数据清洗)

代码实现

头文件 allHead.h

源文件 handleData.cpp

源文件及主函数 naiveBayes.cpp

项目介绍

1-1 通过身高、体重推测性别

1-2 通过身高、体重、肺活量推测性别

1-3 MPI优化

1-4 OpenMP优化

项目1-1的准确度

本项目用了150行C语言代码实现高斯-朴素贝叶斯，虽受限于数据集的特征值少，但预测结果准确度仍有85.9%。

项目可行性：原理简介

高斯-朴素贝叶斯即高斯分布(正态分布)+朴素贝叶斯。

高斯分布和朴素贝叶斯都是高中知识这里就不科普了，具体百度就好。

高斯分布(正态分布)

通过以往的科学研究发现，人类的身高和体重基本满足高斯分布的。

找了许多资料发现，其实人类的身高在几个区间特别集中，图像是有几个峰值的高斯分布曲线。

我们可以通过数据集来获取高斯分布所需要的：平均值(u)、标准差(p)。

//dataSet是数据集

double u=0; //平均值

double p=0; //标准差

//求平均值

for(i=0;i

if(dataSet[i][0]==sex){

sum=sum+1;

u=u+dataSet[i][column];

//printf("%lf\n",dataSet[i][column]);

}

}

u=u/sum;

//求标准差

double numSum=0;

for(i=0;i

if(dataSet[i][0]==sex){

numSum=numSum+pow( dataSet[i][column]-u , 2 );

}

}

p=sqrt(numSum/dataLen);

只要有了平均值和标准差，我们就能通过高斯分布的概率密度函数来获得以下这四个概率:

P(身高=160|性别=女性)、P(体重=50|性别=女性)

P(身高=160|性别=男性)、P(体重=50|性别=男性)

//正态分布 x:随机变量的值 u:样本平均值 p:标准差 y:概率

double gaussianDistribution(double x,double u,double p)

{

double y;

p=pow(p,2);

y = (1 / (2*PI*p)) * exp( -pow((x-u),2) / (2*p) );

return y;

}

得到上述四个概率之后，就可以进行朴素贝叶斯分类了。

朴素贝叶斯

朴素贝叶斯有多朴素呢，基本上就是比较P1和P2 的大小，哪个概率大，结果就是哪个。

P1 = P(身高=160|性别=女性) * P(体重=50|性别=女性)

P2 = P(身高=160|性别=男性) * P(体重=50|性别=男性)

代码流程

获取数据(数据清洗)

读取数据

计算平均值(u)、标准差(p)

通过高斯分布的概率密度函数求出P(特征=xxx|性别=xx)的概率

求出P1、P2

通过P1、P2的大小判断性别

获取数据(数据清洗)

数据集可以在https://www.kaggle.com/找到，搜索关键词height等即可……

本次代码的数据集(韩国大学生的体测数据)如下:

第0列：性别(0表示未知、1代表男性、2代表女性)

第1列：身高(单位:cm)

第2列：体重(单位:kg)

P.S. 不分享该数据集

代码实现

头文件 allHead.h

//allHead.h

#ifndef _STDIO_H_

#define _STDIO_H_

#include

#include

#include

#include

#define PI 3.1415926535898

//单条数据的长度

#define MAX_LINE 100

//数据集的长度

#define DATA_LEN 15300

//特征值的个数

#define EIGEN_NUM 3

/*声明函数*/

//获取数据集

void getData(const char *fileLocation);

//计算概率p(特征列column = x | 性别)

double getProbability(double x,int column,int sex);

//正态分布 x:随机变量的值 u:样本平均值 p:标准差

double gaussianDistribution(double x,double u,double p);

//返回性别字符结果

char* sexResult(float height,float weight);

//返回性别ID结果

int sexIDResult(float height,float weight);

//准确度判断

float precision();

#endif

源文件 handleData.cpp

//handleData.cpp

#include"allHead.h"

/*

basicData.csv:性别、身高、体重 [i][0,1,2]

addVitalCapacityData.csv:性别、身高、体重、肺活量 [i][0,1,2,3]

*/

//性别：0表示未知，1表示男，2表示女

float dataSet[DATA_LEN][EIGEN_NUM];//数据集

int dataLen;//从1开始计算的数据集的长度,数学计算用的长度

char *basicInfo[] = {"性别","身高","体重"};

char *addVitalCapacityInfo[] = {"性别","身高","体重","肺活量"};

//读取数据

void getData(const char *fileLocation)

{

char buf[MAX_LINE];//缓冲区

FILE *fp;//文件指针s

int len;//行字符个数

//读取文件

if((fp = fopen(fileLocation,"r")) == NULL)

{

perror("fail to read");

exit (1) ;

}

//逐行读取及写入数组

char *token;

const char s[2] = ",";

int i=0;

int j=0;

while(fgets(buf,MAX_LINE,fp) != NULL && i< DATA_LEN-1)

{

len = strlen(buf);

buf[len-1] = '\0'; //删去换行符

//分割字符串

token = strtok(buf, s);

//继续分割字符串

j = 0;

while( token != NULL ) {

dataSet[i][j]=atof(token);

//printf("dataSet[%d][%d] = %f\n", i,j, dataSet[i][j] );

//printf( "%f\n", dataRuslt[i][j]);

token = strtok(NULL, s);

j = j+1;

}

i = i + 1;

}

dataLen=i;

}

//计算概率p(特征列column = x | 性别)

double getProbability(double x,int column,int sex)

{

double u=0; //平均值

double p=0; //标准差

double sum=0;

int i=0;

//判断性别

char *gender;

if(sex==1){gender="男性";}

if(sex==2){gender="女性";}

//求平均值

for(i=0;i

if(dataSet[i][0]==sex){

sum=sum+1;

u=u+dataSet[i][column];

//printf("%lf\n",dataSet[i][column]);

}

}

u=u/sum;

//求标准差

double numSum=0;

for(i=0;i

if(dataSet[i][0]==sex){

numSum=numSum+pow( dataSet[i][column]-u , 2 );

}

}

p=sqrt(numSum/dataLen);

//printf("%s%s平均值=%.16lf\n标准差=%.16lf\n",gender,basicInfo[column],u,p);

double Probability;//计算出的概率

Probability = gaussianDistribution(x,u,p);

//printf("p(%s=%lf|性别=%s)=%.16lf\n",basicInfo[column],x,gender,Probability);

return Probability;

}

//正态分布 x:随机变量的值 u:样本平均值 p:标准差 y:概率

double gaussianDistribution(double x,double u,double p)

{

double y;

p=pow(p,2);

y = (1 / (2*PI*p)) * exp( -pow((x-u),2) / (2*p) );

return y;

}

//返回性别字符结果

char* sexResult(float height,float weight)

{

double maleP;//男性概率

double femaleP;//女性概率

double a=0.5; //男女比例各50%

maleP = a * getProbability(height,1,1) * getProbability(weight,2,1);

//printf("\n");

femaleP = a * getProbability(height,1,2) * getProbability(weight,2,2);

//printf("\n");

if(maleP > femaleP){return "男性";}

if(maleP < femaleP){return "女性";}

if(maleP == femaleP){return "未知";}

}

//返回性别ID结果

int sexIDResult(float height,float weight)

{

double maleP;//男性概率

double femaleP;//女性概率

double a=0.5; //男女比例各50%

maleP = a * getProbability(height,1,1) * getProbability(weight,2,1);

//printf("\n");

femaleP = a * getProbability(height,1,2) * getProbability(weight,2,2);

//printf("\n");

if(maleP > femaleP){return 1;}

if(maleP < femaleP){return 2;}

if(maleP == femaleP){return 0;}

}

//准确度判断

float precision()

{

int i;

float preSexID;

float right=0;

float error=0;

for(i=0;i

preSexID=sexIDResult(dataSet[i][1],dataSet[i][2]);

//printf("预测ID:%f 实际ID:%f \n",preSexID,dataSet[i][0]);

if(dataSet[i][0]==preSexID){right=right+1;}

else{error=error+1;}

}

printf("Right:%f\nError:%f\n",right,error);

return right/(error+right);

}

源文件(主函数) naiveBayes.cpp

//naiveBayes.cpp

#include"allHead.h"

extern float dataSet[DATA_LEN][EIGEN_NUM];

extern int dataLen;

/*

性别：0表示未知，1表示男，2表示女

basicData.csv:性别、身高、体重 [i][0,1,2]

*/

void main(){

//给全局变量dataSet赋值

getData("addVitalCapacityData.csv");

/*

//抽查数据集

int i = dataLen-1;

int j = 2;

printf("dataSet[%d][%d] = %f\n", i,j, dataSet[i][j] );

*/

/*

//暴力循环获得准确度

printf("Accuracy：%f\n",precision());

*/

//预测性别

float Height=178;

float Weight=60;

printf("身高:%.2f,体重:%.2f\n根据高斯-朴素贝叶斯算法判断\n此人的性别可能为:%s\n\n",Height,Weight,sexResult(Height,Weight));

}