Author |
Topic: Re: My hobby (second season) (Read 1192 times) |
|
hanzack
Forum Senior Member
lalala...I'm singing a song of arimaa.
Gender: 
Posts: 35
|
 |
Re: My hobby (second season)
« on: Feb 1st, 2011, 12:43am » |
 Quote  Modify
|
I made it! Thanks arimaa babies!
#include <stdio.h>
#include <math.h>
#define N 5
double a[N][N], b[N], x[N];
#define M 250
#define OMEGA (M_PI * 2)
double y[M], c[M][N];
void read_data(char *file)
{
FILE *in;
int k;
double t, w,p[M];
in = fopen(file, "r");
if (in == NULL) { fprintf(stderr, "Read error!\n"); }
for (k = 0; k < M; k++) {
fscanf(in,"%lf %lf",&p[k], &w);
y[k]=w;
}
fclose(in);
}
void print_matrix_a(void)
{
int i, j;
printf("--- matrix a ---\n");
for (i = 0; i < N; i++) {
for (j = 0; j < N; j++) {
printf("%11.6f", a[i][j]);
}
printf("\n");
}
}
void print_vector_b(void)
{
int i;
printf("--- vector b ---\n");
for (i = 0; i < N; i++) {
printf("%10.6f ", b[i]);
}
printf("\n");
}
void print_vector_x(void)
{
int i;
printf("--- vector x ---\n");
for (i = 0; i < N; i++) {
printf("%10.6f ", x[i]);
}
printf("\n");
}
void solve(void)
{
int i, j, k;
double m, sum;
for (k = 0; k < N; k++) {
for (i = k + 1; i < N; i++) {
m = a[i][k] / a[k][k];
for(j=k+1;j<N;j++){
a[i][j]=a[i][j]-m*a[k][j];
}b[i]=b[i]-m*b[k];
}
}
for (k = N - 1; k >= 0; k--) {
sum=0;
for(j=k+1;j<N;j++){
sum=sum+a[k][j]*x[j];
}
x[k] = (b[k] - sum) / a[k][k];
}
}
void make_matrix(void)
{int i,j,k;
double t[M];
a[0][0]=0,b[0]=0;
for(k=0;k<M;k++){
t[k]=0.02*k;
c[k][0]=1;
c[k][1]=cos(OMEGA*t[k]);
c[k][2]=sin(OMEGA*t[k]);
c[k][3]=cos(2*OMEGA*t[k]);
c[k][4]=sin(2*OMEGA*t[k]);
}
for(i=0;i<N;i++){
for(k=0;k<M;k++){
b[i] +=y[k]*c[k][i];
}
}
for(i=0;i<N;i++){
for(j=0;j<N;j++){
for(k=0;k<M;k++){
a[i][j] +=c[k][j]*c[k][i];
}}}
}
void print_rmse(void)
{int j,k;
double suma,sumb,sumc,sumd, RMSE;
suma=0,sumb=0,sumc=0;
for(k=0;k<M;k++){
for(j=0;j<N;j++){
suma +=(c[k][j]*x[j])*(c[k][j]*x[j]);
sumb +=2*y[k]*(c[k][j]*x[j]);
}
}
for(k=0;k<M;k++){
sumc +=y[k]*y[k];
}
sumd = sumc-sumb+suma;
RMSE=sqrt(sumd/M);
printf("RMSE=%10.6f\n ", RMSE);
}
int main(void)
{
read_data("data.txt");
make_matrix();
print_matrix_a();
print_vector_b();
solve();
print_vector_x();
print_rmse();
}
|
|
 IP Logged |
