zengliangfa.cpp

这个是以前开始研究模式识别的关于C均值算法和固定增量法的VC++程序

#include<malloc.h>
#include<stdio.h>
#include<stdlib.h>

int **DArray2(int row,int linea)
{
   int * *m,i;
   if((m=(int * *)malloc(row*sizeof(int *)))==NULL)
   {
     printf("内存空间不够!");
		 exit(0);
   }
   for(i=0;i<row;i++)
	   if((m[i]=(int*)malloc((linea)*sizeof(int)))==NULL)
	   {
	     printf("内存空间不够!");
		 exit(0); 
	   }
     return m;
}

void DiliverDArray2(int * *m,int row)
{
	int i;
	for(i=0;i<row;i++)
		free(m[i]);
	free(m);
}


void main()
{
	int row1,row2,linea;                    
    int i,j,k;
    int * *m;                             
	int * *n;                             
    int * *f;                            
	
	
	printf("input 初始样本集1 row1,linea\n");   
	scanf("%d%d",&row1,&linea);
	printf("input 初始样本集2 row2,linea\n");
	scanf("%d%d",&row2,&linea);	

	m=DArray2(row1,linea+1);  
    n=DArray2(row2,linea+1);  
	
   printf("输入样本集1\n");               
	for(i=0;i<row1;i++) 
	{	
		m[i][0]=1;
		for(j=1;j<=linea;j++)
			scanf("%d",&m[i][j]);
	}
   printf("输入样本集2\n");
    for(i=0;i<row2;i++) 
	{	
		n[i][0]=1;
		for(j=1;j<=linea;j++)
			scanf("%d",&n[i][j]);
	}



   f=DArray2(1,linea+1);
	   for(j=0;j<=linea;j++)
	   {
		   f[0][j]=0;
	   }


   int count=0;                        
   int temp=0;                         
   int flag=1;                         
   while(flag!=0)
   {
	   flag = 0;
   for(i=0;i<row1;i++)
   {  
      count = 0;
    for(j=0;j<=linea;j++)
	{
      count=count+m[i][j]*f[0][j];     
	}
      	  if(count<=0)	   		  
		  {
			  for(k=0;k<=linea;k++)
			  f[0][k]=f[0][k]+m[i][k];
			  flag++;
		  }
   }
for(i=0;i<row2;i++)
   {
   count = 0;
    for(j=0;j<=linea;j++)
	{
      count=count+n[i][j]*f[0][j];     
	}
		  if(count>=0)			  
		  {
			  for(k=0;k<=linea;k++)
			  f[0][k]=f[0][k]-n[i][k];
			  flag++;
		  }
}
    temp++;
}
 
   printf("解向量为:\n");            
       for(i=0;i<=linea;i++)
          printf("%d ",f[0][i]);
   printf("\n");
   printf("temp=%d",temp);

   DiliverDArray2(m,row1);
   DiliverDArray2(n,row2);
   DiliverDArray2(f,1);
}