【模式识别】探秘判别奥秘：Fisher线性判别算法的解密与实战
最编程
2024-02-25 19:14:23
...
// fisher.cpp : Defines the entry point for the console application.
//
// rmbdis.cpp : Defines the entry point for the console application.
//
#include "stdafx.h"
#include "math.h"
#include "conio.h"
#include <fstream>
#include <iomanip>
using namespace std;
#define PNUM 60
unsigned char dat[10][4][8][8][60]={
//0--样本1,1--样本1,......,8--样本9,9--样本10
//0--100,1--50,2--20,3--10
//0--A向,1--B向,2--C向,3--D向,4--新A向,5--新B向,6--新C向,7--新D向
//0--传感1,1--传感2,2--传感3,3--传感4,4--传感5,5--传感6,6--传感7,7--传感8
//
{
#include "样本\\rmb00.txt"
},
{
#include "样本\\rmb01.txt"
},
{
#include "样本\\rmb02.txt"
},
{
#include "样本\\rmb03.txt"
},
{
#include "样本\\rmb04.txt"
},
{
#include "样本\\rmb05.txt"
},
{
#include "样本\\rmb06.txt"
},
{
#include "样本\\rmb07.txt"
},
{
#include "样本\\rmb08.txt"
},
{
#include "样本\\rmb09.txt"
}
};
#define NUM 8
double Eucliden(double x[],double y[],int n)
{
double d;
d=0.0;
for (int i=0;i<n;i++) {
d+=(x[i]-y[i])*(x[i]-y[i]);
}
d=sqrt(d);
return d;
}
double Manhattan(double x[],double y[],int n)
{
double d;
d=0.0;
for (int i=0;i<n;i++) {
d+=fabs(x[i]-y[i]);
}
return d;
}
double Chebyshev(double x[],double y[],int n)
{
double d;
d=0.0;
for (int i=0;i<n;i++) {
if(fabs(x[i]-y[i])>d) d=fabs(x[i]-y[i]);
}
return d;
}
double Minkowski(double x[],double y[],int n,int m)
{
double d;
d=0.0;
for (int i=0;i<n;i++) {
d+=(double)powf((float)(x[i]-y[i]),(float)m);
}
d=(double)powf((float)d,1.0f/m);
return d;
}
double Mahalanobis(double x[],double y[],double matv1[8][8])
{
double dx,dy;
int i,j;
dx=0.0;
for (i=0;i<8;i++) {
dy=0.0;
for (j=0;j<8;j++) {
dy+=matv1[i][j]*(x[j]-y[j]);
}
dx+=dy*(x[i]-y[i]);
}
return dx;
}
void GetMatV(double V[8][8],int k)
{
int i,j,m,n1,n2,n3;
double xm[8],d,x,y;
m=4*8*PNUM;
for (i=0;i<8;i++) {
d=0;
for (n1=0;n1<4;n1++) {
for (n2=0;n2<8;n2++) {
for (n3=0;n3<PNUM;n3++) {
d+=(double)dat[k][n1][n2][i][n3];
}
}
}
d/=m;
xm[i]=d;
}
for (i=0;i<8;i++) {
for (j=0;j<8;j++) {
d=0;
for (n1=0;n1<4;n1++) {
for (n2=0;n2<8;n2++) {
for (n3=0;n3<PNUM;n3++) {
x=(double)dat[k][n1][n2][i][n3]-xm[i];
y=(double)dat[k][n1][n2][j][n3]-xm[j];
d+=x*y;
}
}
}
d/=m-1.0;
V[i][j]=d;
}
}
}
void Gauss_Jordan(double matv[8][8],double matv1[8][8])
{
int n=8;
double mat[8][16],d;
int i,j,l,k;
for (i=0;i<n;i++) {
for (j=0;j<2*n;j++) {
if (j<n)
mat[i][j]=matv[i][j];
else
mat[i][j]=0.0;
}
}
for (i=0;i<n;i++) mat[i][n+i]=1.0;
for (k=0;k<n;k++) {
d=fabs(mat[k][k]);
j=k;
for (i=k+1;i<n;i++) {//选主元
if (fabs(mat[i][k])>d) {
d=fabs(mat[i][k]);
j=i;
}
}
if (j!=k) { //交换
for (l=0;l<2*n;l++) {
d=mat[j][l];
mat[j][l]=mat[k][l];
mat[k][l]=d;
}
}
for (j=k+1;j<2*n;j++) {
mat[k][j]/=mat[k][k];
}
for (i=0;i<n;i++) {
if (i==k) continue;
for (j=k+1;j<2*n;j++) {
mat[i][j]-=mat[i][k]*mat[k][j];
}
}
}
for (i=0;i<n;i++) {
for (j=0;j<n;j++) {
matv1[i][j]=mat[i][j+n];
}
}
}
void getswj(double mats[8][8],double mj[8],unsigned char data[8][60])
{
int i,j,k;
for (i=0;i<8;i++)
{
mj[i]=0.0;
for (k=0;k<PNUM;k++)
{
mj[i]+=(double)data[i][k];
}
mj[i]/=60.0;
}
for (i=0;i<8;i++)
{
for (j=0;j<8;j++)
{
mats[i][j]=0;
for (k=0;k<PNUM;k++)
{
mats[i][j]+=(data[i][k]-mj[i])*(data[j][k]-mj[j]);
}
mats[i][j]/=59.0;
}
}
}
void get4sw(double mats[8][8],double mj[8],unsigned char data[8][8][60])
{
int i,j,k,m;
for (i=0;i<NUM;i++) {
mj[i]=0.0;
for (j=0;j<8;j++) {
for (k=0;k<PNUM;k++) {
mj[i]+=(double)data[j][i][k];
}
}
mj[i]/=8.0*PNUM;
}
for (i=0;i<NUM;i++) {
for (j=0;j<NUM;j++) {
mats[i][j]=0;
for (m=0;m<8;m++) {
for (k=0;k<PNUM;k++) {
mats[i][j]+=(data[m][i][k]-mj[i])*(data[m][j][k]-mj[j]);
}
}
mats[i][j]/=8*PNUM-1;
}
}
}
void getsb(double sb[8][8],double mj[32][8],unsigned char data[4][8][8][60])
{
int i,j,k;
double m[8];
for (i=0;i<8;i++) {
m[i]=0;
for (j=0;j<32;j++) {
for (k=0;k<60;k++) {
m[i]+=data[j/8][j%8][i][k];
}
}
m[i]/=60.0*32.0;
}
for (i=0;i<8;i++) {
for (j=0;j<8;j++) {
sb[i][j]=0;
for (k=0;k<32;k++) {
sb[i][j]+=(mj[k][i]-m[i])*(mj[k][j]-m[j]);
}
sb[i][j]/=32;
}
}
}
void getsw(double swj[32][8][8],double sw[8][8])
{
int i,j,k;
for (i=0;i<8;i++) {
for (j=0;j<8;j++) {
sw[i][j]=0;
for (k=0;k<32;k++) {
sw[i][j]+=swj[k][i][j];
}
sw[i][j]/=32.0;
}
}
}
void MatMul(double mata[8][8],double matb[8][8],double matc[8][8])
{
int i,j,k;
for (i=0;i<NUM;i++) {
for (j=0;j<NUM;j++) {
matc[i][j]=0;
for (k=0;k<NUM;k++) {
matc[i][j]+=mata[i][k]*matb[k][j];
}
}
}
}
void MatAdd(double mata[8][8],double matb[8][8],double matc[8][8])
{
int i,j;
for (i=0;i<NUM;i++) {
for (j=0;j<NUM;j++) {
matc[i][j]=mata[i][j]+matb[i][j];
}
}
}
void MatDec(double mata[8][8],double matb[8][8],double matc[8][8])
{
int i,j;
for (i=0;i<NUM;i++) {
for (j=0;j<NUM;j++) {
matc[i][j]=mata[i][j]-matb[i][j];
}
}
}
void getst(double sw[8][8],double sb[8][8],double st[8][8])
{
MatAdd(sw,sb,st);
}
double MatTrace(double mat[8][8])
{
int i;
double d=0.0;
for(i=0;i<NUM;i++) {
d+=mat[i][i];
}
return d;
}
void OutSw(ofstream outfile,double sw[NUM][NUM])
{
int i,j;
for (i=0;i<NUM;i++) {
for (j=0;j<NUM;j++) {
outfile<<setprecision(5)<<sw[i][j];
if (j<NUM-1) outfile<<",";
else outfile<<endl;
}
}
}
double MulVector(double x[NUM],double y[NUM])
{
int i;
double d;
d=0.0;
for (i=0;i<NUM;i++) {
d+=x[i]*y[i];
}
return d;
}
int main(int argc, char* argv[])
{
double sw[32][8][8];
double mj[32][8];
double sww[8][8];
double sww1[8][8];
int i,j;
// char name[20]="sw100aa.h";
/* for (i=0;i<32;i++) {
getswj(sw[i],mj[i],dat[0][i/8][i%8]);
}
MatAdd(sw[0],sw[8],sww);
Gauss_Jordan(sww,sww1);
ofstream outfile;
outfile.open("sw100ab.txt");
outfile<<"//100A m1: \n";
for (i=0;i<NUM;i++) {
outfile<<setw(5)<<setprecision(3)<<mj[0][i]<<",";
}
outfile<<endl;
outfile<<"//100b m2: \n";
for (i=0;i<NUM;i++) {
outfile<<setw(5)<<setprecision(3)<<mj[8][i]<<",";
}
outfile<<endl;
outfile<<"//100A SW1: \n";
for (i=0;i<NUM;i++) {
for (j=0;j<NUM;j++) {
outfile<<setw(5)<<setprecision(3)<<sw[0][i][j];
if (j<NUM-1) outfile<<",";
else outfile<<endl;
}
}
outfile<<"//100b SW2: \n";
for (i=0;i<NUM;i++) {
for (j=0;j<NUM;j++) {
outfile<<setw(5)<<setprecision(3)<<sw[8][i][j];
if (j<NUM-1) outfile<<",";