|
|
|
hanzack
Forum Senior Member
lalala...I'm singing a song of arimaa.
Gender:
Posts: 35
|
|
Re: My hobby (second season)
« Reply #1 on: Feb 1st, 2011, 12:45am » |
Quote Modify
|
this is the "data.txt"
0.000000 0.148949
0.020000 0.872669
0.040000 0.787320
0.060000 0.542697
0.080000 1.607640
0.100000 0.792586
0.120000 1.213314
0.140000 1.258806
0.160000 1.096778
0.180000 1.413962
0.200000 0.701885
0.220000 0.692180
0.240000 0.467869
0.260000 0.383277
0.280000 0.194057
0.300000 -0.349572
0.320000 0.103754
0.340000 -0.475095
0.360000 -0.762549
0.380000 -1.263494
0.400000 -1.032502
0.420000 -1.409862
0.440000 -1.510348
0.460000 -1.438014
0.480000 -1.941349
0.500000 -1.753632
0.520000 -0.915306
0.540000 -0.872972
0.560000 -0.993482
0.580000 -0.275689
0.600000 -0.731151
0.620000 -0.253896
0.640000 0.034278
0.660000 0.302336
0.680000 0.320775
0.700000 0.453335
0.720000 0.696979
0.740000 0.955246
0.760000 0.996679
0.780000 0.999264
0.800000 0.563160
0.820000 0.954200
0.840000 0.083208
0.860000 0.629852
0.880000 0.206202
0.900000 0.466498
0.920000 -0.175060
0.940000 0.226739
0.960000 0.277326
0.980000 0.694517
1.000000 0.899604
1.020000 0.426493
1.040000 0.638867
1.060000 0.833356
1.080000 1.434481
1.100000 1.057000
1.120000 1.268745
1.140000 1.462938
1.160000 0.820981
1.180000 1.560411
1.200000 1.094464
1.220000 0.613029
1.240000 0.459632
1.260000 0.944514
1.280000 0.443002
1.300000 0.082302
1.320000 -0.454459
1.340000 -0.790384
1.360000 -0.894862
1.380000 -1.321015
1.400000 -1.523030
1.420000 -1.882842
1.440000 -1.703207
1.460000 -1.544685
1.480000 -1.127327
1.500000 -1.159030
1.520000 -1.561615
1.540000 -0.653725
1.560000 -1.207297
1.580000 -0.543838
1.600000 -0.264824
1.620000 0.282943
1.640000 0.319953
1.660000 0.620210
1.680000 0.148273
1.700000 0.565851
1.720000 0.279266
1.740000 0.838175
1.760000 0.313327
1.780000 0.305755
1.800000 0.592538
1.820000 0.663427
1.840000 0.310548
1.860000 0.650993
1.880000 0.476666
1.900000 -0.057217
1.920000 0.277954
1.940000 -0.061010
1.960000 0.438323
1.980000 0.015218
2.000000 0.554333
2.020000 1.096962
2.040000 0.486589
2.060000 1.295227
2.080000 0.998280
2.100000 1.718743
2.120000 1.325217
2.140000 1.380426
2.160000 1.701971
2.180000 1.403068
2.200000 0.743805
2.220000 0.960345
2.240000 0.839129
2.260000 0.879530
2.280000 0.108192
2.300000 0.289297
2.320000 -0.406927
2.340000 -0.542907
2.360000 -0.804938
2.380000 -0.956835
2.400000 -0.969062
2.420000 -1.763540
2.440000 -1.629800
2.460000 -1.763376
2.480000 -1.986885
2.500000 -1.815159
2.520000 -1.304021
2.540000 -1.060269
2.560000 -1.151175
2.580000 -0.125246
2.600000 -0.350666
2.620000 -0.006205
2.640000 0.409014
2.660000 0.382090
2.680000 0.320996
2.700000 0.218711
2.720000 0.702889
2.740000 1.067371
2.760000 0.883676
2.780000 0.610368
2.800000 0.664391
2.820000 0.883116
2.840000 -0.037523
2.860000 0.102343
2.880000 0.631371
2.900000 0.259901
2.920000 -0.063701
2.940000 0.141226
2.960000 0.002918
2.980000 0.763487
3.000000 0.120750
3.020000 1.215525
3.040000 1.354160
3.060000 0.935050
3.080000 0.898151
3.100000 1.726756
3.120000 1.308912
3.140000 1.537892
3.160000 1.303440
3.180000 1.442885
3.200000 1.319862
3.220000 0.475973
3.240000 0.589797
3.260000 0.544649
3.280000 0.385700
3.300000 0.212688
3.320000 -0.088949
3.340000 -0.841775
3.360000 -0.652351
3.380000 -1.068508
3.400000 -0.963318
3.420000 -1.539100
3.440000 -1.521784
3.460000 -1.105703
3.480000 -1.311095
3.500000 -1.552438
3.520000 -1.329230
3.540000 -0.726134
3.560000 -0.686218
3.580000 -0.585860
3.600000 -0.268263
3.620000 0.025169
3.640000 0.066963
3.660000 0.332064
3.680000 -0.007807
3.700000 0.776530
3.720000 0.660876
3.740000 0.722263
3.760000 0.598961
3.780000 0.169825
3.800000 0.359099
3.820000 0.174459
3.840000 0.037562
3.860000 0.547720
3.880000 0.587833
3.900000 0.612495
3.920000 0.305067
3.940000 0.415666
3.960000 0.056643
3.980000 0.284842
4.000000 0.283517
4.020000 0.274994
4.040000 1.099955
4.060000 1.205833
4.080000 1.615293
4.100000 1.148340
4.120000 0.842416
4.140000 1.229825
4.160000 1.059159
4.180000 1.441346
4.200000 1.551181
4.220000 1.314096
4.240000 0.619632
4.260000 0.433915
4.280000 0.173796
4.300000 -0.086280
4.320000 -0.641863
4.340000 -1.095692
4.360000 -1.296427
4.380000 -0.906136
4.400000 -1.657777
4.420000 -1.488468
4.440000 -1.068069
4.460000 -1.725078
4.480000 -1.815086
4.500000 -1.142262
4.520000 -1.596012
4.540000 -0.861357
4.560000 -1.001601
4.580000 -0.798917
4.600000 -0.882130
4.620000 -0.127448
4.640000 -0.086625
4.660000 0.463992
4.680000 0.110359
4.700000 0.340084
4.720000 0.214388
4.740000 0.373932
4.760000 0.854448
4.780000 0.479055
4.800000 1.009229
4.820000 0.661264
4.840000 0.184915
4.860000 0.227685
4.880000 -0.036095
4.900000 0.547944
4.920000 0.686064
4.940000 0.238824
4.960000 0.738175
4.980000 0.021762
|
|
IP Logged |
|
|
|
hanzack
Forum Senior Member
lalala...I'm singing a song of arimaa.
Gender:
Posts: 35
|
|
Re: My hobby (second season)
« Reply #2 on: Feb 1st, 2011, 12:48am » |
Quote Modify
|
here is the function to fit
y(t) = x1 + x2 cos(OMEGA*t) + x3 sin(OMEGA*t) + x4 cos(2*OMEGA*t) + x5 sin(2*OMEGA*t)
yay 
|
|
IP Logged |
|
|
|
|
Home
Help
Search
Members
Login
Register
Reply
Notify of replies
Send Topic
Print