fork download

copy

#include <stdio.h>
 
int main(void){
 
  double sigma[2][2] = { {1, 0}, {0, 1} }, W[2] = {1.0, 0.0};
  //学生番号を入力
  //double n0 = 2.0, n1 = 1.0, n2 = 0.0, n3 = 1.0;
  double n0 = 3.0, n1 = 8.0, n2 = 0.0, n3 = 1.0;
  double x[5][2];
  double sumx[2] = {0.0};
  int i;
 
  x[0][0] = n1;
  x[0][1] = n2 + n3;
  x[1][0] = n0 + n3;
  x[1][1] = n0 + n1;
  x[2][0] = n2 + n3;
  x[2][1] = n1;
  x[3][0] = n1 + n2;
  x[3][1] = n1 + n0;
  x[4][0] = n1 + n0;
  x[4][1] = n2 + n1 + n0;
 
  puts("x");
  for(i = 0; i < 5; i++)
    printf("%d %6.4lf %6.4lf\n", i + 1, x[i][0], x[i][1]);
 
  for(i = 0; i < 5; i++){
    sumx[0] += x[i][0];
    sumx[1] += x[i][1];
  }
 
  printf("\nμ\n%5.3lf %5.3lf\n\n", sumx[0]/5.0, sumx[1]/5.0);
 
  puts("x'");
  for(i = 0; i < 5; i++)
    printf("%d %6.4lf %6.4lf\n", i + 1, x[i][0] - sumx[0]/5.0, x[i][1] - sumx[1]/5.0);
 
  double mu_zx[5] = {0.0}, zz_t[5] = {0.0};
 
  printf("\n<z>n\n");
  for(i = 0; i < 5; i++){
    mu_zx[i] = (x[i][0] - sumx[0] / 5.0) * 0.5;
    printf("%d %6.4lf\n", i + 1, mu_zx[i]);
  }
 
  printf("\n<zz>n\n");
 
  for(i = 0; i < 5; i++){
    zz_t[i] = 0.5 + mu_zx[i] * mu_zx[i];
    printf("%d %6.4lf\n", i + 1, zz_t[i]);
  }
 
  double xx[2][2], zz, xz[2], Sigma[2][2];
 
  for(i = 0; i < 5; i++){
    xx[0][0] += (x[i][0] - sumx[0]/5.0) * (x[i][0] - sumx[0]/5.0);
    xx[0][1] += (x[i][0] - sumx[0]/5.0) * (x[i][1] - sumx[1]/5.0);
    xx[1][0] += (x[i][1] - sumx[1]/5.0) * (x[i][0] - sumx[0]/5.0);
    xx[1][1] += (x[i][1] - sumx[1]/5.0) * (x[i][1] - sumx[1]/5.0);
  }
 
  printf("\nx'x\n%6.4lf  %6.4lf\n%6.4lf  %6.4lf\n", xx[0][0],xx[0][1],xx[1][0],xx[1][1]);
 
  for(i = 0; i < 5; i++){
    zz += zz_t[i];
  }
  printf("\nzz : %6.4lf\n", zz);
 
  for(i = 0; i < 5; i++){
    xz[0] += (x[i][0] - sumx[0] / 5.0) * mu_zx[i];
    xz[1] += (x[i][1] - sumx[1] / 5.0) * mu_zx[i];
  }
 
  printf("\nx'z : %6.4lf %6.4lf\n", xz[0], xz[1]);
 
  printf("\nW = %6.4lf %6.4lf\n", xz[0] * (1.0 / zz), xz[1] * (1.0 / zz) );
 
  Sigma[0][0] = (xx[0][0] - xz[0] * xz[0] * (1.0 / zz));
  Sigma[0][1] = (xx[0][1] - xz[0] * xz[1] * (1.0 / zz));
  Sigma[1][0] = (xx[1][0] - xz[1] * xz[0] * (1.0 / zz));
  Sigma[1][1] = (xx[1][1] - xz[1] * xz[1] * (1.0 / zz));
 
  //この答えをwolframに突っ込んで対角化行列を求め、5で割るとΣ
  printf("\nΣ\n%6.4lf  %6.4lf\n%6.4lf  %6.4lf\n", Sigma[0][0],Sigma[0][1],Sigma[1][0],Sigma[1][1]);
 
  double after[2][2];
 
  after[0][0] = xz[0] * (1.0 / zz) * xz[0] * (1.0 / zz); //対角成分を加算
  after[0][1] = xz[0] * (1.0 / zz) * xz[1] * (1.0 / zz);
  after[1][0] = xz[1] * (1.0 / zz) * xz[0] * (1.0 / zz);
  after[1][1] = xz[1] * (1.0 / zz) * xz[1] * (1.0 / zz); //対角成分を加算
 
  printf("\nAftrer\n%6.4lf  %6.4lf\n%6.4lf  %6.4lf\n", after[0][0],after[0][1],after[1][0],after[1][1]);
 
  return 0;
 
}
 Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
Unshare from #datamodeling
fool Jul 10th at 11:20 PM 

prog.c: In function ‘main’:
prog.c:5:44: warning: unused variable ‘W’ [-Wunused-variable]
   double sigma[2][2] = { {1, 0}, {0, 1} }, W[2] = {1.0, 0.0};
                                            ^
prog.c:5:10: warning: unused variable ‘sigma’ [-Wunused-variable]
   double sigma[2][2] = { {1, 0}, {0, 1} }, W[2] = {1.0, 0.0};
          ^~~~~
prog.c: At top level:
prog.c:99:8: error: stray ‘\342’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
        ^
prog.c:99:9: error: stray ‘\200’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
         ^
prog.c:99:10: error: stray ‘\242’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
          ^
prog.c:99:2: error: unknown type name ‘Wrap’
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
  ^~~~
prog.c:99:16: error: invalid suffix "th" on integer constant
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                ^~~~
prog.c:99:16: error: expected ‘=’, ‘,’, ‘;’, ‘asm’ or ‘__attribute__’ before numeric constant
prog.c:99:33: error: stray ‘\342’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                 ^
prog.c:99:34: error: stray ‘\200’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                  ^
prog.c:99:35: error: stray ‘\242’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                   ^
prog.c:99:37: error: invalid suffix "kB" on integer constant
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                     ^~~
prog.c:99:51: error: stray ‘\342’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                                   ^
prog.c:99:52: error: stray ‘\200’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                                    ^
prog.c:99:53: error: stray ‘\242’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                                     ^
prog.c:99:77: error: stray ‘\302’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                                                             ^
prog.c:99:78: error: stray ‘\240’ in program
  Wrap  • Jul 10th at 11:20 PM • 3kB C snippet • Shared in datamodeling 
                                                                              ^
prog.c:100:14: error: stray ‘#’ in program
 Unshare from #datamodeling
              ^
prog.c:101:10: error: invalid suffix "th" on integer constant
 fool Jul 10th at 11:20 PM
          ^~~~

Standard output is empty