/////test_curve.c
/////「2つのガウシアン+ノイズ」のデータ作成
/////2013/12/25 by Takashi Kojima

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

#define N 30  // N + (N-1) = 2N-1個のデータが出来る。

int main(){

  double d =0.5;  //2つのガウシアンはx=0とx=dに中心がある
  double sigma = 0.1;  //ノイズの大きさ
  int    i,j;
  double F[N];
  double G[N];
  double s[N];
  double t[N];
  int    num=N;
  double x;
  double kankaku=0.1; //データ取得の間隔

  FILE *fp;
  char *filename = "s0.1_d0.5.dat";
  

  //観測データの作成
  for (i=1; i<num ; i++){
    t[i]=-i*kankaku;
    x=t[i];
    G[i]=exp(-x*x)+exp(-(x-d)*(x-d)) + (rand()%10)*0.1*sigma;
  }                               //マイナス側

  for (i=0; i<num ; i++){
    s[i]=i*kankaku;
    x=s[i];
    F[i]=exp(-x*x)+exp(-(x-d)*(x-d)) + (rand()%10)*0.1*sigma;
  }                               //プラス側

  //観測データの書き出し
  fp = fopen(filename, "w");
  if( fp==NULL ){
    printf("Can't open file %s\n", filename);
    exit(1);
  }

  for (i=1; i<num ; i++){
    j=num-i;
    fprintf(fp, "%lf   %lf\n" ,t[j], G[j]);//マイナス方向から
  }

  for (i=0; i<num ; i++){
    fprintf(fp, "%lf   %lf\n" ,s[i], F[i]);//プラス方向へ
  }

  printf("Finished saving.\n");
  fclose(fp);
  
  return 0;
}