void Look_traj() { //make_plot_crd_diff 110 1605 0 vector<double> X; Load_traj(110, 1605, 0, X); int szX = X.size(); if (szX == 0) return; vector<double> diffX(szX-1); for (int i = 0; i < szX-1; ++i) { diffX[i] = X[i+1] - X[i]; } { //------------------------------------------------------------------ char filter[] = TEXT("Plot Text File (*.txt)\0*.txt\0") TEXT("All Files (*.*)\0*.*\0"); char outfile_name[1048]; sprintf(outfile_name, "\0"); DWORD nFilterIndex; if (SaveFileDlg(0, outfile_name, filter, nFilterIndex) == S_OK) { cout << "now saving file " << outfile_name << endl; ofstream ofile; ofile.open(outfile_name, ios::out); for (int i = 0; i < szX-1; ++i) { ofile << X[i] << "," << diffX[i] << endl; } ofile.close(); } } for (int N = 10; N < 100; N += 5) Statistika (diffX, "diffX", N, calcMathOzh); for (int N = 10; N < 100; N += 5) Statistika (X, "linearTrendParameterA", N, calcLinearTrendParameterA); }
int main() { int pil; gotoxy(20,0); printf("Data Mahasiswa JTK Polban Kelas 1-B 2011\n\n"); printf("1. Input data mahasiswa\n"); printf("2. Tampilkan data mahasiwa\n"); printf("3. Daftar Mahasiswa yang tinggal di Bandung\n"); printf("4. Daftar Mahasiswa laki-laki yang tingginya lebih dari 170cm\n"); printf("5. Statistika Mahasiswa berdasarkan kota\n\n"); printf("Masukkan pilihan anda = "); scanf("%d",&pil); switch(pil) { case 1 : Input(); break; case 2 : Tampil(); break; case 3 : Bandung(); break; case 4 : CariMhs(); break; case 5 : Statistika(); break; default : printf("Nomor yang anda masukkan salah!\n"); } return 0; }
int Look_txt() { TCHAR filter[] = TEXT("Ghemical MD results File (*.txt)\0*.txt\0") TEXT("All Files (*.*)\0*.*\0"); TCHAR fpath[1024]; TCHAR filename[1024]; sprintf(filename, "\0"); { DWORD nFilterIndex; vector<string> names; vector<string> *pnames = &names; vector<vector<double> > vectors; vectors.reserve(2000000); while (OpenFileDlg(0, filter, fpath, nFilterIndex) == S_OK) { ReadDatFile(NULL, fpath, filename, &vectors, pnames); pnames = NULL; printf("\nfilename %s\n\n", filename); int cols = names.size(); int rows = vectors.size(); #if WRITE_LOCKED_FORCES int cMom = 4 - 1; int cVx = 5 - 1; int cFxup = 14 - 1; int cFxdw = 17 - 1; int cVxup = 8 - 1; int cVxdw = 11 - 1; #endif #if WRITE_WORKED_FORCES int cMom = 4 - 1; int cVx = 5 - 1; int cVxup = 14 - 1; int cVxdw = 17 - 1; int cVx_wk_up = 8 - 1; int cVx_wk_dw = 11 - 1; int cFx_wk_up = 20 - 1; int cFx_wk_dw = 23 - 1; #endif vector<double> means(cols, 0.0); printf("vectors.size() = %d\n",rows); printf("names.size() = %d\n", cols); for (vector<vector<double> >::iterator it = vectors.begin(); it != vectors.end(); it++) { for (int c = 0; c < cols; c++) { means[c] += (*it).operator [](c); } } for (int c = 0; c < cols; c++) { means[c] /= rows; printf("mean(%s) = %f\n", names[c].c_str(), means[c]); } #if WRITE_LOCKED_FORCES || WRITE_WORKED_FORCES int r0 = 0; cout << "enter r0\n"; cin >> r0; #endif #if WRITE_LOCKED_FORCES vector<double> dF(rows-r0); for (int r = r0; r < rows; r++) { dF[r-r0] = vectors[r][cFxup] - vectors[r][cFxdw]; } Statistika (dF, "dF"); vector<double> Mom(rows-r0); for (r = r0; r < rows; r++) { Mom[r-r0] = vectors[r][cMom]; } Statistika (Mom, "Mom"); vector<double> dV(rows-r0); for (r = r0; r < rows; r++) { dV[r-r0] = vectors[r][cVxup] - vectors[r][cVxdw]; } Statistika (dV, "dV"); vector<double> Vx(rows-r0); for (r = r0; r < rows; r++) { Vx[r-r0] = vectors[r][cVx]; } Statistika (Vx, "Vx"); #endif #if WRITE_WORKED_FORCES vector<double> dF_wk(rows-r0); for (int r = r0; r < rows; r++) { dF_wk[r-r0] = vectors[r][cFx_wk_up] - vectors[r][cFx_wk_dw]; } Statistika (dF_wk, "dF_wk"); vector<double> dV_wk(rows-r0); for (r = r0; r < rows; r++) { dV_wk[r-r0] = vectors[r][cVx_wk_up] - vectors[r][cVx_wk_dw]; } Statistika (dV_wk, "dV_wk"); //if (!worked[n1]) vector<double> Mom(rows-r0); for (r = r0; r < rows; r++) { Mom[r-r0] = vectors[r][cMom]; } Statistika (Mom, "Mom"); vector<double> dV(rows-r0); for (r = r0; r < rows; r++) { dV[r-r0] = vectors[r][cVxup] - vectors[r][cVxdw]; } Statistika (dV, "dV"); vector<double> Vx(rows-r0); for (r = r0; r < rows; r++) { Vx[r-r0] = vectors[r][cVx]; } Statistika (Vx, "Vx"); #endif } } /*else { DWORD nFilterIndex; if (SaveFileDlg(0, filename, filter, nFilterIndex) == S_OK) { SetDlgItemText(ref->hDlg,IDC_EDIT_TRAJFILE2, filename); } }*/ printf("Hello World!\n"); return 0; }