数据处理（蛋白质序列）

#include "stdafx.h"

int _tmain(int argc, _TCHAR* argv[])
{
 int i,j,h;
 char str[5000];
 FILE *readfile=fopen("F:\\风中追风\\蛋白质数据.txt","r");
 FILE *writefile=fopen("F:\\风中追风\\shujujieshou.txt","w+");
 fprintf(writefile,"A,C,D,E,F,G,H,I,K,L,M,N,P,Q,R,S,T,V,W,Y\n");
 while(!feof(readfile))
 {
  double fre[20];
  for(i=0;i<20;i++)
   fre[i]=0;
  i=0;
  fgets(str,5000,readfile);
  fgets(str,5000,readfile);
  while(str[i]!='\n')
  {
  if(str[i]=='A')
  fre[0]++;
  if(str[i]=='C')
   fre[1]++;
  if(str[i]=='D')
   fre[2]++;
  if(str[i]=='E')
   fre[3]++;
  if(str[i]=='F')
   fre[4]++;
  if(str[i]=='G')
   fre[5]++;
  if(str[i]=='H')
   fre[6]++;
  if(str[i]=='I')
   fre[7]++;
  if(str[i]=='K')
   fre[8]++;
  if(str[i]=='L')
   fre[9]++;
  if(str[i]=='M')
   fre[10]++;
  if(str[i]=='N')
   fre[11]++;
  if(str[i]=='P')
   fre[12]++;
  if(str[i]=='Q')
   fre[13]++;
  if(str[i]=='R')
   fre[14]++;
  if(str[i]=='S')
   fre[15]++;
  if(str[i]=='T')
   fre[16]++;
  if(str[i]=='V')
   fre[17]++;
  if(str[i]=='W')
   fre[18]++;
  if(str[i]=='Y')
   fre[19]++;
  i++;
  }
  for(j=0;j<20;j++)
  {
   fre[j]=fre[j]/(i-1);
   fprintf(writefile,"%.3lf,",fre[j]);
  }
  fprintf(writefile,"\n");
 }

 return 0;
}