C语言方差运算

请编制函数 ReadDat()实现从文件 IN.DAT 中读取 1000 个十进制整数到数组 xx 中：

请编制函数 Compute()分别计算出 xx 中偶数的个数 even，奇数的平均值 ave1，偶数的平均值 ave2 以及方差 totfc 的值，后调用函数 WriteDat()把结果输出到 OUT.DAT 文件中。

计算方差的公式如下：

N_____2

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

i=1 设 N 为偶数的个数，xx[i]为偶数，ave2 为偶数的平均值。

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

#include

#include

#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

{ 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

yy[i]=0;

for(i=0;

i

if(xx[i]%2==0) { yy[even++]=xx[i]; ave2+=xx[i];}

else { odd++; ave1+=xx[i];}

if(odd==0) ave1=0;

else ave1/=odd;

if(even==0) ave2=0;

else ave2/=even;

for(i=0;i

totfc+=(yy[i]-ave2)*(yy[i]-ave2)/even;

}

void main()

{

int i;

for(i=0;i

if(ReadDat()){

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

return;

}

Compute();

printf(“OVEN=%d\nAVE1=%f\nAVER2=%f\nTOTFC=%f\n”,even,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",even,ave1,ave2,totfc);

fclose(fp);

}