题目83（无忧id 34 方差运算题）

 

请编制函数ReadDat()实现从文件IN.DAT中读取1000个十进制整数到数组xx中；请编制函数Compute()分别计算出xx中奇数的个数odd，奇数的平均值ave1，偶数的平均值ave2以及所有奇数的方差totfc的值，最后调用函数WriteDat()把结果输出到OUT.DAT文件中。

计算方差的公式如下：

N 2

totfc=1/N ∑ (xx[i]-ave1)

i=1

设N为奇数的个数，xx[i]为奇数，ave1为奇数的平均值。

原始数据文件存放的格式是：每行存放10个数，并用逗号隔开。（每个数均大于0且小于等于2000）

部分源程序存在文件prog1.c中。

请勿改动主函数main()和输出数据函数writeDat()的内容。

#include <stdio.h>

#include <stdlib.h>

#include <string.h>

#define MAX 1000

 

int xx[MAX],odd=0,even=0;

double ave1=0.0,ave2=0.0,totfc=0.0;

void WriteDat(void);

 

int ReadDat(void)

{

int i;

FILE *fp;

 

if((fp=fopen("IN.DAT","r"))==NULL) return 1;

 

/*********编制函数ReadDat()的部分************/

for(i=0;i<MAX;i++)

{fscanf(fp,"%d,",&xx[i]);

if((i+1)%10==0)

fscanf(fp,"\n");

}

/*******************************************/

 

fclose(fp);

return 0;

}

 

void Compute(void)

{ int i,yy[MAX];

for(i=0;i<MAX;i++)

yy[i]=0;

for(i=0;i<MAX;i++)

if(xx[i]%2) { yy[odd++]=xx[i]; ave1+=xx[i];}

else { even++; ave2+=xx[i];}

if(odd==0) ave1=0;

else ave1/=odd;

if(even==0) ave2=0;

else ave2/=even;

for(i=0;i<odd;i++)

totfc+=(yy[i]-ave1)*(yy[i]-ave1)/odd;

}

 

void main()

{

int i;

 

for(i=0;i<MAX;i++)xx[i]=0;

if(ReadDat()){

printf("数据文件IN.DAT不能打开!\007\n");

return;

}

Compute();

printf("ODD=%d\nAVE1=%f\nAVE2=%f\nTOTFC=%f\n",odd,ave1,ave2,totfc);

WriteDat();

}

 

void WriteDat(void)

{

FILE *fp;

int i;

 

fp=fopen("OUT.DAT","w");

fprintf(fp,"%d\n%f\n%f\n%f\n",odd,ave1,ave2,totfc);

fclose(fp);

}

下一篇：2006年9月全国等级考试三级c语言上机题库(八十五)