lda2.c

LDA程式設計

	gsl_vector_set (class1_v21, 2, LDAFea1[20].mean[2]);
	gsl_vector_set (class1_v21, 3, LDAFea1[20].std[0]);
	gsl_vector_set (class1_v21, 4, LDAFea1[20].std[1]);
	gsl_vector_set (class1_v21, 5, LDAFea1[20].std[2]);
	gsl_vector_set (class1_v21, 6, LDAFea1[20].corrcoef[0]);
	gsl_vector_set (class1_v21, 7, LDAFea1[20].corrcoef[1]);
	gsl_vector_set (class1_v21, 8, LDAFea1[20].corrcoef[2]);
	gsl_vector_set (class1_v21, 9, LDAFea1[20].mad[0]);
	gsl_vector_set (class1_v21, 10, LDAFea1[20].mad[1]);
	gsl_vector_set (class1_v21, 11, LDAFea1[20].mad[2]);
	gsl_vector_set (class1_v21, 12, LDAFea1[20].rms[0]);
	gsl_vector_set (class1_v21, 13, LDAFea1[20].rms[1]);
	gsl_vector_set (class1_v21, 14, LDAFea1[20].rms[2]);
	gsl_vector_set (class1_v21, 15, LDAFea1[20].var[0]);
	gsl_vector_set (class1_v21, 16, LDAFea1[20].var[1]);
	gsl_vector_set (class1_v21, 17, LDAFea1[20].var[2]);

	gsl_vector *class1_v22 = gsl_vector_alloc(18);
	gsl_vector_set (class1_v22, 0, LDAFea1[21].mean[0]);
	gsl_vector_set (class1_v22, 1, LDAFea1[21].mean[1]);
	gsl_vector_set (class1_v22, 2, LDAFea1[21].mean[2]);
	gsl_vector_set (class1_v22, 3, LDAFea1[21].std[0]);
	gsl_vector_set (class1_v22, 4, LDAFea1[21].std[1]);
	gsl_vector_set (class1_v22, 5, LDAFea1[21].std[2]);
	gsl_vector_set (class1_v22, 6, LDAFea1[21].corrcoef[0]);
	gsl_vector_set (class1_v22, 7, LDAFea1[21].corrcoef[1]);
	gsl_vector_set (class1_v22, 8, LDAFea1[21].corrcoef[2]);
	gsl_vector_set (class1_v22, 9, LDAFea1[21].mad[0]);
	gsl_vector_set (class1_v22, 10, LDAFea1[21].mad[1]);
	gsl_vector_set (class1_v22, 11, LDAFea1[21].mad[2]);
	gsl_vector_set (class1_v22, 12, LDAFea1[21].rms[0]);
	gsl_vector_set (class1_v22, 13, LDAFea1[21].rms[1]);
	gsl_vector_set (class1_v22, 14, LDAFea1[21].rms[2]);
	gsl_vector_set (class1_v22, 15, LDAFea1[21].var[0]);
	gsl_vector_set (class1_v22, 16, LDAFea1[21].var[1]);
	gsl_vector_set (class1_v22, 17, LDAFea1[21].var[2]);

	gsl_vector *class1_v23 = gsl_vector_alloc(18);
	gsl_vector_set (class1_v23, 0, LDAFea1[22].mean[0]);
	gsl_vector_set (class1_v23, 1, LDAFea1[22].mean[1]);
	gsl_vector_set (class1_v23, 2, LDAFea1[22].mean[2]);
	gsl_vector_set (class1_v23, 3, LDAFea1[22].std[0]);
	gsl_vector_set (class1_v23, 4, LDAFea1[22].std[1]);
	gsl_vector_set (class1_v23, 5, LDAFea1[22].std[2]);
	gsl_vector_set (class1_v23, 6, LDAFea1[22].corrcoef[0]);
	gsl_vector_set (class1_v23, 7, LDAFea1[22].corrcoef[1]);
	gsl_vector_set (class1_v23, 8, LDAFea1[22].corrcoef[2]);
	gsl_vector_set (class1_v23, 9, LDAFea1[22].mad[0]);
	gsl_vector_set (class1_v23, 10, LDAFea1[22].mad[1]);
	gsl_vector_set (class1_v23, 11, LDAFea1[22].mad[2]);
	gsl_vector_set (class1_v23, 12, LDAFea1[22].rms[0]);
	gsl_vector_set (class1_v23, 13, LDAFea1[22].rms[1]);
	gsl_vector_set (class1_v23, 14, LDAFea1[22].rms[2]);
	gsl_vector_set (class1_v23, 15, LDAFea1[22].var[0]);
	gsl_vector_set (class1_v23, 16, LDAFea1[22].var[1]);
	gsl_vector_set (class1_v23, 17, LDAFea1[22].var[2]);

	gsl_vector *class1_v24 = gsl_vector_alloc(18);
	gsl_vector_set (class1_v24, 0, LDAFea1[23].mean[0]);
	gsl_vector_set (class1_v24, 1, LDAFea1[23].mean[1]);
	gsl_vector_set (class1_v24, 2, LDAFea1[23].mean[2]);
	gsl_vector_set (class1_v24, 3, LDAFea1[23].std[0]);
	gsl_vector_set (class1_v24, 4, LDAFea1[23].std[1]);
	gsl_vector_set (class1_v24, 5, LDAFea1[23].std[2]);
	gsl_vector_set (class1_v24, 6, LDAFea1[23].corrcoef[0]);
	gsl_vector_set (class1_v24, 7, LDAFea1[23].corrcoef[1]);
	gsl_vector_set (class1_v24, 8, LDAFea1[23].corrcoef[2]);
	gsl_vector_set (class1_v24, 9, LDAFea1[23].mad[0]);
	gsl_vector_set (class1_v24, 10, LDAFea1[23].mad[1]);
	gsl_vector_set (class1_v24, 11, LDAFea1[23].mad[2]);
	gsl_vector_set (class1_v24, 12, LDAFea1[23].rms[0]);
	gsl_vector_set (class1_v24, 13, LDAFea1[23].rms[1]);
	gsl_vector_set (class1_v24, 14, LDAFea1[23].rms[2]);
	gsl_vector_set (class1_v24, 15, LDAFea1[23].var[0]);
	gsl_vector_set (class1_v24, 16, LDAFea1[23].var[1]);
	gsl_vector_set (class1_v24, 17, LDAFea1[23].var[2]);

	gsl_vector *class1_v25 = gsl_vector_alloc(18);
	gsl_vector_set (class1_v25, 0, LDAFea1[24].mean[0]);
	gsl_vector_set (class1_v25, 1, LDAFea1[24].mean[1]);
	gsl_vector_set (class1_v25, 2, LDAFea1[24].mean[2]);
	gsl_vector_set (class1_v25, 3, LDAFea1[24].std[0]);
	gsl_vector_set (class1_v25, 4, LDAFea1[24].std[1]);
	gsl_vector_set (class1_v25, 5, LDAFea1[24].std[2]);
	gsl_vector_set (class1_v25, 6, LDAFea1[24].corrcoef[0]);
	gsl_vector_set (class1_v25, 7, LDAFea1[24].corrcoef[1]);
	gsl_vector_set (class1_v25, 8, LDAFea1[24].corrcoef[2]);
	gsl_vector_set (class1_v25, 9, LDAFea1[24].mad[0]);
	gsl_vector_set (class1_v25, 10, LDAFea1[24].mad[1]);
	gsl_vector_set (class1_v25, 11, LDAFea1[24].mad[2]);
	gsl_vector_set (class1_v25, 12, LDAFea1[24].rms[0]);
	gsl_vector_set (class1_v25, 13, LDAFea1[24].rms[1]);
	gsl_vector_set (class1_v25, 14, LDAFea1[24].rms[2]);
	gsl_vector_set (class1_v25, 15, LDAFea1[24].var[0]);
	gsl_vector_set (class1_v25, 16, LDAFea1[24].var[1]);
	gsl_vector_set (class1_v25, 17, LDAFea1[24].var[2]);

