移動平均法 C言語

Question

###前提・実現したいこと

初心者です。
C言語でCSVファイルを読み込み、移動平均法で値を抽出し、また別のCSVファイルに出力するという処理を行いたいのですが、データ処移動平均法を用いて値を抽出する部分が思いつかず困っています。

例
測定値 処理
2 (2+5+7)/3
5 (5+7+5)/3
7 (7+5+4)/3
5 …
4

上記のように移動平均法を用いて値を抽出するプログラムを教えて頂きたいです。

ファイルの入出力の部分は下記のようになっています。
#include <stdio.h>

int main(void){
FILE *fp;
double data[1000][2];
int i=0,num,ret;

fp=fopen("data.csv","r");
if( fp == NULL ) {
printf( "file none\n" );
return -1;
}

while( ret = fscanf(fp,"%lf,%lf",&data[i][0],&data[i][1]) != EOF){
i++;
}

num=i;
for(i=0;i<num; i++){
printf("%.1f,%.1f\n",data[i][0],data[i][1]);
}

fclose(fp);

fp=fopen("kekka.csv","w");
for(i=0;i<num; i++){
fprintf(fp,"%f,%f\n",data[i][0],data[i][1]);
}
fclose(fp);
return 0;
}

もしお力添えを頂ける方がいらっしゃいましたらどうかよろしくお願いいたします。

Comments

[ozwk]

今どこの段階ですか? やりたいことはどう分解できて、何が出来て何が出来ませんか? excel使わない理由はありますか?

[yuba]

例示されている処理は移動平均ではありません。移動平均は過去n回の値の平均ですが、例示の処理は未来の値を平均しています。もちろん、未来の値をプログラムが扱えるわけもなく、それで計算できずにいると言えます。例示の処理は正しいか確認してください。

[yp33ahgr]

ozwk様>ご指摘頂きありがとうございます。質問を編集しました。現在データの入出力の処理の部分までは分かっている状況です。課題で指定されているためexcelを使っていません。

Answer 1

移動平均は、アナログ信号の取り込み等でノイズを取り除く場合によく使います。
取りあえず、質問されている事を実現するコードを書いてみました。

#include <stdio.h>

#define	V_AVE_CNT		3
#define	V_INVALID_DT	0x7FFFFFFF
typedef	struct s_avebuf {
	int		   top;
	int		   cnt;
	int		   Total;
	int		   dt[V_AVE_CNT];
} S_AVEBUF;

void							ini_avebuf(
	S_AVEBUF *					pAveBuf
)
{
	pAveBuf->top = 0;
	pAveBuf->cnt = 0;
	pAveBuf->Total = 0.0;
}

int								add_avebuf(
	S_AVEBUF *					pAveBuf,
	int							NewDt
)
{
	int							result = V_INVALID_DT;
	int							pos;

	if (pAveBuf->cnt < V_AVE_CNT) {
		pAveBuf->Total += NewDt;
		pAveBuf->dt[pAveBuf->cnt++] = NewDt;
	}
	else {
		pos = pAveBuf->top;
		pAveBuf->Total -= pAveBuf->dt[pos];
		pAveBuf->dt[pos] = NewDt;
		pAveBuf->Total += pAveBuf->dt[pos];
		pAveBuf->top = (pAveBuf->top + 1) % V_AVE_CNT;
	}
	if (pAveBuf->cnt >= V_AVE_CNT) {
		result = pAveBuf->Total / V_AVE_CNT;
	}
	return(result);
}

int								main(
	int							argc,
	char *						argv[]
)
{
	S_AVEBUF					ave;
	int							indt[] = {2, 5, 7, 6, 4};
	int							i;
	int							r;


	ini_avebuf(&ave);
	for (i = 0; i < sizeof(indt) / sizeof(indt[0]); i++) {
		r = add_avebuf(&ave, indt[i]);
		if (r != V_INVALID_DT) {
			printf("r = %d\n", r);
		}
	}
	return(0);
}

Answer 2

移動平均法は原価算出にたまに使いますが、今回はご提示の様に3個くくりの例として作ってはみました。

#include<stdio.h>
#include<string.h>

void main(void)
{
  int c[31];
  int i = 0;
  printf( "input value & finish Ctl+d \n" );

  while( fscanf( stdin , "%d" , &c[i]) != EOF ) {
    if ( i > 1 ){
        printf("%f \n", (c[i-2]+c[i-1]+c[i])/3.00);
    }
    if ( i++ >= 30 ) return;
  }
}