/* class2*/
	gsl_vector *class2_v1 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v1, 0, LDAFea2[0].mean[0]);
	gsl_vector_set (class2_v1, 1, LDAFea2[0].mean[1]);
	gsl_vector_set (class2_v1, 2, LDAFea2[0].mean[2]);
	gsl_vector_set (class2_v1, 3, LDAFea2[0].std[0]);
	gsl_vector_set (class2_v1, 4, LDAFea2[0].std[1]);
	gsl_vector_set (class2_v1, 5, LDAFea2[0].std[2]);
	gsl_vector_set (class2_v1, 6, LDAFea2[0].corrcoef[0]);
	gsl_vector_set (class2_v1, 7, LDAFea2[0].corrcoef[1]);
	gsl_vector_set (class2_v1, 8, LDAFea2[0].corrcoef[2]);
	gsl_vector_set (class2_v1, 9, LDAFea2[0].mad[0]);
	gsl_vector_set (class2_v1, 10, LDAFea2[0].mad[1]);
	gsl_vector_set (class2_v1, 11, LDAFea2[0].mad[2]);
	gsl_vector_set (class2_v1, 12, LDAFea2[0].rms[0]);
	gsl_vector_set (class2_v1, 13, LDAFea2[0].rms[1]);
	gsl_vector_set (class2_v1, 14, LDAFea2[0].rms[2]);
	gsl_vector_set (class2_v1, 15, LDAFea2[0].var[0]);
	gsl_vector_set (class2_v1, 16, LDAFea2[0].var[1]);
	gsl_vector_set (class2_v1, 17, LDAFea2[0].var[2]);

	gsl_vector *class2_v2 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v2, 0, LDAFea2[1].mean[0]);
	gsl_vector_set (class2_v2, 1, LDAFea2[1].mean[1]);
	gsl_vector_set (class2_v2, 2, LDAFea2[1].mean[2]);
	gsl_vector_set (class2_v2, 3, LDAFea2[1].std[0]);
	gsl_vector_set (class2_v2, 4, LDAFea2[1].std[1]);
	gsl_vector_set (class2_v2, 5, LDAFea2[1].std[2]);
	gsl_vector_set (class2_v2, 6, LDAFea2[1].corrcoef[0]);
	gsl_vector_set (class2_v2, 7, LDAFea2[1].corrcoef[1]);
	gsl_vector_set (class2_v2, 8, LDAFea2[1].corrcoef[2]);
	gsl_vector_set (class2_v2, 9, LDAFea2[1].mad[0]);
	gsl_vector_set (class2_v2, 10, LDAFea2[1].mad[1]);
	gsl_vector_set (class2_v2, 11, LDAFea2[1].mad[2]);
	gsl_vector_set (class2_v2, 12, LDAFea2[1].rms[0]);
	gsl_vector_set (class2_v2, 13, LDAFea2[1].rms[1]);
	gsl_vector_set (class2_v2, 14, LDAFea2[1].rms[2]);
	gsl_vector_set (class2_v2, 15, LDAFea2[1].var[0]);
	gsl_vector_set (class2_v2, 16, LDAFea2[1].var[1]);
	gsl_vector_set (class2_v2, 17, LDAFea2[1].var[2]);

	gsl_vector *class2_v3 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v3, 0, LDAFea2[2].mean[0]);
	gsl_vector_set (class2_v3, 1, LDAFea2[2].mean[1]);
	gsl_vector_set (class2_v3, 2, LDAFea2[2].mean[2]);
	gsl_vector_set (class2_v3, 3, LDAFea2[2].std[0]);
	gsl_vector_set (class2_v3, 4, LDAFea2[2].std[1]);
	gsl_vector_set (class2_v3, 5, LDAFea2[2].std[2]);
	gsl_vector_set (class2_v3, 6, LDAFea2[2].corrcoef[0]);
	gsl_vector_set (class2_v3, 7, LDAFea2[2].corrcoef[1]);
	gsl_vector_set (class2_v3, 8, LDAFea2[2].corrcoef[2]);
	gsl_vector_set (class2_v3, 9, LDAFea2[2].mad[0]);
	gsl_vector_set (class2_v3, 10, LDAFea2[2].mad[1]);
	gsl_vector_set (class2_v3, 11, LDAFea2[2].mad[2]);
	gsl_vector_set (class2_v3, 12, LDAFea2[2].rms[0]);
	gsl_vector_set (class2_v3, 13, LDAFea2[2].rms[1]);
	gsl_vector_set (class2_v3, 14, LDAFea2[2].rms[2]);
	gsl_vector_set (class2_v3, 15, LDAFea2[2].var[0]);
	gsl_vector_set (class2_v3, 16, LDAFea2[2].var[1]);
	gsl_vector_set (class2_v3, 17, LDAFea2[2].var[2]);

	gsl_vector *class2_v4 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v4, 0, LDAFea2[3].mean[0]);
	gsl_vector_set (class2_v4, 1, LDAFea2[3].mean[1]);
	gsl_vector_set (class2_v4, 2, LDAFea2[3].mean[2]);
	gsl_vector_set (class2_v4, 3, LDAFea2[3].std[0]);
	gsl_vector_set (class2_v4, 4, LDAFea2[3].std[1]);
	gsl_vector_set (class2_v4, 5, LDAFea2[3].std[2]);
	gsl_vector_set (class2_v4, 6, LDAFea2[3].corrcoef[0]);
	gsl_vector_set (class2_v4, 7, LDAFea2[3].corrcoef[1]);
	gsl_vector_set (class2_v4, 8, LDAFea2[3].corrcoef[2]);
	gsl_vector_set (class2_v4, 9, LDAFea2[3].mad[0]);
	gsl_vector_set (class2_v4, 10, LDAFea2[3].mad[1]);
	gsl_vector_set (class2_v4, 11, LDAFea2[3].mad[2]);
	gsl_vector_set (class2_v4, 12, LDAFea2[3].rms[0]);
	gsl_vector_set (class2_v4, 13, LDAFea2[3].rms[1]);
	gsl_vector_set (class2_v4, 14, LDAFea2[3].rms[2]);
	gsl_vector_set (class2_v4, 15, LDAFea2[3].var[0]);
	gsl_vector_set (class2_v4, 16, LDAFea2[3].var[1]);
	gsl_vector_set (class2_v4, 17, LDAFea2[3].var[2]);

	gsl_vector *class2_v5 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v5, 0, LDAFea2[4].mean[0]);
	gsl_vector_set (class2_v5, 1, LDAFea2[4].mean[1]);
	gsl_vector_set (class2_v5, 2, LDAFea2[4].mean[2]);
	gsl_vector_set (class2_v5, 3, LDAFea2[4].std[0]);
	gsl_vector_set (class2_v5, 4, LDAFea2[4].std[1]);
	gsl_vector_set (class2_v5, 5, LDAFea2[4].std[2]);
	gsl_vector_set (class2_v5, 6, LDAFea2[4].corrcoef[0]);
	gsl_vector_set (class2_v5, 7, LDAFea2[4].corrcoef[1]);
	gsl_vector_set (class2_v5, 8, LDAFea2[4].corrcoef[2]);
	gsl_vector_set (class2_v5, 9, LDAFea2[4].mad[0]);
	gsl_vector_set (class2_v5, 10, LDAFea2[4].mad[1]);
	gsl_vector_set (class2_v5, 11, LDAFea2[4].mad[2]);
	gsl_vector_set (class2_v5, 12, LDAFea2[4].rms[0]);
	gsl_vector_set (class2_v5, 13, LDAFea2[4].rms[1]);
	gsl_vector_set (class2_v5, 14, LDAFea2[4].rms[2]);
	gsl_vector_set (class2_v5, 15, LDAFea2[4].var[0]);
	gsl_vector_set (class2_v5, 16, LDAFea2[4].var[1]);
	gsl_vector_set (class2_v5, 17, LDAFea2[4].var[2]);

	gsl_vector *class2_v6 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v6, 0, LDAFea2[5].mean[0]);
	gsl_vector_set (class2_v6, 1, LDAFea2[5].mean[1]);
	gsl_vector_set (class2_v6, 2, LDAFea2[5].mean[2]);
	gsl_vector_set (class2_v6, 3, LDAFea2[5].std[0]);
	gsl_vector_set (class2_v6, 4, LDAFea2[5].std[1]);
	gsl_vector_set (class2_v6, 5, LDAFea2[5].std[2]);
	gsl_vector_set (class2_v6, 6, LDAFea2[5].corrcoef[0]);
	gsl_vector_set (class2_v6, 7, LDAFea2[5].corrcoef[1]);
	gsl_vector_set (class2_v6, 8, LDAFea2[5].corrcoef[2]);
	gsl_vector_set (class2_v6, 9, LDAFea2[5].mad[0]);
	gsl_vector_set (class2_v6, 10, LDAFea2[5].mad[1]);
	gsl_vector_set (class2_v6, 11, LDAFea2[5].mad[2]);
	gsl_vector_set (class2_v6, 12, LDAFea2[5].rms[0]);
	gsl_vector_set (class2_v6, 13, LDAFea2[5].rms[1]);
	gsl_vector_set (class2_v6, 14, LDAFea2[5].rms[2]);
	gsl_vector_set (class2_v6, 15, LDAFea2[5].var[0]);
	gsl_vector_set (class2_v6, 16, LDAFea2[5].var[1]);
	gsl_vector_set (class2_v6, 17, LDAFea2[5].var[2]);

	gsl_vector *class2_v7 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v7, 0, LDAFea2[6].mean[0]);
	gsl_vector_set (class2_v7, 1, LDAFea2[6].mean[1]);
	gsl_vector_set (class2_v7, 2, LDAFea2[6].mean[2]);
	gsl_vector_set (class2_v7, 3, LDAFea2[6].std[0]);
	gsl_vector_set (class2_v7, 4, LDAFea2[6].std[1]);
	gsl_vector_set (class2_v7, 5, LDAFea2[6].std[2]);
	gsl_vector_set (class2_v7, 6, LDAFea2[6].corrcoef[0]);
	gsl_vector_set (class2_v7, 7, LDAFea2[6].corrcoef[1]);
	gsl_vector_set (class2_v7, 8, LDAFea2[6].corrcoef[2]);
	gsl_vector_set (class2_v7, 9, LDAFea2[6].mad[0]);
	gsl_vector_set (class2_v7, 10, LDAFea2[6].mad[1]);
	gsl_vector_set (class2_v7, 11, LDAFea2[6].mad[2]);
	gsl_vector_set (class2_v7, 12, LDAFea2[6].rms[0]);
	gsl_vector_set (class2_v7, 13, LDAFea2[6].rms[1]);
	gsl_vector_set (class2_v7, 14, LDAFea2[6].rms[2]);
	gsl_vector_set (class2_v7, 15, LDAFea2[6].var[0]);
	gsl_vector_set (class2_v7, 16, LDAFea2[6].var[1]);
	gsl_vector_set (class2_v7, 17, LDAFea2[6].var[2]);

	gsl_vector *class2_v8 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v8, 0, LDAFea2[7].mean[0]);
	gsl_vector_set (class2_v8, 1, LDAFea2[7].mean[1]);
	gsl_vector_set (class2_v8, 2, LDAFea2[7].mean[2]);
	gsl_vector_set (class2_v8, 3, LDAFea2[7].std[0]);
	gsl_vector_set (class2_v8, 4, LDAFea2[7].std[1]);
	gsl_vector_set (class2_v8, 5, LDAFea2[7].std[2]);
	gsl_vector_set (class2_v8, 6, LDAFea2[7].corrcoef[0]);
	gsl_vector_set (class2_v8, 7, LDAFea2[7].corrcoef[1]);
	gsl_vector_set (class2_v8, 8, LDAFea2[7].corrcoef[2]);
	gsl_vector_set (class2_v8, 9, LDAFea2[7].mad[0]);
	gsl_vector_set (class2_v8, 10, LDAFea2[7].mad[1]);
	gsl_vector_set (class2_v8, 11, LDAFea2[7].mad[2]);
	gsl_vector_set (class2_v8, 12, LDAFea2[7].rms[0]);
	gsl_vector_set (class2_v8, 13, LDAFea2[7].rms[1]);
	gsl_vector_set (class2_v8, 14, LDAFea2[7].rms[2]);
	gsl_vector_set (class2_v8, 15, LDAFea2[7].var[0]);
	gsl_vector_set (class2_v8, 16, LDAFea2[7].var[1]);
	gsl_vector_set (class2_v8, 17, LDAFea2[7].var[2]);

	gsl_vector *class2_v9 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v9, 0, LDAFea2[8].mean[0]);
	gsl_vector_set (class2_v9, 1, LDAFea2[8].mean[1]);
	gsl_vector_set (class2_v9, 2, LDAFea2[8].mean[2]);
	gsl_vector_set (class2_v9, 3, LDAFea2[8].std[0]);
	gsl_vector_set (class2_v9, 4, LDAFea2[8].std[1]);
	gsl_vector_set (class2_v9, 5, LDAFea2[8].std[2]);
	gsl_vector_set (class2_v9, 6, LDAFea2[8].corrcoef[0]);
	gsl_vector_set (class2_v9, 7, LDAFea2[8].corrcoef[1]);
	gsl_vector_set (class2_v9, 8, LDAFea2[8].corrcoef[2]);
	gsl_vector_set (class2_v9, 9, LDAFea2[8].mad[0]);
	gsl_vector_set (class2_v9, 10, LDAFea2[8].mad[1]);
	gsl_vector_set (class2_v9, 11, LDAFea2[8].mad[2]);
	gsl_vector_set (class2_v9, 12, LDAFea2[8].rms[0]);
	gsl_vector_set (class2_v9, 13, LDAFea2[8].rms[1]);
	gsl_vector_set (class2_v9, 14, LDAFea2[8].rms[2]);
	gsl_vector_set (class2_v9, 15, LDAFea2[8].var[0]);
	gsl_vector_set (class2_v9, 16, LDAFea2[8].var[1]);
	gsl_vector_set (class2_v9, 17, LDAFea2[8].var[2]);

	gsl_vector *class2_v10 = gsl_vector_alloc(18);
	gsl_vector_set (class2_v10, 0, LDAFea2[9].mean[0]);
	gsl_vector_set (class2_v10, 1, LDAFea2[9].mean[1]);
	gsl_vector_set (class2_v10, 2, LDAFea2[9].mean[2]);
	gsl_vector_set (class2_v10, 3, LDAFea2[9].std[0]);
	gsl_vector_set (class2_v10, 4, LDAFea2[9].std[1]);
	gsl_vector_set (class2_v10, 5, LDAFea2[9].std[2]);
	gsl_vector_set (class2_v10, 6, LDAFea2[9].corrcoef[0]);
	gsl_vector_set (class2_v10, 7, LDAFea2[9].corrcoef[1]);
	gsl_vector_set (class2_v10, 8, LDAFea2[9].corrcoef[2]);
	gsl_vector_set (class2_v10, 9, LDAFea2[9].mad[0]);
	gsl_vector_set (class2_v10, 10, LDAFea2[9].mad[1]);
	gsl_vector_set (class2_v10, 11, LDAFea2